@@ -222,6 +222,15 @@ def _create_headmask_t1_pipe(name="headmask_t1_pipe", params={}):
222222 inputnode , ('indiv_params' , parse_key , "t1_head_erode" ),
223223 t1_head_erode , "indiv_params" )
224224
225+ # mesh_t1_skull #######
226+ mesh_t1_head = pe .Node (
227+ IsoSurface (),
228+ name = "mesh_t1_head" )
229+
230+ headmask_t1_pipe .connect (
231+ t1_head_erode , "out_file" ,
232+ mesh_t1_head , "nii_file" )
233+
225234 # t1_hmasked
226235 t1_hmasked = pe .Node (interface = ApplyMask (),
227236 name = "t1_hmasked" )
@@ -325,6 +334,26 @@ def _create_skullmask_t1_pipe(name="skullmask_t1_pipe", params={}):
325334 t1_head_erode_skin , "out_file" ,
326335 t1_head_skin_masked , "mask_file" )
327336
337+ # mesh_t1_rawskull #######
338+ mesh_t1_rawskull = pe .Node (
339+ IsoSurface (),
340+ name = "mesh_t1_rawskull" )
341+
342+ skullmask_t1_pipe .connect (
343+ t1_head_skin_masked , "out_file" ,
344+ mesh_t1_rawskull , "nii_file" )
345+
346+ else :
347+
348+ # mesh_t1_rawskull #######
349+ mesh_t1_rawskull = pe .Node (
350+ IsoSurface (),
351+ name = "mesh_t1_rawskull" )
352+
353+ skullmask_t1_pipe .connect (
354+ t1_skull_mask_binary , "out_file" ,
355+ mesh_t1_rawskull , "nii_file" )
356+
328357 if "t1_skull_gcc_erode" in params and \
329358 "t1_skull_gcc_dilate" in params :
330359
@@ -497,8 +526,9 @@ def create_skull_t1_pipe(name="skull_t1_pipe", params={}):
497526 outputnode = pe .Node (
498527 niu .IdentityInterface (
499528 fields = ["t1_skull_mask" , "t1_skull_stl" ,
529+ "t1_rawskull_mask" , "t1_rawskull_stl" ,
500530 "robustt1_skull_mask" , "robustt1_skull_stl" ,
501- "t1_head_mask" ]),
531+ "t1_head_mask" , "t1_head_stl" , ]),
502532 name = 'outputnode' )
503533
504534 # Creating headmask_t1_pipe
@@ -528,6 +558,10 @@ def create_skull_t1_pipe(name="skull_t1_pipe", params={}):
528558 headmask_t1_pipe , "t1_head_erode.out_file" ,
529559 outputnode , "t1_head_mask" )
530560
561+ skull_t1_pipe .connect (
562+ headmask_t1_pipe , "mesh_t1_head.stl_file" ,
563+ outputnode , "t1_head_stl" )
564+
531565 # Creating skullmask_t1_pipe
532566 if "skullmask_t1_pipe" in params :
533567 skullmask_t1_pipe = _create_skullmask_t1_pipe (
@@ -544,14 +578,30 @@ def create_skull_t1_pipe(name="skull_t1_pipe", params={}):
544578 else :
545579 return skull_t1_pipe
546580
581+ # outputnode
547582 skull_t1_pipe .connect (
548583 skullmask_t1_pipe , "mesh_t1_skull.stl_file" ,
549584 outputnode , "t1_skull_stl" )
550585
551586 skull_t1_pipe .connect (
552- skullmask_t1_pipe , "t1_skull_erode .out_file" ,
587+ skullmask_t1_pipe , "t1_skull_mask_binary .out_file" ,
553588 outputnode , "t1_skull_mask" )
554589
590+ # rawskull t1
591+ skull_t1_pipe .connect (
592+ skullmask_t1_pipe , "mesh_t1_rawskull.stl_file" ,
593+ outputnode , "t1_rawskull_stl" )
594+
595+ if "t1_head_erode_skin" in params .keys ():
596+ skull_t1_pipe .connect (
597+ skullmask_t1_pipe , "t1_head_skin_masked.out_file" ,
598+ outputnode , "t1_rawskull_mask" )
599+
600+ else :
601+ skull_t1_pipe .connect (
602+ skullmask_t1_pipe , "t1_skull_mask_binary.out_file" ,
603+ outputnode , "t1_rawskull_mask" )
604+
555605 if "t1_skull_fov" in params .keys ():
556606 skull_t1_pipe .connect (
557607 skullmask_t1_pipe , "t1_skull_fov.out_roi" ,
0 commit comments