Generating multiple treatment plans for radiation therapy

Abstract

A method for generating multiple treatment plans for radiation therapy is provided, where each treatment plan specifies weights for multiple geometrically defined fluence elements, each weight defining an amount of radiation fluence to provide a radiation dose to a target volume. The method is performed by a treatment planning system and includes generating a first set of treatment plans, determining a subset of fluence elements based on the first set of treatment plans, and generating a second set of at least two treatment plans, where the treatment plans include weights for only the subset of fluence elements.

Description

[Technical Field]

[0001] The present disclosure relates to the field of radiation therapy, and in particular to generating plans while constraining them to a subset of fluence elements. [Background technology]

[0002] Multicriteria optimization (MCO), also called multiobjective optimization, for radiation therapy planning allows the user to, for example, select a criterion through a set of slider bars, each representing a criterion that affects the dose distribution. Navigated This is a development that makes it possible to change the dose distribution. Navigated The dose is a convex combination of the dose distributions of the pre-computed master plan set (a convex combination is a weighted average where the weights are non-negative and sum to one). Navigated The dose distribution is updated in real time based on the current slider position, resulting in a feasible treatment plan, i.e., one that takes into account all delivery system limitations on the planning parameters. Navigated Navigation is directly deliverable if the dose distribution can be accurately reproduced.

[0003] Although the relationship between physical dose and spot weight is linear, the directly deliverable navigation of scanning ions is not trivial. Linearity is a convex combination of the spot weights of the master plan (the convex coefficients are Navigated (the same as that used for dose distribution) Navigated This means defining a treatment plan that accurately reproduces the dose distribution. Navigated The spot weights must satisfy certain constraints: each weight must be either zero or between a lower and upper limit. These limits may be fixed or may depend on the beam energy. In ion delivery systems supporting continuous scanning, the spot weight limits may also depend on the length of the spot segment. Navigated Spot weights are generally not feasible due to constraints, even if all spot weights in the master plan are feasible.

[0004] The delivery of arc-based photon beam radiation therapy, such as tomotherapy and intensity-modulated arc therapy (VMAT), is also governed by the weights of the fluence elements, and each weight must satisfy the constraint that it must be zero or between a lower and upper bound. As with ion beam therapy, even if all fluence weights in the basic plan are feasible, Navigated Fluence weighting is generally not feasible due to trade-offs.

[0005] Prior art techniques sometimes use post-processing to arrive at a deliverable treatment plan. For example, the number of nanometers is infeasible due to spot weight limitations. Visited The spot weights are rounded to the nearest feasible value. This post-processing reduces the dose distribution of the deliverable plan. Navigated There may be deviations from the dose distribution. Navigated The need to compensate for discrepancies between the dose distribution and that of the deliverable plan can make the treatment planning workflow a time-consuming trial and error process.

[0006] A computer-assisted method for customizing a dose distribution plan is disclosed in US20130304503A1. Starting from the initial plan, the user specifies new dose values ​​for local groups of voxels, which may be a small fraction (less than 5%) of the total volume covered by the initial plan. The objective may be to avoid local overdose in risk regions or local underdose in target regions. The initial plan is then converted into a navigation plan with the specified new dose values ​​in substantially the same manner as the initial plan. Because the initial plan is only locally modified, the initial plan is largely preserved. The convex combination of the initial plan and the navigation plan derived from it is visualized together with an input means that allows the user to change the weight of each plan. Summary of the Invention [Problem to be solved by the invention]

[0007] One goal is to improve methods for achieving deliverable therapeutic regimens. [Means for solving the problem]

[0008] According to a first aspect, there is provided a method for generating multiple treatment plans for radiation therapy, each treatment plan specifying weights for multiple geometrically defined fluence elements, each weight defining an amount of radiation fluence to provide a radiation dose to a target volume. The method is performed by a treatment planning system and includes generating a first set of treatment plans, determining a subset of fluence elements based on the first set of treatment plans, and generating a second set of at least two treatment plans, the treatment plans including only weights for the subset of fluence elements.

[0009] Each non-zero weight of the fluence elements of the second set of the treatment plan may be greater than or equal to the minimum weight.

[0010] Generating the second set of two treatment plans may include applying the constraint that fluence elements outside the subset shall be zero.

[0011] The first set of treatment plans may be the result of optimization with respect to a first multi-criteria optimization problem, and the second set of treatment plans may be the result of optimization with respect to a second multi-criteria optimization problem, which may differ from the first multi-criteria optimization problem by the constraint that fluence elements outside of a subset shall be zero.

[0012] The method comprises interpolating dose distributions associated with a second set of treatment plans. Navigated The method may further include using the second set of treatment plans with an operator navigation system, which includes calculating a dose distribution.

[0013] Using the second set of treatment plans with the operator navigation system includes: Navigatedproviding a graphical user interface for visualizing the dose distribution; and a navigation control interface, the navigation control interface allowing an operator to Navigated It allows for the dose distribution to be adjusted.

[0014] The step of determining the subset of fluence elements includes discarding fluence elements having a statistical measure less than a threshold weight, the statistical measure being calculated for each fluence element of all treatment plans in the first set of treatment plans.

[0015] Statistical measures may include means and / or percentiles.

[0016] The step of determining the subset of fluence elements can include ensuring a sufficient density of fluence elements throughout the target volume, the effect of which is to ensure a sufficient dose to the entire target volume, such as a tumor, which can correspond to the complete elimination of all colonizing tumor cells.

[0017] Each treatment plan may be configured to be delivered using a scanned ion beam, with each fluence element associated with a scanning spot of the beam, the scanning spot being defined by a scanning position of the beam and a beam energy.

[0018] Each treatment plan may be configured to be delivered using a radiation beam collimated by a binary multi-leaf collimator (MLC), where each leaf of the MLC can be alternately switched between open and closed positions, and where each fluence element is associated with a particular leaf of the MLC at a particular direction of incidence of the beam relative to the target volume.

[0019] Each treatment plan may be configured to be delivered using a radiation beam collimated by an MLC, where the leaves of the MLC are arranged in opposing leaf pairs, and each leaf can assume any one of several positions between a minimum and maximum position. Such an MLC is hereinafter referred to as a "continuous MLC." In this case, each fluence element is associated with a vixel, and each vixel is a surface element in the cross section of the beam at a particular direction of incidence to the target volume.

[0020] The direction of incidence of each radiation beam relative to the target volume may be determined by either or both of the rotating gantry and the movable couch.

[0021] Each treatment plan may be configured to be delivered with the incident direction of the radiation beam relative to the target volume changed over the course of delivery.

[0022] According to a second aspect, there is provided a treatment planning system for generating multiple treatment plans for radiation therapy, each treatment plan specifying weights for multiple geometrically defined fluence elements, each weight defining an amount of radiation fluence to provide a radiation dose to a target volume. The treatment planning system includes a processor and a memory, the memory storing instructions that, when executed by the processor, cause the treatment planning system to generate a first set of treatment plans, determine a subset of the fluence elements based on the first set of treatment plans, and generate a second set of at least two treatment plans, the treatment plans including only weights for the subset of fluence elements.

[0023] According to a third aspect, there is provided a computer program for generating multiple treatment plans for radiation therapy, each treatment plan specifying weights for a plurality of geometrically defined fluence elements, each weight defining an amount of radiation fluence to provide a radiation dose to a target volume, the computer program comprising: computer program code that, when executed on a treatment planning system, causes the treatment planning system to generate a first set of treatment plans, determine a subset of fluence elements based on the first set of treatment plans, and generate a second set of at least two treatment plans, wherein the treatment plans include weights for only the subset of fluence elements.

[0024] According to a fourth aspect, there is provided a computer program product comprising a computer program according to the third aspect and computer readable means on which the computer program is stored.

[0025] In general, all terms used in the claims should be interpreted according to their ordinary meaning in the art, unless expressly defined otherwise herein. All references to "elements, apparatus, components, means, steps, etc." should be interpreted broadly as referring to at least one example of an element, apparatus, component, means, step, etc., unless expressly stated otherwise. The steps of any method disclosed herein do not have to be performed in the exact order disclosed, unless expressly stated otherwise.

[0026] Aspects and embodiments will now be described, by way of example, with reference to the accompanying drawings, in which: [Brief explanation of the drawings]

[0027] [Figure 1] FIG. 1 is a schematic diagram illustrating an environment in which the embodiments presented herein can be applied. [Figure 2] FIG. 2 is a schematic diagram illustrating functional modules of the treatment planning system of FIG. 1, according to one embodiment. [Figure 3]2 is a schematic diagram illustrating the location of the Bragg peaks of different energy layers in the target volume of FIG. 1. FIG. [Figure 4] 4 is a schematic diagram illustrating the lateral distribution of spots in one of the energy layers of FIG. 3, according to one embodiment. [Figure 5] 1 is a schematic perspective view of a treatment machine illustrating a radiation delivery system employing a continuous MLC. [Figure 6] FIG. 6 is a schematic diagram illustrating the MLC of FIG. 5. [Figure 7] 2 is a flowchart illustrating an embodiment of a method for generating multiple radiation therapy treatment plans performed by the treatment planning system of FIG. 1. [Figure 8] FIG. 2 is a schematic diagram illustrating components of the treatment planning system of FIG. 1, according to one embodiment. [Figure 9] FIG. 1 illustrates an example of a computer program product including computer readable means. DETAILED DESCRIPTION OF THE INVENTION

[0028] Aspects of the present disclosure will now be described more fully hereinafter with reference to the accompanying figures, in which specific embodiments of the invention are shown. However, these aspects may be embodied in many different forms and should not be construed as limiting; rather, these embodiments are provided as examples so that this disclosure will be thorough and complete, and will fully convey the scope of all aspects of the present invention to those skilled in the art. Like numbers refer to like elements throughout the description.

[0029] FIG. 1 is a schematic diagram illustrating an environment in which embodiments presented herein can be applied. A treatment planning system 1 determines how radiation is delivered to a target volume 3. More specifically, the treatment planning system provides a treatment plan 12 to a radiation delivery system 2. The treatment plan 12 specifies weights for multiple geometrically defined fluence elements. Each weight defines the amount of radiation fluence that provides a radiation dose to the target volume 3. Organs-at-risk 5 may be located near the target volume 3. In that case, the treatment plan is determined to balance sufficient dose delivery to the target volume 3 while maintaining low dose delivery to the organs-at-risk 5.

[0030] The manner in which the radiation delivery system 2 generates the beam and delivers the dose varies depending on the treatment modality (e.g., photon, electron, or ion) and geometry, as is well known in the art, but the common goal is to deliver a dose to the target volume 3 (i.e., tumor) as close as possible to the prescribed dose while minimizing the dose to organs at risk 5 depending on where the tumor is located.

[0031] In an ion beam embodiment, described in more detail below with reference to FIGS. 3 and 4, the treatment plan is delivered using a scanned ion beam. In this case, each fluence element is associated with a scanning spot of the beam. The scanning spot is defined by the lateral scanning position of the beam and the beam energy. In an ion delivery system supporting continuous scanning, the spot fluence element is defined as the fluence delivered between two scanning positions. The treatment plan consists of a set of energy layers, each layer having a distribution of scanning spots of ion beam therapy. This is communicated to the ion beam system as a treatment plan 12. Based on the treatment plan 12, the ion beam system generates an ion beam 7 that is scanned spot-by-spot over the patient's target volume 3. Each scanning spot generates a spot dose distribution in the patient's target volume 3. In the coordinate system shown in FIG. 1, depth is represented along the z-axis, and the y-axis is upward. Therefore, the view in FIG. 1 can be considered a side view. The location of the maximum dose (Bragg peak) of the spot dose distribution in the depth direction, i.e., along the z-axis, is controlled by the kinetic energy of the ions, with higher energies resulting in greater doses at greater depths. Furthermore, the lateral position along the y- and x-axes (not shown in FIG. 1) is controlled by deflecting the beam 7 using electromagnets. In this way, a scanning spot can be provided to achieve a dose distribution that covers the target volume 3 in three dimensions.

[0032] In an embodiment of arc-based photon beam radiation therapy, described in more detail below with reference to FIGS. 5 and 6, the treatment plan is delivered during rotational movement of the machine gantry and / or patient couch. Additionally, the patient couch can translate during delivery. The positions of the machine gantry and patient couch define the incidence direction at any given time. In one embodiment, the treatment machine is equipped with a binary MLC. The configuration of the binary MLC can be adjusted during movement by fully opening or fully closing each leaf, i.e., by a binary configuration. Each treatment plan is then delivered using a radiation beam collimated by the binary MLC, such that each leaf of the MLC can be alternately switched between an open and a closed position. In this embodiment, each fluence element is associated with a specific leaf of the MLC at a specific incidence direction of the beam relative to the target volume. In another embodiment, instead of a binary MLC, the MLC can be a continuous MLC, whose leaves can occupy any one of several positions between a maximum position (e.g., fully open) and a minimum position (e.g., fully closed). In this embodiment, each fluence element is associated with a vixel, which is a surface element in the cross section of the beam at a particular incidence direction relative to the target volume. As will be explained in more detail below, each incidence direction of the radiation beam relative to the target volume is determined by either or both of a rotating machine gantry and a movable patient couch.

[0033] FIG. 2 is a schematic diagram illustrating functional modules of the treatment planning system of FIG. 1, according to one embodiment.

[0034] The treatment planning system comprises an optimization module 10 and a navigation module 11. Each of these modules 10, 11 can be implemented in software.

[0035] The optimization module 10 generates several base plans optimized with respect to various criteria such as total dose, local dose, minimum / maximum dose, dose in radiosensitive tissues, number of projections, etc. According to the embodiment presented below, the base plans correspond to a second set of treatment plans.

[0036] The navigation module 11 allows the user to navigate through a set of slider bars. Navigated The dose distribution can be modified, with each slider bar representing a criterion that affects the dose distribution. As is well known per se in the field of graphical user interfaces (GUIs), a slider bar is an element that allows setting, modifying, and / or displaying the current value of a scalar quantity. Furthermore, it is known that a slider bar is only one example of an element with this functionality, and the scope of the present invention is not limited to slider bars, but encompasses any equivalent GUI element. Navigated The dose distribution is a convex combination of the dose distributions of a set of master plans, each of which addresses one of the criteria with special emphasis, which may correspond to one of the MCO's objective functions. Navigated The dose distribution is updated in real time based on the current slider position. The update may involve recalculating the convex combination, but typically does not require a new solution of the MCO. Each slider is associated with a criterion, and increasing the slider corresponds to giving a greater weight in the convex combination to master plans that place particular emphasis on that criterion. According to the embodiments presented herein, Navigated The dose distribution can be directly delivered by the radiation delivery system.

[0037] FIG. 3 is a schematic diagram illustrating the Bragg peak positions of energy layers in the target volume 3 of FIG. 1 when radiation is delivered using a scanned ion beam. FIG. 3 is a side view from the same perspective as the view of FIG. 1. As previously mentioned, the depth of the Bragg peak depends on the energy level. Here, the depths of the Bragg peaks for four energy levels 17a-17d are shown in the target volume 3. The first energy level 17a is illustrated by the line at which the Bragg peak occurs when ions of a first energy amount with different transverse deflections are delivered using the ion beam therapy system of FIG. 1. The second energy level 17b is illustrated by the line at which the Bragg peak occurs when ions of a second energy amount are delivered, and so on. Note that the density of the tissue through which the ion beam passes affects the depth. For example, if the beam passes through bone, the Bragg peak depth will be different from if the beam passes only through soft tissue. As a result, the depth of the Bragg peak for each energy level 17a-17d does not need to be a straight line at a particular depth.

[0038] 4 is a schematic diagram illustrating the lateral distribution of scanning spots in one of the energy layers (see 17a-17d) of FIG. 3, according to one embodiment. The energy layer is shown along the xy plane. Although the energy layer need not be perfectly flat in the target volume of the patient, the energy layer is shown here as a flat layer.

[0039] The scanning spots 14, illustrated as circles in FIG. 4, are located across the entire energy layer to cover the target volume 3 at that energy layer. The weights that can be applied can vary between scanning spots, for example, by controlling the scan time at a particular spot. The weight of each scanning spot must be equal to or greater than the minimum spot weight, which may depend on how quickly the kicker magnet can open and close the beam line of the ion beam delivery system.

[0040] FIG. 4 only discloses the distribution of scanning spots in one energy layer; there is a corresponding distribution of scanning spots for each energy layer used in the target volume.

[0041] 5 is a schematic perspective view of a treatment machine illustrating a radiation delivery system for arc-based radiation therapy. Also shown is a Cartesian coordinate system with dimensions x, y, and z. Note that this coordinate system is different from the coordinate systems of FIGS. 1, 3, and 4.

[0042] The gantry 31 is rotatable about a gantry axis, here parallel to the z-axis. A gantry angle 36 defines the range of rotation of the gantry. It is not important from where the gantry angle 36 is defined, as long as this definition is consistent.

[0043] A couch 30 is provided on which the patient (not shown) lies during treatment. Various fixation mechanisms known per se can be applied to securely fix the patient and the treatment volume in a known position. The couch 30 is rotatable about a couch axis, here parallel to the y-axis. A couch angle 35 defines the range of couch rotation. It is not important from where the couch angle 35 is defined, as long as this definition is consistent. Furthermore, the couch may be movable in the z-direction.

[0044] An MLC 33 is provided mounted on the gantry 31 through which radiation is delivered during treatment. The MLC 33 may be rotatable about a collimator axis, which changes its orientation (in a Cartesian coordinate system) as the gantry 31 rotates. A collimator angle 37 defines the range of rotation of the MLC. It is not important from where the collimator angle 37 is defined, as long as this definition is consistent.

[0045] The combination of values ​​for couch angle 35, gantry angle 36, and optionally collimator angle 37 defines the incidence direction. The incidence direction defines the angle at which the radiation will treat the patient. The beam plane is the plane normal to the beam direction, i.e., the collimator axis.

[0046] Each trajectory occurs in an arc from a start time to an end time and defines the movement between the incident directions. In one embodiment, for helical tomotherapy, the movement is helical, where translation between the couch and gantry is allowed along the z-axis.

[0047] Tomotherapy is a form of photon beam therapy in which a patient is irradiated with a slit beam that continuously rotates around the patient. The rotation is discretized into several projections for planning purposes (typically 51 projections per rotation). The slit width is defined by a pair of movable jaws, typically 0.5–5 cm, and the radiation passing through the slit is collimated by a set of pneumatically driven MLC leaves. The collimation provided by the leaves can be binary, meaning that the leaves can only be fully open or fully closed. In the embodiment presented herein, the weight of a fluence element corresponds to the open time of the MLC leaves at a particular incident direction. The weight, i.e., the leaf open time, must be greater than or equal to a lower bound for all open MLC leaves. The lower bound corresponds to the minimum leaf open time, which may depend on the shortest possible time that the leaves can be in the open configuration due to the finite leaf velocity of the MLC. A weight of zero, corresponding to a closed leaf, is also possible. The weight between successive projections must also satisfy a minimum (non-zero) closed time constraint. A leaf does not need to be closed between two successive projections, so a zero closure time is possible.

[0048] In one embodiment, the motion is performed by VMAT. The arcuate trajectory defines the motion, which is implemented using one or more changes in the couch angle 35, collimator angle 37, and gantry angle 36. The MLC here is in a continuous configuration, and each leaf of the MLC can assume any one of several positions between a fully open position and a fully closed position. The leaf positions of a continuous MLC are typically controlled by mechanical motors. The leaves are also arranged in opposing leaf pairs. In this case, each fluence element is associated with a vixel, which is a surface element in the cross section of the beam (i.e., the beam plane) at a specific incident direction to the target volume. The surface element may correspond to the smallest controllable unit of the MLC, such as a leaf or leaf pair. In the embodiment presented herein, the weight of the fluence element, i.e., the vixel weight, corresponds to the amount of radiation fluence delivered while the vixel is not blocked by the MLC leaf. The vixel weight must be greater than or equal to a lower limit, which may depend on the minimum chip gap between opposing leaves and a finite maximum leaf speed. A pixel weight of zero, corresponding to pixels that are always occluded by the MLC leaf, may also be feasible.

[0049] In one embodiment, the radiation is turned on for the entire duration of each arcuate trajectory. The speed of movement during the arcuate trajectory can be constant or can vary.

[0050] Figure 6 is a schematic diagram illustrating the MLC 33 of Figure 5 when applied in a serial configuration. The MLC 33 comprises pairs of leaves 20a-20b, 21a-21b, ..., 26a-26b. Each leaf is movable in only one dimension.

[0051] Each pair of opposing leaves can be positioned with a space between them. In this way, an aperture 28 can be defined through which radiation can pass. The aperture 28 can be adjusted to cover the target volume 3 while reducing radiation to surrounding tissue. Because the leaves 20a-20b, 21a-21b, ..., 26a-26b are movable only along a single dimension, the possible shapes of the aperture 28 depend on the rotation angle 37 of the MLC 33. When the MLC 33 is in a binary configuration, each leaf can only remain in a fully open or fully closed position. In a binary configuration, there can be pairs of leaves similar to those shown in FIG. 6, or there can be only one leaf per configurable aperture, corresponding to each vertical position in FIG. 6.

[0052] 7 is a flow chart illustrating an embodiment of a method for generating multiple radiation therapy treatment plans. The method is performed by the treatment planning system of FIG. 1. Each treatment plan specifies weights for multiple geometrically defined fluence elements. Further, each weight defines an amount of radiation fluence that provides a radiation dose to a target volume.

[0053] In step 40 of generating a first set of treatment plans, the treatment planning system generates the first set of treatment plans. In generating the first set of treatment plans, these are generated with respect to an initial set of fluence elements, which are selected equally for all plans in the first set of treatment plans. Furthermore, the plans are generated without consideration of the (non-zero) minimum weights of the fluence elements. Nevertheless, consideration of the maximum weights can be included in generating the first set of treatment plans. The first set of treatment plans is the result of optimization with respect to a first multi-criteria optimization problem, and several different objective functions may be used (sequentially). Each treatment plan in the first set of treatment plans may express a particular optimization criterion or a particular weighting of more than one criterion. A particular optimization criterion may be expressed in terms of the objective function used to generate that treatment plan and / or in terms of the constraints applied.

[0054] In a fluence element determination step 42, the treatment planning system determines a subset of fluence elements based on the first set of treatment plans. This may include discarding fluence elements having a statistical measure below a threshold weight. The statistical measure is calculated for each fluence element of all treatment plans in the first set of treatment plans. For example, the statistical measure may include an average value or a percentile value. In this way, for example, fluence elements with too small a weight (measured as an average value or a percentile value) in the first set of treatment plans may be discarded because their contribution to the treatment plan is too small, potentially below the minimum deliverable weight.

[0055] Optionally, this step includes ensuring that there is a sufficient density of fluence elements throughout the target volume. Such optional assurance can be applied after the comparison with the threshold weights but before decision-making, thereby excluding from discard those fluence elements that are important for maintaining sufficient density. Thus, fluence elements that would have been discarded if the aforementioned comparison with the threshold weights had been applied exclusively are effectively not discarded, given the optional further goal of maintaining a set of fluence elements that provides sufficient coverage of the target volume.

[0056] A possible process for implementing step 42 of determining the fluence element will now be described. In a first substep, a reference coverage is calculated for each treatment plan of the first set of treatment plans or for all treatment plans of the first set of treatment plans together. The reference coverage is the volume-at-dose, i.e., a specific area V of the target volume V that receives the minimum dose d. d Volume μ(V d ) can be calculated as:

number

[0057] In a second substep, the treatment planning system determines a subset F of fluence elements as described above.

[0058] In a third substep of step 42, it is determined whether the reference coverage can be achieved using the subset F of fluence elements. F V outside of d If there are too many points or voxels in region V, then reference coverage is not considered achievable. d can be assessed by applying the following quantitative criteria, which require that at least a portion of γ remains irradiatable:

number

[0059] Further, under the third sub-step of step 42, region V d An alternative response to the determination that some points or voxels of region V cannot be delivered with a subset of fluence elements is as follows: the treatment planning system determines whether some points or voxels of region V are delivered with a subset of fluence elements that were used in some of the first set of treatment plans and that would be delivered with a subset of fluence elements ... d, and search for additional fluence elements that will allow re-irradiation of a large portion of the target volume. The search may be constrained to fluence elements on or adjacent to the boundary of the subset of fluence elements. Such a boundary may be a point set (or a discretized point set) of a two-dimensional representation of the fluence elements. In that case, the final output of step 42 may include one or more additional fluence elements from this set. In other words, these additional fluence elements are not effectively discarded when the optional additional objective of maintaining a set of fluence elements that provides sufficient coverage of the target volume is taken into account. In summary, the behavior of this optional implementation of step 42, which involves ensuring a sufficient density of fluence elements throughout the target volume, is primarily controlled by the values ​​assigned to the parameters d and γ and, if applicable, by the method for calculating the dose-receiving volume collectively.

[0060] In a generate second set of treatment plans step 44, the treatment planning system generates a second set of at least two treatment plans. These treatment plans include weights for only a subset of fluence elements, with the weights for each fluence element in the subset of fluence elements being constrained to satisfy minimum and maximum weight requirements. Each treatment plan in the second set may correspond to a treatment plan in the first set. Such correspondence may include using the same objective function or using only an objective function modified to include the minimum weight requirement. The second set of treatment plans may be the result of optimization with respect to a second multi-criteria optimization problem.

[0061] The second multi-criteria optimization problem may differ from the first multi-criteria optimization problem in that it includes a constraint that fluence elements outside the subsets shall be zero. As is well known to those skilled in the art, such a constraint may be expressed in terms of modifying the objective function, for example, by adding a barrier function term. In one embodiment, the zero constraint is applied to all fluence elements in the target volume except for the subsets. In another embodiment, the zero constraint is applied to fluence elements that (i) have a non-zero weight in any of the first set of treatment plans and (ii) are outside the determined subsets.

[0062] In one embodiment, each non-zero weight of the fluence elements of the second set of treatment plans is greater than or equal to a minimum weight and less than or equal to a maximum weight. The minimum and maximum weights can be determined based on the physical limitations of the radiation delivery system. For ion beam embodiments, the fluence elements are scanning spots. For tomotherapy embodiments, the fluence elements correspond to specific leaves of the MLC at a specific direction of incidence of the beam relative to the target volume. For VMAT embodiments, the fluence elements correspond to pixels, which are surface elements in a plane perpendicular to a specific direction of incidence of the beam relative to the target volume.

[0063] In a use plans for navigation step 46, the treatment planning system uses the second set of treatment plans in an operator navigation system, such as the navigation module shown in Figure 2 and previously described, by interpolating the dose distributions associated with the second set of treatment plans. Navigated The interpolation can be performed, for example, by forming a convex combination of the dose distributions associated with the second set of treatment plans. Navigated The method may include providing a graphical user interface for visualizing the dose distribution, in which case a navigation control interface is also provided, which allows the operator to navigate, for example using a slider bar. Navigated It allows for the dose distribution to be adjusted.

[0064] The first set of treatment plans can be the result of optimization with respect to a first multi-criteria optimization problem, while the second set of treatment plans can be the result of optimization with respect to a second multi-criteria optimization problem, in other words, the optimization problems can be different for the first and second sets of treatment plans.

[0065] As mentioned above, in one embodiment, each treatment plan is configured to be delivered using a scanned ion beam (see FIGS. 3 and 4). In such a case, each fluence element is associated with a scanning spot of the ion beam.

[0066] Alternatively, each treatment plan can be configured to be delivered using a radiation beam collimated by an MLC, which can be in the form of a binary or continuous configuration. Each treatment plan can be configured to be delivered with the direction of incidence of the radiation beam relative to the target volume changing over the course of delivery.

[0067] Using the embodiments presented herein, a convex combination of treatment plans from the second set is directly deliverable, i.e., can be used directly by the treatment machine, which simplifies clinical decision-making. Furthermore, post-processing after MCO navigation is not required to make the combination of plans deliverable. This is extremely valuable because post-processing can be error-prone and time-consuming to perform.

[0068] 8 is a schematic diagram illustrating components of the treatment planning system of FIG. 1 , according to one embodiment. Processor 60 may be provided using any combination of one or more suitable central processing units (CPUs), multiprocessors, microcontrollers, digital signal processors (DSPs), application specific integrated circuits, etc., capable of executing software instructions 67 stored in memory 64, and thus may be a computer program product. Processor 60 may be configured to perform the method described with reference to FIG. 7 above.

[0069] The memory 64 may be any combination of random access memory (RAM) and read-only memory (ROM). The memory 64 also includes persistent storage, which may be, for example, any single or combination of magnetic memory, optical memory, solid-state memory, or even remotely located memory.

[0070] A data memory 66 is also provided for reading and / or storing data during execution of software instructions in processor 60. Data memory 66 may be any combination of random access memory (RAM) and read only memory (ROM).

[0071] The treatment planning system 1 further comprises an I / O interface 62 for communicating with other external entities. Optionally, the I / O interface 62 also includes a user interface.

[0072] Other components of the treatment planning system 1 have been omitted so as not to obscure the concepts presented herein.

[0073] FIG. 9 illustrates an example of a computer program product including computer-readable means. The computer-readable means may store a computer program 91, which may cause a processor to perform a method according to embodiments described herein. In this example, the computer program product is an optical disc, such as a CD (compact disc), a DVD (digital versatile disc), or a Blu-ray disc. As noted above, the computer program product may also be embodied in the memory of a device, such as computer program product 64 of FIG. 8. While computer program 91 is shown here schematically as a track on the illustrated optical disc, the computer program may be stored in any manner suitable for a computer program product, such as removable solid-state memory, e.g., a Universal Serial Bus (USB) drive.

[0074] Aspects of the present disclosure have been described above primarily with reference to certain embodiments. However, as will be readily apparent to those skilled in the art, other embodiments besides those disclosed above are equally possible within the scope of the present invention, as defined by the appended claims. Thus, while various aspects and embodiments have been disclosed herein, other aspects and embodiments will be apparent to those skilled in the art. The various aspects and embodiments disclosed herein are for purposes of illustration and not limitation, with the true scope and spirit being indicated by the following claims.

Claims

1. A method for calculating a dose distribution based on a plurality of treatment plans for radiation therapy, each treatment plan specifying weights for a plurality of geometrically defined fluence elements, each weight defining an amount of radiation fluence that provides a radiation dose to a target volume, the method being performed by a treatment planning system, generating, by the treatment planning system, a first treatment plan set including the plurality of treatment plans, the first treatment plan set including the initial set of fluence elements; the treatment planning system determining a subset based on the first treatment plan set by discarding the fluence elements having a statistical measure below a threshold weight; generating a second set of treatment plans including at least two treatment plans that include the fluence elements in the subset that have weights equal to or greater than the threshold; calculating a dose distribution by mutually interpolating each treatment plan in the second treatment plan set in an operator navigation system, wherein the treatment planning system includes a graphical user interface and a navigation control interface; A method comprising:

2. The method of claim 1 , wherein non-zero weights of fluence elements in the second treatment plan set are greater than or equal to the threshold weight.

3. 10. The method of claim 1, wherein generating the second treatment plan set includes applying a constraint that fluence elements not included in the determined subset of fluence elements shall be zero.

4. 2. The method of claim 1, wherein a first multi-criteria optimization problem used to optimize the first treatment plan set based on a dose distribution associated with a dose of radiation fluence defined by the first treatment plan set and a second multi-criteria optimization problem used to optimize the second treatment plan set based on a dose distribution associated with a dose of radiation fluence defined by the second treatment plan set are different.

5. 5. The method of claim 4, wherein the second multi-criteria optimization problem differs from the first multi-criteria optimization problem by the constraint that fluence elements not included in the determined subset of fluence elements shall be zero.

6. The method of claim 1 , wherein the graphical user interface is capable of visualizing a dose distribution of the navigation, and the navigation control interface allows an operator to adjust the dose distribution of the navigation.

7. The method of claim 1 , wherein the statistical measure is calculated for a weight of each fluence element of all treatment plans in the first set of treatment plans.

8. The method of claim 7 , wherein the statistical measure comprises a mean or percentile.

9. The method of claim 1 , wherein determining the subset of fluence elements comprises ensuring that there is a sufficient density of geometrically defined fluence elements in the target volume.

10. 10. The method of claim 1, wherein each treatment plan is configured to deliver a radiation fluence defined by each treatment plan to the target volume using a scanned ion beam, each fluence element being associated with a scanning spot of the beam, the scanning spot being defined by a scanning position of the beam and a beam energy.

11. 2. The method of claim 1, wherein each treatment plan is configured to deliver a radiation fluence defined by each treatment plan to the target volume using a radiation beam collimated by a binary multi-leaf collimator MLC, each leaf of the MLC being alternately switchable between an open position and a closed position, and each fluence element is associated with a particular leaf of the MLC at a particular incident direction of the beam relative to the target volume.

12. 12. The method of claim 10 or 11, wherein each incident direction of the radiation beam with respect to the target volume is determined based on the position of either or both of a rotating gantry and a movable couch.

13. 12. The method of claim 11, wherein each treatment plan is configured to deliver a radiation fluence defined by each treatment plan to the target volume, with the incident direction of the radiation beam relative to the target volume being changed over the course of delivery.

14. 2. The method of claim 1, wherein each treatment plan is configured to deliver a radiation fluence defined by each treatment plan to the target volume using a radiation beam collimated by a binary multi-leaf collimator MLC, the leaves of the MLC being arranged in opposing leaf pairs, each leaf being capable of occupying any one of a number of positions between a minimum and a maximum position, each fluence element being associated with a vixel, each vixel being a surface element in a cross section of the beam at a particular direction of incidence relative to the target volume.

15. A treatment planning system for calculating a dose distribution based on a plurality of treatment plans for radiation therapy, each treatment plan specifying weights for a plurality of geometrically defined fluence elements, each weight defining an amount of radiation fluence that provides a radiation dose to a target volume; a processor; a memory for storing instructions; the instructions, when executed by the processor, cause the treatment planning system to: generating a first treatment plan set including the plurality of treatment plans, the first treatment plan set including the initial set of fluence elements; determining a subset based on the first treatment plan set by discarding fluence elements having a statistical measure below a threshold weight; generating a second set of treatment plans including at least two treatment plans that include the fluence elements in the subset having weights equal to or greater than the threshold; calculating a dose distribution by mutually interpolating each treatment plan in the second treatment plan set in an operator navigation system having a graphical user interface and a navigation control interface; Treatment planning system.

16. 1. A computer-readable storage medium having stored thereon a computer program for calculating a dose distribution based on a plurality of radiation therapy treatment plans, each treatment plan specifying weights for a plurality of geometrically defined fluence elements, each weight defining an amount of radiation fluence that provides a radiation dose to a target volume, the computer program, when executed on a treatment planning system, causing the treatment planning system to: generating a first treatment plan set including the plurality of treatment plans, the first treatment plan set including the initial set of fluence elements; determining a subset based on a first treatment plan set by discarding the fluence elements having a statistical measure below a threshold weight; generating a second set of treatment plans including at least two treatment plans that include the fluence elements in the subset having weights equal to or greater than the threshold; calculating a dose distribution by mutually interpolating each treatment plan in the second treatment plan set in an operator navigation system having a graphical user interface and a navigation control interface; A computer-readable storage medium that causes the computer to perform the following:

Images