IP Library Granted Patent US 10,940,333
Granted Patent B2
US 10,940,333 · App. 15/994,077 · Granted Mar 9, 2021

Systems, methods, and devices for radiation beam alignment and radiation beam measurements using electronic portal imaging devices

Inventors: Stefan Thieme (Windisch, CH); Mathias Lehmann (Zurich, CH)
Assignee: VARIAN MEDICAL SYSTEMS INTERNATIONAL AG
A61N5/1075A61B6/4241A61B6/4258A61B6/5258A61N5/1045A61N5/1049A61B6/06A61B6/4233A61N2005/1054A61N2005/1074A61N2005/1089
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 10,940,333
App. No.
15/994,077
Granted
Mar 9, 2021
Kind
B2
Abstract

Systems and methods for using electronic portal imaging devices (EPIDs) as absolute radiation beam measuring devices and as radiation beam alignment devices without implementation of elaborate and complex calibration procedures.

Claims (60)

1. A method of determining radiation beam characteristics using an imaging device, the method comprising:

acquiring one or more images using the imaging device, the imaging device including a plurality of pixels;

determining one or more parameters from the one or more images; and

determining one or more characteristics of the radiation beam from the determined one or more parameters,

wherein the one or more characteristics comprises number of detected photons, radiation beam energy change, radiation beam tilt relative to a collimator axis of rotation, radiation beam symmetry change, radiation beam flatness change, and radiation beam center change,

wherein the determining of the one or more characteristics is independent of variations in pixel gain values.

2. The method of claim 1 , wherein determining the number of detected photons comprises:

generating a plurality of images;

measuring pixel values and one or more noise values of a pixel across the plurality of images;

calculating a mean pixel value of the pixel from the measured pixel values; and

determining a value corresponding to detected photons per pixel based on the mean pixel value and the one or more noise values for the pixel.

3. The method of claim 2 wherein a value corresponding to detected photons is determined for each pixel of the plurality of pixels.

4. The method of claim 2 , wherein the value corresponding to detected photons per pixel is determined using: N=(p/σ) 2 , where N is the number of photons per pixel, p is the mean pixel value, and σ is a standard deviation of the pixel across the plurality of images.

5. The method of claim 2 , wherein the mean pixel value is determined for a plurality of pixels located in a region of interest (ROI) of an image.

6. The method of claim 5 , further comprising normalizing the pixel value using the mean pixel value of the ROI pixels.

7. The method of claim 1 , wherein the imaging device is an electronic portal dose imaging device (EPID).

8. The method of claim 1 , wherein determining the number of detected photons comprises:

determining a region of interest (ROI) in an image, the ROI including a plurality of pixels;

determining gain compensated pixel values for the plurality of pixels;

measuring noise values for the plurality of pixels;

determining a mean pixel value of the gain compensated pixel values;

determining a value corresponding to detected photons per pixel based on the mean pixel value and the noise values.

9. The method of claim 8 , wherein the determining of the gain compensated pixel values for the plurality of pixels includes determining a pixel gain map and correcting each pixel value of the plurality of pixels using the pixel gain map.

10. The method of claim 8 , wherein the value corresponding to detected photons per pixel is determined using: N=(p/σ) 2 , where N is the number of photons per pixel, p is the mean pixel value of the gain compensated pixel values, and σ is a standard deviation noise value for the pixels in the ROI.

11. The method of claim 1 , wherein determining the beam energy change comprises:

generating a first and a second image;

determine a first pixel value for a pixel by measuring pixel signal in the first image, and determine a second pixel value of the pixel by measuring a pixel signal in the second image;

calculate a ratio of the first pixel value to the second pixel value; and

determine a value of a parameter that provides a best fit between the calculated ratio and a pixel location in a pixel panel of the imaging device,

wherein the parameter relates to the beam energy change.

12. The method of claim 11 , wherein the best fit value is determined using a least-square fit algorithm.

13. The method of claim 1 , wherein determining the beam flatness change and the beam symmetry change comprises:

generating a first and a second image;

determine a first pixel value for a pixel by measuring a pixel signal in the first image, and determine a second pixel value of the pixel by measuring a pixel signal in the second image;

calculate a ratio of the first pixel value to the second pixel value;

determine a plurality of parameter values that provide a best fit between the calculated ratio and a pixel location in a pixel panel of the imaging device; and

determine the beam flatness change and the beam symmetry change based on the determined parameter values.

14. The method of claim 13 , wherein a best fit value is determined using a least-square fit algorithm applied to: ρ=d 0 ·(x 2 +y 2 )+d 1 ·x+d 2 ·y+d 3 , where ρ is the log of the ratio, (x, y) represents the pixel location along the x and y imaging panel axis, and d 0 , d 1 , d 2 , and d 3 are first, second, third and fourth parameters, respectively, wherein applying the least-square fit algorithm includes discarding d 3 ,

the method further comprising converting a first parameter value into flatness change, and converting the second and third parameters into symmetry changes along the x and y axis.

15. The method of claim 13 , further comprising correlating the beam flatness change with the beam energy change, and the beam symmetry change with the beam center change.

16. The method of claim 15 , further comprising generating a look-up table including correlation values between the beam flatness change and the beam energy change, and between the beam symmetry change and the beam center change.

17. The method of claim 1 , wherein determining the radiation beam tilt relative to the collimator axis of rotation comprises:

using the imaging device positioned at a predetermined distance (SDD) from a radiation source, acquiring a first image at a first beam energy, and a second image at a second beam energy;

determining first pixel values by measuring pixel signals for pixels in the first image, and second pixel values by measuring pixel signals for corresponding pixels in the second image;

calculating ratios of the first pixel values to the second pixel values;

generating a ratio image from the calculated ratio pixel values, the ratio image depicting a beam shape;

determining a center of the beam shape;

calculating a distance d between the center of the beam shape and a projection of the radiation source on the imaging device along the collimator axis of rotation; and

determining the radiation beam tilt based on the calculated distance d and the predetermined distance (SDD).

18. The method of claim 17 , further comprising calculating the radiation beam tilt relative to a collimator axis of rotation includes calculating a beam tilt angle relative to the collimator axis of rotation using:

sin

α

=

d

SDD

.

19. The method of claim 18 , further comprising translating the calculated beam tilt angle into photon fluence symmetry.

20. The method of claim 1 , further including calibrating a radiation treatment device including the imaging device based on the determined one or more radiation beam characteristics.

Assignments (2)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Jan 10, 2023
From: LEHMANN, MATHIAS; THIEME, STEFAN
To: VARIAN MEDICAL SYSTEMS INTERNATIONAL AG
Reel/Frame 062325/0824 →
CHANGE OF NAME Recorded Jan 10, 2023
From: VARIAN MEDICAL SYSTEMS INTERNATIONAL AG
To: SIEMENS HEALTHINEERS INTERNATIONAL AG
Reel/Frame 062338/0726 →
Continuity (2)
Division 15017280 · Feb 5, 2016
Related Publication 20180272155A1 · Sep 27, 2018
Cited By (2)
US 12,194,315 US 12,317,406