Skip to content

Commit 8f7c0bf

Browse files
authored
ENH: Add expected steps for FreeSurfer 7 recon-all (#3389)
* ENH: Copy FS6 steps for FS7 * ENH: Update expected outputs for FS7 * STY: black
1 parent c0d450e commit 8f7c0bf

File tree

1 file changed

+120
-1
lines changed

1 file changed

+120
-1
lines changed

nipype/interfaces/freesurfer/preprocess.py

Lines changed: 120 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1081,6 +1081,7 @@ class ReconAll(CommandLine):
10811081
#
10821082
# [0] https://surfer.nmr.mgh.harvard.edu/fswiki/ReconAllTableStableV5.3
10831083
# [1] https://surfer.nmr.mgh.harvard.edu/fswiki/ReconAllTableStableV6.0
1084+
# [2] https://surfer.nmr.mgh.harvard.edu/fswiki/ReconAllTableStableV6.0#ReconAllTableStable7.1.1
10841085
_autorecon1_steps = [
10851086
("motioncor", ["mri/rawavg.mgz", "mri/orig.mgz"], []),
10861087
(
@@ -1210,7 +1211,7 @@ class ReconAll(CommandLine):
12101211
("wmparc", ["mri/wmparc.mgz", "stats/wmparc.stats"], []),
12111212
("balabels", ["label/BA.ctab", "label/BA.thresh.ctab"], []),
12121213
]
1213-
else:
1214+
elif Info.looseversion() < LooseVersion("7.0.0"):
12141215
_autorecon2_volonly_steps = [
12151216
("gcareg", ["mri/transforms/talairach.lta"], []),
12161217
("canorm", ["mri/norm.mgz"], []),
@@ -1322,6 +1323,124 @@ class ReconAll(CommandLine):
13221323
[],
13231324
),
13241325
]
1326+
else:
1327+
_autorecon2_volonly_steps = [
1328+
("gcareg", ["mri/transforms/talairach.lta"], []),
1329+
("canorm", ["mri/norm.mgz"], []),
1330+
("careg", ["mri/transforms/talairach.m3z"], []),
1331+
(
1332+
"calabel",
1333+
[
1334+
"mri/aseg.auto_noCCseg.mgz",
1335+
"mri/aseg.auto.mgz",
1336+
"mri/aseg.presurf.mgz",
1337+
],
1338+
[],
1339+
),
1340+
("normalization2", ["mri/brain.mgz"], []),
1341+
("maskbfs", ["mri/brain.finalsurfs.mgz"], []),
1342+
(
1343+
"segmentation",
1344+
["mri/wm.seg.mgz", "mri/wm.asegedit.mgz", "mri/wm.mgz"],
1345+
[],
1346+
),
1347+
(
1348+
"fill",
1349+
[
1350+
"mri/filled.mgz",
1351+
# 'scripts/ponscc.cut.log',
1352+
],
1353+
[],
1354+
),
1355+
]
1356+
_autorecon2_lh_steps = [
1357+
("tessellate", ["surf/lh.orig.nofix"], []),
1358+
("smooth1", ["surf/lh.smoothwm.nofix"], []),
1359+
("inflate1", ["surf/lh.inflated.nofix"], []),
1360+
("qsphere", ["surf/lh.qsphere.nofix"], []),
1361+
("fix", ["surf/lh.inflated", "surf/lh.orig"], []),
1362+
(
1363+
"white",
1364+
[
1365+
"surf/lh.white.preaparc",
1366+
"surf/lh.curv",
1367+
"surf/lh.area",
1368+
"label/lh.cortex.label",
1369+
],
1370+
[],
1371+
),
1372+
("smooth2", ["surf/lh.smoothwm"], []),
1373+
("inflate2", ["surf/lh.inflated", "surf/lh.sulc"], []),
1374+
(
1375+
"curvHK",
1376+
[
1377+
"surf/lh.white.H",
1378+
"surf/lh.white.K",
1379+
"surf/lh.inflated.H",
1380+
"surf/lh.inflated.K",
1381+
],
1382+
[],
1383+
),
1384+
("curvstats", ["stats/lh.curv.stats"], []),
1385+
]
1386+
_autorecon3_lh_steps = [
1387+
("sphere", ["surf/lh.sphere"], []),
1388+
("surfreg", ["surf/lh.sphere.reg"], []),
1389+
("jacobian_white", ["surf/lh.jacobian_white"], []),
1390+
("avgcurv", ["surf/lh.avg_curv"], []),
1391+
("cortparc", ["label/lh.aparc.annot"], []),
1392+
(
1393+
"pial",
1394+
[
1395+
"surf/lh.pial",
1396+
"surf/lh.curv.pial",
1397+
"surf/lh.area.pial",
1398+
"surf/lh.thickness",
1399+
"surf/lh.white",
1400+
],
1401+
[],
1402+
),
1403+
("parcstats", ["stats/lh.aparc.stats"], []),
1404+
("cortparc2", ["label/lh.aparc.a2009s.annot"], []),
1405+
("parcstats2", ["stats/lh.aparc.a2009s.stats"], []),
1406+
("cortparc3", ["label/lh.aparc.DKTatlas.annot"], []),
1407+
("parcstats3", ["stats/lh.aparc.DKTatlas.stats"], []),
1408+
("pctsurfcon", ["surf/lh.w-g.pct.mgh", "stats/lh.w-g.pct.stats"], []),
1409+
]
1410+
_autorecon3_added_steps = [
1411+
(
1412+
"cortribbon",
1413+
["mri/lh.ribbon.mgz", "mri/rh.ribbon.mgz", "mri/ribbon.mgz"],
1414+
[],
1415+
),
1416+
("hyporelabel", ["mri/aseg.presurf.hypos.mgz"], []),
1417+
(
1418+
"aparc2aseg",
1419+
[
1420+
"mri/aparc+aseg.mgz",
1421+
"mri/aparc.a2009s+aseg.mgz",
1422+
"mri/aparc.DKTatlas+aseg.mgz",
1423+
],
1424+
[],
1425+
),
1426+
("apas2aseg", ["mri/aseg.mgz"], ["mri/aparc+aseg.mgz"]),
1427+
("segstats", ["stats/aseg.stats"], []),
1428+
("wmparc", ["mri/wmparc.mgz", "stats/wmparc.stats"], []),
1429+
# Note that this is a very incomplete list; however the ctab
1430+
# files are last to be touched, so this should be reasonable
1431+
(
1432+
"balabels",
1433+
[
1434+
"label/BA_exvivo.ctab",
1435+
"label/BA_exvivo.thresh.ctab",
1436+
"label/lh.entorhinal_exvivo.label",
1437+
"label/rh.entorhinal_exvivo.label",
1438+
"label/lh.perirhinal_exvivo.label",
1439+
"label/rh.perirhinal_exvivo.label",
1440+
],
1441+
[],
1442+
),
1443+
]
13251444

13261445
# Fill out autorecon2 steps
13271446
_autorecon2_rh_steps = [

0 commit comments

Comments
 (0)