SYSTEMS, DEVICES, AND METHODS FOR QUALITY ASSURANCE OF RADIATION THERAPY

• US 20140105355A1
• Filed: 10/12/2012
• Published: 04/17/2014
• Est. Priority Date: 10/12/2012
• Status: Active Grant

• Current Assignee: Siemens Healthineers International AG
• Original Assignee: Varian Medical Systems International AG (Siemens AG)
• Inventors: NORD, Janne, TOIMELA, Lasse Heikki
• Granted Patent: US 9,381,376 B2
• Time in Patent Office: Days
• US Class Current: 378/41
• CPC Class Codes: A61N 2005/1054 using a portal imaging system
• CPC Class Codes: A61N 2005/1072 taking into account movement of the target
• CPC Class Codes: A61N 5/103 Treatment planning systems
• CPC Class Codes: A61N 5/1031 using a specific method of dose optimization
• CPC Class Codes: A61N 5/1064 for adjusting radiation treatment in response to monitoring
• CPC Class Codes: A61N 5/1071 for verifying the dose delivered by the treatment plan

Abstract

Systems, devices, and methods for quality assurance for verification of radiation dose delivery in arc-based radiation therapy devices using a 3D gamma evaluation method.

Claims

1. A method for verifying quantities of interest of a radiation beam in an arc-based radiation therapy device including a gantry, comprising:

generating radiation dose distribution images at predetermined radiation fields, each radiation field corresponding to a predetermined arc segment of the gantry; and

comparing the generated radiation dose distribution images with corresponding predicted radiation dose distribution images using a three-dimensional (3D) gamma evaluation method comprising calculating one of or in combination of dose differences, spatial differences, and angular differences between points in the generated images and corresponding points in the predicted images as parameters.

2. The method of claim 1, further comprising storing all generated radiation dose distribution images in a data structure where positional information of the points in the generated images are stored together with corresponding values of the points.

3. The method of claim 2, wherein the storing includes storing all generated radiation dose distribution images as a first data set mapped as a first 3D array including 3D positional information of the points in the generated images.

4. The method of claim 3, wherein the 3D positional information includes spatial and angular locations of the delivered beams (X, Y, θ

).

5. The method of claim 3, further comprising storing all predicted radiation dose distribution images as a second data set mapped as a second 3D array including 3D position information of the points in the predicted images.

6. The method of claim 5, wherein the 3D positional information includes spatial and angular locations of the predicted beams (X, Y, θ

).

7. The method of claim 5, wherein the comparing includes comparing the first and second data sets using the 3D gamma evaluation method.

8. The method of claim 1, wherein the comparing includes computing a 3D gamma error matrix based on a gamma error function:

γ

(r_p)=min{Γ

(r_m,r_p,θ

)}∀

{r_m},

wherein each gamma value γ

in the 3D matrix is a minimum value of a distance metric Γ

indicating a measure of acceptability of a quantity of interest.

9. The method of claim 8, further comprising calculating distance metric values for each point in the first data set using:

Γ

((x_m,y_m,θ

_m),(x_p,y_p,θ

_p))=√

{square root over (s_xy(Δ

x²+Δ

y²)+s_d(Δ

D)²+s_θ

Δ

θ

²)}{square root over (s_xy(Δ

x²+Δ

y²)+s_d(Δ

D)²+s_θ

Δ

θ

²)}

wherein Δ

D is the dosimetric difference between measured and predicted dose values, Δ

x=x_m−

x_p and Δ

x=y_m−

y_p are axial differences between corresponding points in the first and second data sets; and

Δ

θ

=θ

_m−

θ

_p are angular differences between corresponding points in the first and second data sets.

10. The method of claim 9, further comprising reducing the 3D gamma error matrix to a 2D gamma error matrix.

11. The method of claim 1, wherein quantities of interest include radiation dose.

12. The method of claim 11, further including determining dose delivery errors based on the comparison.

13. A system for verifying quantities of interest of a radiation beam in an arc-based radiation therapy device including a gantry, comprising:

a portal dose imaging device adapted to measure incident radiation dose from each predetermined radiation field, each radiation field corresponding to a predetermined arc segment of the gantry, the portal dose imaging device being further configured to generate two-dimensional (2D) portal images for each arc segment; and

a processing device operably connected to the portal dose imaging device and being configured to store a plurality of measured 2D portal dose images in a first 3D array having spatial and angular locations of the delivered beams (X, Y, θ

) as dimensions, the processing device being further configured to store a plurality of predicted 2D portal dose images in a second 3D array having spatial and angular locations of the predicted beams (X, Y, θ

) as dimensions;

the processing device being further configured to compare points in the first 3D array with corresponding points in the second 3D array using a 3D gamma evaluation method, the gamma evaluation method including dose differences, spatial differences, and angular differences between corresponding points in the first and second 3D arrays as parameters,

wherein errors in the quantities of interest are determined based on the comparison.

14. The system of claim 13, wherein the arc-field radiation therapy device is an intensity modulated radiation therapy device.

15. The system of claim 14, wherein portal imaging device is an electronic portal imaging device (EPID).

16. The system of claim 15, wherein the processing device is further configured to generate a 3D gamma error matrix based on a gamma error function:

γ

(r_p)=min{Γ

(r_m,r_p,θ

)}∀

{r_m},

wherein each gamma value γ

in the 3D matrix is a minimum value of a distance metric Γ

indicating a measure of acceptability of the quantity of interest.

17. The system of claim 16, wherein the processing device is further configured to calculate distance metric values for each point in the first 3D array using:

Γ

((x_m,y_m,θ

_m),(x_p,y_p,θ

_p))=√

{square root over (s_xy(Δ

x²+Δ

y²)+s_d(Δ

D)²+s_θ

Δ

θ

²)}{square root over (s_xy(Δ

x²+Δ

y²)+s_d(Δ

D)²+s_θ

Δ

θ

²)}

wherein Δ

D is the dosimetric difference between measured and predicted dose values, Δ

x=x_m−

x_p and Δ

x=y_m−

y_p are spatial differences between corresponding points in the first and second 3D arrays; and

Δ

θ

=θ

_m−

θ

_p are angular differences between corresponding points in the first and second 3D arrays.

18. The system of claim 17, wherein the processing device is further configured to reduce the 3D gamma error matrix to a 2D gamma error matrix.

19. A non-transitory computer-readable storage medium upon which is embodied a sequence of programmed instructions for verifying quantities of interest of a radiation beam in a system including an arc-based radiation therapy device including a gantry and a computer processing system which executes the sequence of programmed instructions embodied on the computer-readable storage medium to cause the computer processing system to perform the steps of:

generating radiation dose distribution images at predetermined radiation fields, each radiation field corresponding to a predetermined arc segment of the gantry; and

comparing the generated radiation dose distribution images with corresponding predicted radiation dose distribution images using a three-dimensional (3D) gamma evaluation method including dose differences, spatial differences, and angular differences between points in the generated images and corresponding points in the predicted images as parameters.

20. The non-transitory computer-readable medium of claim 19, further comprising storing all generated radiation dose distribution images in a data structure where positional information of the points in the generated images are stored together with corresponding values of the points.

21. The non-transitory computer-readable medium of claim 20, wherein the storing includes storing all generated radiation dose distribution images as a first data set mapped as a first 3D array including 3D positional information of the points in the generated images.

22. The non-transitory computer-readable medium of claim 21, wherein the 3D positional information includes spatial and angular locations of the delivered beams (X, Y, θ

).

23. The non-transitory computer-readable medium of claim 21, further comprising storing all predicted radiation dose distribution images as a second data set mapped as a second 3D array including 3D position information of the points in the predicted images.

24. The non-transitory computer-readable medium of claim 23, wherein the 3D positional information includes spatial and angular locations of the predicted beams (X, Y, θ

).

25. The non-transitory computer-readable medium of claim 19, wherein the comparing includes comparing the first and second data sets using the 3D gamma evaluation method.

26. The non-transitory computer-readable medium of claim 25, wherein the comparing includes computing a 3D gamma error matrix based on a gamma error function:

γ

(r_p)=min{Γ

(r_m,r_p,θ

)}∀

{r_m}

wherein each gamma value γ

in the 3D matrix is a minimum value of a distance metric Γ

indicating a measure of acceptability of a quantity of interest.

27. The non-transitory computer-readable medium of claim 26, further comprising calculating distance metric values for each point in the first data set using:

Γ

((X_m,y_m,θ

_m),(x_p,y_p,θ

_p))=√

{square root over (s_xy(Δ

x²+Δ

y²)+s_d(Δ

D)²+s_θ

Δ

θ

²)}{square root over (s_xy(Δ

x²+Δ

y²)+s_d(Δ

D)²+s_θ

Δ

θ

²)}

wherein Δ

D is the dosimetric difference between measured and predicted dose values, Δ

x=x_m−

x_p and Δ

x=y_m−

y_p are axial differences between corresponding points in the first and second data sets; and

Δ

θ

=θ

_m−

θ

_p are angular differences between corresponding points in the first and second data sets.

28. The non-transitory computer-readable medium of claim 27, further comprising reducing the 3D gamma error matrix to a 2D gamma error matrix.

29. The non-transitory computer-readable medium of claim 19, wherein quantities of interest include radiation dose.

30. The non-transitory computer-readable medium of claim 29, further including determining dose delivery errors based on the comparison.

31. A method for verifying quantities of interest of a radiation beam in an arc-based radiation therapy device including a gantry, comprising:

delivering radiation at a plurality of radiation fields corresponding to a predetermined arc segment of the gantry;

measuring transmitted or absorbed radiation from each arc segment using an electronic portal dose imaging device;

generating 2D images based on the measured radiation;

converting the 2D images into corresponding 2D portal dose images (PDIs) using a portal dose reconstruction model; and

comparing the obtained 2D portal dose images with corresponding predicted radiation dose images using a three-dimensional (3D) gamma evaluation method including dose differences, spatial differences, and angular differences between points in the generated images and corresponding points in the predicted images as parameters.