Treatment plan optimization method for radiation therapy

• US 7,801,270 B2
• Filed: 06/19/2008
• Issued: 09/21/2010
• Est. Priority Date: 06/19/2008
• Status: Active Grant

• Current Assignee: Siemens Healthineers International AG
• Original Assignee: Varian Medical Systems International AG (Siemens AG)
• Inventors: Peltola, Jarkko Y., Nord, Janne Ilmari
• Primary Examiner(s): KAO, CHIH CHENG G
• Application Number:
  US12/142,181
• Publication Number:
  US 20090316858A1
• Time in Patent Office: 824 Days
• Field of Search: 378/65
• US Class Current: 378/65
• CPC Class Codes: A61N 5/1031 using a specific method of dose optimization

Abstract

A treatment planning method and system for optimizing a treatment plan used to irradiate a treatment volume including a target volume, such as a tumor, is disclosed. In accordance with the present invention, gradients of a cost function that defines a treatment volume are back projected to a surface of interest. The method and system of preferred embodiments of the present invention calculate the gradients with respect to the dose received by the treatment volume. Machine parameters that are associated with the surface of interest may then be optimized based on the back projected gradients.

Claims

1. A method of optimizing a treatment plan for irradiating a treatment volume within a patient, comprising:

establishing at least one treatment parameter;

establishing a cost function corresponding to a treatment volume, wherein the treatment volume includes a target volume intended to receive a therapeutic dose of radiation and a non-target volume intended to receive a smaller dose of radiation, wherein the cost function is responsive to the radiation dose received by the treatment volume;

back-projecting a plurality of gradients of said cost function to an optimization surface, wherein said gradients are responsive to changes in said at least one treatment parameter and wherein said optimization surface is defined by said at least one treatment parameter; and

changing said at least one treatment parameter, wherein said change in said at least one treatment parameter causes a corresponding change in said back projected gradients;

wherein at least one of establishing at least one treatment parameter, back projecting a plurality of gradients, or changing said at least one treatment parameter is performed with a treatment planning computer.

2. The method of claim 1 further comprising normalizing said back-projected gradients within said optimization surface.

3. The method of claim 2 wherein normalizing said back projection gradients comprises grouping said back projection gradients into at least two groups, wherein said groups are responsive to changes in said at least one treatment parameter.

4. The method of claim 3, wherein said at least one treatment parameter is changed to optimize said groups, wherein optimizing said groups comprises minimizing a sum of said back projection gradients within said groups.

5. The method of claim 4 wherein the sum of said back projected gradients further include an additional factor proportional to a magnitude of the change in said at least one treatment parameter, wherein said factor increases the sum of said back projection gradients, whereby minimizing the sum of said back projection gradients within said groups is favored by smaller changes in said at least one treatment parameter.

6. The method of claim 5 wherein the steps are iterated a plurality of times.

7. The method of claim 6 wherein said treatment plan comprises Intensity Modulated Radiation Therapy.

8. The method of claim 6 wherein said treatment plan comprises Arc treatment.

9. The method of claim 5 wherein the at least one treatment parameter is a treatment machine parameter and said gradients of said cost function are not calculated with respect to said treatment machine parameter.

10. The method of claim 5 wherein the at least one treatment parameter is a treatment machine parameter and said gradients of said cost function are calculated with respect to said radiation dose received by said treatment volume.

11. A method for generating a radiation therapy treatment plan comprising:

establishing at least one treatment machine parameter;

establishing a cost function corresponding to a treatment volume, wherein the treatment volume includes a target volume intended to receive a therapeutic dose of radiation and a non-target volume intended to receive a smaller dose of radiation, wherein the cost function is responsive to the radiation dose received by the treatment volume;

back projecting gradients of said cost function to a surface, wherein said surface is defined by said at least one treatment machine parameter, wherein said gradients are calculated with respect to said radiation dose, and wherein said gradients are responsive to changes in said at least one treatment machine parameter; and

changing said at least one treatment machine parameter, wherein said change in said at least one treatment machine parameter causes a corresponding change in said back projected gradients;

wherein at least one of establishing at least one treatment machine parameter, back projecting gradients of said cost function, or changing said at least one treatment machine parameter is performed with a treatment planning computer.

12. The method of claim 11 further comprising:

segmenting said surface into a matrix, wherein said matrix has at least two elements;

aggregating said back projection gradients that project into each of said at least two matrix elements;

grouping said matrix elements into at least two groups, wherein said groups contain matrix elements that compete for changes in said treatment machine parameter.

13. The method of claim 12 wherein said at least one treatment machine parameter is changed to minimize the sum of said back projected gradients within said groups.

14. The method of claim 13 wherein the sum of said back projected gradients further include an additional factor proportional to a magnitude of the change in said at least one treatment machine parameter, wherein said factor increases the sum of said back projection gradients, whereby minimizing the sum of said back projection gradients within said groups is favored by smaller changes in said at least one treatment machine parameter.

15. The method of claim 14 wherein the steps are iterated a plurality of times.

16. A method for optimizing leaf positions of a multi-leaf collimator used to deliver a radiation treatment plan to a patient comprising:

establishing an initial position for a plurality of leaves of a multi-leaf collimator;

establishing a cost function corresponding to a treatment volume, wherein the treatment volume includes a target volume intended to receive a therapeutic dose of radiation and a non-target volume intended to receive a smaller dose of radiation, wherein the cost function is responsive to the radiation dose received by the treatment volume;

back projecting gradients of said cost function to a plane established by the leaves of said multi-leaf collimator, wherein said gradients are calculated with respect to said radiation dose and wherein said gradients are responsive to changes in said leaf positions; and

changing positions of said leaves in said multi leaf collimator, wherein said change in said positions of said leaves of said multi leaf collimator causes a corresponding change in said back projected gradients.

17. The method of claim 16 further comprising:

establishing a matrix in said plane, wherein each element of said matrix corresponds to a position of a leaf in said multi leaf collimator;

summing all back projected gradients of said cost function that back project to the same matrix element;

grouping said matrix elements into at least two groups, wherein said groups contain matrix elements that compete for changes in said positions of said leaves in said multi leaf collimator.

18. The method of claim 17 wherein said position of said leaves in said multi leaf collimator is changed to minimize the sum of said back projected gradients within said groups.

19. The method of claim 18 wherein the sum of said back projected gradients further include an additional factor proportional to a distance traveled by said leaves of said multi leaf collimator, wherein said factor increases the sum of said back projection gradients, whereby minimizing the sum of said back projection gradients within said groups is favored by smaller distances traveled by said leaves of said multi leaf collimator.

20. The method of claim 19 wherein the steps are iterated a plurality of times.

21. A non-transitory computer readable medium storing thereon a program for execution by a computer, the program having code for performing the steps of:

establishing at least one treatment parameter;

establishing a cost function corresponding to a treatment volume, wherein the treatment volume includes a target volume intended to receive a therapeutic dose of radiation and a non-target volume intended to receive a smaller dose of radiation, wherein the cost function is responsive to the radiation dose received by the treatment volume;

back-projecting a plurality of gradients of said cost function to an optimization surface, wherein said gradients are responsive to changes in said at least one treatment parameter and wherein said optimization surface is defined by said at least one treatment parameter; and

changing said at least one treatment parameter, wherein said change in said at least one treatment parameter causes a corresponding change in said back projected gradients.

22. The non-transitory computer readable medium of claim 21 further comprising:

normalizing said back-projected gradients within said optimization surface, wherein normalizing said back projected gradients comprises grouping said back projection gradients into at least two groups, wherein said groups are responsive to changes in said at least one treatment parameter.

23. The non-transitory computer readable medium of claim 22, wherein said at least one treatment parameter is changed to optimize said groups, wherein optimizing said groups comprises minimizing a sum of said back projection gradients within said groups.

24. The non-transitory computer readable medium of claim 23, wherein the sum of said back projected gradients further include an additional factor proportional to a magnitude of the change in said at least one treatment parameter, wherein said factor increases the sum of said back projection gradients, whereby minimizing the sum of said back projection gradients within said groups is favored by smaller changes in said at least one treatment parameter.