Respiratory Motion Estimation from Slowly Rotating X-Ray Projections

As radiotherapy has become increasingly conformal, geometric uncertainties caused by breathing and organ motion have become an important issue. Accurate motion estimates may lead to improved treatment planning and dose calculation in radiation therapy. However, respiratory motion is difficult to study by conventional X-ray CT imaging since object motion causes inconsistent projection views leading to artifacts in reconstructed images. We propose to estimate the parameters of a nonrigid motion model from a set of projection views of the thorax that are acquired using a slowly rotating cone-beam CT scanner, such as a radiotherapy simulator. We use a conventionally reconstructed 3D thorax image, acquired by breath-hold CT, as a reference volume. We represent respiratory motion using a flexible parametric nonrigid motion model based on B-splines. The motion parameters are estimated by optimizing a regularized cost function that includes the squared error between the measured projection views and the reprojections of the deformed reference image. Preliminary 2D simulation results show that there is good agreement between the estimated motion and the true motion.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85878/1/Fessler197.pd

Estimating 3-D Respiratory Motion From Orbiting Views by Tomographic Image Registration

Respiratory motion remains a significant source of errors in treatment planning for the thorax and upper abdomen. Recently, we proposed a method to estimate two-dimensional (2-D) object motion from a sequence of slowly rotating X-ray projection views, which we called deformation from orbiting views (DOVs). In this method, we model the motion as a time varying deformation of a static prior of the anatomy. We then optimize the parameters of the motion model by maximizing the similarity between the modeled and actual projection views. This paper extends the method to full three-dimensional (3-D) motion and cone-beam projection views. We address several practical issues for using a cone-beam computed tomography (CBCT) scanner that is integrated in a radiotherapy system, such as the effects of Compton scatter and the limited gantry rotation for one breathing cycle. We also present simulation and phantom results to illustrate the performance of this method.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85995/1/Fessler38.pd

Estimating 3D Respiratory Motion from Orbiting Views

This paper describes a method for estimating 3D respiratory motion so as to characterize tumor motion. This method uses two sets of measurements. One is a reference thorax volume obtained from a conventional fast CT scanner under breath-hold condition. The other is a sequence of projection views of the same patient (acquired at treatment time) using a slowly rotating cone-beam system (1 minute per rotation) during free breathing. We named this method deformation from orbiting views (DOV). Breathing motion over the entire acquisition period is estimated by deforming the reference volume through time so that its projections best match the measured projection views. The nonrigid breathing motion is described by a B-spline based deformation model. The parameters of this model are estimated by minimizing a regularized squared error cost function, using a conjugate gradient descent algorithm. Performance of this approach was evaluated by simulation. Results showed good agreement between the estimated and synthesized motion, with a mean absolute error of 1.63 mm. Relatively larger errors tended to occur in uniform regions, which would not have significant effects on generating deformed volumes based on the estimated motion. The results indicate that it is feasible to estimate realistic nonrigid motion from a sequence of slowly rotating cone beam projection views.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85996/1/Fessler214.pd

Reconstruction of implanted marker trajectories from cone-beam CT projection images using interdimensional correlation modeling

PURPOSE: Cone-beam CT (CBCT) is a widely used imaging modality for image-guided radiotherapy. Most vendors provide CBCT systems that are mounted on a linac gantry. Thus, CBCT can be used to estimate the actual 3-dimensional (3D) position of moving respiratory targets in the thoracic/abdominal region using 2D projection images. The authors have developed a method for estimating the 3D trajectory of respiratory-induced target motion from CBCT projection images using interdimensional correlation modeling. METHODS: Because the superior-inferior (SI) motion of a target can be easily analyzed on projection images of a gantry-mounted CBCT system, the authors investigated the interdimensional correlation of the SI motion with left-right and anterior-posterior (AP) movements while the gantry is rotating. A simple linear model and a state-augmented model were implemented and applied to the interdimensional correlation analysis, and their performance was compared. The parameters of the interdimensional correlation models were determined by least-square estimation of the 2D error between the actual and estimated projected target position. The method was validated using 160 3D tumor trajectories from 46 thoracic/abdominal cancer patients obtained during CyberKnife treatment. The authors' simulations assumed two application scenarios: (1) retrospective estimation for the purpose of moving tumor setup used just after volumetric matching with CBCT; and (2) on-the-fly estimation for the purpose of real-time target position estimation during gating or tracking delivery, either for full-rotation volumetric-modulated arc therapy (VMAT) in 60 s or a stationary six-field intensity-modulated radiation therapy (IMRT) with a beam delivery time of 20 s. RESULTS: For the retrospective CBCT simulations, the mean 3D root-mean-square error (RMSE) for all 4893 trajectory segments was 0.41 mm (simple linear model) and 0.35 mm (state-augmented model). In the on-the-fly simulations, prior projections over more than 60° appear to be necessary for reliable estimations. The mean 3D RMSE during beam delivery after the simple linear model had established with a prior 90° projection data was 0.42 mm for VMAT and 0.45 mm for IMRT. CONCLUSIONS: The proposed method does not require any internal/external correlation or statistical modeling to estimate the target trajectory and can be used for both retrospective image-guided radiotherapy with CBCT projection images and real-time target position monitoring for respiratory gating or tracking.NHMRC, National Research Foundation of Kore