I am a user of 3Dslicer Mac version 5.2.1, with MarkUpsToModel, Sandbox and SegmentEditorExtraEffects extensions.
I used the threshold tool on a paralumbar muscle after obtaining its volume and extraction but it did not succeed.
My goal is to obtain a muscle volume or a muscle area in binary color: black for the pure muscle and white for the fat in the muscle, via a predefined cutoff, like this exemple:
This is what I get when using threshold, and it doesn’t work:
Thank you for your help.
Indeed, I tried different threshold tools (SegmentEditor → Local threshold // Volume → threshold // Filtering → threshold scalar volume), but not this one.
Once the threshold is done and saved, how to get a percentage of black voxel and white voxel in a table? Currently, I get this table in Quantification → SegmentStatistics but it doesn’t separate the 2 colors I want to get.
Within a fixed region, If you set one threshold to define muscle + fat and set it as a segment, then create a second segment and threshold specific to the muscle only, you can do what I said above, which is
(number of voxels in segment1 - number of voxels in segment2) / number of voxels in segment1 *100
which will give you the percentage of fat voxels in the defined area.
But you have to keep the area exactly the same in both segments.
I tried to apply your method but there is an error:
If I explode it in voxels, my volume of interest is 621650 voxels.
With the same thresholding at 75 :
If I choose to take the fat (white) in my volume of interest, it makes 188201 voxels.
If I choose to take the pure muscle (black) in my volume of interest, it makes 431195 voxels.
If I sum the two, I should in theory find the voxels of my total volume, that is 621650 voxels. Now, 188201+431195=619396.
→ Do you know where the lost voxels are?
Moreover, in the table created with the Segment statistics tool, they describe volume (1) (columns B C D), volume (2) (columns E F G), volume (3) (columns N O) cf image. What do the (1) (2) (3) correspond to?
Actually I found my mistake do it this way and you will get correct values
Download MRhead sample data
create three empty segments
In the first segment using the paint brush in 3D mode create a spherical region (I used 10%) in the center of the brain
Copy this region to segment_2 and Segment_3 using Logical Operators
Switch to Segment_2, go to threshold tool and enter these threshold values 0 - 34.99 and then set the Masking fields to: Editable area to Inside Segment_2 and then modify other segments: Overwrite Visible (your segment_3 should be visible, that important). Then hit apply.
Turn the visibility on for segment_1, go to segment editor and calculate values.
You will see that sum of voxel counts of segment_2 and segment_3 are now identical to the segment_1. This procedure split the original region into two none overlapping segmentations. You can adopt this procedure to your own dataset.
I was able to apply your threshold method on MRhead and then on my MRI. This technique is really easier, thank you.
Concerning the volumes and the “Segment Statistics” table, I understood that (1) (2) and (3) represent different volumes, but according to you, which volume represents better the reality (i.e. the real volume) between Labelmap statistics, scalar volume statistics or closed surface statistics?
A topic has already been opened concerning this question but there is no answer ( Differences between closed surface volume, labelmap volume and scalar volume ).