An investigation of 4D cone-beam CT algorithms for slowly rotating scanners.

PURPOSE: To evaluate several algorithms for 4D cone-beam computed tomography (4D CBCT) with slow rotating devices. 4D CBCT is used to perform phase-correlated (PC) reconstructions of moving objects, such as breathing patients, for example. Such motion phase-dependent reconstructions are especially useful for updating treatment plans in radiation therapy. The treatment plan can be registered more precisely to the motion of the tumor and, in consequence, the irradiation margins for the treatment, the so-called planning target volume, can be reduced significantly METHODS: In the study, several algorithms were evaluated for kilovoltage cone-beam CT units attached to linear particle accelerators. The reconstruction algorithms were the conventional PC reconstruction, the McKinnon-Bates (MKB) algorithm, the prior image constrained compressed sensing (PICCS) approach, a total variation minimization (ASD-POCS) algorithm, and the auto-adaptive phase correlation (AAPC) algorithm. For each algorithm, the same motion-affected raw data were used, i.e., one simulated and one measured data set. The reconstruction results from the authors' implementation of these algorithms were evaluated regarding their noise and artifact levels, their residual motion blur, and their computational complexity and convergence. RESULTS: In general, it turned out that the residual motion blur was lowest in those algorithms which exclusively use data from a single motion phase. Algorithms using the image from the full data set as initialization or as a reference for the reconstruction were not capable of fully removing the motion blurring. The iterative algorithms, especially approaches based on total variation minimization, showed lower noise and artifact levels but were computationally complex. The conventional methods based on a single filtered backprojection were computationally inexpensive but suffered from higher noise and streak artifacts which limit the usability. In contrast, these methods showed to be less demanding and more predictable in their outcome than the total variation minimization based approaches. CONCLUSIONS: The reconstruction algorithms including at least one iterative step can reduce the 4 CBCT specific artifacts. Nevertheless, the algorithms that use the full data set, at least for initialization, such as MKB and PICCS in the authors' implementation, are only a trade-off and may not fully achieve the optimal temporal resolution. A predictable image quality as seen in conventional reconstruction methods, i.e., without total variation minimization, is a desirable property for reconstruction algorithms. Fast, iterative approaches such as the MKB can therefore be seen as a suitable tradeoff.

Autoadaptive phase-correlated (AAPC) reconstruction for 4D CBCT.

PURPOSE: Kilovoltage cone-beam computed tomography (CBCT) is widely used in image-guided radiation therapy for exact patient positioning prior to the treatment. However, producing time series of volumetric images (4D CBCT) of moving anatomical structures remains challenging. The presented work introduces a novel method, combining high temporal resolution inside anatomical regions with strong motion and image quality improvement in regions with little motion. METHODS: In the proposed method, the projections are divided into regions that are subject to motion and regions at rest. The latter ones will be shared among phase bins, leading thus to an overall reduction in artifacts and noise. An algorithm based on the concept of optical flow was developed to analyze motion-induced changes between projections. The technique was optimized to distinguish patient motion and motion deriving from gantry rotation. The effectiveness of the method is shown in numerical simulations and patient data. RESULTS: The images reconstructed from the presented method yield an almost the same temporal resolution in the moving volume segments as a conventional phase-correlated reconstruction, while reducing the noise in the motionless regions down to the level of a standard reconstruction without phase correlation. The proposed simple motion segmentation scheme is yet limited to rotation speeds of less than 3 degrees/s. CONCLUSIONS: The method reduces the noise in the reconstruction and increases the image quality. More data are introduced for each phase-correlated reconstruction, and therefore the applied dose is used more efficiently.

A modified McKinnon-Bates (MKB) algorithm for improved 4D cone-beam computed tomography (CBCT) of the lung.

PURPOSE: Four-dimensional (4D) cone-beam computed tomography (CBCT) of the lung is an effective tool for motion management in radiotherapy but presents a challenge because of slow gantry rotation times. Sorting the individual projections by breathing phase and using an established technique such as Feldkamp-Davis-Kress (FDK) to generate corresponding phase-correlated (PC) three-dimensional (3D) images results in reconstructions (FDK-PC) that often contain severe streaking artifacts due to the sparse angular sampling distributions. These can be reduced by further slowing down the gantry at the expense of incurring unwanted increases in scan times and dose. A computationally efficient alternative is the McKinnon-Bates (MKB) reconstruction algorithm that has shown promise in reducing view aliasing-induced streaking but can produce ghosting artifacts that reduce contrast and impede the determination of motion trajectories. The purpose of this work was to identify and correct shortcomings in the MKB algorithm. METHODS: In the general MKB approach, a time-averaged 3D prior image is first reconstructed. The prior is then forward-projected at the same angles as the original projection data creating time-averaged reprojections. These reprojections are subsequently subtracted from the original (unblurred) projections to create motion-encoded difference projections. The difference projections are reconstructed into PC difference images that are added to the well-sampled 3D prior to create the higher quality 4D image. The cause of the ghosting in the traditional 4D MKB images was studied and traced to motion-induced streaking in the prior that, when reprojected, has the undesirable effect of re-encoding for motion in what should be a purely time-averaged reprojection. A new method, designated as the modified McKinnon-Bates (mMKB) algorithm, was developed based on destreaking the prior. This was coupled with a postprocessing 4D bilateral filter for noise suppression and edge preservation (mMKBbf ). The algorithms were tested with the 4D XCAT phantom using four simulated scan times (57, 60, 120, 180 s) and with two in vivo thorax studies (acquisition time of 60 and 90 s). Contrast-to-noise ratios (CNRs) of the target lesions and overall visual quality of the images were assessed. RESULTS: Prior destreaking (mMKB algorithm) reduced ghosting artifacts and increased CNRs for all cases, with the biggest impacts seen in the end inhale (EI) and end exhale (EE) phases of the respiratory cycle. For the XCAT phantom, mMKB lesion CNR was 44% higher than the MKB lesion CNR and was 81% higher than the FDK-PC lesion CNR (EI and EE phases). The bilateral filter provided a further average CNR improvement of 87% with the highest increases associated with longer scan times. Across all phases and scan times, the maximum mMKBbf -to-FDK-PC CNR improvement was over 300%. In vivo results agreed with XCAT results. Significantly less ghosting was observed throughout the mMKB images including near the lesions-of-interest and the diaphragm allowing for, in one case, visualization of a small tumor with nearly 30 mm of motion. The maximum FDK-PC-to-MKBbf CNR improvement for Patient 1's lesion was 261% and for Patient 2's lesion was 318%. CONCLUSIONS: The 4D mMKB algorithm yields good quality coronal and sagittal images in the thorax that may provide sufficient information for patient verification.

Empirical dual energy calibration (EDEC) for cone-beam computed tomography.

Material-selective imaging using dual energy CT (DECT) relies heavily on well-calibrated material decomposition functions. These require the precise knowledge of the detected x-ray spectra, and even if they are exactly known the reliability of DECT will suffer from scattered radiation. We propose an empirical method to determine the proper decomposition function. In contrast to other decomposition algorithms our empirical dual energy calibration (EDEC) technique requires neither knowledge of the spectra nor of the attenuation coefficients. The desired material-selective raw data p1 and p2 are obtained as functions of the measured attenuation data q1 and q2 (one DECT scan = two raw data sets) by passing them through a polynomial function. The polynomial's coefficients are determined using a general least squares fit based on thresholded images of a calibration phantom. The calibration phantom's dimension should be of the same order of magnitude as the test object, but other than that no assumptions on its exact size or positioning are made. Once the decomposition coefficients are determined DECT raw data can be decomposed by simply passing them through the polynomial. To demonstrate EDEC simulations of an oval CTDI phantom, a lung phantom, a thorax phantom and a mouse phantom were carried out. The method was further verified by measuring a physical mouse phantom, a half-and-half-cylinder phantom and a Yin-Yang phantom with a dedicated in vivo dual source micro-CT scanner. The raw data were decomposed into their components, reconstructed, and the pixel values obtained were compared to the theoretical values. The determination of the calibration coefficients with EDEC is very robust and depends only slightly on the type of calibration phantom used. The images of the test phantoms (simulations and measurements) show a nearly perfect agreement with the theoretical micro values and density values. Since EDEC is an empirical technique it inherently compensates for scatter components. The empirical dual energy calibration technique is a pragmatic, simple, and reliable calibration approach that produces highly quantitative DECT images.