You can generally do this with the Slicer RT extension, but you need to also output the geometry, I.e. “region” or “cell” or however your code calls it, as a voxel type labelmap. Then you can use the DVH functionality of Slicer RT.
Note: Uncheck “show dose volumes only” because even though your volume is a dose volume, it’s probably not RT dose volume so it won’t immediately show up. The results will be displayed as “intensity volume histogram” but it’s equivalent to DVH.
For the first option: As far as I know, the dose volume (RT dose) and segmentation (RT structure) files are requered for the DVH to calculate and plot the DVHs, and these files can only be obtained from DICOM files.
Please, tell me how can i plot the DVH with 3D dose distribution file.
Thanks in advance.
You usually have to run a separate tally, or calculation mode to output the geometry from your MC program. I haven’t used the PENELOPE code, but for example in many particle transport
codes, there is an option to sample the material number (integer) , or region ID (integer) number in the voxel mesh (instead, or at the same time as, the dose tally) – there should be a description of how to do this within
PENELOPE’s documentation. Not sure if this helps, but how I usually make DVHs in Slicer is to:
a) run the dose simulation with the MC code to get 3D dose distribution which is imported into Slicer as a ‘volume’ (a reader may need to be made for it as cpinter mentioned)
b) run the geometry option with the MC code to get a voxel representation of the geometry, and import that into Slicer as a segmentation, or as a ‘labelmap’ which can be immediately converted to segmentation within Slicer
c) then calculate DVH within the DVH module of the slicer RT add-on