segment statistics for each slice of segmented volume

to calculate the radiation dose received by the fetus during a CT exam, I need to know the mean Hounsfield Unit for each slice of my segmented volume (the uterus). However, I am only able to obtain statistics for the whole segmented volume. I need to know the same information but for each slice.
Is there any way to do so?

Thank you in advance.

You can use the SegmentGeometry extension to compute average intensity in a segmentation along a chosen axis.

1 Like

Thank you for your response.
I tried using SegmentGeometry, however, when I compared the mean slice brightness to the mean Hounsfield Unit of a single slice, the results were not the same.
If only I could get the mean Hounsfield Unit of each slice rather than the mean Hounsfield Unit of the whole segmented volume.

Do you need the average intensity of the entire slice or only of the part of the slice that is inside the segmented structure? Maybe you could add a screenshot to explain it.

I have uploaded two images. The first one is showing the mean HU of the whole volume, and the second one shows the mean intensity of each slice of the volume.
What I would like is to have the mean HU of each slice and not of the whole volume.

What you show above looks correct. Mean brightness values in the table are higher because the irrelevant air voxels (with value of -1024) are excluded from the statistics.

The air voxels must be excluded because otherwise the computed mean intensity would depend on the bore diameter of the CT, the chosen padding value (voxel values outside the reconstructed image), patient support, table appearance (if the image acquisition software does not remove it), cropping settings, etc.

Ok, I see what you’re saying. As a result, getting the mean HU by slice for a volume is not an option.
Thank you for your responses; I appreciate your assistance.

You can easily get the mean HU for the entire image slice, it is just not a meaningful value in general.

Ah, the reason the values are different is because there was an error in the code. I forgot intensity values could be negative when working with Hounsfield units and was filtering out anything zero or below. I have fixed the issue and the update will be available tomorrow. If you have the current stable release, you should be able to update Segment Geometry through the extension manager. If you have the preview release, you’ll need to download and install tomorrow’s build.

1 Like

Thanks a lot for the update Jonathan! I will try again tomorrow then.
Kind Regards

1 Like