SYSTEM AND METHOD FOR PROTON THERAPY TREATMENT PLANNING WITH PROTON ENERGY AND SPOT OPTIMIZATION

Abstract

Embodiments of the present invention disclose methods and systems for proton therapy planning that includes proton energy and spot optimization that discretizes layers and spots using an optimization algorithm to produce an optimal distribution of layer energies and spots with a relatively smooth dose distribution. The treatment planning algorithms disclosed herein can freely choose the number of spots and the energy levels of the spots. In this way, each spot can be treated as its own layer and is not constrained by the requirements of other spots/layers. Thereafter, the spots defined by the algorithm can be sorted in a list according to energy levels/depth, and the spots can be grouped into blocks according to intensity and location. The blocks can be assigned energy levels based on the corresponding spots, such as an average of all the spots associated with the block. The blocks then are used as the energy layers applied by the proton therapy treatment system.

Claims

1. A system for proton therapy treatment, comprising: a gantry comprising a nozzle configured to emit a controllable proton beam; a proton therapy treatment system operable to automatically control the gantry according to a treatment plan; and a treatment planning system comprising: a memory; and a processor operable to perform a method of generating the treatment plan, the method comprising: receiving treatment planning parameters; determining a number of spots and energy levels of the spots according to the treatment planning parameters; calculating intensities for the spots according to the energy levels; discretizing the spots into a plurality of blocks according to locations of the spots and intensities of the spots, wherein the treatment plan comprises the plurality of blocks, and wherein the blocks represent an energy layer applied by the controllable proton beam; and outputting the treatment plan.

2. The system as described in claim 1, wherein the method further comprises treating a patient using the gantry according to the treatment plan.

3. The system as described in claim 1, wherein said discretizing the spots into a plurality of blocks is performed iteratively until a clinical goal is achieved.

4. The system as described in claim 1, wherein the method further comprises determining a plurality of angles for performing proton therapy treatment and wherein the determining a number of spots and energy levels of the spots is further based on the plurality of angles.

5. The system as described in claim 1, wherein the method further comprises receiving a number of layers as user input, where the determining a number of spots and energy levels of the spots is further based on the number of layers.

6. The system as described in claim 1, wherein the method further comprises sorting the spots according to depth before the discretizing the spots.

7. The system as described in claim 1, wherein the method further comprises sorting the spots according to energy level before discretizing the spots.

8. The system as described in claim 1, wherein the energy layers comprise a fixed spacing.

9. The system as described in claim 1, wherein the energy layers have a fixed total intensity.

10. A method of proton therapy treatment, the method comprising: receiving treatment planning parameters; determining a number of spots and energy levels of the spots according to the treatment planning parameters; sorting the spots based on a depth of a respective spot; calculating intensities for the spots according to the energy levels; discretizing the spots into a plurality of blocks according to locations of the spots and intensities of the spots; and outputting a treatment plan, wherein the treatment plan comprises the plurality of blocks, and wherein the blocks represent an energy layer applied by a controllable proton beam.

11. The method as described in claim 10, further comprising controlling a gantry structure according to the treatment plan to perform particle therapy on a patient.

12. The method as described in claim 10, wherein said discretizing the spots into a plurality of blocks is performed iteratively until a clinical goal is achieved.

13. The method as described in claim 10, further comprising determining a plurality of angles for performing proton therapy treatment and wherein the determining a number of spots and energy levels of the spots is further based on the plurality of angles.

14. The method as described in claim 10, further comprising receiving a number of layers as user input, where the determining a number of spots and energy levels of the spots is further based on the number of layers.

15. The method as described in claim 10, wherein the energy layers comprise a fixed spacing.

16. The method as described in claim 10, wherein the energy layers comprise a fixed total intensity.

17. A non-transitory computer-readable storage medium having embedded therein program instructions, which when executed by one or more processors of a device, causes the device to execute a method of proton therapy treatment, the method comprising: receiving treatment planning parameters; determining a number of spots and energy levels of the spots according to the treatment planning parameters; sorting the spots based on a depth of respective spots; calculating intensities for the spots according to associated energy levels thereof; discretizing the spots into a plurality of blocks according to locations of the spots and intensities of the spots, wherein the blocks represent an energy layer applied during proton therapy treatment; and outputting a treatment plan, wherein the treatment plan comprises the plurality of blocks, and wherein the blocks represent an energy layer applied by a controllable proton beam.

18. The non-transitory computer-readable storage medium as described in claim 17, wherein the method further comprising automatically controlling a gantry structure according to the treatment plan to perform particle therapy on a patient.

19. The non-transitory computer-readable storage medium as described in claim 17, wherein said discretizing the spots into a plurality of blocks is performed iteratively until a clinical goal is achieved.

20. The non-transitory computer-readable storage medium as described in claim 17, wherein the method further comprising determining a plurality of angles for performing proton therapy treatment and wherein the determining a number of spots and energy levels of the spots is further based on the plurality of angles.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0019] The accompanying drawings, which are incorporated in and form a part of this specification, illustrate embodiments of the invention and, together with the description, serve to explain the principles of the invention:

[0020] FIG. 1 shows a block diagram of an example of a computing system upon which the embodiments described herein may be implemented.

[0021] FIG. 2 is a block diagram showing selected components of a radiation treatment system upon which embodiments according to the present invention can be implemented.

[0022] FIG. 3 illustrates elements of a radiation treatment system in accordance with embodiments according to the invention.

[0023] FIG. 4 is a block diagram illustrating components in a process for creating an optimized proton therapy treatment plan in embodiments according to the present invention.

[0024] FIG. 5 is a flow-chart depicting an exemplary sequence steps for automatically creating an optimized proton therapy treatment plan according to embodiments of the present invention.

[0025] FIG. 6 is a diagram of an exemplary proton therapy treatment plan with predetermined energy layers.

[0026] FIG. 7 is a diagram of an exemplary proton therapy treatment plan optimized according to a spot distribution using customized energy layers having a fixed intensity per layer.

DETAILED DESCRIPTION

[0027] Reference will now be made in detail to several embodiments. While the subject matter will be described in conjunction with the alternative embodiments, it will be understood that they are not intended to limit the claimed subject matter to these embodiments. On the contrary, the claimed subject matter is intended to cover alternative, modifications, and equivalents, which may be included within the spirit and scope of the claimed subject matter as defined by the appended claims.

[0028] Furthermore, in the following detailed description, numerous specific details are set forth in order to provide a thorough understanding of the claimed subject matter. However, it will be recognized by one skilled in the art that embodiments may be practiced without these specific details or with equivalents thereof. In other instances, well-known methods, procedures, components, and circuits have not been described in detail as not to unnecessarily obscure aspects and features of the subject matter.

[0029] Portions of the detailed description that follows are presented and discussed in terms of a method. Although steps and sequencing thereof are disclosed in a figure herein describing the operations of this method, such steps and sequencing are exemplary. Embodiments are well suited to performing various other steps or variations of the steps recited in the flowchart (e.g., FIG. 5) of the figures herein, and in a sequence other than that depicted and described herein.

[0030] Some portions of the detailed description are presented in terms of procedures, steps, logic blocks, processing, and other symbolic representations of operations on data bits that can be performed on computer memory. These descriptions and representations are the means used by those skilled in the data processing arts to most effectively convey the substance of their work to others skilled in the art. A procedure, computer-executed step, logic block, process, etc., is here, and generally, conceived to be a self-consistent sequence of steps or instructions leading to a desired result. The steps are those requiring physical manipulations of physical quantities. Usually, though not necessarily, these quantities take the form of electrical or magnetic signals capable of being stored, transferred, combined, compared, and otherwise manipulated in a computer system. It has proven convenient at times, principally for reasons of common usage, to refer to these signals as bits, values, elements, symbols, characters, terms, numbers, or the like.

[0031] It should be borne in mind, however, that all of these and similar terms are to be associated with the appropriate physical quantities and are merely convenient labels applied to these quantities. Unless specifically stated otherwise as apparent from the following discussions, it is appreciated that throughout, discussions utilizing terms such as “generating,” “writing,” “including,” “storing,” “transmitting,” “traversing,” “associating,” “identifying,” “optimizing” or the like, refer to the action and processes of a computer system, or similar electronic computing device, that manipulates and transforms data represented as physical (electronic) quantities within the computer system's registers and memories into other data similarly represented as physical quantities within the computer system memories or registers or other such information storage, transmission or display devices.

[0032] Some embodiments may be described in the general context of computer-executable instructions, such as program modules, executed by one or more computers or other devices. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. Typically the functionality of the program modules may be combined or distributed as desired in various embodiments.

Energy and Spot Optimization for Proton Therapy Treatment

[0033] The following description is presented to enable a person skilled in the art to make and use the embodiments of this invention; it is presented in the context of a particular application and its requirements. Various modifications to the disclosed embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be applied to other embodiments and applications without departing from the spirit and scope of the present disclosure. Thus, the present invention is not limited to the embodiments shown, but is to be accorded the widest scope consistent with the principles and features disclosed herein.

[0034] Methods and systems for proton therapy planning including proton energy and spot optimization that discretizes layers and spots using an optimization process are described herein. The optimization process produces an optimal distribution of layer energies and spots with a relatively smooth dose distribution, and can freely choose the number of spots and the energy levels of the spots. In this way, each spot can be treated as its own layer and is not constrained by the requirements of other spots/layers. Thereafter, the spots defined by the process can be sorted into a list according to energy levels/depths, and the spots can be grouped into blocks. The blocks can be assigned energy levels based on the corresponding spots, for example, using an average of all of the spots associated with the block. The blocks are then used to define the energy levels applied by the proton therapy treatment system.

[0035] Embodiments of the present invention can thereby produce fewer layers and potentially fewer spots, thereby advantageously reducing the time required to treat a patient. The optimization processes, e.g., algorithms, described herein can split individual spots among different layers to avoid creating spots that are below the minimum energy threshold deliverable by the treatment system. The spots can be divided into layers such that the layers receive substantially equal intensities. For example, the spots can be divided into layers so that each layer receives a substantially equal number of protons delivered by the proton therapy treatment system. In general, a treatment plan that is optimized using fewer layers leads to a faster delivery time. On the other hand, embodiments of the present invention can also optimize proton treatment plans for accuracy by using more spots and/or layers to better fit the target volume. The spots can be distributed at fixed intervals (e.g., 3 mm spacing) or distributed based on variations such as beam energy and/or beam range. In all cases, the general principle holds that the intensity of a layer is the sum of the intensities of all spots on a given layer.

[0036] FIG. 1 shows a block diagram of an example of a computing system 100 upon which the embodiments described herein may be implemented. In a basic configuration, the system 100 includes at least one processing unit 102 and a memory 104. This most basic configuration is illustrated in FIG. 1 by dashed line 106. The system 100 may also have additional features and/or functionality. For example, the system 100 may also include additional storage (removable and/or non-removable) including, but not limited to, magnetic or optical disks or tape. Such additional storage is illustrated in FIG. 1 by removable storage 108 and non-removable storage 120. The system 100 may also contain communications connection(s) 122 that allow the device to communicate with other devices, e.g., in a networked environment using logical connections to one or more remote computers.

[0037] The system 100 of FIG. 1 also includes input device(s) 124 such as keyboard, mouse, pen, voice input device, touch input device, etc. Output device(s) 126 such as a display device, speakers, printer, etc., are also included.

[0038] In the example of FIG. 1, the memory 104 includes computer-readable instructions, data structures, program modules, and the like. Depending on how it is to be used, the system 100—by executing the appropriate instructions or the like—can be used to implement a planning system used that discretizes layers and spots using an optimization algorithm (comprising computer readable instructions) to produce an optimal distribution of layer energies and spots with a relatively smooth dose distribution. The treatment planning algorithms executed by CPU 102 can freely choose the number of spots and the energy levels of the spots so that each spot can be treated as its own layer unconstrained by the requirements of other spots/layers. Thereafter, the spots defined by the algorithm can be sorted in a list according to energy levels/depth using CPU 102, and the spots can be grouped into blocks. The blocks can be assigned energy levels based on corresponding spots, for example, using an average of all spots associated with the block. More generally, system 100 can be used to generate and/or optimize proton therapy treatment plans in accordance with the present invention.

[0039] FIG. 2 is a block diagram showing selected components of a radiation treatment system 200 upon which embodiments according to the present invention can be implemented. In the example of FIG. 2, the system 200 includes an accelerator and beam transport system 204 that generate and/or accelerate a beam 201. Embodiments according to the invention can generate and deliver beams of various types including, for instance, proton beams, electron beams, neutron beams, photon beams, ion beams, or atomic nuclei beams (e.g., using elements such as carbon, helium, or lithium) and the like. The operations and parameters of the accelerator and beam transport system 204 are controlled so that the intensity, energy, size, and/or shape of the beam are dynamically modulated or controlled during treatment of a patient according to an optimized radiation treatment plan produced by system 100 as discussed above.

[0040] A recent radiobiology study has demonstrated the effectiveness of delivering an entire, relatively high therapeutic radiation dose to a target within a single, short period of time. This type of treatment is referred to generally herein as FLASH radiation therapy (FLASH RT). Evidence to date suggests that FLASH RT advantageously spares normal, healthy tissue from damage when that tissue is exposed to only a single irradiation for only a very short period of time. For FLASH RT, the accelerator and beam transport system 204 can generate beams that can deliver at least four (4) grays (Gy) in less than one second, and may deliver as much as 20 Gy or 50 Gy, or more, in less than one second. The control system 210 can execute a treatment plan for FLASH RT, and the plan can be generated or optimized by system 100 executing an optimization algorithm in accordance with embodiments of the present invention.

[0041] The nozzle 206 is used to aim the beam toward various locations (a target) within a patient supported on the patient support device 208 (e.g., a chair, couch, or table) located in a treatment room. A target may be an organ, a portion of an organ (e.g., a volume or region within the organ), a tumor, diseased tissue, or a patient outline, for instance.

[0042] The nozzle 206 may be mounted on or may be a part of a gantry (FIG. 3) that can be moved relative to the patient support device 208, which may also be moveable. In embodiments, the accelerator and beam transport system 204 are also mounted on, or are a part of, the gantry structure; in another embodiment, the accelerator and beam transport system are separate from (but in communication with) the gantry structure.

[0043] The control system 210 of FIG. 2 receives and implements a prescribed treatment plan generated and/or optimized according to embodiments of the present invention. In embodiments, the control system 210 includes a computing system having a processor, memory, an input device (e.g., a keyboard), and perhaps optionally a display; the system 100 of FIG. 1 is an example of such a platform for the control system 210. The control system 210 can receive data regarding the operation of the system 200. The control system 210 can control parameters of the accelerator and beam transport system 204, nozzle 206, and patient support device 208, including parameters such as the energy, intensity, size, and/or shape of the beam, direction of the nozzle, and position of the patient support device (and the patient) relative to the nozzle, according to data the control system 210 receives and according to the radiation treatment plan.

[0044] FIG. 3 illustrates elements of a radiation treatment system 300 for treating a patient 304 in embodiments according to the present invention. The system 300 is an example of an implementation of the radiation treatment system 200 of FIG. 2, for example. In embodiments, the gantry 302 and nozzle 306 can be moved up and down the length of the patient 304 and/or around the patient, and the gantry and nozzle can move independently of one another. While the patient 304 is supine in the example of FIG. 3, the invention is not so limited to this orientation. For example, the patient 304 can instead be seated in a chair or standing. The gantry 302 can be controlled by a treatment system using an optimized treatment plan generated according to embodiments of the present invention.

[0045] With regard to FIG. 4, an exemplary proton therapy system 400 for imaging and treating a patient 304 is depicted according to embodiments of the present invention. In the example of FIG. 4, patient 304 is imaged using an image system 402 that uses, for example, x-rays, magnetic resonance imaging (MM), and computed tomography (CT). When CT or MM imagery, for example, is used, a series of two-dimensional (2D) images are taken from a 3D volume. Each 2D image is an image of a cross-sectional “slice” of the 3D volume. The resulting collection of 2D cross-sectional slices can be combined to create a 3D model or reconstruction of the patient's anatomy (e.g., internal organs). The 3D model will contain organs of interest, which may be referred to as structures of interest herein. Those organs of interest include the organ targeted for radiation therapy (a target), as well as other organs that may be at risk of radiation exposure during treatment. According to some embodiments, the imaging process is a separate process from the treatment planning process, and the treatment planning process can include receiving imaging data from a prior imaging session, for example.

[0046] One purpose of the 3D model is the preparation of a radiation treatment plan. To develop a patient-specific radiation treatment plan, information is extracted from the 3D model to determine parameters such as organ shape, organ volume, tumor shape, tumor location in the organ, and the position or orientation of several other structures of interest as they relate to the affected organ and any tumor. The radiation treatment plan can specify, for example, how many radiation beams to use and from which angle each of the beams will be delivered.

[0047] In embodiments according to the present invention, the images from the image system 402 are input to a planning system 404. In embodiments, the planning system 404 includes a computing system having a processor, memory, an input device (e.g., a keyboard), and a display. The system 100 of FIG. 1 is an example of a platform for the planning system 404.

[0048] Continuing with reference to FIG. 4, the planning system 404 executes software that is capable of producing an optimized treatment plan for treating patient 304. The treatment planning system 404 can receive imagery data generated by image system 402 to determine constraints regarding patient anatomy along with treatment planning parameters, such as other patient-related data and machine parameters. The treatment planning parameters are discretized into spots with defined energy levels, and the spots are sorted according to their depth/energy level. The planning system 404 can optimize the treatment plan according to an optimization algorithm that determines an optimal distribution of layers with a relatively smooth dose distribution. The treatment plan can be optimized in accordance with various prescribed results. For instance, the treatment plan can be optimized to use fewer layers leading to a faster delivery time, or can be optimized for accuracy or 3D dose homogeneity by using more layers to conform to the target volume. The planning system 404 can also receive user input that defines the number of layers used in the optimized plan 408, for example.

[0049] The treatment planning system 404 outputs an optimized plan 408 according to an iterative optimizing algorithm 406 as described herein according to embodiments of the present invention. The spots can be distributed at fixed intervals (e.g., 3 mm spacing) or based on variations such as beam energy and/or beam range, and the spots are added to a list of spots that is sorted according to energy level or depth of the spot. Based on the energy levels of the spots and/or the depth/location of a spot, the spots are assigned intensities, where the intensity of a layer is the sum of the intensities of all spots of the layer. The spots are then discretized into blocks according to the energies and intensities, where each block represents a distinct layer. Discretizing the spots into blocks can be an iterative process that is performed repeatedly until a clinical goal or optimization goal is achieved. The optimized plan 408 is then used to configure treatment system 300 for performing proton therapy treatment on patient 304 using gantry 302, for example.

[0050] With regard to FIG. 5, an exemplary sequence of computer implemented steps for automatically generating a proton therapy treatment plan is depicted according to embodiments of the present invention. The treatment plan receives treatment planning parameters at step 501, and the parameters may include patient anatomy data and machine parameters. The treatment planning parameters are discretized into spots with defined energy levels at step 502. At this point each spot may have a separate energy level. At step 503, the spots are added to a list of spots and the list is sorted according to energy level or depth of the spot. An optimization algorithm is executed at step 504 to produce an optimal distribution of layer energies and spots with a relatively smooth dose distribution. The spots can be distributed at fixed intervals (e.g., 3 mm spacing) or based on variations such as beam energy and/or beam range, and specifically, the distribution can be based on the energy levels of the spots and/or the depth/location of a spot, the spots are assigned intensities, where the intensity of a layer is the sum of the intensities of all spots of the layer. The spots are then discretized into blocks according to the energies, intensities, and the location of the spots, and each block represents a distinct energy layer to be applied by a proton therapy system. Discretizing the spots into blocks can be an iterative process that is performed repeatedly until a clinical goal or optimization goal is achieved.

[0051] At step 505, the optimized treatment plan is output, for example, as a computer readable data file. The data file can be stored in memory and used by a proton therapy treatment system to emit a controllable proton beam according to the treatment plan, using a gantry, for example.

[0052] According to some embodiments, the optimization algorithm also determines optimal beam angles for a clinical goal or optimization goal and the spots are generated and discretized according to the beam angles.

[0053] With regard to FIG. 6, an exemplary proton therapy treatment plan 600 applied to a target volume 602 (e.g., a tumor or organ) surrounded by normal tissue 601 for proton therapy treatment is depicted. The treatment plan is generated using limited optimization based on dose volume objectives for the target volume 602 and organs at risk (e.g., normal tissue 601), and machine parameters. Dose distributions are calculated according to the constraints of beam characteristics of proton beam 604 and a machine specific dose calibration. However, the treatment plan depicted in FIG. 6 may not use the full system/machine capability to make use of the system in the most efficient and reliable manner. Moreover, the application of treatment plan 600 at the machine may fail or not achieve the optimal delivery efficiency as requested by plan objectives during treatment planning due to machine specific capability limitations of combined plan parameters that are not taken into account.

[0054] As depicted in FIG. 6, the treatment plan is constrained because it is limited to fixed energy layers for discretizing spots produced by proton beam 604 for treating the target volume 602. The energy levels of the predefined layers are defined based on limited considerations of the patient anatomy and characteristics of the treatment device, and a scanning pattern using layers 603 is applied by beam 604 based on these characteristics. Furthermore, using this approach, the deepest layer of the target 602 receives the highest intensity, and layers closer to the surface receive significantly lower intensity, making it difficult to optimize for the surface layers.

[0055] Moreover, using the approach depicted in FIG. 6, no less than a certain minimum intensity can be applied to a single spot. Therefore, after optimization, if the intensity of a spot is less than the minimum intensity, the spot must be removed from the plan to generate a plan that is deliverable by the machine. If all spots of an energy layer are below the minimum intensity, the entire layer may need to be removed. A simple removal of spots/layers can change the treatment plan such that the dose distribution is no longer optimal. Moreover, failure to remove these spots or layers from the treatment plan during the planning phase may lead to a failure of the treatment, in which case the treatment must be stopped and restarted after the layer is removed. A more efficient approach involves customizing the number of layers and positioning them as necessary to achieve the clinical objectives, as depicted in FIGS. 7 and 8 below.

[0056] With regard to FIG. 7, an exemplary optimized proton therapy treatment plan 700 applied to a target volume 702 (e.g., a tumor or organ) surrounded by normal tissue 701 for proton therapy treatment is depicted. The treatment plan 700 is initially discretized to layers as if there is infinite number of layers, and an iterative optimization process distributes spots using the layers to achieve the clinical goal of treatment plan 700 and ensures that the plan is actually deliverable by the treatment device (e.g., gantry 302). According to some embodiments, an input number of layers is received as user input (for instance). The optimization algorithm then places the spots in a list sorted by range/energy level, and calculates intensities for each spot (e.g., based on energy level/depth) to determine an resultant (e.g., optimal) distribution of layers for achieving the clinical goal of treatment plan 700 using proton beam 704. As depicted in FIG. 7, the treatment plan 700 includes 6 layers 703, and each layer applies ⅙ of the total intensity (number of protons) of treatment plan 700. In this way, the intensities of the layers are fixed, and the spacing of the layers in unconstrained. Splitting the spots among the layers using the principle of equal layer intensities avoids creating spots and layers that are below the minimum intensity deliverable by the treatment system. Note that the principle of equal layer intensities is used here as an example only and other principles could be used to distribute the spots into energy layers. For example, an arbitrary energy-dependent function could be designed to distribute the intensity of the energy layers unevenly. In general, a treatment plan that is optimized using fewer layers leads to a faster delivery time. Moreover, optimizing the treatment plan for accuracy tends to use more layers to better fit the target volume. The spots can be distributed at fixed intervals (e.g., 3 mm spacing) using proton beam 704, or can be distributed based on variations such as beam energy and/or beam range. In either case, the intensity of a layer is the sum of the intensities of all spots on a given layer.

[0057] According to some embodiments, the distribution of the layers and/or the energy levels thereof are manually defined by a user. The user-defined energy levels can include layers with an intensity equal to 0 such that no protons are delivered at certain energies, and these layers can be automatically removed by the treatment planning system after the optimization process. The layers can be defined by the user at fixed distances, or custom distances, for example.

[0058] Embodiments of the present invention are thus described. While the present invention has been described in particular embodiments, it should be appreciated that the present invention should not be construed as limited by such embodiments, but rather construed according to the following claims.