Hi there,
I have been wondering about how to compare different volreg methods to see which produces the best alignment. I have been looking in the GUI between and within blocks/ overlaying them etc. But I wanted a more quantifiable method
I was wondering if this approach, related to running volreg twice and comparing enorms makes sense??
This is what I did
1. Run my volreg1 using chosen method (generates enorm and dfile 1D files)
2. Rerun registration (volreg2) on pb01.volreg blocks to generate emorn and dfiles for this ‘fixed’ pb01 block – to see how much redisual movement is left after you did your volreg1
3. Plot enorms for volreg2 against each other on the same plot (to see how much movement remains after volreg1) and compare enorm means between volreg 1 and 2
(example code below)
Could someone let me know if this is a legitimate way to do this?
Or, is there a better/ recommended way for quantifying volreg success?
Thank you for your help
Harriet
-----------
1.
afni_proc.py -subj_id $subj \
-dsets $data_dir/P007_FR_RUN* \
-tcat_remove_first_trs 0 \
-blocks volreg \
-volreg_align_to MIN_OUTLIER \
-volreg_interp -Fourier \
-volreg_warp_dxyz 0.8 \
-volreg_motsim
-----------------------
2.
set cap = 1
# Re run volreg on volreg'd files to get motion parameters after correction
foreach run ( 01 02 03 04 05 06 )
3dvolreg -verbose -zpad 1 \
-base vr_base_min_outlier+orig \
-1Dfile dfile.r{$run}_2.1D \
-prefix pb01.P007_FR_{$cap}.r{$run}.volreg_2 \
-Fourier \
pb01.P007_FR_{$cap}.r{$run}.volreg+orig
end
# combine new dfiles
cat dfile.r01_2.1D dfile.r02_2.1D dfile.r03_2.1D dfile.r04_2.1D dfile.r05_2.1D dfile.r06_2.1D > dfile_rall_2.1D
----
3.
# Calc enorm
1d_tool.py -infile dfile_rall_2.1D -set_nruns 6 -derivative -collapse_cols euclidean_norm -write motion_P007_FR_{$cap}_enorm_2.1D
# Plot
1dplot -one motion_P007_FR_{$cap}_enorm.1D motion_P007_FR_{$cap}_enorm_2.1D "1D: 2000@0.3" &