Hi, I’m having some trouble loading certain DICOM segmentation’s into the Slicer scene. Importing the file is no problem (imported using DICOM database), but the error log says “could not load: [filename] as a DICOMSegmentation”. Does anyone know why this might be? Thanks
If you only have this problem with certain segmentation files then we need one of those files so that we can investigate. You can upload the files anywhere and post the link here (make sure no patient health information is in the files).
If you cannot share the files, it may also be helpful to share the full error log (Help > Report a bug).
Are you able to load some DICOM segmentations, but not others, or you cannot load any?
Hello, here is an example of a segmentation file that I am unable to upload. The full error message is just “could not load as a DICOMSegmentation”. Thanks everyone for your help.
Thank you, I can reproduce the problem.
I believe the segmentation series here is the one with SeriesDescription “VOI PE Segmentation thresh=70” from this study available on IDC: https://viewer.imaging.datacommons.cancer.gov/viewer/188.8.131.52.4.1.145184.108.40.206.7695.1700.186130315530827043916665138479.
There are several issues:
- the MR series corresponding to this segmentation is a 4D series
- the segmentation is encoded as FRACTIONAL
It looks like this combination is not handled well at the moment by either Slicer or IDC/OHIF Viewer.
Loading 4D MRI series should work well in Slicer already. Slicer can also represent FRACTIONAL segmentations, although right now fractional segmentations are still just loaded as binary. So, all segmentations should be possible to load.
The problem with this file is that it contains an empty segmentation. It has a single slice, with all zero values. Probably it is just a placeholder.
While checking this, I’ve found a bug in the segmentation importer for loading single-frame segmentations, and I’ve submitted a fix.
Andras, thank you for the fix. I merged it. I don’t think it will solve the underlying issue for the user though. We need to add a check to the SEG DICOM plugin to account for the situation when the segment is empty. I admit I didn’t even think of a possibility of someone saving a segment that is blank. Thank you for investigating this so quickly!
In this case, there was nothing in the segmentation, so the user got all that was available (nothing).
But yes, the importer has to be prepared for blank segments. For example, if a segmentation has blank and non-blank segments, probably import stops due to the failed import of the first blank segment and the rest of the segments are not imported. At least we should file submit a bug report to QuantitativeReporting to remember to fix this at some point.
I’ve just added a quick fix here: https://github.com/QIICR/QuantitativeReporting/commit/1ca31ce1681f9e04dc3129fdca42abe815c31e06
I confirmed with this fix there is no error, segmentation node is created, but has no segments, and there is a warning in the error log. I think now the behavior of the application makes sense.
It could be nice to create the segment anyway and set its name, color, terminology, etc. instead of skipping it completely (if it is easy to implement).
Yes, I agree it would be nice, and it’s an interesting “hack” to work around the lack of an equivalent of a “pick list” in DICOM SEG/SR.
I believe at the moment the behavior of empty segment import is such that no empty segment is created. The corresponding line is highlighted below. I think this would need to be implemented in the segmentations logic. I do not know how difficult it would be to implement.
It would be hard to tell what and how many empty segments
segmentationsLogic.ImportLabelmapToSegmentationNode should create. We could add an option that would create an empty segment for each entry defined in the color node associated with the input labelmap. But this would be quite a special case, so I’m not sure if it would worth the time to implement.
Instead DICOMSegmentationPlugin could be updated so that if
_importSegmentAndRemoveLabel returns an empty segmentation then it adds a new segment (using
This actually would be very easy, since when dcmqi converts DICOM SEG, each segment is saved into a separate file. So when that import function is called on a label map, the label map has one and only segment.
In any case, I created a ticket to keep track of this idea, if/when we get to it: https://github.com/QIICR/QuantitativeReporting/issues/253.