Computer-implemented method for radiation treatment planning, computer program product and computer system for performing the method

Abstract

A computer-based method for optimizing a radiation treatment plan for a patient is proposed, wherein the complete treatment, including both external beam radiation therapy and brachytherapy, is optimized in one step using an optimization problem with an objective function designed to optimize the total dose distribution as a combination of a first dose distribution delivered by a first radiation set and a second dose distribution delivered by a second radiation set, one of which is external beam radiation therapy and the other is brachytherapy. The optimization is based on the desired total dose for the entire treatment, images of the patient before treatment begins, and estimated images of the patient after the first radiation set is delivered.

Description

[Technical Field]

[0001] The present invention relates to computer-implemented methods, computer program products and apparatus for performing such methods for radiation treatment planning, and in particular to treatment planning for radiation therapy that delivers both external beam radiation therapy and brachytherapy for the same patient. [Background technology]

[0002] Most radiation treatments are given as external radiation delivered to the patient from an external radiation source, known as external beam radiation therapy (i.e., EBRT). This is usually delivered in many fractions, e.g., 15, 30, or more. Alternatively, radiation may be delivered from a radiation source placed inside the patient's body. This is known as brachytherapy, and involves placing one or more needles or other types of devices into a target inside the patient's body to expose it to radiation from within. This is usually done in fewer fractions, e.g., one or three.

[0003] For EBRT treatments, planning is currently primarily performed as an inverse treatment planning process in an optimization treatment planning system. The goal is usually to achieve a minimal or uniform dose throughout the target. Brachytherapy (i.e., BT) plans require the placement of a device to deliver radiation within the target, either in the body or in a virtual environment, and the dose is usually developed in a forward planning process, where the device's location and resulting dose are determined. In low-dose-rate BT, one or more radiation sources, known as seeds, are implanted into the target and typically remain there for a period of time. On the other hand, in high-dose-rate and pulsed-dose-rate BT, radiation is delivered by a radiation source that moves through a hollow channel in an implanted device, such as a needle, catheter, or applicator. Typically, BT dose distribution is not uniform but instead is concentrated around the implanted device.

[0004] Combining EBRT and brachytherapy in the same patient is known in the art. Typically, the EBRT fraction is delivered first, followed by brachytherapy, although the reverse order is also possible. Patient geometry usually changes during the first treatment. Also, the devices inserted into the patient for brachytherapy alter the shape of the target and surrounding tissue. For both of these reasons, the two modalities are delivered to different patient geometries. When both brachytherapy and external beam radiation therapy are used on the same patient, traditionally, two separate treatment plans are developed: one for the EBRT portion of the treatment and one for the brachytherapy portion of the treatment. These are typically based on the total dose and a predetermined division of the dose between the different types of treatment. For example, a total dose of 80 Gy may be set, with the EBRT treatment providing 60 Gy and the brachytherapy providing the remaining 20 Gy. D1 discloses a radiation treatment planning method that allows for the combination of two or more modalities, which may include a combination of brachytherapy and external beam radiation therapy. According to D1, the plan is iteratively optimized, first selecting and optimizing the most promising modality, then adding another modality and optimizing the contribution of this other modality.

[0005] It is an object of the present invention to provide an improved planning method for treatment planning including both EBRT and brachytherapy. Summary of the Invention

[0006] The present invention provides the following a. obtaining input data including a first image of a patient; b. obtaining an optimization problem including an objective function designed to optimize a total dose distribution based on input data as a combination of a first dose distribution provided by a first radiation set and a second dose distribution provided by a second radiation set based on a dose criterion for the total dose distribution, wherein one of the radiation sets is external beam therapy and the other is brachytherapy; c. optimizing the treatment plan as a combination of external beam therapy and brachytherapy using an optimization problem; The present invention relates to a computer-based method for optimizing a radiation treatment plan for a patient, including:

[0007] Thus, in accordance with the present invention, the treatment parameters of the two radiation sets are simultaneously optimized using a joint optimization problem, resulting in a better overall treatment plan than the separate plans of the prior art. The first and second doses may be delivered in any suitable order, or alternately. The input data may also include a second image of the patient, which may be an estimated image of the patient after delivery of a portion of the treatment plan. The second image may then be used to provide information about the patient's anatomy after delivery of the portion. Typically, this portion is the portion delivered by the first radiation set.

[0008] Preferably, the method includes the further steps of determining a dose delivered to the patient after delivery of a first portion of the treatment plan, providing at least one updated image of the patient, and re-optimizing the remaining treatment plan using a second optimization problem taking into account the delivered dose and the updated image, thereby allowing one or more remaining portions of the treatment plan to be refined based on the actual results of delivery of the first portion, where the first portion typically corresponds to a dose delivered by a first radiation set and the remaining portions typically correspond to a dose delivered by a second radiation set.

[0009] In some embodiments, the first portion of the treatment plan is the portion delivered as an external beam therapy. In this case, the at least one updated image preferably includes an updated image of the patient and an updated image of the patient with the brachytherapy device applied to take into account the actual patient geometry including the brachytherapy device during delivery. The remaining treatment plan is updated based on the updated image.

[0010] In other embodiments, the first portion of the treatment plan is the portion to be delivered as brachytherapy, and preferably in such cases, step b is further performed based on a current image of the patient with the brachytherapy device applied, to take into account the actual patient geometry during delivery.

[0011] The at least one updated image may include at least one image of the patient taken after delivery of the first portion, which provides the most accurate information about the actual patient geometry.

[0012] Alternatively or additionally, the at least one updated image may include at least one simulated image based on an estimate of the patient's geometry after delivery of the portion, which may be preferable if for some reason it is not feasible to take a new image of the patient after delivery of the portion or with the brachytherapy device inserted.

[0013] In a preferred embodiment, the optimization includes transforming at least one of the doses delivered by the first and second radiation sets, respectively, to obtain a common geometry for the treatment portion and accumulating them using a biological model, and the objective function is a set of penalties for the accumulated dose and the radiation set-specific dose.

[0014] Preferably, robust planning is used to account for uncertainties in brachytherapy device placement, EBRT delivery, and / or the determined delivered dose. Methods for performing robust planning are known in the art.

[0015] The present invention also relates to a computer program product comprising computer readable code means configured to, when run on a computer, cause the computer to carry out the method according to any of the above embodiments. The computer program product may be stored on any suitable type of non-transitory storage medium.

[0016] The invention also relates to a computer system comprising a processor and at least one program memory, characterized in that the program memory holds a computer program as defined above.

[0017] In a preferred embodiment, the present invention involves planning a brachytherapy plan that also takes into account the cumulative dose already delivered to the patient by several other modalities, which means that deviations from the planned dose to the target can be compensated for, and also excessive doses to organs at risk can be compensated for by altering the brachytherapy treatment plan accordingly.

[0018] The invention will be explained in more detail below, by way of example and with reference to the accompanying drawings, in which: [Brief explanation of the drawings]

[0019] [Figure 1a] 1 is a cross-section of a medical image of a patient. [Figure 1b] 1 is a cross-section of a medical image of a patient. [Figure 1c] 1 is a cross-section of a medical image of a patient. [Figure 2] 1 is a flow chart of a general embodiment of the method. [Figure 3] 1 is a flow chart of a more specific embodiment of the method including BT after EBRT. [Figure 4] 1 is a flow chart of a second, more specific embodiment of the method including EBRT after BT. [Figure 5] 1 is a schematic overview of a computer system on which embodiments of the present invention may be implemented. DETAILED DESCRIPTION OF THE INVENTION

[0020] In external beam radiation therapy (i.e., EBRT), radiation is administered to a patient in the form of externally delivered beams. The radiation may be any type of radiation, such as photons, electrons, protons, or other ions. In brachytherapy (i.e., BT), some type of device is inserted into a target inside the patient's body and used to deliver radiation from one or more points within the target. The device may include many small needles and / or catheters, one or more larger applicators, one or more seeds, or any combination of different types of devices. Depending on the number and size of the devices, the geometry of the target and surrounding patient is modified.

[0021] Due to the different nature of the two radiation sets, EBRT and BT plans use different treatment parameters. Treatment parameters for EBRT treatments include beam and beam-limiting device configurations. Treatment parameters for BT treatments include variables such as device position and dwell time. Each radiation set typically requires radiation to be delivered in one or more fractions, a greater number for EBRT than for BT, which can typically, but not necessarily, even be delivered in one single fraction.

[0022] As discussed above, inverse treatment planning using optimization problems is common for EBRT planning, but this has not traditionally been used for brachytherapy.

number

number

[0023] The present invention relates to the simultaneous optimization of treatment parameters for EBRT and BT treatments, which can be solved by solving the optimization problem in equation (3):

number

[0024] Typically, optimization involves transforming doses into a common geometric shape and accumulating them using a biological model, with the objective function being a set of accumulated doses and penalties for radiation set-specific doses.

[0025] FIGS. 1a, 1b, and 1c are simplified examples of medical images taken at different times during a procedure of the present invention, as discussed in more detail below in connection with FIG. 2. FIG. 1a is a cross-section 11 of a conceptual medical image of a patient's abdomen, focusing on a target 13 and organs at risk 15, used for treatment planning according to an embodiment of the present invention. FIG. 1b is a corresponding cross-section 11′ of a medical image of the same patient after a first type of treatment, illustrating schematically the changes that may occur in the patient's geometry due to this treatment. Naturally, the target 13 has shrunk due to the treatment, which is typically the desired result. FIG. 1c is a corresponding cross-section 11″ of a medical image of the same patient with needles inserted into the target to deliver brachytherapy to the patient. These needles are shown as small dots 17 within the target. As shown, this also changes the geometry of the target 13″ and the area of ​​the patient surrounding the target.

[0026] FIG. 2 is a flowchart of an overall method according to one embodiment of the present invention.

[0027] In a first step S21, images of the relevant portion of the patient are obtained, as discussed in connection with Figures 1a-1c. In step S22, dose criteria are determined for the total dose distribution to be delivered as a combined plan including both EBRT and brachytherapy.

[0028] In step S23, an optimization problem is defined based on one or more images and dose criteria for the desired total dose. The dose criteria are established as is common in the art. They typically include a minimum dose for all voxels of the target and often a maximum dose for one or more organs at risk. For example, the dose criteria may specify a total dose of at least 60 Gy in each target voxel, and a maximum of 30% of the organs at risk to receive a total dose of more than 40 Gy. The dose criteria may also include a partial or complete dose distribution. In S24, the optimization problem is used to optimize the treatment plan. The optimization problem includes an objective function, such as function (3) above.

[0029] In step S25, a portion of the treatment plan is delivered to the patient, and in step S26, the cumulative dose delivered to the patient from that portion of the treatment plan is estimated. The cumulative dose may be determined in any suitable manner. Methods for doing this are well known in the art and are typically based on at least one medical image, for example, many fraction images taken throughout the delivery of the first portion of the treatment plan.

[0030] In step S27, a new image of the same portion of the patient is obtained to confirm the new patient geometry after the partial delivery in step S25. If applicable, other modifications, such as the insertion of a brachytherapy device, may be made, and an image may be taken that reflects the resulting geometry. The new image may be an image of the patient taken at this stage or a synthetic image based on an estimate of the patient's new geometry.

[0031] In step S28, the remainder of the treatment plan is re-optimized using an inverse treatment planning method and based on an optimization function, discussed in more detail below, that takes into account the cumulative dose of the previous partial deliveries. The optimization problem must then include an objective function related to:

number

number

number

[0032] where g is another objective function, which may or may not be equal to f, and the dose delivered from the first radiation set over all fractions (measured or estimated) is used as a fixed background dose for planning for the second radiation set.

[0033] In step S28, the reoptimized remainder of the plan is delivered to the patient.

[0034] 3 is a flowchart of one embodiment of the method of the present invention, in which the overall treatment plan includes first an EBRT treatment and then a BT treatment. Input data S31 to this plan includes current medical images of the patient, dose criteria for the desired dose distribution, and a predicted model of the patient's geometry after the EBRT treatment. The predicted model may be a calibrated or synthetic medical image. The input data may also include a predicted model of the patient's geometry after the EBRT treatment with the EB device included. The medical images may be CT images or any other suitable modality images, such as MR or ultrasound images.

[0035] In method step S32, an optimization problem is obtained based on the dose criteria and input data. The optimization problem includes an objective function according to equation (3) based on the total dose of both the EBRT and BT fractions of the treatment and, optionally, the dose delivered for each radiation set. This is typically the dose d EBRT and d BT The method includes transforming at least one of the images into a common geometry and accumulating them using a suitable biological model. This may include establishing a combination of accumulated dose and a penalty for a particular dose. Models for establishing the common geometry, typically involving non-rigid registration of the images, are known. Models for determining accumulated dose are also known to those skilled in the art. For example, the biological concept EQD2 may be applied to provide an estimate of the total effective dose.

[0036] Subsequent method step S33 involves optimization based on a common geometry and accumulated dose. The optimization problem includes an objective function according to equation (3) above. Of course, the objective function can be extended to depend on treatment parameters, and the optimization problem may also include constraints depending on the dose or treatment parameters. The output S34 from optimization step S33 is an overall treatment plan that includes one fraction for each radiation set: one EBRT fraction and one BT fraction. Each fraction includes the portion of the dose delivered by the corresponding radiation set and the number of fractions for delivering it.

[0037] The EBRT subportion of the treatment plan is then delivered to the patient in step S35, and the actual delivered dose is determined or estimated from this delivery in step S36. Preferably, the situation after EBRT delivery is evaluated and used to refine the BT subportion of the treatment plan as outlined in the following steps.

[0038] In step S37, updated images of the patient after EBRT treatment are obtained. These include new images of the patient to compensate for geometric changes that occurred during EBRT treatment. The images also include images of the patient with the BT device inserted, since the BT device causes some deformation of the target and surrounding patient geometry depending on the type of device.

[0039] In a subsequent step S38, the BT fraction of the treatment is re-optimized, taking into account the delivered dose from step S35 and the new image obtained in step S37. The optimization problem in this case is the following equation (2):

number

number

[0040] The output from step S38 is a new optimized BT treatment plan S39, which is preferably delivered to the patient.

[0041] As discussed with respect to FIG. 2, steps S36-S39 may be performed without first performing the preceding steps, i.e., the BT plan may take into account prior EBRT treatments even if the two initial combined plans do not exist.

[0042] 4 is a flowchart of a method in which the first radiation set delivered is BT. Input data S41 for this plan includes current medical images of the patient, medical images of the patient with the BT device inserted, and dose criteria for the desired dose distribution. Preferably, this input data also includes a predicted model of the patient's geometry after BT treatment. The medical images may be CT images or any other suitable modality images, such as MR or ultrasound images.

[0043] In method step S42, the optimization problem is defined in a manner similar to step S32. Images of the patient with the BT device inserted when the first BT dose portion is delivered are already available as input. These images can be non-rigidly registered to provide a geometric correspondence between the treatment geometries. As in step S32, the dose d EBRT and d BT are transformed into a common geometry and accumulated using a suitable biological model. This may involve establishing a combination of accumulated dose and a penalty for a particular dose. Models for establishing the common geometry are known, typically involving non-rigid registration of the images. Models for determining accumulated dose are also known to those skilled in the art.

[0044] In a subsequent method step S43, optimization is performed based on a common geometry and deposited dose. The optimization problem includes an objective function according to equation (3) above. Of course, the optimization problem can be extended to depend on treatment parameters and include other objective functions and / or constraints. The output S44 from the optimization step S43 is an overall treatment plan including one sub-fraction for each radiation set, i.e., one BT sub-fraction and one EBRT sub-fraction.

[0045] The BT subportion of the treatment plan is then delivered to the patient in step S45, and the actual delivered dose is determined or estimated from this delivery in step S46. Preferably, the post-BT delivery situation is evaluated and used to refine the EBRT subportion of the treatment plan, as outlined in the steps below.

[0046] In step S47, an updated image of the patient after BT treatment is obtained to compensate for the geometric changes that occurred during the BT treatment.

[0047] Thereafter, in step S48, the EBRT fraction of the treatment is re-optimized taking into account the delivered dose from step S45 and the new images obtained in step S47. The optimization problem in this case is the following equation (2):

number

number

[0048] The output from step S48 is a new optimized EBRT treatment plan S49, which is preferably delivered to the patient.

[0049] As discussed with respect to Figures 2 and 3, steps S46-S49 may be performed without first performing the preceding steps, i.e., the EBRT plan may take into account previously delivered BT treatments even in the absence of these two initial combined plans.

[0050] It is also possible to create a plan in which the BT and EBRT fractions are not given as two consecutive fractions, but instead the BT fraction is distributed between the EBRT fractions. In this type of treatment, one or more fractions of the treatment that have not yet been delivered can be reoptimized taking into account the delivered dose. Both the delivered and non-delivered fractions of the treatment are typically a combination of BT and EBRT. The optimization is then calculated using the following equation (7):

number

[0051] In all of the above methods, joint optimization must be performed with some care to ensure that the individual doses of each radiation set are still satisfactory on their own. One possible adverse effect of joint optimization is that the EBRT dose has cold spots at the target that are later filled by the BT dose. This can be mitigated, for example, by incorporating robustness to device positioning uncertainties and deformation effects into the model. Treatment-specific objective functions are also possible (similar to the current beam-set-specific objective functions).

[0052] As with any radiation therapy plan, there are sources of uncertainty, such as patient placement, BT machine positioning, and estimated delivered dose. Robust planning may be used to compensate for this. In particular, deformations between images result in an approximate deposited dose, the quality of which depends on the accuracy of the non-rigid registration. To avoid over-optimization for a deposited dose that differs from the actual delivered dose, methods for robust planning can be used that go beyond the indication of uncertainty. Various degrees of refinement can be used, for example: Margins as ITV across the predicted image or just a smearing of the area to be treated can be applied. Robust planning may be applied using scenarios generated as rigid body shifts of the patient geometry independently for each radiation set. Robust planning using scenarios generated by multiple non-rigid registrations can be applied. If an EBRT partial dose is delivered before a BT partial dose, this requires the use of multiple predictions. If a BT partial dose is delivered before an EBRT partial dose, this requires the use of registration perturbations between the acquired images. These deformations can also result from anatomical deformations during EBRT, e.g., tumor shrinkage.

[0053] Methods according to embodiments of the present invention may also be combined with multi-criteria optimization, where navigation can be performed in several dimensions with several trade-off objectives targeting either the total dose or the individual treatment doses.

[0054] Figure 5 is a schematic diagram of a computer system capable of carrying out the method of the present invention. The computer 31 comprises a processor 33 and a program memory 36 connected to first and second data memories 34, 35. One or more user input means 38, 39 are also present, preferably in the form of a keyboard, mouse, joystick, voice recognition means or any other available user input means. The user input means may also be arranged to receive data from an external memory unit.

[0055] A first data memory 34 contains data necessary to perform the method, such as necessary images. A second data memory 35 holds data relating to one or more current patients for whom a treatment plan is to be developed. A program memory 36 holds a computer program configured to cause a computer to perform the method steps, for example as discussed in connection with any of Figures 2, 3 and 4.

[0056] It will be appreciated that the data memories 34, 35 and the program memory 36 are shown and discussed schematically. There may be several data memory units, each holding one or more different types of data, or there may be one data memory that holds all data in a suitably structured manner, which holds the program memory. One or more of the components may be present in the cloud environment, as long as they can communicate with each other.

Claims

1. 1. A computer-implemented method for optimizing a radiation treatment plan for a patient, comprising the steps of: a. a computer receiving input data including a first image of the patient and dose criteria for a desired total dose distribution (S21; S31; S41); b) obtaining, based on the input data, an optimization problem (S23; S32; S42) including an objective function designed to optimize a total dose distribution as a combination of a first dose distribution delivered by a first ray set of a first radiation therapy and a second dose distribution delivered by a second ray set of a second radiation therapy, wherein one of the first radiation therapy and the second radiation therapy is external beam therapy and the other is brachytherapy; c. The computer optimizes a treatment plan as a combination of external beam therapy and brachytherapy using the optimization problem (S24; S33; S43); Including, optimizing includes transforming the first dose distribution and the second dose distribution into a common geometric shape and accumulating dose using a biological model, and the objective function includes a set of accumulated dose and radiation set-specific dose penalties; method.

2. The method of claim 1 , wherein the input data further comprises a second image of the patient.

3. The method of claim 2 , wherein the second image is an estimated result image of the patient after the first radiation therapy.

4. d. After delivery of radiation by a first radiation therapy of the treatment plan, the computer determines a dose delivered to the patient by the first radiation therapy (S26; S36; S46), the computer receives at least one updated image of the patient (S27; S37; S47), and the computer re-optimizes the second radiation therapy treatment plan using a second optimization problem taking into account the delivered dose and the updated image (S28; S38; S48). The method of any one of claims 1 to 3, further comprising:

5. 5. The method of claim 4, wherein a first radiation therapy of the treatment plan is external beam therapy, the at least one updated image includes an updated image of the patient after the external beam therapy and an updated image of the patient with a subsequent application of a brachytherapy device, and wherein the treatment plan for the second radiation therapy, brachytherapy, is optimized based on the updated images.

6. 5. The method of claim 4, wherein a first radiation therapy of the treatment plan is brachytherapy, and the input data further includes an image of the patient with a brachytherapy device applied.

7. The method of any one of claims 4 to 6, wherein the at least one updated image comprises at least one image of the patient taken after delivery of radiation from the first radiation therapy.

8. 7. The method of claim 4, wherein the at least one updated image comprises at least one simulated image based on an estimate of a patient's geometry after delivery of radiation from the first radiation therapy.

9. 9. The method of any one of claims 1 to 8, wherein robust planning is used to account for uncertainties in the brachytherapy, the external beam therapy, and / or the determined delivered dose.

10. A computer program comprising computer readable code means adapted to, when run on a computer, cause said computer to carry out the method of any one of claims 1 to 9.

11. A computer system (31) comprising a processor (33) and at least one program memory (36), characterized in that said program memory holds a computer program according to claim 10.

Images