@@ -659,82 +659,86 @@ class ReconAll(CommandLine):
659
659
_steps = [
660
660
# autorecon1
661
661
('motioncor' , ['mri/rawavg.mgz' , 'mri/orig.mgz' ]),
662
- ('talairach' , ['mri/transforms/talairach.auto.xfm' ,
663
- 'mri/transforms/talairach.xfm' ]),
662
+ ('talairach' , ['mri/orig_nu.mgz' ,
663
+ 'mri/transforms/talairach.auto.xfm' ,
664
+ 'mri/transforms/talairach.xfm' ,
665
+ # 'mri/transforms/talairach_avi.log',
666
+ ]),
664
667
('nuintensitycor' , ['mri/nu.mgz' ]),
665
668
('normalization' , ['mri/T1.mgz' ]),
666
- ('skullstrip' ,
667
- [ 'mri/brainmask.auto.mgz' ,
668
- 'mri/brainmask.mgz' ]),
669
+ ('skullstrip' , [ 'mri/talairach_with_skull.lta' ,
670
+ 'mri/brainmask.auto.mgz' ,
671
+ 'mri/brainmask.mgz' ]),
669
672
# autorecon2
670
673
('gcareg' , ['mri/transforms/talairach.lta' ]),
671
674
('canorm' , ['mri/norm.mgz' ]),
672
675
('careg' , ['mri/transforms/talairach.m3z' ]),
673
- ('careginv' , ['mri/transforms/talairach.m3z.inv.x.mgz' ,
674
- 'mri/transforms/talairach.m3z.inv.y.mgz' ,
675
- 'mri/transforms/talairach.m3z.inv.z.mgz' ]),
676
- ('rmneck' , ['mri/nu_noneck.mgz' ]),
677
- ('skull-lta' , ['mri/transforms/talairach_with_skull_2.lta' ]),
678
- ('calabel' ,
679
- ['mri/aseg.auto_noCCseg.mgz' , 'mri/aseg.auto.mgz' , 'mri/aseg.mgz' ]),
676
+ ('calabel' , ['mri/aseg.auto_noCCseg.mgz' ,
677
+ 'mri/aseg.auto.mgz' ,
678
+ 'mri/aseg.mgz' ]),
680
679
('normalization2' , ['mri/brain.mgz' ]),
681
680
('maskbfs' , ['mri/brain.finalsurfs.mgz' ]),
682
- ('segmentation' , ['mri/wm.asegedit.mgz' , 'mri/wm.mgz' ]),
683
- ('fill' , ['mri/filled.mgz' ]),
681
+ ('segmentation' , ['mri/wm.seg.mgz' ,
682
+ 'mri/wm.asegedit.mgz' ,
683
+ 'mri/wm.mgz' ]),
684
+ ('fill' , ['mri/filled.mgz' ,
685
+ # 'scripts/ponscc.cut.log',
686
+ ]),
684
687
('tessellate' , ['surf/lh.orig.nofix' , 'surf/rh.orig.nofix' ]),
685
688
('smooth1' , ['surf/lh.smoothwm.nofix' , 'surf/rh.smoothwm.nofix' ]),
686
689
('inflate1' , ['surf/lh.inflated.nofix' , 'surf/rh.inflated.nofix' ]),
687
690
('qsphere' , ['surf/lh.qsphere.nofix' , 'surf/rh.qsphere.nofix' ]),
688
691
('fix' , ['surf/lh.orig' , 'surf/rh.orig' ]),
689
- ('white' ,
690
- ['surf/lh.white' ,
691
- 'surf/rh.white' ,
692
- 'surf/lh.curv' ,
693
- 'surf/rh.curv' ,
694
- 'surf/lh.area' ,
695
- 'surf/rh.area' ,
696
- 'label/lh.cortex.label' ,
697
- 'label/rh.cortex.label' ]),
692
+ ('white' , ['surf/lh.white.preaparc' , 'surf/rh.white.preaparc' ,
693
+ 'surf/lh.curv' , 'surf/rh.curv' ,
694
+ 'surf/lh.area' , 'surf/rh.area' ,
695
+ 'label/lh.cortex.label' , 'label/rh.cortex.label' ]),
698
696
('smooth2' , ['surf/lh.smoothwm' , 'surf/rh.smoothwm' ]),
699
- ('inflate2' ,
700
- ['surf/lh.inflated' ,
701
- 'surf/rh.inflated' ,
702
- 'surf/lh.sulc' ,
703
- 'surf/rh.sulc' ,
704
- 'surf/lh.inflated.H' ,
705
- 'surf/rh.inflated.H' ,
706
- 'surf/lh.inflated.K' ,
707
- 'surf/rh.inflated.K' ]),
697
+ ('inflate2' , ['surf/lh.inflated' , 'surf/rh.inflated' ,
698
+ 'surf/lh.sulc' , 'surf/rh.sulc' ]),
699
+ ('curvHK' , ['surf/lh.white.H' , 'surf/rh.white.H' ,
700
+ 'surf/lh.white.K' , 'surf/rh.white.K' ,
701
+ 'surf/lh.inflated.H' , 'surf/rh.inflated.H' ,
702
+ 'surf/lh.inflated.K' , 'surf/rh.inflated.K' ]),
703
+ ('curvstats' , ['stats/lh.curv.stats' , 'stats/rh.curv.stats' ]),
708
704
# autorecon3
709
705
('sphere' , ['surf/lh.sphere' , 'surf/rh.sphere' ]),
710
706
('surfreg' , ['surf/lh.sphere.reg' , 'surf/rh.sphere.reg' ]),
711
707
('jacobian_white' , ['surf/lh.jacobian_white' ,
712
708
'surf/rh.jacobian_white' ]),
713
709
('avgcurv' , ['surf/lh.avg_curv' , 'surf/rh.avg_curv' ]),
714
710
('cortparc' , ['label/lh.aparc.annot' , 'label/rh.aparc.annot' ]),
715
- ('pial' ,
716
- ['surf/lh.pial' ,
717
- 'surf/rh.pial' ,
718
- 'surf/lh.curv.pial' ,
719
- 'surf/rh.curv.pial' ,
720
- 'surf/lh.area.pial' ,
721
- 'surf/rh.area.pial' ,
722
- 'surf/lh.thickness' ,
723
- 'surf/rh.thickness' ]),
724
- ('cortparc2' , ['label/lh.aparc.a2009s.annot' ,
725
- 'label/rh.aparc.a2009s.annot' ]),
726
- ('parcstats2' ,
727
- ['stats/lh.aparc.a2009s.stats' ,
728
- 'stats/rh.aparc.a2009s.stats' ,
729
- 'stats/aparc.annot.a2009s.ctab' ]),
711
+ ('pial' , ['surf/lh.pial' , 'surf/rh.pial' ,
712
+ 'surf/lh.curv.pial' , 'surf/rh.curv.pial' ,
713
+ 'surf/lh.area.pial' , 'surf/rh.area.pial' ,
714
+ 'surf/lh.thickness' , 'surf/rh.thickness' ]),
715
+ # TODO: Optional -T2pial / -FLAIRpial
730
716
('cortribbon' , ['mri/lh.ribbon.mgz' , 'mri/rh.ribbon.mgz' ,
731
717
'mri/ribbon.mgz' ]),
718
+ ('parcstats' , ['stats/lh.aparc.astats' , 'stats/rh.aparc.stats' ,
719
+ 'stats/aparc.annot.ctab' ]),
720
+ ('cortparc2' , ['label/lh.aparc.a2009s.annot' ,
721
+ 'label/rh.aparc.a2009s.annot' ]),
722
+ ('parcstats2' , ['stats/lh.aparc.a2009s.stats' ,
723
+ 'stats/rh.aparc.a2009s.stats' ,
724
+ 'stats/aparc.annot.a2009s.ctab' ]),
725
+ ('cortparc3' , ['label/lh.aparc.DKTatlas.annot' ,
726
+ 'label/rh.aparc.DKTatlas.annot' ]),
727
+ ('parcstats3' , ['stats/lh.aparc.DKTatlas.stats' ,
728
+ 'stats/rh.aparc.DKTatlas.stats' ,
729
+ 'stats/aparc.annot.DKTatlas.ctab' ]),
730
+ ('pctsurfcon' , ['surf/lh.w-g.pct.mgh' , 'surf/rh.w-g.pct.mgh' ]),
731
+ ('hyporelabel' , ['mri/aseg.presurf.hypos.mgz' ]),
732
+ ('aparc2aseg' , ['mri/aparc+aseg.mgz' ,
733
+ 'mri/aparc.a2009s+aseg.mgz' ,
734
+ 'mri/aparc.DKTatlas+aseg.mgz' ]),
735
+ ('apas2aseg' , ['mri/aseg.mgz' ]), # XXX: Will not run because of calabel
732
736
('segstats' , ['stats/aseg.stats' ]),
733
- ('aparc2aseg' , ['mri/aparc+aseg.mgz' , 'mri/aparc.a2009s+aseg.mgz' ]),
734
737
('wmparc' , ['mri/wmparc.mgz' , 'stats/wmparc.stats' ]),
735
- ('balabels' , ['BA.ctab' , 'BA.thresh.ctab' ]),
736
- ('label-exvivo-ec' , ['label/lh.entorhinal_exvivo.label' ,
737
- 'label/rh.entorhinal_exvivo.label' ])]
738
+ ('balabels' , ['BA.ctab' , 'BA.thresh.ctab' ,
739
+ 'label/lh.entorhinal_exvivo.label' ,
740
+ 'label/rh.entorhinal_exvivo.label' ]),
741
+ ]
738
742
739
743
def _gen_subjects_dir (self ):
740
744
return os .getcwd ()
0 commit comments