Name

mri_watershed - strip skull and other outer non-brain tissue

Arguments

Positional Arguments

invol	input volume
outvol	output volume

Optional Flagged Arguments

-w weight	preweight the input image using atlas information (no_w -> 1)
-no_wta	don't use the preweighting for the template deformation
-b proba_merging	use the basins merging using atlas information (no_b -> 1)
-copy	just copy input to output, ignore other options
-atlas	use the atlas information to correct the segmentation.	When the segmented brain is not correct, this option might help you.
-surf [surfname]	save the BEM surfaces.	In order to get the surfaces consistent with tkmedit, you have to use the option -useSRAS.
-useSRAS	use the surface RAS coordinates (not the scanner RAS) for surfaces.
-noT1	don't do T1 analysis. (Useful when running out of memory)
-less	shrink the surface
-more	expand the surface
-wat	use only the watershed algorithm
-T1	specify T1 input volume (T1 grey value = 110)
-wat+temp	watershed algo and first template smoothing
-first_temp	use only the first template smoothing + local matching
-surf_debug	visualize the surfaces onto the output volume
-brainsurf [surfname]	save the brain surface
-shk_br_surf [int_h surfname]	to save the brain surface shrank inward of int_h mm
-s [int_i int_j int_k]	add a seed point
-c [int_i int_j int_k]	specify the center of the brain (in voxel unit)
-r int_r	specify the radius of the brain (in voxel unit)
-t int_threshold	change the threshold in the watershed analyze process
-h int_hpf	precize the preflooding height (in percent)
-n	not use the watershed analyze process
-LABEL	labelize the output volume into scalp, skull, csf, gray and white
-man [int_csf int_trn int_gray]	to change the different parameters csf_max, transition_intensity and GM_intensity
-mask	mask a volume with the brain mask
--help	show usage message
--version	show the current version

Outputs

brainvol	skull stripped brain volume
BEMsurfaces	when you specify the option -brainsurf surfname

Produce the brain volume from T1 volume or the scanned volume.

Description

This program strips skull and other non-brain tissue and produces the brain volume from T1 volume or the scanned volume.

If you used the -atlas option, then the template was used to correct the surface.

The default parameters are: -w 0.82 -b 0.32 -h 10 -seedpt -ta -wta

You can use one of the first five flags to change these default parameters.

Examples

Example 1

mri-watershed -atlas T1 brain

where T1 is the T1 volume and brain is the output brain volume. When the cerebellum is cut-off from the brain or getting the left/right asymmetric brain, you should first try this -atlas option.

Example 2

mri-watershed T1 brain

The same as the first example, but no correction is applied to the intermediate result.

Bugs

None

Links

FreeSurfer, FsFast

Methods Description

The "watershed" segmentation algorithm was used to determine the intensity values for white matter, grey matter, and CSF.
A force field was then used to fit a spherical surface to the brain. The shape of the surface fit was then evaluated against a previously derived template.

The finely grained sphere was fit to the brain.

(Segonne 2004)

"-surf" Description:
The BEM surfaces are generated after the initial mri_watershed procedure has run and produced a "brain" surface file. This brain surface file is then grown by 3mm outward to generate the "inner skull" surface. This growth is meant to account for the CSF located between the pial matter and the inner skull (see Figure 17 in Segonne et al. 2004 for a description). Next a fifth order icosahedral surface is fit arround the outer edge of the Volume and smoothed making the "outer skin" surface. Finally, this surface is shrunk by 3mm to attempt to account for the SCALP (Skin, Connective tissue, Aponeurosis, Loose areolar connective tissue, and the Pericranium). This surface is saved as the "outer skull" surface.

The -surf tool is intended as an added convenience and not as a high quality estimation of these tissue boundaries (one would not use an MPRAGE type volume for defining these tissue boundaries). For a tool designed to do that please see mri_make_bem_surfaces.

Reporting Bugs

Report bugs to < freesurfer@nmr.mgh.harvard.edu > YasunariTosa

-  ⇤ ← Revision 14 as of 2004-08-18 18:07:49 → 
  Size: 3437
  Editor: YasunariTosa
  Comment:
+   ← Revision 29 as of 2018-01-04 10:17:45 → ⇥
  Size: 5169
  Editor: MorganFogarty
  Comment:
-Deletions are marked like this.
+Additions are marked like this.
 Line 1:
+## Note: This page was created with the CommandTemplate
##
## If you're modifying this page please take a look at the
## latest version of CommandTemplate to ensure that you're
## using the latest version of the CommandTemplate
##
-Line 2:
+Line 8:
-##
## Note: If you are modifying this page please refer to the
## CommandTemplate page and confirm that this page has the latest
## CommandTemplate formatting
+<<Navigation(children)>> '''Index''' <<TableOfContents>>
-Line 8:
+Line 11:
-mri_watershed - strips skull and other outer non-brain voxels from an acquired volume (usually T1).
+mri_watershed - strip skull and other outer non-brain tissue
-Line 11:
+Line 13:
-= Synopsis =
mri_watershed  [<options>] <input volume> <brain volume>
+= Arguments =
== Positional Arguments ==
||invol ||input volume ||
||outvol ||output volume ||
-Line 14:
+Line 18:
-= Positional Arguments =
-Line 16:
+Line 19:
-= Required Flagged Arguments =
+== Optional Flagged Arguments ==
||-w weight || preweight the input image using atlas information (no_w -> 1) || ||
||-no_wta || don't use the preweighting for the template deformation || ||
||-b proba_merging || use the basins merging using atlas information (no_b -> 1) || ||
||-copy || just copy input to output, ignore other options || ||
||-atlas ||use the atlas information to correct the segmentation. ||When the segmented brain is not correct, this option might help you. ||
||-surf [surfname] ||save the BEM surfaces. ||In order to get the surfaces consistent with tkmedit, you have to use the option -useSRAS. ||
||-useSRAS ||use the surface RAS coordinates (not the scanner RAS) for surfaces. || ||
||-noT1 ||don't do T1 analysis.  (Useful when running out of memory) || ||
||-less ||shrink the surface || ||
||-more ||expand the surface || ||
||-wat ||use only the watershed algorithm || ||
||-T1 ||specify T1 input volume (T1 grey value = 110) || ||
||-wat+temp ||watershed algo and first template smoothing || ||
||-first_temp ||use only the first template smoothing + local matching || ||
||-surf_debug ||visualize the surfaces onto the output volume || ||
||-brainsurf [surfname] ||save the brain surface || ||
||-shk_br_surf [int_h surfname] ||to save the brain surface shrank inward of int_h mm || ||
||-s [int_i int_j int_k] ||add a seed point || ||
||-c [int_i int_j int_k] ||specify the center of the brain (in voxel unit) || ||
||-r int_r ||specify the radius of the brain (in voxel unit) || ||
||-t int_threshold ||change the threshold in the watershed analyze process || ||
||-h int_hpf ||precize the preflooding height (in percent) || ||
||-n ||not use the watershed analyze process || ||
||-LABEL ||labelize the output volume into scalp, skull, csf, gray and white || ||
||-man [int_csf int_trn int_gray] ||to change the different parameters csf_max, transition_intensity and GM_intensity || ||
||-mask ||mask a volume with the brain mask || ||
||--help ||show usage message || ||
||--version ||show the current version || ||
-Line 18:
+Line 49:
-= Optional Flagged Arguments =
|| -atlas || use the atlas information to correct the segmentation.  When the segmented brain is not correct, this option might help you. ||
|| -surf [surfname] || save the BEM surfaces.  In order to get the surfaces consistent with tkmedit, you have to use the option -useSRAS.  ||
|| -useSRAS || use the surface RAS coordinates (not the scanner RAS) for surfaces. ||
|| -noT1    || don't do T1 analysis.  (Useful when running out of memory) ||
|| -less    || shrink the surface ||
|| -more    || expand the surface ||
|| -wat     || use only the watershed algorithm ||
|| -T1      || specify T1 input volume (T1 grey value = 110) ||
|| -wat+temp  || watershed algo and first template smoothing ||
|| -first_temp || use only the first template smoothing + local matching ||
|| -surf_debug || visualize the surfaces onto the output volume ||
|| -brainsurf surfname  || save the brain surface ||
|| -shk_br_surf int_h surfname || to save the brain surface shrank inward of int_h mm ||
|| -s int_i int_j int_k || add a seed point ||
|| -c int_i int_j int_k || specify the center of the brain (in voxel unit) ||
|| -r int_r             || specify the radius of the brain (in voxel unit) ||
|| -t int_threshold     || change the threshold in the watershed analyze process ||
|| -h int_hpf           || precize the preflooding height (in percent) ||
|| -n                   || not use the watershed analyze process ||
|| -LABEL               || labelize the output volume into scalp, skull, csf, gray and white ||
|| -man int_csf int_trn int_gray ||  to change the different parameters csf_max, transition_intensity and GM_intensity  ||
|| -mask                || mask a volume with the brain mask ||
|| --help               || show this usage message ||
|| --version            || show the current version ||
-Line 45:
+Line 51:
-|| <brain volume> || skull stripped brain volume ||
|| [BEM surfaces] || when you specify the option -brainsurf surfname ||
+||brainvol ||skull stripped brain volume ||
||BEMsurfaces ||when you specify the option -brainsurf surfname ||


Produce the brain volume from T1 volume or the scanned volume.
-Line 49:
+Line 58:
-See Florent Segonne et al. paper for details.
+This program strips skull and other non-brain tissue and produces the brain volume from T1 volume or the scanned volume.
-Line 51:
+Line 60:
-First uses "watershed" process to figure out the grey scale values for white matter, grey matter, and CSF.
Second uses the force fields to shrink the rough sphere onto the brain.
Third evaluate the shape using the template ($FREESURFER_HOME/average/rigidly_aligned_brain_template.tiff)
for correctness.  If you used -atlas option, then use the template to correct the shape.   
Fourth use the finely grained sphere to fit onto the brain.
+If you used the -atlas option, then the template was used to correct the surface.
-Line 57:
+Line 62:
-= Example 1 =
mri-watershed foo -i f -o out
+The default parameters are: -w 0.82 -b 0.32 -h 10 -seedpt -ta -wta
-Line 60:
+Line 64:
-description
+You can use one of the first five flags to change these default parameters.
-Line 62:
+Line 66:
-= Example 2 =
mri-watershed foo -i f -o out -f fvalue
+= Examples =
== Example 1 ==
{{{
mri-watershed -atlas T1 brain
}}}
where T1 is the T1 volume and brain is the output brain volume.    When the cerebellum is cut-off from the brain or getting the left/right asymmetric brain, you should first try this -atlas  option.
-Line 65:
+Line 73:
-description
+== Example 2 ==
{{{
mri-watershed T1 brain
}}}
The same as the first example, but no correction is applied to the intermediate result.

= Bugs =
None
-Line 68:
+Line 83:
-["othercommand1"], ["othercommand2"]
+[[mri_normalize]]
-Line 75:
+Line 90:
-description
description
+The "watershed" segmentation algorithm was used to determine the intensity values for white matter, grey matter, and CSF.
A force field was then used to fit a spherical surface to the brain. The shape of the surface fit was then evaluated against a previously derived template.

The finely grained sphere was fit to the brain.

(Segonne 2004)

"-surf" Description:
The BEM surfaces are generated after the initial mri_watershed procedure has run and produced a "brain" surface file. This brain surface file is then grown by 3mm outward to generate the "inner skull" surface. This growth is meant to account for the CSF located between the pial matter and the inner skull (see Figure 17 in Segonne et al. 2004 for a description). Next a fifth order icosahedral surface is fit arround the outer edge of the Volume and smoothed making the "outer skin" surface. Finally, this surface is shrunk by 3mm to attempt to account for the SCALP (Skin, Connective tissue, Aponeurosis, Loose areolar connective tissue, and the Pericranium). This surface is saved as the "outer skull" surface.
-Line 78:
+Line 100:
-= References =
["References/Lastname###"]
+The -surf tool is intended as an added convenience and not as a high quality estimation of these tissue boundaries (one would not use an MPRAGE type volume for defining these tissue boundaries). For a tool designed to do that please see [[mri_make_bem_surfaces]].
-Line 83:
+Line 103:
-Report bugs to <analysis-bugs@nmr.mgh.harvard.edu>

= Author/s =
JaneSmith
+Report bugs to < freesurfer@nmr.mgh.harvard.edu > YasunariTosa