Pascal Paysan

Self‐adapting cyclic registration for motion‐compensated cone‐beam CT in image‐guided radiation therapy

PURPOSE In image-guided radiation therapy an additional kV imaging system next to the linear particle accelerator provides information for an accurate patient positioning. However, the acquisition time of the system is much longer than the patients breathing cycle due to the low gantry rotation speed. Our purpose is a cyclic registration in the context of motion-compensated image reconstruction that provides high quality respiratory-correlated 4D volumes for on-board flat panel detector cone-beam CT scans. METHODS Based on the small motion assumption, widely used within registration algorithms, a strategy is developed for motion estimation. In this strategy temporal restrictions are incorporated, for example, the cyclic motion patterns of respiration. The resultant cyclic registration method is to show less sensitivity on image artifacts, in particular on artifacts due to projection data sparsification. Using a new cyclic registration method a motion estimation is performed on respiratory-correlated reconstructions, and the obtained motion vector fields are used for motion compensation. RESULTS The proposed cyclic registration is evaluated in the context of motion-compensated image reconstruction using simulated data and patient data. Motion artifacts of 3D standard reconstructions can be significantly reduced by the resulting cyclic motion compensation. The method outperforms the respiratory-correlated reconstructions regarding sparse-view artifacts and maintains the high temporal resolution at the same time. Image artifacts show only minor to almost no effect on the motion estimation using the cyclic registration. CONCLUSIONS The cyclic motion compensation approach provides respiratory-correlated volumes with high image quality. The cyclic motion estimation is of such low sensitivity to sparse-view artifacts, that it is capable to determine high quality motion vector fields based on reconstructions of low sampled data.

Artifact‐resistant motion estimation with a patient‐specific artifact model for motion‐compensated cone‐beam CT

PURPOSE In image-guided radiation therapy (IGRT) valuable information for patient positioning, dose verification, and adaptive treatment planning is provided by an additional kV imaging unit. However, due to the limited gantry rotation speed during treatment the typical acquisition time is quite long. Tomographic images of the thorax suffer from motion blurring or, if a gated 4D reconstruction is performed, from significant streak artifacts. Our purpose is to provide a method that reliably estimates respiratory motion in presence of severe artifacts. The estimated motion vector fields are then used for motion-compensated image reconstruction to provide high quality respiratory-correlated 4D volumes for on-board cone-beam CT (CBCT) scans. METHODS The proposed motion estimation method consists of a model that explicitly addresses image artifacts because in presence of severe artifacts state-of-the-art registration methods tend to register artifacts rather than anatomy. Our artifact model, e.g., generates streak artifacts very similar to those included in the gated 4D CBCT images, but it does not include respiratory motion. In combination with a registration strategy, the model gives an error estimate that is used to compensate the corresponding errors of the motion vector fields that are estimated from the gated 4D CBCT images. The algorithm is tested in combination with a cyclic registration approach using temporal constraints and with a standard 3D-3D registration approach. A qualitative and quantitative evaluation of the motion-compensated results was performed using simulated rawdata created on basis of clinical CT data. Further evaluation includes patient data which were scanned with an on-board CBCT system. RESULTS The model-based motion estimation method is nearly insensitive to image artifacts of gated 4D reconstructions as they are caused by angular undersampling. The motion is accurately estimated and our motion-compensated image reconstruction algorithm can correct for it. Motion artifacts of 3D standard reconstruction are significantly reduced, while almost no new artifacts are introduced. CONCLUSIONS Using the artifact model allows to accurately estimate and compensate for patient motion, even if the initial reconstructions are of very low image quality. Using our approach together with a cyclic registration algorithm yields a combination which shows almost no sensitivity to sparse-view artifacts and thus ensures both high spatial and high temporal resolution.

Acuros CTS: A fast, linear Boltzmann transport equation solver for computed tomography scatter – Part II: System modeling, scatter correction, and optimization

PURPOSE To correct for scatter in kV cone-beam CT (CBCT) projection data on a clinical system using a new tool, Acuros® CTS, that estimates scatter images rapidly and accurately by deterministically solving the linear Boltzmann transport equation. METHODS Phantom and patient CBCT scans were acquired on TrueBeam® radiotherapy machines. A first-pass reconstruction was used to create water and bone density maps of the imaged object, which was updated to include a more accurate representation of the patient couch. The imaging system model accounted for the TrueBeam x-ray source (polychromatic spectrum, beam filtration, bowtie filter, and collimation hardware) and x-ray detection system (antiscatter grid, flat-panel imager). Acuros CTS then used the system and object models to estimate the scatter component of each projection image, which was subtracted from the measured projections. The corrected projections were then reconstructed to produce the final result. We examined the tradeoff between run time and accuracy using a Pareto optimization of key parameters, including the voxel size of the down-sampled object model, the number of pixels in the down-sampled detector, and the number of scatter images (angular down-sampling). All computations and reconstructions were performed on a research workstation containing two graphics processing units (GPUs). In addition, we established a method for selecting a subset of projections for which scatter images were calculated. The projections were selected to minimize interpolation errors in the remaining projections. Image quality improvement was assessed by measuring the accuracy of the reconstructed phantom and patient images. RESULTS The Pareto optimization yielded a set of parameters with an average run time of 26 seconds for scatter correction while maintaining high accuracy of scatter estimation. This was achieved in part by means of optimizing the projection angles that were processed, thus favoring the use of more angles in the lateral (i.e., horizontal) direction and fewer angles in the AP direction. In a 40 cm solid water phantom reconstruction, nonuniformities were decreased from 217 HU without scatter correction to 51 HU with conventional (kernel-based) scatter correction to 17 HU with Acuros CTS-based scatter correction. In clinical pelvis scans, nonuniformities in the bladder were reduced from 85 HU with conventional scatter correction to 14 HU with Acuros CTS. CONCLUSIONS Acuros CTS is a promising new tool for fast and accurate scatter correction for CBCT imaging. By carefully modeling the imaging chain and optimizing several parameters, we achieved high correction accuracies with computation times compatible with the clinical workflow. The improvement in image quality enables better soft-tissue visualization and potentially enables applications such as adaptive radiotherapy.

Five-dimensional motion compensation for respiratory and cardiac motion with cone-beam CT of the thorax region

We propose an adapted method of our previously published five-dimensional (5D) motion compensation (MoCo) algorithm1, developed for micro-CT imaging of small animals, to provide for the first time motion artifact-free 5D cone-beam CT (CBCT) images from a conventional flat detector-based CBCT scan of clinical patients. Image quality of retrospectively respiratory- and cardiac-gated volumes from flat detector CBCT scans is deteriorated by severe sparse projection artifacts. These artifacts further complicate motion estimation, as it is required for MoCo image reconstruction. For high quality 5D CBCT images at the same x-ray dose and the same number of projections as todays 3D CBCT we developed a double MoCo approach based on motion vector fields (MVFs) for respiratory and cardiac motion. In a first step our already published four-dimensional (4D) artifact-specific cyclic motion-compensation (acMoCo) approach is applied to compensate for the respiratory patient motion. With this information a cyclic phase-gated deformable heart registration algorithm is applied to the respiratory motion-compensated 4D CBCT data, thus resulting in cardiac MVFs. We apply these MVFs on double-gated images and thereby respiratory and cardiac motion-compensated 5D CBCT images are obtained. Our 5D MoCo approach processing patient data acquired with the TrueBeam 4D CBCT system (Varian Medical Systems). Our double MoCo approach turned out to be very efficient and removed nearly all streak artifacts due to making use of 100% of the projection data for each reconstructed frame. The 5D MoCo patient data show fine details and no motion blurring, even in regions close to the heart where motion is fastest.

WE-AB-207A-08: BEST IN PHYSICS (IMAGING): Advanced Scatter Correction and Iterative Reconstruction for Improved Cone-Beam CT Imaging On the TrueBeam Radiotherapy Machine

PURPOSE To improve CBCT image quality for image-guided radiotherapy by applying advanced reconstruction algorithms to overcome scatter, noise, and artifact limitations METHODS: CBCT is used extensively for patient setup in radiotherapy. However, image quality generally falls short of diagnostic CT, limiting soft-tissue based positioning and potential applications such as adaptive radiotherapy. The conventional TrueBeam CBCT reconstructor uses a basic scatter correction and FDK reconstruction, resulting in residual scatter artifacts, suboptimal image noise characteristics, and other artifacts like cone-beam artifacts. We have developed an advanced scatter correction that uses a finite-element solver (AcurosCTS) to model the behavior of photons as they pass (and scatter) through the object. Furthermore, iterative reconstruction is applied to the scatter-corrected projections, enforcing data consistency with statistical weighting and applying an edge-preserving image regularizer to reduce image noise. The combined algorithms have been implemented on a GPU. CBCT projections from clinically operating TrueBeam systems have been used to compare image quality between the conventional and improved reconstruction methods. Planning CT images of the same patients have also been compared. RESULTS The advanced scatter correction removes shading and inhomogeneity artifacts, reducing the scatter artifact from 99.5 HU to 13.7 HU in a typical pelvis case. Iterative reconstruction provides further benefit by reducing image noise and eliminating streak artifacts, thereby improving soft-tissue visualization. In a clinical head and pelvis CBCT, the noise was reduced by 43% and 48%, respectively, with no change in spatial resolution (assessed visually). Additional benefits include reduction of cone-beam artifacts and reduction of metal artifacts due to intrinsic downweighting of corrupted rays. CONCLUSION The combination of an advanced scatter correction with iterative reconstruction substantially improves CBCT image quality. It is anticipated that clinically acceptable reconstruction times will result from a multi-GPU implementation (the algorithms are under active development and not yet commercially available). All authors are employees of and (may) own stock of Varian Medical Systems.

Motion vector field upsampling for improved 4D cone-beam CT motion compensation of the thorax

To improve the accuracy of motion vector fields (MVFs) required for respiratory motion compensated (MoCo) CT image reconstruction without increasing the computational complexity of the MVF estimation approach, we propose a MVF upsampling method that is able to reduce the motion blurring in reconstructed 4D images. While respiratory gating improves the temporal resolution, it leads to sparse view sampling artifacts. MoCo image reconstruction has the potential to remove all motion artifacts while simultaneously making use of 100% of the rawdata. However the MVF accuracy is still below the temporal resolution of the CBCT data acquisition. Increasing the number of motion bins would increase reconstruction time and amplify sparse view artifacts, but not necessarily the accuracy of MVF. Therefore we propose a new method to upsample estimated MVFs and use those for MoCo. To estimate the MVFs, a modified version of the Demons algorithm is used. Our proposed method is able to interpolate the original MVFs up to a factor that each projection has its own individual MVF. To validate the method we use an artificially deformed clinical CT scan, with a breathing pattern of a real patient, and patient data acquired with a TrueBeamTM4D CBCT system (Varian Medical Systems). We evaluate our method for different numbers of respiratory bins, each again with different upsampling factors. Employing our upsampling method, motion blurring in the reconstructed 4D images, induced by irregular breathing and the limited temporal resolution of phase–correlated images, is substantially reduced.

Motion vector field phase-to-amplitude resampling for 4D motion-compensated cone-beam CT

We propose a phase-to-amplitude resampling (PTAR) method to reduce motion blurring in motion-compensated (MoCo) 4D cone-beam CT (CBCT) image reconstruction, without increasing the computational complexity of the motion vector field (MVF) estimation approach. PTAR is able to improve the image quality in reconstructed 4D volumes, including both regular and irregular respiration patterns. The PTAR approach starts with a robust phase-gating procedure for the initial MVF estimation and then switches to a phase-adapted amplitude gating method. The switch implies an MVF-resampling, which makes them amplitude-specific. PTAR ensures that the MVFs, which have been estimated on phase-gated reconstructions, are still valid for all amplitude-gated reconstructions. To validate the method, we use an artificially deformed clinical CT scan with a realistic breathing pattern and several patient data sets acquired with a TrueBeamTM integrated imaging system (Varian Medical Systems, Palo Alto, CA, USA). Motion blurring, which still occurs around the area of the diaphragm or at small vessels above the diaphragm in artifact-specific cyclic motion compensation (acMoCo) images based on phase-gating, is significantly reduced by PTAR. Also, small lung structures appear sharper in the images. This is demonstrated both for simulated and real patient data. A quantification of the sharpness of the diaphragm confirms these findings. PTAR improves the image quality of 4D MoCo reconstructions compared to conventional phase-gated MoCo images, in particular for irregular breathing patterns. Thus, PTAR increases the robustness of MoCo reconstructions for CBCT. Because PTAR does not require any additional steps for the MVF estimation, it is computationally efficient. Our method is not restricted to CBCT but could rather be applied to other image modalities.

Optimization-based reconstruction for reduction of CBCT artifact in IGRT

Kilo-voltage cone-beam computed tomography (CBCT) plays an important role in image guided radiation therapy (IGRT) by providing 3D spatial information of tumor potentially useful for optimizing treatment planning. In current IGRT CBCT system, reconstructed images obtained with analytic algorithms, such as FDK algorithm and its variants, may contain artifacts. In an attempt to compensate for the artifacts, we investigate optimization-based reconstruction algorithms such as the ASD-POCS algorithm for potentially reducing arti- facts in IGRT CBCT images. In this study, using data acquired with a physical phantom and a patient subject, we demonstrate that the ASD-POCS reconstruction can significantly reduce artifacts observed in clinical re- constructions. Moreover, patient images reconstructed by use of the ASD-POCS algorithm indicate a contrast level of soft-tissue improved over that of the clinical reconstruction. We have also performed reconstructions from sparse-view data, and observe that, for current clinical imaging conditions, ASD-POCS reconstructions from data collected at one half of the current clinical projection views appear to show image quality, in terms of spatial and soft-tissue-contrast resolution, higher than that of the corresponding clinical reconstructions.

Iterative motion-compensated reconstruction for image-guided radiation therapy

In image-guided radiation therapy (IGRT) an additional kV imaging system next to the linear particle accelerator provides information for an accurate patient positioning. However, due to the limited gantry rotation speed during treatment the typical acquisition time is much longer than the patients breathing cycle resulting in low image quality. In particular, respiratory motion causes severe artifacts such as blurring and streaks in tomographic images. Compensating for motion is an interesting option and capable of providing high quality respiratory-correlated 4D volumes. The particular challenge is to do this without knowledge from prior scans and without specific requirements on the acquisition. State-of-the-art methods for estimation of the motion vector fields suffer from image artifacts that appear in intermediated image frames. In applications like ours conventional registration algorithms tend to register artifacts rather than anatomy. Our proposed registration method addresses the image artifact problem by additionally using temporal constraints within the registration process like the cyclic motion patterns of respiration. This cyclic regularization avoids the algorithm being sensitive to the above-mentioned streak artifacts. Hence, this registration algorithm consists of two components, the spatial and the temporal part, and one of them might be dominating. For this very reason the cyclic registration method is supplemented by an iterative refinement to avoid a potential underestimation of motion and simultaneously to comply with the temporal constraints. The iterative and cyclic registration method is verified using simulated rawdata and patient data. It shows low sensitivity on image artifacts and deals with potential underestimation of motion at the same time. In this way, the motion is accurately estimated and a motion compensation corrects for it. Our iterative and cyclic registration method is capable of high quality motion estimation only on basis of respiratory-correlated reconstructions. Thus, motion-compensated image reconstruction without knowledge from prior scans and without specific requirements on the acquisition becomes feasible in image-guided radiation therapy.