Patent Application
20250065155
The present invention relates to radiation therapy delivery systems for
scanned proton beams, and in particular to quality assurance of treatment delivered from such systems.
A goal of radiation therapy, or radiotherapy is typically to deliver sufficient dose to a target such as a tumor, to eliminate, or at least reduce the size of the tumor. Radiation therapy delivery systems for delivering scanned ion beams to patients are known in the art. The ions are typically protons, but other types of ions may be used as well.
Radiotherapy treatment always involves delivering dose outside of the target, to healthy tissue or organs at risk, and therefore there is always a risk that healthy organs or tissue are damaged by the radiation. One emerging treatment method that appears to involve less undesired damage is FLASH therapy, which involves treatment at a much higher dose rate than conventional therapy, for example, 70 Gy/s. In the literature related to FLASH, various lower limits to the dose rates are suggested, such as at least 40 Gy/s or 50 Gy/s. If, for example, a dose of 20 Gy is to be delivered with a dose rate of 70 Gy/s, the whole dose will be delivered in 0.29 s. In contrast, conventional radiotherapy treatment is delivered at much lower dose rates, a typical dose rate for conventional radiotherapy treatment being a few Gy per minute.
It has been found that with FLASH therapy the damage done to healthy tissue by a particular dose is lower than with conventional therapy, while the effect on the target, that is, the tumorous tissue response, remains the same. Therefore, FLASH therapy can have significant advantage over conventional radiotherapy treatment, since they can affect the tumor with less damage done to the healthy tissue.
The FLASH effect is highly dependent on the time structure of the radiation delivery, to ensure that the dose rate is high enough to activate the FLASH effect.
Therefore, to ensure that a treatment plan has the intended effect on the patient, it is vital to know the time structure of the beam delivery. The rate of energy deposition in a voxel depends on a number of factors, including
When delivering FLASH treatments by scanned proton beams the dose is the resulting accumulation of energy depositions. The energy deposition over time is non-trivial and depends on the rate of energy deposition in a voxel, which is a function of the scanning characteristics of the delivery system, the beam current, the treatment plan and the beam quality.
Normally, a prediction of the system characteristics may be provided by the system vendors, in the form of nominal data or a model of the system. This may not always be the case, and even if it is, there may still be a desire to verify the data to ensure that treatment delivered to a patient will have the intended properties. As explained above, this is particularly important in the case of FLASH therapy, where the time structure is essential for achieving the desired effect.
Accordingly, it is an object of the present disclosure to ensure that radiotherapy treatment is delivered as intended by the delivery apparatus.
The invention relates to a computer-based method for use in a radiotherapy treatment delivery system arranged to deliver ion radiotherapy treatment according to a plan, comprising the steps of
This enables the assessment or verification of system characteristics, by comparing a simulation of a plan based on the model of the system to the actual outcome of delivery of the plan. As explained above, the system characteristics become particularly important in treatment plans such as FLASH treatment plans, that are highly dependent on the time structure of the delivery, which means that if the model is not correct, the effective dose can change considerably. The plan may be a treatment plan arranged for delivery to a patient, or a test plan, which is normally simpler than a treatment plan and used to test aspects of the delivery system. Adapting the model may involve changing the model itself or adjusting one or more parameters of the model.
The step of comparing typically includes comparing the time-resolved measurement to the time-resolved simulation.
According to embodiments, the method further comprises the steps of
The Relative Biological Effectiveness (RBE) weighted dose means the effective dose to the patient, which may differ from the physical dose, influenced by the dose rate, type of radiation and other factors as is known to the skilled person. The FLASH dose is a special case of RBE weighted dose, in which very high dose rates causes even greater changes between the physical dose and the RBE weighted dose.
In some preferred embodiments, the time-resolved measurement is obtained from a dry-run of the plan. This enables the adjustment of model parameters before actual delivery of a treatment plan to the patient.
In other embodiments, the time-resolved measurement is obtained during delivery of the plan to a patient. In this case, data obtained from the measurement may be used to adjust the remaining fractions of the plan based on any deviations between the dose according to the plan and the dose actually delivered.
As is common in the art, the treatment plan may be designed to be delivered to the patient in a number of fractions, wherein the time-resolved simulation is obtained for one or more of these fractions. In these cases, simulations may be obtained at different stages of the treatment, typically shortly before the delivery of one of the fractions, to account for the fact that changes in environmental factors such as temperature and atmospheric pressure may affect the simulation.
The time-resolved measurement may be obtained during delivery of the same one or more of the fractions as the time-resolved simulation, wherein the step of taking action includes adapting delivery of a subsequent fraction. This is particularly useful if the measurement is obtained as a dry run of the system delivering the subsequent fraction. In this case, the dry run may be used to obtain information about the current function of the system and what will be the result of delivery of the subsequent fraction, which may be used to adapt the subsequent fraction to the current function of the system.
Preferably, the time-resolved measurement is obtained during delivery of the same one or more of the fractions as the time-resolved simulation.
The method may include the steps of obtaining delivered FLASH dose information from the time-resolved measurement and simulated FLASH dose information from the time-resolved simulation and the step of comparing includes comparing the delivered FLASH dose information to the simulated FLASH dose information, and the step of taking action includes amending the plan for subsequent fractions.
The step of comparing may include comparing accumulated dose information from one fraction or from a number of fractions. If the time-resolved measurement is obtained from a dry-run before the delivery of a particular fraction, the result of the comparison may be used to adapt that particular fraction in view of the difference.
The step of taking action may include adjusting the model based on the result of the comparison.
The step of taking action may include assessing the function of the delivery system based on the result of the comparison.
The measuring may be performed in any suitable manner and using any suitable detector. For example, a detector such as a radiation detector array or a laterally integrating detector may be used.
The model of the system may include parameter values that are based on
nominal data provided for the delivery system. Such parameter values are normally provided by the vendor of the system. Alternatively, the model includes data obtained by an experimental characteristic of the system. This means that during parameter values are determined based on testing of the delivery system. In both cases, the method enables verification and/or adjustment of the parameters of the model. In both cases, the method may also be used to adjust the model itself.
The disclosure also relates to a computer program product comprising computer readable code means which, when run in a computer, will cause the computer to perform the method according to any one of the embodiments discussed in this document. The computer program product may include non-transitory storage medium on which the code means is stored.
The disclosure also relates to a computer system comprising a processor, a program memory and optionally a data memory, wherein the program memory comprises such a computer program product.
The invention will be described in more detail in the following, by way of examples and with reference to the appended drawings.
A typical dose rate and dose accumulation time trace of a voxel in a PBS beam is show in
According to the present disclosure, it is desired to obtain such time trace data, related to dose rate or dose accumulation, or both, to be used as input data when modelling the effect of the treatment plan, in particular of a FLASH treatment plan. Such time trace data may be acquired in a number of different ways. For example, in a proton beam scanning delivery system, the plan includes a number of energy levels and a number of spots within each energy level, as the skilled person is aware.
In step S25, the simulated signal and the measured signal are compared. Step S26 is a decision step for deciding which action to take in dependence of the result of the comparison. The outcome may be that some action needs to be taken S27. As will be discussed in more detail below, the actions may include adjusting the model parameters or improving the model itself, adapting the plan, or determining that the delivery system needs to be adjusted. The result of the comparison may also be that no action needs to be taken S28.
In each of the embodiments in this document, the model S22, S32, S42, S52 of the treatment delivery system may be a system model including nominal parameters provided by the vendor of the delivery system. It may also be based on an experimental characterization of the system. The model includes parameter values that likewise may be provided by the vendor of the delivery system or based on an experimental characterization of the system. The simulated signal in S23, S33, S43, S53 therefore represents a simulation of the treatment delivery system, which may deviate from the actual behaviour of the treatment delivery system.
The measuring may be performed using any suitable form of detector, typically a type of ion chamber detector. Such detectors are available as detector arrays and as laterally integrating detectors arranged to measure a total signal across the beam. The measured quantity may include the beam current, and/or one or more other parameters that change with time may be used in the delivery data, for example, one or more of the following:
In each of the embodiments in this document, the time resolved measurement S24, S34, S44, S54 may be performed as a dry-run of the plan using the treatment delivery system. A dry-run is performed without a patient present, and may be performed with or without a phantom. Another way of obtaining such data is during delivery of a plan to a patient. In this case a measuring device must be provided near the patient to measure the actual radiation provided, for example beam current as a function of time. The measurements must have a sufficiently high time resolution to enable the desired radiobiological modeling, preferably at least to one millisecond. Even more preferred would be a fraction of a millisecond, for example 1/10 of a millisecond. In both cases, the time-resolved measurement represents the actual beam-delivery sequence from the treatment delivery system.
In both the above cases, the measurement can be made by methods that are known per se, for example, by picking up the signal from one of the existing ion chamber monitors or by inserting a dedicated and independent detector in the beam or in a phantom being irradiated. The measurements may be made over a suitable time period matching the time of radiation delivery. The detector may typically comprise an array of detectors, each arranged to detect the time when it receives dose, and the magnitude of the dose received. Alternatively, a single ion chamber detector may be used, such as a broad plane parallel ion chamber that measures the current across the whole beam. Both radiation detector arrays and single ion chamber detectors are available to the skilled person. By matching the plan data (spot data and energy layer data) with the time trace an accurate plan delivery timing characterization will be possible. The data will have essentially the same structure as shown in
As is common in the art, the treatment plan may be divided into a number of treatment fractions, to be delivered to the patient at certain time intervals, for example at the rate of one fraction per day. The simulation in step S23, S33, S43, S53 and the measurement in step S24, S34, S44, S54 may be performed for a complete plan including one or more fractions. Alternatively, the simulation and measuremen steps may be performed for one or more fractions.
As will be discussed below, the comparison may be performed by comparing the simulated time-resolved signal to the time-resolved measured signal, or the signals may be processed in a suitable way before comparison. For example, as will be discussed in more detail in connection with
For each of the embodiments discussed in connection with
As before, the input data include a treatment plan S51 and a model S52 of the treatment delivery system. A simulated signal that is a time-resolved measured signal of a quantity representative of a beam delivery sequence for delivery of one or more of the fractions of the plan is obtained in step S53. A measured signal that is a time-resolved measured signal of a quantity representative of a beam delivery sequence for delivery of the one or more of the fractions of the plan is obtained in step S54. Optionally, steps S53′ and S54′ may be included to calculate simulated dose information and measured dose information based on the simulated signal and the measured signal, respectively, in the same way as in steps S43′ and S44′. In step S55, the simulated signal and the measured signal, and/or the simulated and measured doses, are compared. Step S56 is a decision step in which the need for further action is determined.
In a first implementation of the third embodiment, the measured signal is obtained in step S54 during a dry run of a considered fraction that has not yet been delivered to the patient, preferably a fraction that is to be delivered shortly after the measurement, for example on the same day. In this case, the result of the comparison in step S55 will be an indicator of how the system will perform the delivery and what will be the actual dose resulting from the considered fraction. The dry run may also consider more than one fraction. In this implementation, the result of the comparison can be used to change the considered fraction, and possibly other remaining fractions, to adjust to the system's current properties. It can also be used to determine that a whole new plan needs to be obtained, if the simulated signal and the measured signal are too different.
In a second implementation of the third embodiment, the measured signal is obtained in step S54 may have been obtained during delivery of one or more previous fractions, or as a dry run of the one or more previous fraction or fractions. In step S54′, therefore, actual accumulated dose information from the previous fraction or fractions is obtained, based either on a dry run of the system or on actual delivery of the fractions to the patient. The simulated accumulated dose is calculated in step S53′ based on the simulated signal. The information regarding the actual accumulated dose may be used to recalculate the remaining fractions, to take into account the difference between the actual accumulated dose and the simulated accumulated dose.
The data memory 64 holds data that may be used in the method, such as the nominal system data provided by the vendor, one or more treatment plans, and/or simulation results. As will be understood, the data memory 64 is only shown schematically. There may be several data memory units, each holding one or more different types of data.
The program memory 65 holds a computer program arranged to control the processor to perform the procedure according to any of the embodiments discussed above. Like the data memory 64, the program memory may also be implemented as one or several units as is seen fit.
| Number | Date | Country | Kind |
|---|---|---|---|
| 22152263.4 | Jan 2022 | EP | regional |
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/EP2022/085199 | 12/9/2022 | WO |
