Inverse treatment planning method

• US 8,183,541 B2
• Filed: 02/02/2007
• Issued: 05/22/2012
• Est. Priority Date: 02/09/2006
• Status: Active Grant

• Current Assignee: Deutsches Krebsforschungszentrum Stiftung Des Oeffentlichen Rechts
• Original Assignee: Deutsches Krebsforschungszentrum Des Oeffentlichen Rechts
• Inventors: Wilkens, Jan, Oelfke, Uwe
• Primary Examiner(s): Vanore, David A
• Assistant Examiner(s): Ippolito, Nicole
• Application Number:
  US12/223,712
• Publication Number:
  US 20090060130A1
• Time in Patent Office: 1,936 Days
• Field of Search: 600/1, 600/2, 607/1, 250/491.1, 250/492.1, 250/492.3
• US Class Current: 250/492.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/1042 with spatial modulation of the radiation beam within the treatment head

Abstract

The invention relates to a radiation treatment apparatus and an inverse treatment planning method for intensity modulated particle therapy for the treatment of a target within a biological system using at least two fields, each field comprising a plurality of Bragg peaks, the at least two fields being planned to place a defined number of beam spots j from different directions with certain weights w_j within the target. The inverse treatment planning method optimizes the weights w_j of the beam spots for the at least two fields simultaneously in order to produce a prescribed biological effect within the target by minimizing an objective function based on biological effects ε, the biological effects ε being treated in a linear-quadratic model, which describes the biological effects in the target by two parameters α and β, where ε=αD+βD², D denoting a dose, and wherein the two parameters α and β for each voxel i of the target (α_i and β_i) are calculated as dose-averaged mean values of α_i,j- and √{square root over (β_i,j)}-components, which relate to all beam spots j contributing to a total dose D_i in the voxel i.

Claims

1. An inverse treatment planning method for intensity modulated particle therapy for the treatment of a target within a biological system using at least two fields, each field comprising a plurality of Bragg peaks, the at least two fields being planned to place a defined number of beam spots j from different directions with certain weights w_j within the target, wherein the inverse treatment planning method optimizes the weights w_j of the beam spots for the at least two fields simultaneously in order to produce a prescribed biological effect within the target by minimizing an objective function based on biological effects ε

, the biological effects ε

being treated in a linear-quadratic model which describes the biological effects in the target by two parameters α and

β

, where ε

=α

D+β

D², D denoting a dose, and wherein the two parameters α and

β

for each voxel i of the target α

_i and β

_i are calculated as dose-averaged mean values of α

_i,j- and √

{square root over (β

_i,j)}-components, which relate to all beam spots j contributing to a total dose D_i in the voxel i,

wherein the two parameters α and

β

for each voxel i of the target (α

_i and β

_i) are calculated as dose-averaged mean values of the form

2. The inverse treatment planning method according to claim 1, wherein an objective function is minimized, which is based on the formula

3. The inverse treatment planning method according to claim 1, wherein an objective function is minimized, which is based on the formula

4. The inverse treatment planning method according to claim 1, wherein an objective function is minimized, which includes constraints relating to minimum and maximum biological effect levels for the target and maximum effect levels for organs at risk or normal tissue within the biological system.

5. The inverse treatment planning method according to claim 1, wherein the α

_i,j- and √

{square root over (β

_i,j)}-components are derived from a radiobiological model or from measured data for α and

β

as a function of depth in a single Bragg peak.

6. The inverse treatment planning method according to claim 1, wherein a constant value β

_i,j=β

_i is set, which depends only on a tissue type at voxel i.

7. The inverse treatment planning method according to claim 1, wherein the α

_i,j-components are assumed to be laterally constant and are extracted from a single curve as a function of residual range, independent of an initial beam energy.

8. Use of an inverse treatment planning method for intensity modulated particle therapy for the treatment of a target within a biological system using at least two fields, each field comprising a plurality of Bragg peaks, the at least two fields being planned to place a defined number of beam spots j from different directions with certain weights w_j within the target, wherein the inverse treatment planning method optimizes the weights w_j of the beam spots for the at least two fields simultaneously in order to produce a prescribed biological effect within the target by minimizing an objective function based on biological effects ε

, the biological effects ε

being treated in a linear-quadratic model which describes the biological effects in the target by two parameters α and

β

, where ε

=α

D+β

D², D denoting a dose, and wherein the two parameters α and

β

for each voxel i of the target α

_i and β

_i are calculated as dose-averaged mean values of α

_i,j- and √

{square root over (β

_i,j)}-components, which relate to all beam spots j contributing to a total dose D_i in the voxel i for planning the treatment of a tumor with an ion beam containing ions with an atomic number ≧

2.

9. A radiation treatment apparatus for delivering a radiation treatment to a target within a biological system comprising

a field unit for generating at least two fields, each field comprising a plurality of Bragg peaks, the at least two fields being provided to place a defined number of beam spots j from different directions with certain weights w_j within the target and

a unit for inverse treatment planning being provided for optimizing the weights w_j of the beam spots for the at least two fields simultaneously before the radiation treatment in order to produce a prescribed biological effect within the target by minimizing an objective function based on biological effects ε

, the biological effects ε

being treated in a linear-quadratic model, which describes the biological effects in the target by two parameters α and

β

, where ε

=α

D+β

D², D denoting a dose, and wherein the two parameters α and

β

for each voxel i of the target (α

_i, β

_i) are calculated as dose-averaged mean values of α

_i,j- and √

{square root over (β

_i,j)}-components, which relate to all beam spots j contributing to a total dose D_i in the voxel i.