Hi everyone, I am a beginner at using 3D Slicer. Before that, I used to work with Mimics. I was wondering is it possible to extract/export the whole set of Hounsfield Unit values from a segmentation performed via threshold tool? I want to use them for calculating bone density. There is a possibilit…

Yes, of course, getting intensities of voxels in a segmented region is very simple, and the best is that you don’t even need to export the file, then load it into your own Python script for analysis and display, but you can do everything using a small Python code snippet within Slicer (copy-paste th…

Dear Andras, Thanks for your answer. I did not work with Python, so I am a very beginner in this area. I tried to copy-paste the code that you sent. I renamed the “masterVolumeNode”, “segmentationNode” and “labelmapVolumeNode” but it did not work. I was wondering could you please edit the code ba…

You can set segmentationNode like this: segmentationNode = getNode('P3-L4')

Thank you for your answer. I tried this code and it ran without any error, but I did not see any results: segmentationNode = getNode(‘P3-L4’) labelmapVolumeNode = getNode(‘P3-L4’) masterVolumeNode = getNode(‘P3-whole’) # Get segmentation as labelmap volume node labelmapVolumeNode = slicer.mrml…

There are a few small issues in your script, for example: remove labelmapVolumeNode = getNode('P3-L4') - you set the variable a few lines below *sgmentVoxels should be segmentVoxels Always put code between triple-ticks: ```python your code is here just another line ``` Otherwise, the code loo…

I have the same issue and followed your post for a few days. Today I tried your code but the same as your result. Also I was prompted that I can try this plugin. Thanks to the experts who wrote it. Still can’t get the whole set of Hounsfield Unit values from a segmentation. It was just enough for m…

I’ve tested this example script and it works perfectly for both latest Slicer Stable Release and Slicer Preview Release. [image] All you need to do is to replace the first part (that generates input data) with setting your data as masterVolumeNode, segmentationNode, labelValue.

Thanks for your answer. The reason that I said the table does not make sense: As I mentioned before, the table which I need should contain 4 columns, like this: [image] And considering the CT-scans that I used (0.3 mm) I expected around 250,000 data for this sample (something like the data that…

Since you have the volume and labels as numpy arrays, you already have all these tables. You just need to specify what you need. For example, to get coordinates and HU at all segment voxel positions: coordinates = np.where(labelArray==labelValue) hu = volumeArray[coordinates] You can put together…

Thank you for your answer dear Andras. I still can not get the table that I want, but I am working on it. :frowning: I need to get density for each element in the sample, so I wrote a code in Matlab to obtain the density. One of the input data for that code is the Hu value based on CT images which…

Export Hounsfield Unit values as a table from each segmentation

Support

Topic		Replies	Views
Extract HounsField Units of Segmentation and Save Support segmentation	19	1424	March 22, 2022
Export Hounsfield Unit values from segmentation Support segmentation	1	335	April 16, 2021
Histogram with the Hounsfield scale of segmented Support segmentation , python	7	2014	January 30, 2024
How to export voxel intensity Support	19	3569	December 5, 2022
Extracting Hounsfield units of segmentation Support segmentation , dicom , stl	0	657	January 14, 2021

Export Hounsfield Unit values as a table from each segmentation

Related topics