@@ -155,7 +155,7 @@ fsleyes mt1.nii.gz -cm greyscale -a 100.0 label/template/PAM50_t2.nii.gz -cm gre
155155
156156
157157
158- # Computing MTR
158+ # Computing MTR using MT0/MT1 coregistration
159159# ======================================================================================================================
160160
161161# Register mt0->mt1 using z-regularized slicewise translations (algo=slicereg)
@@ -169,6 +169,27 @@ sct_compute_mtr -mt0 mt0_reg.nii.gz -mt1 mt1.nii.gz
169169
170170
171171
172+ # Contrast-agnostic registration
173+ # ======================================================================================================================
174+
175+ # 1. T2w preprocessing (cropping around spinal cord)
176+ cd ../t2
177+ sct_deepseg_sc -i t2.nii.gz -c t2 -qc ~ /qc_singleSubj
178+ sct_create_mask -i t2.nii.gz -p centerline,t2_seg.nii.gz -size 35mm -f cylinder -o mask_t2.nii.gz
179+ sct_crop_image -i t2.nii.gz -m mask_t2.nii.gz
180+
181+ # 2. T1w preprocessing (cropping around spinal cord)
182+ cd ../t1
183+ sct_deepseg_sc -i t1.nii.gz -c t1 -qc ~ /qc_singleSubj
184+ sct_create_mask -i t1.nii.gz -p centerline,t1_seg.nii.gz -size 35mm -f cylinder -o mask_t1.nii.gz
185+ sct_crop_image -i t1.nii.gz -m mask_t1.nii.gz
186+
187+ # 3. Perform registration
188+ # NB: `-dseg` is not necessary for registration, but is provided for the `-qc` reporting to help with spinal cord visualization
189+ sct_register_multimodal -i t1_crop.nii.gz -d ../t2/t2_crop.nii.gz -param step=1,type=im,algo=dl -qc ~ /qc_singleSubj -dseg ../t2/t2_seg.nii.gz
190+
191+
192+
172193# Gray/white matter: Segmentation
173194# ======================================================================================================================
174195
0 commit comments