Differences between revisions 15 and 16
Deletions are marked like this. Additions are marked like this.
Line 48: Line 48:
 10. '''`recon-all -stage2 -subjid `''your_subject_name'''''  10. After you have made edits to your wm.mgz volume, run the next step of recon-all:
   
'''`recon-all -stage2 -subjid `''your_subject_name'''''
Line 50: Line 51:
 11. (Possibly make more edits to wm.mgz volume)

 12. '''`recon-all -stage3 -stage4a -stage4b ?!!!!-avgtif -gcs!!!!? -subjid `''your_subject_name'''''
 11. Check your inflated surfaces and wm.mgz volumes to make sure they are accurate.
  11.1 If more edits are needed, run the following command after making edits:
    '''`recon-all -stage2 -subjid `''your_subject_name'''''
  11.2 If no more edits are needed, continue on to the next step of recon-all:
    '''`recon-all -stage3 -stage4a -stage4b ?!!!!-avgtif -gcs!!!!? -subjid `''your_subject_name'''''
Line 65: Line 68:

== Recon-all steps and stages ==

||<rowbgcolor='#80FF80'>'''recon-all step'''||'''Individual Flag'''||'''Input'''||'''Command Line'''||'''Output'''||
||<bgcolor='#FFFFE0'>none||<bgcolor='#FFFFE0'>none||<bgcolor='#E0E0FF'>file.dcm||<bgcolor='#E0E0FF'>["mri_convert"] file.dcm orig/001.mgz||<bgcolor='#E0E0FF'>orig/001.mgz||
||<|17(bgcolor='#FFFFE0'>["recon-all"] -stage1 -subjid subj||
||<|3(bgcolor='#FFFFE0'>-motioncor||<bgcolor='#E0E0FF'>orig/001.mgz||<|2(bgcolor='#E0E0FF'>["mri_motion_correct2"] -i orig/001.mgz -i orig/002.mgz -o rawavg.mgz||<|2(bgcolor='#E0E0FF'>rawavg.mgz||
||<rowbgcolor='#E0E0FF'>orig/002.mgz||
||<rowbgcolor='#E0E0FF'>rawavg.mgz||["mri_convert"] rawavg.mgz orig.mgz --conform||orig.mgz||
||<|3(bgcolor='#FFFFE0'>-nuintensitycor||<bgcolor='#E0E0FF'>orig.mgz||<bgcolor='#E0E0FF'>["mri_convert"] orig.mgz orig.mnc||<bgcolor='#E0E0FF'>orig.mnc||
||<rowbgcolor='#E0E0FF'>orig.mnc||(4 iterations of) ["nu_correct"] -clobber nu0.mnc nu1.mnc||nu4.mnc||
||<rowbgcolor='#E0E0FF'>nu4.mnc||["mri_convert"] nu4.mnc nu.mgz||nu.mgz||
||<bgcolor='#FFFFE0'>-talairach||<bgcolor='#E0E0FF'>nu.mgz||<bgcolor='#E0E0FF'>["talairach2"] subjid -mgz||<bgcolor='#E0E0FF'>transforms/talairach.xfm||
||<bgcolor='#FFFFE0'>-normalization||<bgcolor='#E0E0FF'>nu.mgz||<bgcolor='#E0E0FF'>["mri_normalize"] nu.mgz T1.mgz||<bgcolor='#E0E0FF'>T1.mgz||
||<bgcolor='#FFFFE0'>-skullstrip||<bgcolor='#E0E0FF'>T1.mgz||<bgcolor='#E0E0FF'>["mri_watershed"] T1.mgz brain.mgz||<bgcolor='#E0E0FF'>brain.mgz||
||<bgcolor='#FFFFE0'>-segmentation||<bgcolor='#E0E0FF'>brain.mgz||<bgcolor='#E0E0FF'>["mri_segment"] brain.mgz wm.mgz||<bgcolor='#E0E0FF'>wm.mgz||
||<|2(bgcolor='#FFFFE0'>-fill||<bgcolor='#E0E0FF'>wm.mgz||<|2(bgcolor='#E0E0FF'>["mri_fill"] -a ../scripts/ponscc.cut.log -xform transforms/talairach.xfm -segmentation aseg.mgz wm.mgz filled.mgz||<bgcolor='#E0E0FF'>filled.mgz||
||<rowbgcolor='#E0E0FF'>aseg.mgz||../scripts/ponscc.cut.log||
||<|2((bgcolor='#FFFFE0'>-tessellate||<|2(bgcolor='#E0E0FF'>filled.mgz||<bgcolor='#E0E0FF'>["mri_tessellate"] filled.mgz 255 ../surf/lh.orig||<bgcolor='#E0E0FF'>lh.orig||
||<rowbgcolor='#E0E0FF'>["mri_tessellate"] filled.mgz 128 ../surf/rh.orig||rh.orig||
||<bgcolor='#FFFFE0'>-smooth1||<bgcolor='#E0E0FF'>?h.orig||<bgcolor='#E0E0FF'>["mris_smooth"] ?h.orig ?h.smoothwm||<bgcolor='#E0E0FF'>?h.smoothwm||
||<bgcolor='#FFFFE0'>-inflate1||<bgcolor='#E0E0FF'>?h.smoothwm||<bgcolor='#E0E0FF'>["mris_inflate"] ?h.smoothwm ?h.inflated||<bgcolor='#E0E0FF'>?h.inflated||
||<-5(rowbgcolor='#FF8080'>Check skullstrip (brain.mgz), talairach (transforms/talairach.xfm), and normalization (brain.mgz or T1.mgz - mean wm voxel value = 110). Perform manual edits to wm.mgz volume.||
||<|7(bgcolor='#FFFFE0'>["recon-all"] –stage2 -subjid subj||
||<|2(bgcolor='#FFFFE0'>-fill||<bgcolor='#E0E0FF'>wm.mgz||<|2(bgcolor='#E0E0FF'>["mri_fill"] -a ../scripts/ponscc.cut.log -xform transforms/talairach.xfm -segmentation aseg.mgz wm.mgz filled.mgz||<bgcolor='#E0E0FF'>filled.mgz||
||<rowbgcolor='#E0E0FF'>aseg.mgz||../scripts/ponscc.cut.log||
||<|2((bgcolor='#FFFFE0'>-tessellate||<|2(bgcolor='#E0E0FF'>filled.mgz||<bgcolor='#E0E0FF'>["mri_tessellate"] filled.mgz 255 ../surf/lh.orig||<bgcolor='#E0E0FF'>lh.orig||
||<rowbgcolor='#E0E0FF'>["mri_tessellate"] filled.mgz 128 ../surf/rh.orig||rh.orig||
||<bgcolor='#FFFFE0'>-smooth1||<bgcolor='#E0E0FF'>?h.orig||<bgcolor='#E0E0FF'>["mris_smooth"] ?h.orig ?h.smoothwm||<bgcolor='#E0E0FF'>?h.smoothwm||
||<bgcolor='#FFFFE0'>-inflate1||<bgcolor='#E0E0FF'>?h.smoothwm||<bgcolor='#E0E0FF'>["mris_inflate"] ?h.smoothwm ?h.inflated||<bgcolor='#E0E0FF'>?h.inflated||
||<-5(rowbgcolor='#FF8080'>Check inflated surface and wm volume, make more edits if necessary||
||<|12(bgcolor='#FFFFE0'>["recon-all"] –stage3 -subjid subj||
||<|2(bgcolor='#FFFFE0'>-fill||<bgcolor='#E0E0FF'>wm.mgz||<|2(bgcolor='#E0E0FF'>["mri_fill"] -a ../scripts/ponscc.cut.log -xform transforms/talairach.xfm -segmentation aseg.mgz wm.mgz filled.mgz||<bgcolor='#E0E0FF'>filled.mgz||
||<rowbgcolor='#E0E0FF'>aseg.mgz||../scripts/ponscc.cut.log||
||<|2((bgcolor='#FFFFE0'>-tessellate||<|2(bgcolor='#E0E0FF'>filled.mgz||<bgcolor='#E0E0FF'>["mri_tessellate"] filled.mgz 255 ../surf/lh.orig||<bgcolor='#E0E0FF'>lh.orig||
||<rowbgcolor='#E0E0FF'>["mri_tessellate"] filled.mgz 128 ../surf/rh.orig||rh.orig||
||<bgcolor='#FFFFE0'>-smooth1||<bgcolor='#E0E0FF'>?h.orig||<bgcolor='#E0E0FF'>["mris_smooth"] ?h.orig ?h.smoothwm||<bgcolor='#E0E0FF'>?h.smoothwm||
||<bgcolor='#FFFFE0'>-inflate1||<bgcolor='#E0E0FF'>?h.smoothwm||<bgcolor='#E0E0FF'>["mris_inflate"] ?h.smoothwm ?h.inflated||<bgcolor='#E0E0FF'>?h.inflated||
||<bgcolor='#FFFFE0'>-qsphere||<bgcolor='#E0E0FF'>?h.inflated||<bgcolor='#E0E0FF'>["mris_sphere"] -w 0 -inflate -in 200 -q ?h.inflated ?h.qsphere||<bgcolor='#E0E0FF'>?h.qsphere||
||<bgcolor='#FFFFE0'>-fix||<bgcolor='#E0E0FF'>?h.qsphere||<bgcolor='#E0E0FF'>["mris_fix_topology"] subj ?h||<bgcolor='#E0E0FF'>?h.orig||
||<bgcolor='#FFFFE0'>-euler||<bgcolor='#E0E0FF'>?h.orig||<bgcolor='#E0E0FF'>["mris_euler_number"] ?h.orig||<bgcolor='#E0E0FF'>?h.orig.euler||
||<bgcolor='#FFFFE0'>-smooth2||<bgcolor='#E0E0FF'>?h.orig||<bgcolor='#E0E0FF'>["mris_smooth"] ?h.orig ?h.smoothwm||<bgcolor='#E0E0FF'>?h.smoothwm||
||<bgcolor='#FFFFE0'>-inflate2||<bgcolor='#E0E0FF'>?h.smoothwm||<bgcolor='#E0E0FF'>["mris_inflate"] ?h.smoothwm ?h.inflated||<bgcolor='#E0E0FF'>?h.inflated||
||<|7(bgcolor='#FFFFE0'>["recon-all"] –stage4a -subjid subj||
||<|3(bgcolor='#FFFFE0'>-finalsurfs||<bgcolor='#E0E0FF'>brain.mgz||<|3(bgcolor='#E0E0FF'>["mris_make_surfaces"] -w 0 subj ?h||<bgcolor='#E0E0FF'>?h.white||
||<|2(rowbgcolor='#E0E0FF'>?h.orig||?h.pial||
||<rowbgcolor='#E0E0FF'>?h.thickness||
||<|3(bgcolor='#FFFFE0'>-cortribbon||<bgcolor='#E0E0FF'>orig.mgz||<|3(bgcolor='#E0E0FF'>["mri_surf2vol"] --mkmask --hemi ?h --fillribbon --template orig.mgz --volregidentity subj --outvol ?h.ribbon.mgz||<|3(bgcolor='#E0E0FF'>?h.ribbon.mgz||
||<rowbgcolor='#E0E0FF'>?h.white||
||<rowbgcolor='#E0E0FF'>?h.pial||
||<|7(bgcolor='#FFFFE0'>["recon-all"] –stage4b -subjid subj||
||<bgcolor='#FFFFE0'>-sphere||<bgcolor='#E0E0FF'>?h.inflated||<bgcolor='#E0E0FF'>["mris_sphere"] -w 0 ?h.inflated ?h.sphere||<bgcolor='#E0E0FF'>?h.sphere||
||<bgcolor='#FFFFE0'>-surfreg||<bgcolor='#E0E0FF'>?h.sphere||<bgcolor='#E0E0FF'>["mris_register"] -w 0 -curv ?h.sphere $AvgTif ?h.sphere.reg||<bgcolor='#E0E0FF'>?h.sphere.reg||
||<|2(bgcolor='#FFFFE0'>-contrasurfreg||<bgcolor='#E0E0FF'>rh.sphere||<bgcolor='#E0E0FF'>["mris_register"] -w 0 -curv -reverse rh.sphere $AvgTif rh.lh.sphere.reg||<bgcolor='#E0E0FF'>rh.lh.sphere.reg||
||<rowbgcolor='#E0E0FF'>lh.sphere||["mris_register"] -w 0 -curv -reverse lh.sphere $AvgTif lh.rh.sphere.reg||lh.rh.sphere.reg||
||<bgcolor='#FFFFE0'>-avgcurv||<bgcolor='#E0E0FF'>?h.sphere.reg||<bgcolor='#E0E0FF'>["mrisp_paint"] -a 5 $AvgTif#6 ?h.sphere.reg ?h.avg_curv||<bgcolor='#E0E0FF'>?h.avg_curv||

[wiki:FreeSurferWorkFlows top] | [wiki:TroubleshootingReconstruction previous] | [wiki:SubcorticalSegmentation next]

This workflow uses the older version of FreeSurfer, in which manual editting of the wm volume was required and subcortical segmentations were not a part of the regular processing stream. See HistoricalReconTable for the steps run by this script. Continue using this workflow if you've already begun processing your subjects this way.

FreeSurfer Historical Reconstruction Work Flow

  1. Source the correct version of FreeSurfer:

    • source your_freesurfer_dir/SetUpFreeSurfer.csh

  2. Set your SUBJECTS_DIR variable to your subjects directory:
    • setenv SUBJECTS_DIR your_subjects_dir

  3. Make a subjects directory within your SUBJECTS_DIR:
    • mksubjdirs your_subject_name

  4. Convert the DICOM files to .mgz files in the $SUBJECTS_DIR/your_subjects_name/mri/orig directory:
    • mri_convert your_dicom_file your_subject_name/mri/orig/001.mgz

If you have multiple scans from the same session, convert additional scans like:

  • mri_convert your_other_dicom_file your_subject_name/mri/orig/002.mgz

  1. If you have a single scan, conform the volume to 1 mm^3 and 256x256x256:
    • mri_convert --conform your_subject_name/mri/orig/002.mgz your_subject_name/mri/orig.mgz

  2. Start the recon-all script:
    • recon-all -stage1 -subjid your_subject_name

  3. Check the talairach transform:
    • 7.1 Make sure the talairach subject is linked to $SUBJECTS_DIR

      • if not: ln -s $FREESURFER_HOME/subjects/talairach $SUBJECTS_DIR/.

      7.2 tkregister2 --mgz --s your_subject_name --fstal --surf orig For instructions on how to fix bad talairachs tranforms, refer to:

      • ["FsTutorial/Talairach"]
  4. Check the skull strip:

    tkmedit your_subject_name brain.mgz -aux T1.mgz For intructions on how to fix poor normalizations or skull strips, refer to:

  5. Edit the wm.mgz volume For instructions on editing the wm.mgz, refer to:
  6. After you have made edits to your wm.mgz volume, run the next step of recon-all:
    • recon-all -stage2 -subjid your_subject_name

  7. Check your inflated surfaces and wm.mgz volumes to make sure they are accurate.
    • 11.1 If more edits are needed, run the following command after making edits:
      • recon-all -stage2 -subjid your_subject_name

      11.2 If no more edits are needed, continue on to the next step of recon-all:
      • recon-all -stage3 -stage4a -stage4b ?!!!!-avgtif -gcs!!!!?  -subjid your_subject_name

  8. Check the white and pial surfaces: Note: these Tutorial sections show data processed using Basic workflow, so large defects were autofilled.

HistoricalReconstruction (last edited 2021-09-22 11:45:30 by DevaniCordero)