macapype icon indicating copy to clipboard operation
macapype copied to clipboard

Options when subjects fail in crop_aladin stage

Open julfou81 opened this issue 3 years ago • 1 comments

Inspired from subject sub-032159 in site-mountsinai-S from PRIME-DE, which fail miserably since the first step of macapype.

Here is for instance the image after the first stage: Capture d’écran 2022-09-22 à 18 03 45

  • option 1 -- Prepare the data with FSL-FAST: fast -l 8 -B sub-032159_ses-001_run-1_T1w.nii.gz -- Apply the reg_aladin command used in macapype to the bias corrected image:
reg_aladin -aff sub-032159_ses-001_run-1_T1w_aff.txt -flo sub-032159_ses-001_run-1_T1w_restore.nii.gz -ln 12 -lp 10 -nac -noSym -omp 1  -ref /Users/jsein/Documents/Centre_IRMf/DATA/BIDS/Covoprim/inia19-t1.nii -res sub-032159_ses-001_run-1_T1w_restore_res.nii.gz -rigOnly -smooR 1.000000

-- Apply the calculated transform to the original image

reg_resample -flo sub-032159_ses-001_run-1_T1w.nii.gz -ref /Users/jsein/Documents/Centre_IRMf/DATA/BIDS/Covoprim/inia19-t1.nii  -res test.nii.gz -trans sub-032159_ses-001_run-1_T1w_aff.tx

This is the result: Capture d’écran 2022-09-22 à 18 09 02

  • option 2 (much faster) -- change the parameters of reg_aladin on the original image:
reg_aladin -aff sub-032159_ses-001_run-1_T1w_aff.txt -flo sub-032159_ses-001_run-1_T1w.nii.gz -ln 8 -lp 10 -nac -noSym -omp 1  -ref /Users/jsein/Documents/Centre_IRMf/DATA/BIDS/Covoprim/inia19-t1.nii -res sub-032159_ses-001_run-1_T1w_res.nii.gz -rigOnly -smooR 1.000000

=> same result

julfou81 avatar Sep 22 '22 16:09 julfou81

With a little more exploration, using sites ucdavis and mount sinai-S, I had to choose between two sets of parameters for reg_aladin: (with ref inia19)

  1. reg_aladin -ln 12 -lp 10 -nac -noSym -omp 1 -ref -rigOnly -smooR 1 => works for all subjects of ucdavis and all but two subjects of mount sinai-S

  2. reg_aladin -ln 5 -lp 10 -nac -noSym -omp 1 -ref -rigOnly -smooR 1 => works for all subjects of ucdavis, mount sinai-S ( even for the two subjects of mount sinai-S that weren't working)

For INT macaques, option worked for all the one I tried, except Galileo, which needed to be cropped before and use of option 1.

julfou81 avatar Sep 23 '22 22:09 julfou81