#@# Fix Topology Copy lh Mon May 31 11:15:08 CEST 2021 /mnt/newdata/projects/avoss/HCP/subjects/avoss_01/T1w/avoss_01/scripts cp ../surf/lh.orig.nofix ../surf/lh.orig cp ../surf/lh.inflated.nofix ../surf/lh.inflated #-------------------------------------------- #@# Fix Topology Copy rh Mon May 31 11:15:08 CEST 2021 /mnt/newdata/projects/avoss/HCP/subjects/avoss_01/T1w/avoss_01/scripts cp ../surf/rh.orig.nofix ../surf/rh.orig cp ../surf/rh.inflated.nofix ../surf/rh.inflated #@# Fix Topology lh Mon May 31 11:15:08 CEST 2021 mris_fix_topology -rusage /mnt/newdata/projects/avoss/HCP/subjects/avoss_01/T1w/avoss_01/touch/rusage.mris_fix_topology.lh.dat -mgz -sphere qsphere.nofix -ga -seed 1234 avoss_01 lh reading spherical homeomorphism from 'qsphere.nofix' using genetic algorithm with optimized parameters setting seed for random number genererator to 1234 ************************************************************* Topology Correction Parameters retessellation mode: genetic search number of patches/generation : 10 number of generations : 10 surface mri loglikelihood coefficient : 1.0 volume mri loglikelihood coefficient : 10.0 normal dot loglikelihood coefficient : 1.0 quadratic curvature loglikelihood coefficient : 1.0 volume resolution : 2 eliminate vertices during search : 1 initial patch selection : 1 select all defect vertices : 0 ordering dependant retessellation: 0 use precomputed edge table : 0 smooth retessellated patch : 2 match retessellated patch : 1 verbose mode : 0 ************************************************************* INFO: assuming .mgz format $Id: mris_fix_topology.c,v 1.50.2.1 2016/10/27 22:25:58 zkaufman Exp $ $Id: mrisurf.c,v 1.781.2.6 2016/12/27 16:47:14 zkaufman Exp $ before topology correction, eno=-236 (nv=170046, nf=340564, ne=510846, g=119) using quasi-homeomorphic spherical map to tessellate cortical surface... Correction of the Topology Finding true center and radius of Spherical Surface...done Surface centered at (0,0,0) with radius 100.0 in 10 iterations marking ambiguous vertices... 31240 ambiguous faces found in tessellation segmenting defects... 66 defects found, arbitrating ambiguous regions... analyzing neighboring defects... -merging segment 18 into 11 -merging segment 21 into 15 -merging segment 45 into 44 63 defects to be corrected 0 vertices coincident reading input surface /mnt/newdata/projects/avoss/HCP/subjects/avoss_01/T1w/avoss_01/surf/lh.qsphere.nofix... reading brain volume from brain... reading wm segmentation from wm... Computing Initial Surface Statistics -face loglikelihood: -9.1568 (-4.5784) -vertex loglikelihood: -6.1795 (-3.0897) -normal dot loglikelihood: -3.5667 (-3.5667) -quad curv loglikelihood: -5.9168 (-2.9584) Total Loglikelihood : -24.8197 CORRECTING DEFECT 0 (vertices=65, convex hull=66, v0=6) After retessellation of defect 0 (v0=6), euler #=-47 (152704,454667,301916) : difference with theory (-60) = -13 CORRECTING DEFECT 1 (vertices=186, convex hull=131, v0=90) After retessellation of defect 1 (v0=90), euler #=-48 (152723,454805,302034) : difference with theory (-59) = -11 CORRECTING DEFECT 2 (vertices=132, convex hull=69, v0=531) After retessellation of defect 2 (v0=531), euler #=-47 (152734,454864,302083) : difference with theory (-58) = -11 CORRECTING DEFECT 3 (vertices=12, convex hull=30, v0=1522) After retessellation of defect 3 (v0=1522), euler #=-46 (152735,454878,302097) : difference with theory (-57) = -11 CORRECTING DEFECT 4 (vertices=25, convex hull=61, v0=1811) After retessellation of defect 4 (v0=1811), euler #=-45 (152743,454924,302136) : difference with theory (-56) = -11 CORRECTING DEFECT 5 (vertices=85, convex hull=48, v0=1827) After retessellation of defect 5 (v0=1827), euler #=-44 (152751,454965,302170) : difference with theory (-55) = -11 CORRECTING DEFECT 6 (vertices=29, convex hull=65, v0=2014) After retessellation of defect 6 (v0=2014), euler #=-43 (152770,455047,302234) : difference with theory (-54) = -11 CORRECTING DEFECT 7 (vertices=254, convex hull=93, v0=2636) After retessellation of defect 7 (v0=2636), euler #=-42 (152799,455172,302331) : difference with theory (-53) = -11 CORRECTING DEFECT 8 (vertices=96, convex hull=48, v0=4296) After retessellation of defect 8 (v0=4296), euler #=-41 (152807,455211,302363) : difference with theory (-52) = -11 CORRECTING DEFECT 9 (vertices=51, convex hull=44, v0=4384) After retessellation of defect 9 (v0=4384), euler #=-40 (152812,455242,302390) : difference with theory (-51) = -11 CORRECTING DEFECT 10 (vertices=11, convex hull=35, v0=6265) After retessellation of defect 10 (v0=6265), euler #=-39 (152815,455263,302409) : difference with theory (-50) = -11 CORRECTING DEFECT 11 (vertices=9542, convex hull=2193, v0=6475) XL defect detected... After retessellation of defect 11 (v0=6475), euler #=-48 (153382,458232,304802) : difference with theory (-49) = -1 CORRECTING DEFECT 12 (vertices=10, convex hull=31, v0=7190) After retessellation of defect 12 (v0=7190), euler #=-47 (153384,458248,304817) : difference with theory (-48) = -1 CORRECTING DEFECT 13 (vertices=474, convex hull=226, v0=9704) After retessellation of defect 13 (v0=9704), euler #=-46 (153479,458646,305121) : difference with theory (-47) = -1 CORRECTING DEFECT 14 (vertices=54, convex hull=45, v0=13781) After retessellation of defect 14 (v0=13781), euler #=-45 (153488,458686,305153) : difference with theory (-46) = -1 CORRECTING DEFECT 15 (vertices=408, convex hull=184, v0=15433) After retessellation of defect 15 (v0=15433), euler #=-43 (153500,458813,305270) : difference with theory (-45) = -2 CORRECTING DEFECT 16 (vertices=49, convex hull=28, v0=16848) After retessellation of defect 16 (v0=16848), euler #=-42 (153502,458827,305283) : difference with theory (-44) = -2 CORRECTING DEFECT 17 (vertices=48, convex hull=69, v0=19458) After retessellation of defect 17 (v0=19458), euler #=-41 (153526,458928,305361) : difference with theory (-43) = -2 CORRECTING DEFECT 18 (vertices=126, convex hull=93, v0=20010) After retessellation of defect 18 (v0=20010), euler #=-40 (153550,459043,305453) : difference with theory (-42) = -2 CORRECTING DEFECT 19 (vertices=5, convex hull=19, v0=21093) After retessellation of defect 19 (v0=21093), euler #=-39 (153550,459048,305459) : difference with theory (-41) = -2 CORRECTING DEFECT 20 (vertices=53, convex hull=46, v0=24663) After retessellation of defect 20 (v0=24663), euler #=-38 (153555,459080,305487) : difference with theory (-40) = -2 CORRECTING DEFECT 21 (vertices=756, convex hull=137, v0=32841) After retessellation of defect 21 (v0=32841), euler #=-37 (153571,459194,305586) : difference with theory (-39) = -2 CORRECTING DEFECT 22 (vertices=320, convex hull=59, v0=50960) After retessellation of defect 22 (v0=50960), euler #=-36 (153583,459251,305632) : difference with theory (-38) = -2 CORRECTING DEFECT 23 (vertices=39, convex hull=36, v0=53620) After retessellation of defect 23 (v0=53620), euler #=-35 (153593,459294,305666) : difference with theory (-37) = -2 CORRECTING DEFECT 24 (vertices=38, convex hull=34, v0=54552) After retessellation of defect 24 (v0=54552), euler #=-34 (153600,459326,305692) : difference with theory (-36) = -2 CORRECTING DEFECT 25 (vertices=658, convex hull=222, v0=58122) After retessellation of defect 25 (v0=58122), euler #=-33 (153617,459495,305845) : difference with theory (-35) = -2 CORRECTING DEFECT 26 (vertices=296, convex hull=268, v0=60820) After retessellation of defect 26 (v0=60820), euler #=-33 (153766,460078,306279) : difference with theory (-34) = -1 CORRECTING DEFECT 27 (vertices=26, convex hull=55, v0=61654) After retessellation of defect 27 (v0=61654), euler #=-32 (153778,460137,306327) : difference with theory (-33) = -1 CORRECTING DEFECT 28 (vertices=102, convex hull=46, v0=65747) After retessellation of defect 28 (v0=65747), euler #=-31 (153797,460213,306385) : difference with theory (-32) = -1 CORRECTING DEFECT 29 (vertices=11, convex hull=43, v0=67238) After retessellation of defect 29 (v0=67238), euler #=-30 (153800,460238,306408) : difference with theory (-31) = -1 CORRECTING DEFECT 30 (vertices=38, convex hull=65, v0=68474) After retessellation of defect 30 (v0=68474), euler #=-29 (153821,460328,306478) : difference with theory (-30) = -1 CORRECTING DEFECT 31 (vertices=6, convex hull=32, v0=70278) After retessellation of defect 31 (v0=70278), euler #=-28 (153823,460345,306494) : difference with theory (-29) = -1 CORRECTING DEFECT 32 (vertices=24, convex hull=44, v0=72579) After retessellation of defect 32 (v0=72579), euler #=-27 (153837,460404,306540) : difference with theory (-28) = -1 CORRECTING DEFECT 33 (vertices=14, convex hull=26, v0=74424) After retessellation of defect 33 (v0=74424), euler #=-26 (153838,460417,306553) : difference with theory (-27) = -1 CORRECTING DEFECT 34 (vertices=45, convex hull=35, v0=79490) After retessellation of defect 34 (v0=79490), euler #=-25 (153847,460455,306583) : difference with theory (-26) = -1 CORRECTING DEFECT 35 (vertices=318, convex hull=261, v0=80532) After retessellation of defect 35 (v0=80532), euler #=-25 (153992,461019,307002) : difference with theory (-25) = 0 CORRECTING DEFECT 36 (vertices=11, convex hull=33, v0=82016) After retessellation of defect 36 (v0=82016), euler #=-24 (153996,461046,307026) : difference with theory (-24) = 0 CORRECTING DEFECT 37 (vertices=17, convex hull=30, v0=84129) After retessellation of defect 37 (v0=84129), euler #=-23 (153999,461064,307042) : difference with theory (-23) = 0 CORRECTING DEFECT 38 (vertices=71, convex hull=43, v0=92829) After retessellation of defect 38 (v0=92829), euler #=-22 (154012,461117,307083) : difference with theory (-22) = 0 CORRECTING DEFECT 39 (vertices=42, convex hull=62, v0=93771) After retessellation of defect 39 (v0=93771), euler #=-21 (154033,461210,307156) : difference with theory (-21) = 0 CORRECTING DEFECT 40 (vertices=547, convex hull=396, v0=96735) After retessellation of defect 40 (v0=96735), euler #=-21 (154088,461576,307467) : difference with theory (-20) = 1 CORRECTING DEFECT 41 (vertices=21, convex hull=21, v0=103368) After retessellation of defect 41 (v0=103368), euler #=-20 (154091,461591,307480) : difference with theory (-19) = 1 CORRECTING DEFECT 42 (vertices=233, convex hull=131, v0=107376) After retessellation of defect 42 (v0=107376), euler #=-18 (154159,461859,307682) : difference with theory (-18) = 0 CORRECTING DEFECT 43 (vertices=689, convex hull=316, v0=112629) After retessellation of defect 43 (v0=112629), euler #=-17 (154293,462426,308116) : difference with theory (-17) = 0 CORRECTING DEFECT 44 (vertices=36, convex hull=24, v0=113615) After retessellation of defect 44 (v0=113615), euler #=-16 (154296,462443,308131) : difference with theory (-16) = 0 CORRECTING DEFECT 45 (vertices=94, convex hull=47, v0=116583) After retessellation of defect 45 (v0=116583), euler #=-15 (154305,462488,308168) : difference with theory (-15) = 0 CORRECTING DEFECT 46 (vertices=185, convex hull=114, v0=117540) After retessellation of defect 46 (v0=117540), euler #=-14 (154317,462577,308246) : difference with theory (-14) = 0 CORRECTING DEFECT 47 (vertices=48, convex hull=44, v0=121337) After retessellation of defect 47 (v0=121337), euler #=-13 (154327,462628,308288) : difference with theory (-13) = 0 CORRECTING DEFECT 48 (vertices=22, convex hull=25, v0=121501) After retessellation of defect 48 (v0=121501), euler #=-12 (154334,462656,308310) : difference with theory (-12) = 0 CORRECTING DEFECT 49 (vertices=28, convex hull=60, v0=124924) After retessellation of defect 49 (v0=124924), euler #=-11 (154343,462712,308358) : difference with theory (-11) = 0 CORRECTING DEFECT 50 (vertices=28, convex hull=33, v0=127232) After retessellation of defect 50 (v0=127232), euler #=-10 (154350,462749,308389) : difference with theory (-10) = 0 CORRECTING DEFECT 51 (vertices=80, convex hull=30, v0=128956) After retessellation of defect 51 (v0=128956), euler #=-9 (154354,462772,308409) : difference with theory (-9) = 0 CORRECTING DEFECT 52 (vertices=12, convex hull=32, v0=131341) After retessellation of defect 52 (v0=131341), euler #=-8 (154358,462792,308426) : difference with theory (-8) = 0 CORRECTING DEFECT 53 (vertices=33, convex hull=26, v0=131989) After retessellation of defect 53 (v0=131989), euler #=-7 (154363,462813,308443) : difference with theory (-7) = 0 CORRECTING DEFECT 54 (vertices=36, convex hull=69, v0=132392) After retessellation of defect 54 (v0=132392), euler #=-6 (154372,462866,308488) : difference with theory (-6) = 0 CORRECTING DEFECT 55 (vertices=42, convex hull=68, v0=135908) After retessellation of defect 55 (v0=135908), euler #=-5 (154387,462933,308541) : difference with theory (-5) = 0 CORRECTING DEFECT 56 (vertices=234, convex hull=129, v0=141235) After retessellation of defect 56 (v0=141235), euler #=-4 (154437,463140,308699) : difference with theory (-4) = 0 CORRECTING DEFECT 57 (vertices=283, convex hull=123, v0=145442) After retessellation of defect 57 (v0=145442), euler #=-3 (154479,463323,308841) : difference with theory (-3) = 0 CORRECTING DEFECT 58 (vertices=24, convex hull=55, v0=155621) After retessellation of defect 58 (v0=155621), euler #=-2 (154489,463376,308885) : difference with theory (-2) = 0 CORRECTING DEFECT 59 (vertices=46, convex hull=48, v0=156711) After retessellation of defect 59 (v0=156711), euler #=-1 (154491,463398,308906) : difference with theory (-1) = 0 CORRECTING DEFECT 60 (vertices=21, convex hull=52, v0=162481) After retessellation of defect 60 (v0=162481), euler #=0 (154501,463448,308947) : difference with theory (0) = 0 CORRECTING DEFECT 61 (vertices=25, convex hull=81, v0=164661) After retessellation of defect 61 (v0=164661), euler #=1 (154515,463520,309006) : difference with theory (1) = 0 CORRECTING DEFECT 62 (vertices=20, convex hull=60, v0=166569) After retessellation of defect 62 (v0=166569), euler #=2 (154523,463563,309042) : difference with theory (2) = 0 computing original vertex metric properties... storing new metric properties... computing tessellation statistics... vertex spacing 0.89 +- 0.30 (0.04-->18.04) (max @ vno 17245 --> 29510) face area 0.00 +- 0.00 (0.00-->0.00) performing soap bubble on retessellated vertices for 0 iterations... vertex spacing 0.89 +- 0.30 (0.04-->18.04) (max @ vno 17245 --> 29510) face area 0.00 +- 0.00 (0.00-->0.00) tessellation finished, orienting corrected surface... 234 mutations (33.6%), 463 crossovers (66.4%), 2438 vertices were eliminated building final representation... 15523 vertices and 0 faces have been removed from triangulation after topology correction, eno=2 (nv=154523, nf=309042, ne=463563, g=0) writing corrected surface to /mnt/newdata/projects/avoss/HCP/subjects/avoss_01/T1w/avoss_01/surf/lh.orig... 0.000 % of the vertices (0 vertices) exhibit an orientation change topology fixing took 332.7 minutes 0 defective edges removing intersecting faces 000: 1108 intersecting 001: 39 intersecting 002: 4 intersecting mris_fix_topology utimesec 19956.712609 mris_fix_topology stimesec 14.538138 mris_fix_topology ru_maxrss 978836 mris_fix_topology ru_ixrss 0 mris_fix_topology ru_idrss 0 mris_fix_topology ru_isrss 0 mris_fix_topology ru_minflt 3092953 mris_fix_topology ru_majflt 8 mris_fix_topology ru_nswap 0 mris_fix_topology ru_inblock 1464 mris_fix_topology ru_oublock 14880 mris_fix_topology ru_msgsnd 0 mris_fix_topology ru_msgrcv 0 mris_fix_topology ru_nsignals 0 mris_fix_topology ru_nvcsw 8180 mris_fix_topology ru_nivcsw 11605 FSRUNTIME@ mris_fix_topology lh 5.5451 hours 1 threads #@# Fix Topology rh Mon May 31 16:47:51 CEST 2021 mris_fix_topology -rusage /mnt/newdata/projects/avoss/HCP/subjects/avoss_01/T1w/avoss_01/touch/rusage.mris_fix_topology.rh.dat -mgz -sphere qsphere.nofix -ga -seed 1234 avoss_01 rh reading spherical homeomorphism from 'qsphere.nofix' using genetic algorithm with optimized parameters setting seed for random number genererator to 1234 ************************************************************* Topology Correction Parameters retessellation mode: genetic search number of patches/generation : 10 number of generations : 10 surface mri loglikelihood coefficient : 1.0 volume mri loglikelihood coefficient : 10.0 normal dot loglikelihood coefficient : 1.0 quadratic curvature loglikelihood coefficient : 1.0 volume resolution : 2 eliminate vertices during search : 1 initial patch selection : 1 select all defect vertices : 0 ordering dependant retessellation: 0 use precomputed edge table : 0 smooth retessellated patch : 2 match retessellated patch : 1 verbose mode : 0 ************************************************************* INFO: assuming .mgz format $Id: mris_fix_topology.c,v 1.50.2.1 2016/10/27 22:25:58 zkaufman Exp $ $Id: mrisurf.c,v 1.781.2.6 2016/12/27 16:47:14 zkaufman Exp $ before topology correction, eno=-198 (nv=166804, nf=334004, ne=501006, g=100) using quasi-homeomorphic spherical map to tessellate cortical surface... Correction of the Topology Finding true center and radius of Spherical Surface...done Surface centered at (0,0,0) with radius 100.0 in 12 iterations marking ambiguous vertices... 25542 ambiguous faces found in tessellation segmenting defects... 91 defects found, arbitrating ambiguous regions... analyzing neighboring defects... -merging segment 3 into 1 -merging segment 7 into 1 -merging segment 30 into 14 -merging segment 42 into 39 -merging segment 43 into 39 -merging segment 45 into 39 -merging segment 60 into 57 -merging segment 71 into 70 -merging segment 85 into 82 -merging segment 86 into 82 81 defects to be corrected 0 vertices coincident reading input surface /mnt/newdata/projects/avoss/HCP/subjects/avoss_01/T1w/avoss_01/surf/rh.qsphere.nofix... reading brain volume from brain... reading wm segmentation from wm... Computing Initial Surface Statistics -face loglikelihood: -9.2100 (-4.6050) -vertex loglikelihood: -6.1191 (-3.0596) -normal dot loglikelihood: -3.5662 (-3.5662) -quad curv loglikelihood: -5.9860 (-2.9930) Total Loglikelihood : -24.8813 CORRECTING DEFECT 0 (vertices=5, convex hull=26, v0=0) After retessellation of defect 0 (v0=0), euler #=-80 (152406,453937,301451) : difference with theory (-78) = 2 CORRECTING DEFECT 1 (vertices=230, convex hull=191, v0=43) After retessellation of defect 1 (v0=43), euler #=-77 (152461,454192,301654) : difference with theory (-77) = 0 CORRECTING DEFECT 2 (vertices=173, convex hull=44, v0=49) After retessellation of defect 2 (v0=49), euler #=-76 (152470,454236,301690) : difference with theory (-76) = 0 CORRECTING DEFECT 3 (vertices=233, convex hull=169, v0=1066) After retessellation of defect 3 (v0=1066), euler #=-76 (152522,454480,301882) : difference with theory (-75) = 1 CORRECTING DEFECT 4 (vertices=28, convex hull=64, v0=1566) After retessellation of defect 4 (v0=1566), euler #=-75 (152539,454556,301942) : difference with theory (-74) = 1 CORRECTING DEFECT 5 (vertices=28, convex hull=24, v0=1688) After retessellation of defect 5 (v0=1688), euler #=-74 (152540,454564,301950) : difference with theory (-73) = 1 CORRECTING DEFECT 6 (vertices=52, convex hull=24, v0=2598) After retessellation of defect 6 (v0=2598), euler #=-73 (152541,454572,301958) : difference with theory (-72) = 1 CORRECTING DEFECT 7 (vertices=32, convex hull=29, v0=4188) After retessellation of defect 7 (v0=4188), euler #=-72 (152544,454589,301973) : difference with theory (-71) = 1 CORRECTING DEFECT 8 (vertices=27, convex hull=58, v0=5513) After retessellation of defect 8 (v0=5513), euler #=-71 (152553,454638,302014) : difference with theory (-70) = 1 CORRECTING DEFECT 9 (vertices=12, convex hull=31, v0=6459) After retessellation of defect 9 (v0=6459), euler #=-70 (152554,454648,302024) : difference with theory (-69) = 1 CORRECTING DEFECT 10 (vertices=24, convex hull=34, v0=7952) After retessellation of defect 10 (v0=7952), euler #=-69 (152557,454667,302041) : difference with theory (-68) = 1 CORRECTING DEFECT 11 (vertices=363, convex hull=210, v0=8813) After retessellation of defect 11 (v0=8813), euler #=-68 (152618,454959,302273) : difference with theory (-67) = 1 CORRECTING DEFECT 12 (vertices=7290, convex hull=1501, v0=10262) XL defect detected...