Hi - I would not use the 5th dimension as this is not yet fully supported in most implementations of NIFTI - I would find a workaround for that and re-test. Cheers. On 26 Jul 2007, at 11:11, Alle Meije Wink wrote: > Hi -thanks for the quick reply! > > What I saw (with linear interpolation earlier, I thought sinc would be > better for the real analysis) was that also took very long. I thought > that something was happening outside just resampling. > > In the analysis there are some 4D files of 10 volumes -actually > they are > 5D, because they are not time series so we put everything in dim[5]- > which are in <double> format. So they are quite big... > > But then, I assume that files are resampled on a 3D (volume by volume) > basis? Why would it be necessary then to load the whole series of > volume > into the memory? > > My computer has 2 G of memory, so 91*109*91*8*10 = 72210320 = 69 M > should not be a problem? > > I'll have another test with trilinear interpolation to see what that > does to the speed. I'm afraid about what it will do to the quality > though... > > I have written B-spline interpolation code before (whose interpolants > for degrees >2 are almost sinc functions) and that didn't seem to > suffer > from increasing the degree that much. Are you sure this is what > makes it > so slow? > > Thanks! > Alle Meije > > Steve Smith wrote: >> Hi - I would have expected changing the interpolation to have a large >> effect. It is possible instead that applying flirt to the 4D >> dataset is >> slow because doing that all in one is requiring too much RAM and your >> computer is doing lots of swapping - in which case I would recommend >> using avwsplit to split it down to 3D files, applying the >> transform with >> flirt, and then recombining with avwmerge. >> >> Cheers. > >> On 26 Jul 2007, at 10:41, Alle Meije Wink wrote: >>> Hi, >>> >>> I am using flirt for standard space mapping in my analyses. The >>> way I >>> have set it up is as follows: >>> >>> All single-subject time series (4D) are in NifTI format, the time >>> series >>> averages (3D) are called *NatAverage.nii, and the statistic >>> images (in >>> native space) are all called Nat*.nii. >>> >>> I align the *NatAverage.nii file to the template $TEMPL (which in >>> this >>> case is EPI.nii) and then use that transformation to apply to all >>> the >>> other Nat*.nii files as well. >>> >>> The bash scripting bit that does the jobe looks like this >>> >>> for PREFIX in PREFIXES; do >>> >>> # all Nat's except NatAverage >>> NATS=${PREFIX}Nat[BMNSV]*.nii >>> >>> echo "computing standard space mapping" >>> ${FLIRT} \ >>> -in ${PREFIX}NatAverage.nii -ref ${TEMPL} \ >>> -out ${PREFIX}StdAverage.nii -omat ${PREFIX}Affine.txt \ >>> -bins 256 -cost corratio -dof 12 \ >>> -searchrx -90 90 -searchry -90 90 -searchrz -90 90 \ >>> -interp sinc -sincwidth 7 -sincwindow hanning >>> >>> echo "applying standard space mapping" >>> for NAT in $NATS; do >>> >>> printf "." >>> # use FLIRT to apply transformation to other Nat files >>> ${FLIRT} -in ${NAT} -ref ${TEMPL} -out ${NAT//Nat/Std} \ >>> -applyxfm -init ${PREFIX}Affine.txt \ >>> -interp sinc -sincwidth 7 -sincwindow hanning >>> >>> done >>> printf "\n" >>> >>> done >>> >>> And it works... but it is amazingly slow for the second part (the >>> -applyxfm bit). It almost seems that flirt starts the >>> registration all >>> over again. >>> >>> I assumed that the other *Nat.nii images would just be resampled >>> using >>> the transformation parameters in *Affine.txt. That is correct, is >>> it? >>> >>> Is there any way to speed up the process? Changing the >>> interpolation to >>> linear seems to have only a marginal effect. Resampling a 90x100x90 >>> volume should be done well inside a second? >>> >>> Thanks for any hints >>> Alle Meije ------------------------------------------------------------------------ --- Stephen M. Smith, Professor of Biomedical Engineering Associate Director, Oxford University FMRIB Centre FMRIB, JR Hospital, Headington, Oxford OX3 9DU, UK +44 (0) 1865 222726 (fax 222717) [log in to unmask] http://www.fmrib.ox.ac.uk/~steve ------------------------------------------------------------------------ ---