Measurements on a 3D Region of Interest  

Analyzing and taking measurements on a part of a 3D volume is a common task. ImageMeasurements, extended to work with 3D volumes, provides a large number of properties that can be computed on a region of interest.

Segment a volume to the region of interest.

In[1]:=

X -Graphics-

In[2]:=

X mask = ImageForestingComponents[vol, {{1, 1, 124}, {102, 1, 1}}]; vol2 = ImageMultiply[Image3D[Erosion[mask - 1, 2]], vol]; seg1 = MorphologicalComponents[vol2, 0.075]; roi = Image3D[SelectComponents[seg1, "Count", -1], "Bit"]

Out[2]=

Compute a single measure.

In[3]:=

X ImageMeasurements[vol, "Mean", Masking -> roi]

Out[3]=

Compute a list of properties.

In[4]:=

X props = {"Min", "Max", "Mean", "Median", "StandardDeviation", "Total", "Skew", "IntensityCentroid"}; ImageMeasurements[vol, props, Masking -> roi]

Out[4]=

Out[5]//TraditionalForm=