If you know a software for drawing dose volume histogram (DVH) using numbers in any format other than DICOM format, it can be helpful. thank you.
You're welcome Maryam Zamanian and wish you merry & healthy holidays!
Matlab is very EASY, professional and scientific tool for plotting and programming. I suggest MATLAB, see e.g., https://www.youtube.com/watch?v=mu0NCzgBPaE
https://au.mathworks.com/videos/what-is-medical-imaging-toolbox-1663045797378.html
You can much more via googling further....
Simple practice: https://www.tutorialspoint.com/matlab/matlab_plotting.htm
Muhammad Ali Thank you for the links you sent.
The point of drawing the dose volume histogram (DVH-in Radiotherapy) is that there is no direct relationship between X and Y values at the end of the diagram. At the end of all lines should hit the x-axis without the final value being zero, which is unclear to me.
Okay Maryam Zamanian,
You can plot everything related to graphing or coding.
Simle Histogram ploting - MATLAB:: https://au.mathworks.com/help/matlab/ref/matlab.graphics.chart.primitive.histogram.html
Muhammad Ali Thank you Muhammad Ali. I must say that even though this plot is called a histogram, it is not a histogram, and is most like a linear plot.
the challenge is that the end of the plot is at the highest dose value on the x-axis, even though the corresponding value on the y-axis (volume) is about 80- 90, for correctness from a clinical point of view (in radiotherapy to check the accuracy of the dose) it must hit to the x-axis and have a value of zero (it is drawn with the treatment planning software in bed mode).
Thus, it may be interesting for you to say that drawing it outside the treatment planning software is not compatible with the simple rule of plots.
import matplotlib.pyplot as plt
import numpy as np
# Define the dose values and volume values
dose = np.array([100, 90, 80, 70, 60, 50, 40, 30, 20, 10])
volume = np.array([10, 20, 30, 40, 50, 60, 70, 80, 90, 100])
# Create a figure and a subplot
fig, ax = plt.subplots()
# Plot the DVH using a bar chart
ax.bar(volume, dose, width=5)
# Set the x-axis and y-axis labels
ax.set_xlabel('Volume (% of total volume)')
ax.set_ylabel('Dose (Gy)')
# Show the plot
plt.show()
import matplotlib.pyplot as plt
# Load the data from an Excel file
import pandas as pd
df = pd.read_excel('dvh_data.xlsx')
# Extract the dose and volume columns
dose = df['dose']
volume = df['volume']
# Create a figure and axis
fig, ax = plt.subplots()
# Plot the DVH using a bar chart
ax.bar(dose, volume)
# Add axis labels and a title
ax.set_xlabel('Dose (Gy)')
ax.set_ylabel('Volume (cc)')
ax.set_title('Dose-Volume Histogram')
# Display the plot
plt.show()
Thank you, dear Suganth Velumani, your python code was exact and perfect.
You might be interested to see the code result compared to the expected DVH histogram.
hi,
this page may help you.
DVH Analysis - PlatiPy documentation (pyplati.github.io)
- Badges
- Science topic