Dear Jesper,
Thanks for this very useful information. I was able to run a similar script to the one you suggested and corrected in just the PE we collect in. I will add the -interp argument to the next scripts. I have questions about inspecting the images and ensuring the correction procedure did what we were hoping they would do. Do you have general suggestions about how to inspect the data afterwards? What should I look out for, and sanity checks can I perform?
When I inspect the corrected images alone, it is difficult for me to tell what to look out for. I created a difference map between corrected and uncorrected images, and I see that values have changed across the entire brain, and seem to be more extreme around the eye sockets, brainstem, and the cerebellum area which makes sense to me. I am most concerned about correcting distortion around the vmPFC/OFC areas. I also have questions about the set up file b02b0.cnf. Can you tell me generally how the default values are set up? Is it assumed these should suffice for most datasets, and we shouldn't edit unless doing very specific tweaking?
Thanks very much,
Helen
b02b0.cnf:
# Resolution (knot-spacing) of warps in mm
--warpres=20,16,14,12,10,6,4,4,4
# Subsampling level (a value of 2 indicates that a 2x2x2 neighbourhood is collapsed to 1 voxel)
--subsamp=2,2,2,2,2,1,1,1,1
# FWHM of gaussian smoothing
--fwhm=8,6,4,3,3,2,1,0,0
# Maximum number of iterations
--miter=5,5,5,5,5,10,10,20,20
# Relative weight of regularisation
--lambda=0.005,0.001,0.0001,0.000015,0.000005,0.0000005,0.00000005,0.0000000005,0.00000000001
# If set to 1 lambda is multiplied by the current average squared difference
--ssqlambda=1
# Regularisation model
--regmod=bending_energy
# If set to 1 movements are estimated along with the field
--estmov=1,1,1,1,1,0,0,0,0
# 0=Levenberg-Marquardt, 1=Scaled Conjugate Gradient
--minmet=0,0,0,0,0,1,1,1,1
# Quadratic or cubic splines
--splineorder=3
# Precision for calculation and storage of Hessian
--numprec=double
# Linear or spline interpolation
--interp=spline
# If set to 1 the images are individually scaled to a common mean intensity
--scale=1
|