Patent Grant
8849373
1. Field of the Invention
The present invention generally relates to radiotherapy, and in particular to methods for estimating internal target position for use by systems able to continuously align a radiation beam with the target.
2. Background Description
Real-time internal target position estimation is of high interest in radiotherapy, particularly with the recent development of robotic, linear accelerator, DMLC and couch-based systems which can continuously align the radiation beam with the target. It would make an important contribution to effective use of these systems to develop an algorithm to combine x-ray and external respiratory monitoring sources for estimating real-time target position. Currently, there is no system in radiotherapy to estimate internal position from a single x-ray imager and an optical system. What is needed is a methodology for combining x-ray and external respiratory sources for estimating real-time target position.
It is therefore an object of the invention to provide a methodology for combining x-ray and external respiratory sources for estimating target position in real-time.
An aspect of the invention is a method for real-time position estimation of an oncology target in a patient, the method comprising continuous monitoring of a respiratory signal of the patient; imaging of the target using a single imaging device, the imaging being done periodically to limit radiation dosage administered by the single imaging device; applying a model to estimate a centroid position of the target from a combination of the respiratory signal and target position data from the imaging; determining and updating parameters of the model by using position displacement data from the imaging; and using the updated model in real-time to apply a treatment program.
A further aspect of the invention is the optional implanting of a marker in the patient and then using the location information of the marker obtained in the periodic imaging to improve accuracy of the target position data. It is also an aspect of the invention to adjust the periodicity of imaging to implement a tradeoff between position estimation accuracy and said radiation dosage limitation.
Yet another aspect of the invention is to make an initial determination of the model parameters in a pre-treatment session by using four dimensional computed tomography (4D CT) or four dimensional cone beam computed tomography (4D CBCT). In a preferred embodiment the imaging device is an X-ray machine, and more particularly a linear accelerator.
The foregoing and other objects, aspects and advantages will be better understood from the following detailed description of a preferred embodiment of the invention with reference to the drawings, in which:
Table 1 is a list of linear accelerator vendors, their configurations and comments on the applicability of using a single x-ray imager and respiratory monitoring system to determine target position.
Table 2 is a list of terms with their symbols and descriptions.
Real-time internal target position estimation is of high interest in radiotherapy, particularly with the recent development of robotic, linear accelerator, DMLC (Dynamic Multi-Leaf Collimation) and couch-based systems which can continuously align the radiation beam with the target. To estimate the internal target position, one interesting approach is to combine x-ray images with respiratory monitoring. X-ray images give internal anatomic information, particularly if radio-opaque objects such as implanted markers are used. However, this procedure also gives additional radiation dose to the patient.
Respiratory monitoring gives a motion signal without additional ionizing radiation but does not yield target position information. The current invention focuses on the use of a single x-ray imager with 4D CT or 4D CBCT prior information. Note that a methodology for this approach with orthogonal simultaneous x-ray imagers and without 4D CT information has been developed and is commercially available. Simultaneous imaging from multiple directions can be considered an extension of the current algorithm where the time difference between acquired images at different angles is small or zero.
The reason for choosing a single imager and respiratory signal configuration is its high accuracy (an overall accuracy from the 160 3D tumor motion data in seen to be within an average rms error of 1 mm for the entire data set and 1.5 mm for datasets in which the average tumor motion is over 5 mm per respiratory cycle), compromise between radiation dose and position monitoring accuracy and widespread applicability to commercially available linear accelerator (and potentially proton accelerator) designs, as explained in Table 1. Also, there are many commercially available sources of respiratory signal, including optical monitoring, strain gauge, pressure sensor, nasal thermometer and spirometer.
This invention describes an algorithm to combine x-ray and external respiratory monitoring sources for estimating real-time target position.
In the scenario contemplated by the invention, there is a quasi-continuous external position measurement, and occasional internal position measurements. This approach can be generalized to cover multiple sources of the internal position. The objective is accurate and precise estimates of internal position quasi-continuously.
It is assumed that internal and external motion have the same basis function, F, but differences in position, motion ratio and time lag/lead (phase differences) have been observed. The form of F will be different for different patients, and we do not wish to characterize this as a given function, e.g. sinusoidal. However, the form of F will be characterized by a period (which can change), because this is necessary for estimates as described below.
Thus we can formulate:
I(t)=xI+aIF(t+Δt)
E(t)=xE+aEF(t)
Using the discrete measurements of Ĩ and {tilde over (E)} at n discrete time points, iε{1:n} spanning an integral number of respiratory cycles (e.g. 4D CT or fluoroscopy data), we wish to build a model for Î(t) based on E(t).
Thus
xI and xE can be found from the mean of Ĩ and {tilde over (E)} respectively.
can be estimated from the dispersion of Ĩ/{tilde over (E)}. Δt is estimated by maximizing the correlation between Î(t) and Ĩ over the same sampling period which also gives an estimate of the lower bound of the uncertainty.
Generalizing to different directions with different Δt values will result in hysteresis.
The present invention describes an algorithm to combine x-ray and external respiratory monitoring sources for estimating real-time target position. The target position with respect to the linac isocenter, {circumflex over (T)}(t), can be calculated using the general equation {circumflex over (T)}(t)=Î(t)+{circumflex over (M)}θ(t)R(t), where Î(t) is the estimate of the internal/external displacement at time t, {circumflex over (M)}θ(t) is the estimated motion relationship between internal motion and external motion as a function of phase at time t and R(t) is the external respiratory signal. The coefficients of Î(t) and {circumflex over (M)}θ(t) are updated periodically using the 2D target position from the x-ray projection, {tilde over (T)}(t)2Dg,c. The initial coefficient values are given from 4D CT (Four Dimensional Computed Tomography) or 4D CBCT (Cone Beam Computed Tomography). An estimate of the optimal operating characteristics of the invention may be made using a cost-benefit analysis of the trade-off of accuracy versus imaging dose.
Method Description
For clarity a definition of terms is given in Table 2. A schematic of the estimation process is given in
The aim of position monitoring is to estimate the target centroid position with respect to the linac (linear accelerator) isocenter, {circumflex over (T)}(x,y,z,t)≡{circumflex over (T)}(t), throughout a course of therapy. In principle, position monitoring can be used to generate higher dimensionality information rather than just an estimated target position. However such monitoring will be the subject of future research.
General Formalism
In order to reduce dose the x-ray imager is triggered periodically (rather than continuously) to determine and update the model parameters used to estimate {circumflex over (T)}(t) based on the measured respiratory signal, {tilde over (R)}(t) The target position with respect to the linac isocenter, {circumflex over (T)}(t), is calculated using the general equation:
{circumflex over (T)}(t)=Î(t)+{circumflex over (M)}θ(t)R(t) (1)
Note that Î(t) could be incorporated directly into {circumflex over (M)}θ(t), however given the potentially large variation between Ĩp and Î(t) (e.g. due to a change in the day to day position of the respiratory monitor) this displacement was included separately.
Initial Estimates
The initial estimates of the parameters for equation (1) are acquired from the 4D CT (or 4D CBCT) scan. For each phase of the 4D CT scan the external signal position and internal position are computed. Thus Ĩp, {tilde over (M)}p(θ) and {tilde over (T)}p(θ) are determined, and an initial discretized form of {tilde over (M)}p(θ) is obtained. However, as the phase is continually changing with time, an appropriate curvilinear relationship is obtained to ensure Mp(θ) is continuous over [0,2π).
Periodic Updating
It is assumed that the internal/external relationship changes with time (otherwise the initial information obtained from the 4D CT scan would be sufficient and no further imaging would be necessary). In general, a single planar x-ray image will not be able to resolve the marker positions along the beam axis, so there will be some residual uncertainty in the estimation of the marker position.
The coefficients of Î(t) and {circumflex over (M)}θ(t) are updated with information from the projection images, through {tilde over (T)}(t)2Dg,c where the update of the coefficients are weighted by the resolving power, i.e. for an anterior-posterior (AP) aligned imager the superior-inferior coefficients can be updated however the AP coefficients remain the same. As the gantry rotates, all of the coefficients will be updated. Differences between the estimated position {circumflex over (T)}(t), projected along the same gantry angle as the imager and couch angle as the patient, and {tilde over (T)}(t)2Dg,c, are used to update the correlation model.
Any feedback system has a finite response time. The foregoing system and method can be implemented with motion prediction to account for this response time, depending on the time for computation and response device used for these calculations.
While the invention has been described in terms of preferred embodiments, those skilled in the art will recognize that the invention can be practiced with modification within the spirit and scope of the appended claims.
This application claims the benefit of U.S. Provisional Application No. 60/917,637 of the same title filed on May 11, 2007.
This invention was made with Government support under contract CA093626 awarded by the National Institutes of Health. The Government has certain rights in this invention.
| Number | Name | Date | Kind |
|---|---|---|---|
| 6144875 | Schweikard et al. | Nov 2000 | A |
| 6501981 | Schweikard et al. | Dec 2002 | B1 |
| 7260426 | Schweikard et al. | Aug 2007 | B2 |
| 7570738 | Khamene et al. | Aug 2009 | B2 |
| 7620444 | Le et al. | Nov 2009 | B2 |
| 20080021300 | Allison | Jan 2008 | A1 |
| 20080081991 | West et al. | Apr 2008 | A1 |
| Number | Date | Country | |
|---|---|---|---|
| 20080281192 A1 | Nov 2008 | US |
| Number | Date | Country | |
|---|---|---|---|
| 60917637 | May 2007 | US |
