Intensity-modulated proton therapy

• US 9,220,920 B2
• Filed: 02/03/2014
• Issued: 12/29/2015
• Est. Priority Date: 12/06/2011
• Status: Active Grant

• Current Assignee: Haifa University Economic Corporation Ltd, LOMA Linda University Medical Center (Mission Healthcare Services LLC)
• Original Assignee: LOMA Linda University Medical Center (Mission Healthcare Services LLC)
• Inventors: Schulte, Reinhard W., Censor, Yair
• Primary Examiner(s): Bayat, Ali
• Application Number:
  US14/171,490
• Publication Number:
  US 20140350322A1
• Time in Patent Office: 694 Days
• Field of Search: 382/100
• US Class Current: 1/1
• CPC Class Codes: A61N 2005/1087 Ions; Protons
• CPC Class Codes: A61N 5/1031 using a specific method of dose optimization
• CPC Class Codes: A61N 5/1039 using functional images, e.g. PET or MRI
• CPC Class Codes: A61N 5/1077 Beam delivery systems
• CPC Class Codes: A61N 5/1084 for delivering multiple intersecting beams at the same time, e.g. gamma knives

Abstract

The therapeutic treatment of a patient using intensity-modulated proton therapy is described. In one example, a method of creating a proton treatment plan is presented that divides volumes of interest into sub-volumes, applies dose constraints to the sub-volumes, finds one or more feasible configurations of a proton therapy system, and selects a proton beam configuration that improves or optimizes one or more aspects of proton therapy. In some implementations, the method of dividing volumes into sub-volumes includes creating fractional sub-volumes based at least in part on proximity to a target volume boundary. In some implementations, the method of finding an improved or optimal proton beam configuration from a set of feasible configurations includes finding a minimum of a cost function that utilizes weighting factors associated with treatment sites.

Claims

1. A method for performing intensity-modulated ion therapy, the method comprising:

one or more hardware processors configured to perform the steps of;

generating a representation of a patient by;

for a plurality of ions, estimating a path taken by each ion within the patient based at least in part on trajectories determined from measured data;

determining a system of equations of the form Ax=b where x represents a distribution parameter associated with the patient, b represents interaction quantities of the plurality of ions resulting from interactions within the patient, and A represents an operator that operates on x to yield b, the operating having information about the estimated paths of the plurality of ions in the patient;

determining one or more feasible solutions to the system of equations by perturbing an estimated initial solution to determine whether the perturbation yields a superior reconstruction than a non-perturbed solution; and

calculating the distribution parameter associated with the patient;

identifying a volume of interest from the representation of the patient, the volume of interest divided into a total number of voxels ordered according to distance to a feature of interest;

defining a plurality of sub-volumes wherein a sub-volume comprises a number of consecutive ordered voxels, wherein a ratio of the number consecutive ordered voxels comprising the sub-volume to the total number of voxels is approximately equal to a target fractional value for the sub-volume; and

generating an ion treatment plan that satisfies dose constraints for each of the plurality of sub-volumes.

2. The method of claim 1, wherein the ions include protons.

3. The method of claim 1, wherein the ions include carbon ions.

4. The method of claim 1 further comprising delivering ions to the patient based on the ion treatment plan.

5. The method of claim 1 further comprising setting a dose constraint for each of the plurality of sub-volumes.

6. The method of claim 1, wherein generating the ion treatment plan comprises using an inverse problem solver to determine ion beam characteristics predicted to deliver radiation doses to the plurality of sub-volumes within the dose constraints.

7. The method of claim 1, wherein estimating a path taken by each ion within the patient comprises calculating an estimated path corresponding to an actual path taken by each ion within the patient, the estimated path providing an estimate of the actual path.

8. The method of claim 1, wherein the system of equations of the form Ax=b is a system of linear equations wherein A is a matrix, x is a vector, and b is a vector.

9. The method of claim 8, wherein an element of the matrix A corresponds to an estimated intersection length of a selected ion in a corresponding voxel, the estimated intersection length being calculated as a straight-line approximation of the estimated path of the selected ion in the corresponding voxel.

10. The method of claim 1, wherein the interaction quantity comprise an energy loss of a corresponding ion resulting from its passage through the patient.

11. The method of claim 1, wherein the distribution parameter comprises a distribution of an electron density-based quantity.

12. The method of claim 11, wherein the electron density-based quantity comprises relative proton stopping power with respect to a substantially uniform material.

13. The method of claim 12, wherein the substantially uniform material comprises water.

14. An intensity-modulated ion therapy system comprising:

an ion delivery system configured to deliver a plurality of ions according to an ion treatment plan;

a computed tomography system configured to generate a representation of a patient, the computed tomography system comprising a hardware processor configured to;

for a plurality of ions, estimate a path taken by each ion within the patient based at least in part on trajectories determined from measured data;

determine a system of equations of the form Ax=b where x represents a distribution parameter associated with the patient, b represents interaction quantities of the plurality of ions resulting from interactions within the patient, and A represents an operator that operates on x to yield b, the operating having information about the estimated paths of the plurality of ions in the patient;

determine one or more feasible solutions to the system of equations by perturbing an estimated initial solution to determine whether the perturbation yields a superior reconstruction than a non-perturbed solution; and

calculate the distribution parameter associated with the patient; and

a treatment plan system configured to generate the ion treatment plan, the treatment plan system comprising a hardware processor configured to;

identify a volume of interest from the representation of the patient, the volume of interest divided into a total number of voxels ordered according to distance to a feature of interest;

define a plurality of sub-volumes wherein a sub-volume comprises a number of consecutive ordered voxels, wherein a ratio of the number consecutive ordered voxels comprising the sub-volume to the total number of voxels is approximately equal to a target fractional value for the sub-volume; and

generate the ion treatment plan that satisfies dose constraints for each of the plurality of sub-volumes.

15. The system of claim 14, wherein the plurality of ions include protons.

16. The system of claim 14, wherein the plurality of ions include carbon ions.

17. The system of claim 14 further comprising an inverse problem solver module within the treatment plan system, the inverse problem solver module configured to determine ion beam characteristics predicted to deliver radiation doses to the plurality of sub-volumes within the dose constraints.

18. The system of claim 14 further comprising a forward problem solver module within the treatment plan system, the forward problem solver module configured to determine predicted doses to the patient based at least in part on ion beam characteristics.

19. The system of claim 14, wherein the computed tomography system is configured to generate a three-dimensional representation of the patient.

20. The system of claim 14 further comprising a detector system configured to measure trajectories of ions before and after passing through a portion of the patient and to measure energies of ions after passing through the portion of the patient.