Streaking artifacts reduction in four-dimensional cone-beam computed tomography

A review on 4D cone-beam CT (4D-CBCT) in radiation therapy: Technical advances and clinical applications

Cone-beam CT (CBCT) is the most commonly used onboard imaging technique for target localization in radiation therapy. Conventional 3D CBCT acquires x-ray cone-beam projections at multiple angles around the patient to reconstruct 3D images of the patient in the treatment room. However, despite its wide usage, 3D CBCT is limited in imaging disease sites affected by respiratory motions or other dynamic changes within the body, as it lacks time-resolved information. To overcome this limitation, 4D-CBCT was developed to incorporate a time dimension in the imaging to account for the patient's motion during the acquisitions. For example, respiration-correlated 4D-CBCT divides the breathing cycles into different phase bins and reconstructs 3D images for each phase bin, ultimately generating a complete set of 4D images. 4D-CBCT is valuable for localizing tumors in the thoracic and abdominal regions where the localization accuracy is affected by respiratory motions. This is especially important for hypofractionated stereotactic body radiation therapy (SBRT), which delivers much higher fractional doses in fewer fractions than conventional fractionated treatments. Nonetheless, 4D-CBCT does face certain limitations, including long scanning times, high imaging doses, and compromised image quality due to the necessity of acquiring sufficient x-ray projections for each respiratory phase. In order to address these challenges, numerous methods have been developed to achieve fast, low-dose, and high-quality 4D-CBCT. This paper aims to review the technical developments surrounding 4D-CBCT comprehensively. It will explore conventional algorithms and recent deep learning-based approaches, delving into their capabilities and limitations. Additionally, the paper will discuss the potential clinical applications of 4D-CBCT and outline a future roadmap, highlighting areas for further research and development. Through this exploration, the readers will better understand 4D-CBCT's capabilities and potential to enhance radiation therapy.

Surrogate-driven respiratory motion model for projection-resolved motion estimation and motion compensated cone-beam CT reconstruction from unsorted projection data

Objective.As the most common solution to motion artefact for cone-beam CT (CBCT) in radiotherapy, 4DCBCT suffers from long acquisition time and phase sorting error. This issue could be addressed if the motion at each projection could be known, which is a severely ill-posed problem. This study aims to obtain the motion at each time point and motion-free image simultaneously from unsorted projection data of a standard 3DCBCT scan.Approach.Respiration surrogate signals were extracted by the Intensity Analysis method. A general framework was then deployed to fit a surrogate-driven motion model that characterized the relation between the motion and surrogate signals at each time point. Motion model fitting and motion compensated reconstruction were alternatively and iteratively performed. Stochastic subset gradient based method was used to significantly reduce the computation time. The performance of our method was comprehensively evaluated through digital phantom simulation and also validated on clinical scans from six patients.Results.For digital phantom experiments, motion models fitted with ground-truth or extracted surrogate signals both achieved a much lower motion estimation error and higher image quality, compared with non motion-compensated results.For the public SPARE Challenge datasets, more clear lung tissues and less blurry diaphragm could be seen in the motion compensated reconstruction, comparable to the benchmark 4DCBCT images but with a higher temporal resolution. Similar results were observed for two real clinical 3DCBCT scans.Significance.The motion compensated reconstructions and motion models produced by our method will have direct clinical benefit by providing more accurate estimates of the delivered dose and ultimately facilitating more accurate radiotherapy treatments for lung cancer patients.

A deep unsupervised learning framework for the 4D CBCT artifact correction

Objective. Four-dimensional cone-beam computed tomography (4D CBCT) has unique advantages in moving target localization, tracking and therapeutic dose accumulation in adaptive radiotherapy. However, the severe fringe artifacts and noise degradation caused by 4D CBCT reconstruction restrict its clinical application. We propose a novel deep unsupervised learning model to generate the high-quality 4D CBCT from the poor-quality 4D CBCT. Approach. The proposed model uses a contrastive loss function to preserve the anatomical structure in the corrected image. To preserve the relationship between the input and output image, we use a multilayer, patch-based method rather than operate on entire images. Furthermore, we draw negatives from within the input 4D CBCT rather than from the rest of the dataset. Main results. The results showed that the streak and motion artifacts were significantly suppressed. The spatial resolution of the pulmonary vessels and microstructure were also improved. To demonstrate the results in the different directions, we make the animation to show the different views of the predicted correction image in the supplementary animation. Significance. The proposed method can be integrated into any 4D CBCT reconstruction method and maybe a practical way to enhance the image quality of the 4D CBCT.

Spatiotemporal structure-aware dictionary learning-based 4D CBCT reconstruction

PURPOSE

Four-dimensional cone-beam computed tomography (4D CBCT) is developed to reconstruct a sequence of phase-resolved images, which could assist in verifying the patient's position and offering information for cancer treatment planning. However, 4D CBCT images suffer from severe streaking artifacts and noise due to the extreme sparse-view CT reconstruction problem for each phase. As a result, it would cause inaccuracy of treatment estimation. The purpose of this paper was to develop a new 4D CBCT reconstruction method to generate a series of high spatiotemporal 4D CBCT images.

METHODS

Considering the advantage of (DL) on representing structural features and correlation between neighboring pixels effectively, we construct a novel DL-based method for the 4D CBCT reconstruction. In this study, both a motion-aware dictionary and a spatially structural 2D dictionary are trained for 4D CBCT by excavating the spatiotemporal correlation among ten phase-resolved images and the spatial information in each image, respectively. Specifically, two reconstruction models are produced in this study. The first one is the motion-aware dictionary learning-based 4D CBCT algorithm, called motion-aware DL based 4D CBCT (MaDL). The second one is the MaDL equipped with a prior knowledge constraint, called pMaDL. Qualitative and quantitative evaluations are performed using a 4D extended cardiac torso (XCAT) phantom, simulated patient data, and two sets of patient data sets. Several state-of-the-art 4D CBCT algorithms, such as the McKinnon-Bates (MKB) algorithm, prior image constrained compressed sensing (PICCS), and the high-quality initial image-guided 4D CBCT reconstruction method (HQI-4DCBCT) are applied for comparison to validate the performance of the proposed MaDL and prior constraint MaDL (pMaDL) pmadl reconstruction frameworks.

RESULTS

Experimental results validate that the proposed MaDL can output the reconstructions with few streaking artifacts but some structural information such as tumors and blood vessels, may still be missed. Meanwhile, the results of the proposed pMaDL demonstrate an improved spatiotemporal resolution of the reconstructed 4D CBCT images. In these improved 4D CBCT reconstructions, streaking artifacts are suppressed primarily and detailed structures are also restored. Regarding the XCAT phantom, quantitative evaluations indicate that an average of 58.70%, 45.25%, and 40.10% decrease in terms of root-mean-square error (RMSE) and an average of 2.10, 1.37, and 1.37 times in terms of structural similarity index (SSIM) are achieved by the proposed pMaDL method when compared with piccs, PICCS, MaDL(2D), and MaDL(2D), respectively. Moreover the proposed pMaDL achieves a comparable performance with HQI-4DCBCT algorithm in terms of RMSE and SSIM metrics. However, pMaDL has a better ability to suppress streaking artifacts than HQI-4DCBCT.

CONCLUSIONS

The proposed algorithm could reconstruct a set of 4D CBCT images with both high spatiotemporal resolution and detailed features preservation. Moreover the proposed pMaDL can effectively suppress the streaking artifacts in the resultant reconstructions, while achieving an overall improved spatiotemporal resolution by incorporating the motion-aware dictionary with a prior constraint into the proposed 4D CBCT iterative framework.

Accurate and robust sparse-view angle CT image reconstruction using deep learning and prior image constrained compressed sensing (DL-PICCS)

BACKGROUND

Sparse-view CT image reconstruction problems encountered in dynamic CT acquisitions are technically challenging. Recently, many deep learning strategies have been proposed to reconstruct CT images from sparse-view angle acquisitions showing promising results. However, two fundamental problems with these deep learning reconstruction methods remain to be addressed: (1) limited reconstruction accuracy for individual patients and (2) limited generalizability for patient statistical cohorts.

PURPOSE

The purpose of this work is to address the previously mentioned challenges in current deep learning methods.

METHODS

A method that combines a deep learning strategy with prior image constrained compressed sensing (PICCS) was developed to address these two problems. In this method, the sparse-view CT data were reconstructed by the conventional filtered backprojection (FBP) method first, and then processed by the trained deep neural network to eliminate streaking artifacts. The outputs of the deep learning architecture were then used as the needed prior image in PICCS to reconstruct the image. If the noise level from the PICCS reconstruction is not satisfactory, another light duty deep neural network can then be used to reduce noise level. Both extensive numerical simulation data and human subject data have been used to quantitatively and qualitatively assess the performance of the proposed DL-PICCS method in terms of reconstruction accuracy and generalizability.

RESULTS

Extensive evaluation studies have demonstrated that: (1) quantitative reconstruction accuracy of DL-PICCS for individual patient is improved when it is compared with the deep learning methods and CS-based methods; (2) the false-positive lesion-like structures and false negative missing anatomical structures in the deep learning approaches can be effectively eliminated in the DL-PICCS reconstructed images; and (3) DL-PICCS enables a deep learning scheme to relax its working conditions to enhance its generalizability.

CONCLUSIONS

DL-PICCS offers a promising opportunity to achieve personalized reconstruction with improved reconstruction accuracy and enhanced generalizability.

Fast and effective single-scan dual-energy cone-beam CT reconstruction and decomposition denoising based on dual-energy vectorization

PURPOSE

Flat-panel detector (FPD) based dual-energy cone-beam computed tomography (DE-CBCT) is a promising imaging technique for dedicated clinical applications. In this paper, we proposed a fully analytical method for fast and effective single-scan DE-CBCT image reconstruction and decomposition.

METHODS

A rotatable Mo filter was inserted between an x-ray source and imaged object to alternately produce low and high-energy x-ray spectra. First, filtered-backprojection (FBP) method was applied on down-sampled projections to reconstruct low and high-energy images. Then, the two images were converted into a vectorized form represented with an amplitude and an argument image. Using amplitude image as a guide, a joint bilateral filter was applied to denoise the argument image. Then, high-quality dual-energy images were recovered from the amplitude image and the denoised argument image. Finally, the recovered dual-energy images were further used for low-noise material decomposition and electron density synthesis. Imaging was conducted on a Catphan^® 600 phantom and an anthropomorphic head phantom. The proposed method was evaluated via comparison with the traditional two-scan method and a commonly used filtering method (HYPR-LR).

RESULTS

On the Catphan^® 600 phantom, the proposed method successfully reduced streaking artifacts and preserved spatial resolution and noise-power-spectrum (NPS) pattern. In the electron density image, the proposed method increased contrast-to-noise ratio (CNR) by more than 2.5 times and achieved <1.2% error for electron density values. On the anthropomorphic head phantom, the proposed method greatly improved the soft-tissue contrast and the fine detail differentiation ability. In the selected ROIs on different human tissues, the differences between the CT number obtained by the proposed method and that by the two-scan method were less than 4 HU. In the material images, the proposed method suppressed noise by over 75.5% compared with two-scan results, and by over 40.4% compared with HYPR-LR results. Implementation of the whole algorithm took 44.5 s for volumetric imaging, including projection preprocessing, FBP reconstruction, joint bilateral filtering, and material decomposition.

CONCLUSIONS

Using down-sampled projections in single-scan DE-CBCT, the proposed method could effectively and efficiently produce high-quality DE-CBCT images and low-noise material decomposition images. This method demonstrated superior performance on spatial resolution enhancement, NPS preservation, noise reduction, and electron density accuracy, indicating better prospect in material differentiation and dose calculation.

Four-dimensional inverse-geometry computed tomography: a preliminary study

This study introduces and evaluates respiratory-correlated four-dimensional (4D) inverse geometry computed tomography (IGCT). The projection data of the IGCT were acquired in a single gantry rotation over 120 s. Three virtual phantoms-static Defrise, 4D Shepp-Logan, and 4D extended cardiac-torso (XCAT)-were used to obtain projection data for the IGCT and cone-beam computed tomography (CBCT). The projection acquisition parameters were determined to eliminate vacancies in the Radon space for an accurate rebinning process. Phase-based sorting was conducted within 10 phase bins, and the sorted projection data were binned into a cone beam geometry. Finally, Feldkamp-Davis-Kress reconstruction was conducted independently at each phase. The reconstructed images were compared using the structural similarity index measure (SSIM) and root mean square error (RMSE). The vertical profile of the Defrise reconstruction image was uniform, and the cone beam artefact was reduced in the IGCT image. Under an ideal projection acquisition condition, the mean coronal plane SSIMs of the Shepp-Logan and 4D XCAT phantoms were 0.899 and 0.706, respectively, which were higher than those of the CBCT (0.784 and 0.623, respectively). Similarly, the mean RMSEs of the coronal plane IGCT (0.036 and 0.158) exhibited an improvement over those of the CBCT (0.165 and 0.261, respectively). The mean standard deviations of the SSIM and RMSE were lower for IGCT than for CBCT. In particular, the SSIM and RMSE of the sagittal and coronal planes of the Shepp-Logan IGCT images were stable in all phase bins; however, those of the CBCT changed depending on the phase bins. Poor image quality was observed for IGCT under inappropriate conditions. This was caused by a vacancy in the Radon space, owing to an inappropriate scan setting. Overall, the proposed 4D IGCT exhibited better image quality than conventional CBCT.

[Four-dimensional cone-beam CT reconstruction based on motion-compensated robust principal component analysis]

OBJECTIVE

To propose a motion compensation reconstruction method based on robust principal component analysis (RPCA) to reduce the influence of streak artifacts on accurate estimation of interphase motion deformation fields.

OBJECTIVE

We propose a RPCA motion compensation reconstruction algorithm to improve the estimation of motion deformation fields based on the traditional MC-FDK algorithm. RPCA was used to decompose the cone-beam computed tomography (CBCT) images into low-rank and sparse components, and the motion deformation fields between different phase images were then estimated using Horn and Schunck optical flow method from the low-rank images to reduce the influence of striping artifacts on the accuracy of estimation of interphase motion deformation fields. The performance of the algorithm was evaluated using simulation data and real data. The simulation phantom data was obtained by back-projection of 4D-CT images acquired from Philips 16-slice spiral CT using MATLAB software programming according to the scanning geometry of Varian Edge accelerator. The real patient data were obtained using the Elekta Synergy system of CBCT scanning system with half-fan mode CB projection data from lung cancer patients.

OBJECTIVE

Compared with images reconstructed using the traditional MC-FDK algorithm, the reconstructed image using the proposed method had clearer tissue boundaries with reduced motion artifact was reduced. The results of phantom data reconstruction showed that compared with the MC- FDK algorithm, the proposed algorithms resulted in improvements of PSNR by 25.4% and SSIM by 7.6%; compared with the FDK algorithm, PSNR was improved by 37.9% and SSIM by 17.6%.

OBJECTIVE

The proposed algorithm can achieve accurate estimation of inter-phase motion deformation fields and improve the quality of the reconstructed CBCT images.

Fully Automatic Sliding Motion Compensated and Simultaneous 4D-CBCT via Bilateral Filtering

Purpose

To incorporate the bilateral filtering into the Deformable Vector Field (DVF) based 4D-CBCT reconstruction for realizing a fully automatic sliding motion compensated 4D-CBCT.

Materials and Methods

Initially, a motion compensated simultaneous algebraic reconstruction technique (mSART) is used to generate a high quality reference phase (e.g. 0% phase) by using all phase projections together with the initial 4D-DVFs. The initial 4D-DVF were generated via Demons registration between 0% phase and each other phase image. The 4D-DVF will then kept updating by matching the forward projection of the deformed high quality 0% phase with the measured projection of the target phase. The loss function during this optimization contains an projection intensity difference matching criterion plus a DVF smoothing constrain term. We introduce a bilateral filtering kernel into the DVF constrain term to estimate the sliding motion automatically. The bilateral filtering kernel contains three sub-kernels: 1) an spatial domain Guassian kernel; 2) an image intensity domain Guassian kernel; and 3) a DVF domain Guassian kernel. By choosing suitable kernel variances, the sliding motion can be extracted. A non-linear conjugate gradient optimizer was used. We validated the algorithm on a non-uniform rotational B-spline based cardiac-torso (NCAT) phantom and four anonymous patient data. For quantification, we used: 1) the Root-Mean-Square-Error (RMSE) together with the Maximum-Error (MaxE); 2) the Dice coefficient of the extracted lung contour from the final reconstructed images and 3) the relative reconstruction error (RE) to evaluate the algorithm's performance.

Results

For NCAT phantom, the motion trajectory's RMSE/MaxE are 0.796/1.02 mm for bilateral filtering reconstruction; and 2.704/4.08 mm for original reconstruction. For patient pilot study, the 4D-Dice coefficient obtained with bilateral filtering are consistently higher than that without bilateral filtering. Meantime several image content such as the rib position, the heart edge definition, the fibrous structures all has been better corrected with bilateral filtering.

Conclusion

We developed a bilateral filtering based fully automatic sliding motion compensated 4D-CBCT scheme. Both digital phantom and initial patient pilot studies confirmed the improved motion estimation and image reconstruction ability. It can be used as a 4D-CBCT image guidance tool for lung SBRT treatment.

Evaluation of a Motion Correction Algorithm for C-Arm Computed Tomography Acquired During Transarterial Chemoembolization

Purpose

The aim of this retrospective study was to evaluate the feasibility of a motion correction 3D reconstruction prototype technique for C-arm computed tomography (CACT).

Material and Methods

We included 65 consecutive CACTs acquired during transarterial chemoembolization of 54 patients (47 m,7f; 67 ± 11.3 years). All original raw datasets (CACT_Org) underwent reconstruction with and without volume punching of high-contrast objects using a 3D image reconstruction software to compensate for motion (CACT_{MC_bone};CACT_{MC_no bone}). Subsequently, the effect on image quality (IQ) was evaluated using objective (image sharpness metric) and subjective criteria. Subjective criteria were defined by vessel geometry, overall IQ, delineation of tumor feeders, the presence of foreign material-induced artifacts and need for additional imaging, assessed by two independent readers on a 3-(vessel geometry and overall IQ) or 2-point scale, respectively. Friedman rank-sum test and post hoc analysis in form of pairwise Wilcoxon signed-rank test were computed and inter-observer agreement analyzed using kappa test.

Results

Objective IQ as defined by an image sharpness metric, increased from 273.5 ± 28 (CACT_Org) to 328.5 ± 55.1 (CACT_{MC_bone}) and 331 ± 57.8 (CACT_{MC_no bone}; all p < 0.0001). These results could largely be confirmed by the subjective analysis, which demonstrated predominantly good and moderate inter-observer agreement, with best agreement for CACT_{MC_no bone} in all categories (e.g., vessel geometry: CACT_Org: κ = 0.51, CACT_{MC_bone}: κ = 0.42, CACT_{MC_no bone}: κ = 0.69).

Conclusion

The application of a motion correction algorithm was feasible for all data sets and led to an increase in both objective and subjective IQ parameters.

Level of Evidence

3

A modern review of the uncertainties in volumetric imaging of respiratory-induced target motion in lung radiotherapy

Radiotherapy has become a critical component for the treatment of all stages and types of lung cancer, often times being the primary gateway to a cure. However, given that radiation can cause harmful side effects depending on how much surrounding healthy tissue is exposed, treatment of the lung can be particularly challenging due to the presence of moving targets. Careful implementation of every step in the radiotherapy process is absolutely integral for attaining optimal clinical outcomes. With the advent and now widespread use of stereotactic body radiation therapy (SBRT), where extremely large doses are delivered, accurate, and precise dose targeting is especially vital to achieve an optimal risk to benefit ratio. This has largely become possible due to the rapid development of image-guided technology. Although imaging is critical to the success of radiotherapy, it can often be plagued with uncertainties due to respiratory-induced target motion. There has and continues to be an immense research effort aimed at acknowledging and addressing these uncertainties to further our abilities to more precisely target radiation treatment. Thus, the goal of this article is to provide a detailed review of the prevailing uncertainties that remain to be investigated across the different imaging modalities, as well as to highlight the more modern solutions to imaging motion and their role in addressing the current challenges.

SPARE: Sparse-view reconstruction challenge for 4D cone-beam CT from a 1-min scan

PURPOSE

Currently, four-dimensional (4D) cone-beam computed tomography (CBCT) requires a 3-4 min full-fan scan to ensure usable image quality. Recent advancements in sparse-view 4D-CBCT reconstruction have opened the possibility to reduce scan time and dose. The aim of this study is to provide a common framework for systematically evaluating algorithms for 4D-CBCT reconstruction from a 1-min scan. Using this framework, the AAPM-sponsored SPARE Challenge was conducted in 2018 to identify and compare state-of-the-art algorithms.

METHODS

A clinically realistic CBCT dataset was simulated using patient CT volumes from the 4D-Lung database. The selected patients had multiple 4D-CT sessions, where the first 4D-CT was used as the prior CT, and the rest were used as the ground truth volumes for simulating CBCT projections. A GPU-based Monte Carlo tool was used to simulate the primary, scatter, and quantum noise signals. A total of 32 CBCT scans of nine patients were generated. Additional qualitative analysis was performed on a clinical Varian and clinical Elekta dataset to validate the simulation study. Participants were blinded from the ground truth, and were given 3 months to apply their reconstruction algorithms to the projection data. The submitted reconstructions were analyzed in terms of root-mean-squared-error (RMSE) and structural similarity index (SSIM) with the ground truth within four different region-of-interests (ROI) - patient body, lungs, planning target volume (PTV), and bony anatomy. Geometric accuracy was quantified as the alignment error of the PTV.

RESULTS

Twenty teams participated in the challenge, with five teams completing the challenge. Techniques involved in the five methods included iterative optimization, motion-compensation, and deformation of the prior 4D-CT. All five methods rendered significant reduction in noise and streaking artifacts when compared to the conventional Feldkamp-Davis-Kress (FDK) algorithm. The RMS of the three-dimensional (3D) target registration error of the five methods ranged from 1.79 to 3.00 mm. Qualitative observations from the Varian and Elekta datasets mostly concur with those from the simulation dataset. Each of the methods was found to have its own strengths and weaknesses. Overall, the MA-ROOSTER method, which utilizes a 4D-CT motion model for temporal regularization, had the best and most consistent image quality and accuracy.

CONCLUSION

The SPARE Challenge represents the first framework for systematically evaluating state-of-the-art algorithms for 4D-CBCT reconstruction from a 1-min scan. Results suggest the potential for reducing scan time and dose for 4D-CBCT. The challenge dataset and analysis framework are publicly available for benchmarking future reconstruction algorithms.

Quantitative evaluation of 4D Cone beam CT scans with reduced scan time in lung cancer patients

•

Fast (2 min) 4D CBCT can be simulated accurately from long (4 min) scans.

•

Registration was accurate for 96.6% of simulated 2 min scans.

•

Acquired 2 min scan registration was accurate in 6/8 patients.

•

2 min 4D CBCT produces sufficient image quality for IGRT in lung cancer patients.

Purpose

Image guided radiotherapy (IGRT) based on respiration correlated cone-beam CT (4D-CBCT) provides accurate tumour localisation in lung cancer patients by taking into account respiratory motion when deriving setup correction. However, 4D-CBCT scan times are typically longer than for acquisition of 3D-CBCT scans, e.g. 4 min. This work aims to quantitatively evaluate the effect of reduced scan times on 4D-CBCT image quality and registration accuracy in lung cancer patients.

Methods and materials

Scan times down to 1 min were simulated by retaining only projection images corresponding to every second, third or fourth respiratory cycle in forty-four 4D-CBCTs from 15 lung cancer patients. In addition twenty 2-minute scans were acquired for 12 lung cancer patients. Image quality was quantified by assessing registration accuracy in the shorter scan times, comparing to the 4-minute scan registration result where available as reference.

Results

Use of 2-minute scans had little impact on registration accuracy or ability to detect tumour motion: automatic registration accuracy was within 2 mm in 6/8 scans analysed with 2-minute acquisitions, and 96.6% of registration discrepancies were within 2 mm for the simulated scans. When the scan time simulated was below 2 min, automatic registration results still agreed within 2 mm for 84.7% of scans, however visual image quality was considerably degraded.

Conclusion

A 4D-CBCT acquisition time of 2 min produces scans of sufficient image quality for IGRT in most lung cancer patients, as demonstrated quantitatively by assessing the impact on automatic registration accuracy in simulated and real acquisitions.

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.

A fast 4D cone beam CT reconstruction method based on the OSC-TV algorithm

BACKGROUND

Four-dimensional cone beam computed tomography allows for temporally resolved imaging with useful applications in radiotherapy, but raises particular challenges in terms of image quality and computation time.

OBJECTIVE

The purpose of this work is to develop a fast and accurate 4D algorithm by adapting a GPU-accelerated ordered subsets convex algorithm (OSC), combined with the total variation minimization regularization technique (TV).

METHODS

Different initialization schemes were studied to adapt the OSC-TV algorithm to 4D reconstruction: each respiratory phase was initialized either with a 3D reconstruction or a blank image. Reconstruction algorithms were tested on a dynamic numerical phantom and on a clinical dataset. 4D iterations were implemented for a cluster of 8 GPUs.

RESULTS

All developed methods allowed for an adequate visualization of the respiratory movement and compared favorably to the McKinnon-Bates and adaptive steepest descent projection onto convex sets algorithms, while the 4D reconstructions initialized from a prior 3D reconstruction led to better overall image quality.

CONCLUSION

The most suitable adaptation of OSC-TV to 4D CBCT was found to be a combination of a prior FDK reconstruction and a 4D OSC-TV reconstruction with a reconstruction time of 4.5 minutes. This relatively short reconstruction time could facilitate a clinical use.

Motion guided Spatiotemporal Sparsity for high quality 4D-CBCT reconstruction

Conventional cone-beam computed tomography is often deteriorated by respiratory motion blur, which negatively affects target delineation. On the other side, the four dimensional cone-beam computed tomography (4D-CBCT) can be considered to describe tumor and organ motion. But for current on-board CBCT imaging system, the slow rotation speed limits the projection number at each phase, and the associated reconstructions are contaminated by noise and streak artifacts using the conventional algorithm. To address the problem, we propose a novel framework to reconstruct 4D-CBCT from the under-sampled measurements—Motion guided Spatiotemporal Sparsity (MgSS). In this algorithm, we try to divide the CBCT images at each phase into cubes (3D blocks) and track the cubes with estimated motion field vectors through phase, then apply regional spatiotemporal sparsity on the tracked cubes. Specifically, we recast the tracked cubes into four-dimensional matrix, and use the higher order singular value decomposition (HOSVD) technique to analyze the regional spatiotemporal sparsity. Subsequently, the blocky spatiotemporal sparsity is incorporated into a cost function for the image reconstruction. The phantom simulation and real patient data are used to evaluate this algorithm. Results show that the MgSS algorithm achieved improved 4D-CBCT image quality with less noise and artifacts compared to the conventional algorithms.

A Novel Method of Cone Beam CT Projection Binning Based on Image Registration

Accurate sorting of beam projections is important in 4D cone beam computed tomography (4D CBCT) to improve the quality of the reconstructed 4D CBCT image by removing motion-induced artifacts. We propose image registration-based projection binning (IRPB), a novel marker-less binning method for 4D CBCT projections, which combines intensity-based feature point detection and trajectory tracking using random sample consensus. IRPB extracts breathing motion and phases by analyzing tissue feature point trajectories. We conducted experiments with two phantom and six patient datasets, including both regular and irregular respirations. In experiments, we compared the performance of the proposed IRPB, Amsterdam Shroud method (AS), Fourier transform-based method (FT), and local intensity feature tracking method (LIFT). The results showed that the average absolute phase shift of IRPB was 3.74 projections and 0.48 projections less than that of FT and LIFT, respectively. AS lost the most breathing cycles in the respiration extraction for the five patient datasets, so we could not compare the average absolute phase shift between IRPB and AS. Based on the peak signal-to-noise ratio (PSNR) of the reconstructed 4D CBCT images, IRPB had 5.08, 1.05, and 2.90 dB larger PSNR than AS, FT, and LIFT, respectively. The average Structure SIMilarity Index (SSIM) of the 4D CBCT image reconstructed by IRPB, AS, and LIFT were 0.87, 0.74, 0.84, and 0.70, respectively. These results demonstrated that IRPB has superior performance to the other standard methods.

High quality 4D cone-beam CT reconstruction using motion-compensated total variation regularization

Four dimensional cone-beam computed tomography (4D-CBCT) has great potential clinical value because of its ability to describe tumor and organ motion. But the challenge in 4D-CBCT reconstruction is the limited number of projections at each phase, which result in a reconstruction full of noise and streak artifacts with the conventional analytical algorithms. To address this problem, in this paper, we propose a motion compensated total variation regularization approach which tries to fully explore the temporal coherence of the spatial structures among the 4D-CBCT phases. In this work, we additionally conduct motion estimation/motion compensation (ME/MC) on the 4D-CBCT volume by using inter-phase deformation vector fields (DVFs). The motion compensated 4D-CBCT volume is then viewed as a pseudo-static sequence, of which the regularization function was imposed on. The regularization used in this work is the 3D spatial total variation minimization combined with 1D temporal total variation minimization. We subsequently construct a cost function for a reconstruction pass, and minimize this cost function using a variable splitting algorithm. Simulation and real patient data were used to evaluate the proposed algorithm. Results show that the introduction of additional temporal correlation along the phase direction can improve the 4D-CBCT image quality.

Evaluation of three presets for four-dimensional cone beam CT in lung radiotherapy verification by visual grading analysis

OBJECTIVE

To evaluate three image acquisition presets for four-dimensional cone beam CT (CBCT) to identify an optimal preset for lung tumour image quality while minimizing dose and acquisition time.

METHODS

Nine patients undergoing radical conventionally fractionated radiotherapy for lung cancer had verification CBCTs acquired using three presets: Preset 1 on Day 1 (11 mGy dose, 240 s acquisition time), Preset 2 on Day 2 (9 mGy dose, 133 s acquisition time) and Preset 3 on Day 3 (9 mGy dose, 67 s acquisition time). The clarity of the tumour and other thoracic structures, and the acceptability of the match, were retrospectively graded by visual grading analysis (VGA). Logistic regression was used to identify the most appropriate preset and any factors that might influence the result.

RESULTS

Presets 1 and 2 met a clinical requirement of 75% of structures to be rated "Clear" or above and 75% of matches to be rated "Acceptable" or above. Clarity is significantly affected by preset, patient, observer and structure. Match acceptability is significantly affected by preset.

CONCLUSION

The application of VGA in this initial study enabled a provisional selection of an optimal preset (Preset 2) to be made.

ADVANCES IN KNOWLEDGE

This was the first application of VGA to the investigation of presets for CBCT.

Measuring interfraction and intrafraction lung function changes during radiation therapy using four-dimensional cone beam CT ventilation imaging

PURPOSE

Adaptive ventilation guided radiation therapy could minimize the irradiation of healthy lung based on repeat lung ventilation imaging (VI) during treatment. However the efficacy of adaptive ventilation guidance requires that interfraction (e.g., week-to-week), ventilation changes are not washed out by intrafraction (e.g., pre- and postfraction) changes, for example, due to patient breathing variability. The authors hypothesize that patients undergoing lung cancer radiation therapy exhibit larger interfraction ventilation changes compared to intrafraction function changes. To test this, the authors perform the first comparison of interfraction and intrafraction lung VI pairs using four-dimensional cone beam CT ventilation imaging (4D-CBCT VI), a novel technique for functional lung imaging.

METHODS

The authors analyzed a total of 215 4D-CBCT scans acquired for 19 locally advanced non-small cell lung cancer (LA-NSCLC) patients over 4-6 weeks of radiation therapy. This set of 215 scans was sorted into 56 interfraction pairs (including first day scans and each of treatment weeks 2, 4, and 6) and 78 intrafraction pairs (including pre/postfraction scans on the same-day), with some scans appearing in both sets. VIs were obtained from the Jacobian determinant of the transform between the 4D-CBCT end-exhale and end-inhale images after deformable image registration. All VIs were deformably registered to their corresponding planning CT and normalized to account for differences in breathing effort, thus facilitating image comparison in terms of (i) voxelwise Spearman correlations, (ii) mean image differences, and (iii) gamma pass rates for all interfraction and intrafraction VI pairs. For the side of the lung ipsilateral to the tumor, we applied two-sided t-tests to determine whether interfraction VI pairs were more different than intrafraction VI pairs.

RESULTS

The (mean ± standard deviation) Spearman correlation for interfraction VI pairs was r̄(Inter)=0.52±0.25, which was significantly lower than for intrafraction pairs (r̄(Intra)=0.67±0.20, p = 0.0002). Conversely, mean absolute ventilation differences were larger for interfraction pairs than for intrafraction pairs, with |ΔV̄(Inter)|=0.42±0.65 and |ΔV̄(Intra)|=0.32±0.53, respectively (p < 10(-15)). Applying a gamma analysis with ventilation/distance tolerance of 25%/10 mm, we observed mean pass rate of (69% ± 20%) for interfraction VIs, which was significantly lower compared to intrafraction pairs (80% ± 15%, with p ∼ 0.0003). Compared to the first day scans, all patients experienced at least one subsequent change in median ipsilateral ventilation ≥10%. Patients experienced both positive and negative ventilation changes throughout treatment, with the maximum change occurring at different weeks for different patients.

CONCLUSIONS

The authors' data support the hypothesis that interfraction ventilation changes are larger than intrafraction ventilation changes for LA-NSCLC patients over a course of conventional lung cancer radiation therapy. Longitudinal ventilation changes are observed to be highly patient-dependent, supporting a possible role for adaptive ventilation guidance based on repeat 4D-CBCT VIs. We anticipate that future improvement of 4D-CBCT image reconstruction algorithms will improve the capability of 4D-CBCT VI to resolve interfraction ventilation changes.

Generation of hybrid sinograms for the recovery of kV-CT images with metal artifacts for helical tomotherapy

PURPOSE

The overall goal of this study is to restore kilovoltage computed tomography (kV-CT) images which are disfigured by patients' metal prostheses. By generating a hybrid sinogram that is a combination of kV and megavoltage (MV) projection data, the authors suggest a novel metal artifact-reduction (MAR) method that retains the image quality to match that of kV-CT and simultaneously restores the information of metal prostheses lost due to photon starvation.

METHODS

CT projection data contain information about attenuation coefficients and the total length of the attenuation. By normalizing raw kV projections with their own total lengths of attenuation, mean attenuation projections were obtained. In the same manner, mean density projections of MV-CT were obtained by the normalization of MV projections resulting from the forward projection of density-calibrated MV-CT images with the geometric parameters of the kV-CT device. To generate the hybrid sinogram, metal-affected signals of the kV sinogram were identified and replaced by the corresponding signals of the MV sinogram following a density calibration step with kV data. Filtered backprojection was implemented to reconstruct the hybrid CT image. To validate the authors' approach, they simulated four different scenarios for three heads and one pelvis using metallic rod inserts within a cylindrical phantom. Five inserts describing human body elements were also included in the phantom. The authors compared the image qualities among the kV, MV, and hybrid CT images by measuring the contrast-to-noise ratio (CNR), the signal-to-noise ratio (SNR), the densities of all inserts, and the spatial resolution. In addition, the MAR performance was compared among three existing MAR methods and the authors' hybrid method. Finally, for clinical trials, the authors produced hybrid images of three patients having dental metal prostheses to compare their MAR performances with those of the kV, MV, and three existing MAR methods.

RESULTS

The authors compared the image quality and MAR performance of the hybrid method with those of other imaging modalities and the three MAR methods, respectively. The total measured mean of the CNR (SNR) values for the nonmetal inserts was determined to be 14.3 (35.3), 15.3 (37.8), and 25.5 (64.3) for the kV, MV, and hybrid images, respectively, and the spatial resolutions of the hybrid images were similar to those of the kV images. The measured densities of the metal and nonmetal inserts in the hybrid images were in good agreement with their true densities, except in cases of extremely low densities, such as air and lung. Using the hybrid method, major streak artifacts were suitably removed and no secondary artifacts were introduced in the resultant image. In clinical trials, the authors verified that kV and MV projections were successfully combined and turned into the resultant hybrid image with high image contrast, accurate metal information, and few metal artifacts. The hybrid method also outperformed the three existing MAR methods with regard to metal information restoration and secondary artifact prevention.

CONCLUSIONS

The authors have shown that the hybrid method can restore the overall image quality of kV-CT disfigured by severe metal artifacts and restore the information of metal prostheses lost due to photon starvation. The hybrid images may allow for the improved delineation of structures of interest and accurate dose calculations for radiation treatment planning for patients with metal prostheses.

A new scheme for real-time high-contrast imaging in lung cancer radiotherapy: a proof-of-concept study

Visualization of anatomy in real time is of critical importance for motion management in lung cancer radiotherapy. To achieve real-time, and high-contrast in-treatment imaging, we propose a novel scheme based on the measurement of Compton scatter photons. In our method, a slit x-ray beam along the superior-inferior direction is directed to the patient, (intersecting the lung region at a 2D plane) containing most of the tumor motion trajectory. X-ray photons are scattered off this plane primarily due to the Compton interaction. An imager with a pinhole or a slat collimator is placed at one side of the plane to capture the scattered photons. The resulting image, after correcting for incoming fluence inhomogeneity, x-ray attenuation, scatter angle variation, and outgoing beam geometry, represents the linear attenuation coefficient of Compton scattering. This allows the visualization of the anatomy on this plane. We performed Monte Carlo simulation studies both on a phantom and a patient for proof-of-principle purposes. In the phantom case, a small tumor-like structure could be clearly visualized. The contrast-resolution calculated using tumor/lung as foreground/background for kV fluoroscopy, cone beam computed tomography (CBCT), and scattering image were 0.037, 0.70, and 0.54, respectively. In the patient case, tumor motion could be clearly observed in the scatter images. Imaging dose to the voxels directly exposed by the slit beam was ~0.4 times of that under a single CBCT projection. These studies demonstrated the potential feasibility of the proposed imaging scheme to capture the instantaneous anatomy of a patient on a 2D plane with a high image contrast. Clear visualization of the tumor motion may facilitate marker-less tumor tracking.

Cardiorespiratory motion-compensated micro-CT image reconstruction using an artifact model-based motion estimation

PURPOSE

Cardiac in vivo micro-CT imaging of small animals typically requires double gating due to long scan times and high respiratory rates. The simultaneous respiratory and cardiac gating can either be done prospectively or retrospectively. In any case, for true 5D imaging, i.e., three spatial dimensions plus one respiratory-temporal dimension plus one cardiac temporal dimension, the amount of information corresponding to a given respiratory and cardiac phase is orders of magnitude lower than the total amount of information acquired. Achieving similar image quality for 5D than for usual 3D investigations would require increasing the amount of data and thus the applied dose to the animal. Therefore, the goal is phase-correlated imaging with high image quality but without increasing the dose level.

METHODS

To achieve this, the authors propose a new image reconstruction algorithm that makes use of all available projection data, also of that corresponding to other motion windows. In particular, the authors apply a motion-compensated image reconstruction approach that sequentially compensates for respiratory and cardiac motion to decrease the impact of sparsification. In that process, all projection data are used no matter which motion phase they were acquired in. Respiratory and cardiac motion are compensated for by using motion vector fields. These motion vector fields are estimated from initial phase-correlated reconstructions based on a deformable registration approach. To decrease the sensitivity of the registration to sparse-view artifacts, an artifact model-based approach is used including a cyclic consistent nonrigid registration algorithm.

RESULTS

The preliminary results indicate that the authors' approach removes the sparse-view artifacts of conventional phase-correlated reconstructions while maintaining temporal resolution. In addition, it achieves noise levels and spatial resolution comparable to that of nongated reconstructions due to the improved dose usage. By using the proposed motion estimation, no sensitivity to streaking artifacts has been observed.

CONCLUSIONS

Using sequential double gating combined with artifact model-based motion estimation allows to accurately estimate respiratory and cardiac motion from highly undersampled data. No sensitivity to streaking artifacts introduced by sparse angular sampling has been observed for the investigated dose levels. The motion-compensated image reconstruction was able to correct for both, respiratory and cardiac motion, by applying the estimated motion vector fields. The administered dose per animal can thus be reduced for 5D imaging allowing for longitudinal studies at the highest image quality.

Common-mask guided image reconstruction (c-MGIR) for enhanced 4D cone-beam computed tomography

Compared to 3D cone beam computed tomography (3D CBCT), the image quality of commercially available four-dimensional (4D) CBCT is severely impaired due to the insufficient amount of projection data available for each phase. Since the traditional Feldkamp-Davis-Kress (FDK)-based algorithm is infeasible for reconstructing high quality 4D CBCT images with limited projections, investigators had developed several compress-sensing (CS) based algorithms to improve image quality. The aim of this study is to develop a novel algorithm which can provide better image quality than the FDK and other CS based algorithms with limited projections. We named this algorithm 'the common mask guided image reconstruction' (c-MGIR).In c-MGIR, the unknown CBCT volume is mathematically modeled as a combination of phase-specific motion vectors and phase-independent static vectors. The common-mask matrix, which is the key concept behind the c-MGIR algorithm, separates the common static part across all phase images from the possible moving part in each phase image. The moving part and the static part of the volumes were then alternatively updated by solving two sub-minimization problems iteratively. As the novel mathematical transformation allows the static volume and moving volumes to be updated (during each iteration) with global projections and 'well' solved static volume respectively, the algorithm was able to reduce the noise and under-sampling artifact (an issue faced by other algorithms) to the maximum extent. To evaluate the performance of our proposed c-MGIR, we utilized imaging data from both numerical phantoms and a lung cancer patient. The qualities of the images reconstructed with c-MGIR were compared with (1) standard FDK algorithm, (2) conventional total variation (CTV) based algorithm, (3) prior image constrained compressed sensing (PICCS) algorithm, and (4) motion-map constrained image reconstruction (MCIR) algorithm, respectively. To improve the efficiency of the algorithm, the code was implemented with a graphic processing unit for parallel processing purposes.Root mean square error (RMSE) between the ground truth and reconstructed volumes of the numerical phantom were in the descending order of FDK, CTV, PICCS, MCIR, and c-MGIR for all phases. Specifically, the means and the standard deviations of the RMSE of FDK, CTV, PICCS, MCIR and c-MGIR for all phases were 42.64  ±  6.5%, 3.63  ±  0.83%, 1.31%  ±  0.09%, 0.86%  ±  0.11% and 0.52 %  ±  0.02%, respectively. The image quality of the patient case also indicated the superiority of c-MGIR compared to other algorithms.The results indicated that clinically viable 4D CBCT images can be reconstructed while requiring no more projection data than a typical clinical 3D CBCT scan. This makes c-MGIR a potential online reconstruction algorithm for 4D CBCT, which can provide much better image quality than other available algorithms, while requiring less dose and potentially less scanning time.

Impact of scanning parameters and breathing patterns on image quality and accuracy of tumor motion reconstruction in 4D CBCT: a phantom study

Four-dimensional, cone-beam CT (4D CBCT) substantially reduces respiration-induced motion blurring artifacts in three-dimension (3D) CBCT. However, the image quality of 4D CBCT is significantly degraded which may affect its accuracy in localizing a mobile tumor for high-precision, image-guided radiation therapy (IGRT). The purpose of this study was to investigate the impact of scanning parameters hereinafter collectively referred to as scanning sequence) and breathing patterns on the image quality and the accuracy of computed tumor trajectory for a commercial 4D CBCT system, in preparation for its clinical implementation. We simulated a series of periodic and aperiodic sinusoidal breathing patterns with a respiratory motion phantom. The aperiodic pattern was created by varying the period or amplitude of individual sinusoidal breathing cycles. 4D CBCT scans of the phantom were acquired with a manufacturer-supplied scanning sequence (4D-S-slow) and two in-house modified scanning sequences (4D-M-slow and 4D-M-fast). While 4D-S-slow used small field of view (FOV), partial rotation (200°), and no imaging filter, 4D-M-slow and 4D-M-fast used medium FOV, full rotation, and the F1 filter. The scanning speed was doubled in 4D-M-fast (100°/min gantry rotation). The image quality of the 4D CBCT scans was evaluated using contrast-to-noise ratio (CNR), signal-to-noise ratio (SNR), and motion blurring ratio (MBR). The trajectory of the moving target was reconstructed by registering each phase of the 4D CBCT with a reference CT. The root-mean-squared-error (RMSE) analysis was used to quantify its accuracy. Significant decrease in CNR and SNR from 3D CBCT to 4D CBCT was observed. The 4D-S-slow and 4D-M-fast scans had comparable image quality, while the 4D-M-slow scans had better performance due to doubled projections. Both CNR and SNR decreased slightly as the breathing period increased, while no dependence on the amplitude was observed. The difference of both CNR and SNR between periodic and aperiodic breathing patterns was insignificant (p > 0.48). At end-exhale phases, the motion blurring was negligible for both periodic and aperiodic breathing patterns; at mid-inhale phase, the motion blurring increased as the period, the amplitude or the amount of cycle-to-cycle variation on amplitude increased. Overall, the accuracy of localizing the moving target in 4D CBCT was within 2 mm under all studied cases. No difference in the RMSEs was noticed among the three scanning sequences. The 4D-M-fast scans, free of volume truncation artifacts, exhibited comparable image quality and accuracy in tumor motion reconstruction as the 4D-S-slow scans with reduced imaging dose (0.60 cGy vs. 0.99 cGy) due to the use of faster gantry rotation and the F1 filter, suggesting its suitability for clinical use.

Artifacts: The downturn of CBCT image

Cone-beam computed tomography (CBCT) has been accepted as a useful tool for diagnosis and treatment planning in dentistry. Despite a growing trend of CBCT in dentistry, it has some disadvantages like artifacts. Artifacts are discrepancies between the reconstructed visual image and the actual content of the subject which degrade the quality of CBCT images, making them diagnostically unusable. Additionally, structures that do not exist in the subject may appear within images. Such structures can occur because of patient motion, the image capture and reconstruction process. To optimize image quality, it is necessary to understand the types of artifacts. This article aims to throw light on the various types of artifacts associated with CBCT images.

Imaging Findings of the Distal Radio-Ulnar Joint in Trauma

Traumatic lesions of the distal radio-ulnar joint (DRUJ) occur frequently in conjunction with fractures of the distal radius. They are a common cause of pain and limited range of motion after distal radial fractures. Due to the complex anatomy they are however often ignored or underappreciated. Distal radial fractures and luxations of the DRUJ often disturb the normal curvature of the radial notch and cause damage to the cartilage of this joint. The growth of the radius may be disrupted, resulting in a positive ulnar variance, and possibly give rise to complications such as ulnar abutment and motion restriction.

Ulnar styloid fractures – sometimes barely visible on plain film – may give rise to symptomatic bony pseudarthrosis, dislocation and laceration of the tendon of the m. extensor carpi ulnaris and a rare posttraumatic deformity of the ulnar epiphysis. Also the possibility of lesions at the adjacent triangular fibrocartilage complex and the joint capsule should be kept in mind.

This paper presents a pictorial review of the complex functional anatomy and pathologic conditions of this joint and emphasises why the DRUJ should be evaluated independently and thoroughly. The merit of each imaging modality is mentioned.

A correction article relating to Fig. 2 and Fig. 27 can be found here: http://dx.doi.org/10.5334/jbr-btr.966

Direct comparison of conventional radiography and cone-beam CT in small bone and joint trauma

OBJECTIVE

To compare the diagnostic value of cone-beam computed tomography (CBCT) and conventional radiography (CR) after acute small bone or joint trauma.

MATERIALS AND METHODS

Between March 2013 and January 2014, 231 patients with recent small bone or joint trauma underwent CR and subsequent CBCT. CR and CBCT examinations were independently assessed by two readers, blinded to the result of the other modality. The total number of fractures as well as the number of complex fractures were compared, and inter- and intraobserver agreement for CBCT was calculated. In addition, radiation doses and evaluation times for both modalities were noted and statistically compared.

RESULTS

Fracture detection on CBCT increased by 35% and 37% for reader 1 and reader 2, respectively, and identification of complex fractures increased by 236% and 185%. Interobserver agreement for CBCT was almost perfect, as was intraobserver agreement for reader 1. The intraobserver agreement for reader 2 was substantial. Radiation doses and evaluation time were significantly higher for CBCT.

CONCLUSION

CBCT detects significantly more small bone and joint fractures, in particular complex fractures, than CR. In the majority of cases, the clinical implication of the additionally detected fractures is limited, but in some patients (e.g., fracture-dislocations), the management is significantly influenced by these findings. As the radiation dose for CBCT substantially exceeds that of CR, we suggest adhering to CR as the first-line examination after small bone and joint trauma and keeping CBCT for patients with clinical-radiographic discordance or suspected complex fractures in need of further (preoperative) assessment.

A method for volumetric imaging in radiotherapy using single x-ray projection

PURPOSE

It is an intriguing problem to generate an instantaneous volumetric image based on the corresponding x-ray projection. The purpose of this study is to develop a new method to achieve this goal via a sparse learning approach.

METHODS

To extract motion information hidden in projection images, the authors partitioned a projection image into small rectangular patches. The authors utilized a sparse learning method to automatically select patches that have a high correlation with principal component analysis (PCA) coefficients of a lung motion model. A model that maps the patch intensity to the PCA coefficients was built along with the patch selection process. Based on this model, a measured projection can be used to predict the PCA coefficients, which are then further used to generate a motion vector field and hence a volumetric image. The authors have also proposed an intensity baseline correction method based on the partitioned projection, in which the first and the second moments of pixel intensities at a patch in a simulated projection image are matched with those in a measured one via a linear transformation. The proposed method has been validated in both simulated data and real phantom data.

RESULTS

The algorithm is able to identify patches that contain relevant motion information such as the diaphragm region. It is found that an intensity baseline correction step is important to remove the systematic error in the motion prediction. For the simulation case, the sparse learning model reduced the prediction error for the first PCA coefficient to 5%, compared to the 10% error when sparse learning was not used, and the 95th percentile error for the predicted motion vector was reduced from 2.40 to 0.92 mm. In the phantom case with a regular tumor motion, the predicted tumor trajectory was successfully reconstructed with a 0.82 mm error for tumor center localization compared to a 1.66 mm error without using the sparse learning method. When the tumor motion was driven by a real patient breathing signal with irregular periods and amplitudes, the average tumor center error was 0.6 mm. The algorithm robustness with respect to sparsity level, patch size, and presence or absence of diaphragm, as well as computation time, has also been studied.

CONCLUSIONS

The authors have developed a new method that automatically identifies motion information from an x-ray projection, based on which a volumetric image is generated.

Investigation of different sparsity transforms for the PICCS algorithm in small-animal respiratory gated CT

Respiratory gating helps to overcome the problem of breathing motion in cardiothoracic small-animal imaging by acquiring multiple images for each projection angle and then assigning projections to different phases. When this approach is used with a dose similar to that of a static acquisition, a low number of noisy projections are available for the reconstruction of each respiratory phase, thus leading to streak artifacts in the reconstructed images. This problem can be alleviated using a prior image constrained compressed sensing (PICCS) algorithm, which enables accurate reconstruction of highly undersampled data when a prior image is available. We compared variants of the PICCS algorithm with different transforms in the prior penalty function: gradient, unitary, and wavelet transform. In all cases the problem was solved using the Split Bregman approach, which is efficient for convex constrained optimization. The algorithms were evaluated using simulations generated from data previously acquired on a micro-CT scanner following a high-dose protocol (four times the dose of a standard static protocol). The resulting data were used to simulate scenarios with different dose levels and numbers of projections. All compressed sensing methods performed very similarly in terms of noise, spatiotemporal resolution, and streak reduction, and filtered back-projection was greatly improved. Nevertheless, the wavelet domain was found to be less prone to patchy cartoon-like artifacts than the commonly used gradient domain.

A hybrid reconstruction algorithm for fast and accurate 4D cone-beam CT imaging

PURPOSE

4D cone beam CT (4D-CBCT) has been utilized in radiation therapy to provide 4D image guidance in lung and upper abdomen area. However, clinical application of 4D-CBCT is currently limited due to the long scan time and low image quality. The purpose of this paper is to develop a new 4D-CBCT reconstruction method that restores volumetric images based on the 1-min scan data acquired with a standard 3D-CBCT protocol.

METHODS

The model optimizes a deformation vector field that deforms a patient-specific planning CT (p-CT), so that the calculated 4D-CBCT projections match measurements. A forward-backward splitting (FBS) method is invented to solve the optimization problem. It splits the original problem into two well-studied subproblems, i.e., image reconstruction and deformable image registration. By iteratively solving the two subproblems, FBS gradually yields correct deformation information, while maintaining high image quality. The whole workflow is implemented on a graphic-processing-unit to improve efficiency. Comprehensive evaluations have been conducted on a moving phantom and three real patient cases regarding the accuracy and quality of the reconstructed images, as well as the algorithm robustness and efficiency.

RESULTS

The proposed algorithm reconstructs 4D-CBCT images from highly under-sampled projection data acquired with 1-min scans. Regarding the anatomical structure location accuracy, 0.204 mm average differences and 0.484 mm maximum difference are found for the phantom case, and the maximum differences of 0.3-0.5 mm for patients 1-3 are observed. As for the image quality, intensity errors below 5 and 20 HU compared to the planning CT are achieved for the phantom and the patient cases, respectively. Signal-noise-ratio values are improved by 12.74 and 5.12 times compared to results from FDK algorithm using the 1-min data and 4-min data, respectively. The computation time of the algorithm on a NVIDIA GTX590 card is 1-1.5 min per phase.

CONCLUSIONS

High-quality 4D-CBCT imaging based on the clinically standard 1-min 3D CBCT scanning protocol is feasible via the proposed hybrid reconstruction algorithm.

Improving thoracic four-dimensional cone-beam CT reconstruction with anatomical-adaptive image regularization (AAIR)

Total-variation (TV) minimization reconstructions can significantly reduce noise and streaks in thoracic four-dimensional cone-beam computed tomography (4D CBCT) images compared to the Feldkamp-Davis-Kress (FDK) algorithm currently used in practice. TV minimization reconstructions are, however, prone to over-smoothing anatomical details and are also computationally inefficient. The aim of this study is to demonstrate a proof of concept that these disadvantages can be overcome by incorporating the general knowledge of the thoracic anatomy via anatomy segmentation into the reconstruction. The proposed method, referred as the anatomical-adaptive image regularization (AAIR) method, utilizes the adaptive-steepest-descent projection-onto-convex-sets (ASD-POCS) framework, but introduces an additional anatomy segmentation step in every iteration. The anatomy segmentation information is implemented in the reconstruction using a heuristic approach to adaptively suppress over-smoothing at anatomical structures of interest. The performance of AAIR depends on parameters describing the weighting of the anatomy segmentation prior and segmentation threshold values. A sensitivity study revealed that the reconstruction outcome is not sensitive to these parameters as long as they are chosen within a suitable range. AAIR was validated using a digital phantom and a patient scan and was compared to FDK, ASD-POCS and the prior image constrained compressed sensing (PICCS) method. For the phantom case, AAIR reconstruction was quantitatively shown to be the most accurate as indicated by the mean absolute difference and the structural similarity index. For the patient case, AAIR resulted in the highest signal-to-noise ratio (i.e. the lowest level of noise and streaking) and the highest contrast-to-noise ratios for the tumor and the bony anatomy (i.e. the best visibility of anatomical details). Overall, AAIR was much less prone to over-smoothing anatomical details compared to ASD-POCS and did not suffer from residual noise/streaking and motion blur migrated from the prior image as in PICCS. AAIR was also found to be more computationally efficient than both ASD-POCS and PICCS, with a reduction in computation time of over 50% compared to ASD-POCS. The use of anatomy segmentation was, for the first time, demonstrated to significantly improve image quality and computational efficiency for thoracic 4D CBCT reconstruction. Further developments are required to facilitate AAIR for practical use.

Metal artifact correction for x-ray computed tomography using kV and selective MV imaging

PURPOSE

The overall goal of this work is to improve the computed tomography (CT) image quality for patients with metal implants or fillings by completing the missing kilovoltage (kV) projection data with selectively acquired megavoltage (MV) data that do not suffer from photon starvation. When both of these imaging systems, which are available on current radiotherapy devices, are used, metal streak artifacts are avoided, and the soft-tissue contrast is restored, even for regions in which the kV data cannot contribute any information.

METHODS

Three image-reconstruction methods, including two filtered back-projection (FBP)-based analytic methods and one iterative method, for combining kV and MV projection data from the two on-board imaging systems of a radiotherapy device are presented in this work. The analytic reconstruction methods modify the MV data based on the information in the projection or image domains and then patch the data onto the kV projections for a FBP reconstruction. In the iterative reconstruction, the authors used dual-energy (DE) penalized weighted least-squares (PWLS) methods to simultaneously combine the kV/MV data and perform the reconstruction.

RESULTS

The authors compared kV/MV reconstructions to kV-only reconstructions using a dental phantom with fillings and a hip-implant numerical phantom. Simulation results indicated that dual-energy sinogram patch FBP and the modified dual-energy PWLS method can successfully suppress metal streak artifacts and restore information lost due to photon starvation in the kV projections. The root-mean-square errors of soft-tissue patterns obtained using combined kV/MV data are 10-15 Hounsfield units smaller than those of the kV-only images, and the structural similarity index measure also indicates a 5%-10% improvement in the image quality. The added dose from the MV scan is much less than the dose from the kV scan if a high efficiency MV detector is assumed.

CONCLUSIONS

The authors have shown that it is possible to improve the image quality of kV CTs for patients with metal implants or fillings by completing the missing kV projection data with selectively acquired MV data that do not suffer from photon starvation. Numerical simulations demonstrated that dual-energy sinogram patch FBP and a modified kV/MV PWLS method can successfully suppress metal streak artifacts and restore information lost due to photon starvation in kV projections. Combined kV/MV images may permit the improved delineation of structures of interest in CT images for patients with metal implants or fillings.

Image quality in thoracic 4D cone-beam CT: a sensitivity analysis of respiratory signal, binning method, reconstruction algorithm, and projection angular spacing

PURPOSE

Respiratory signal, binning method, and reconstruction algorithm are three major controllable factors affecting image quality in thoracic 4D cone-beam CT (4D-CBCT), which is widely used in image guided radiotherapy (IGRT). Previous studies have investigated each of these factors individually, but no integrated sensitivity analysis has been performed. In addition, projection angular spacing is also a key factor in reconstruction, but how it affects image quality is not obvious. An investigation of the impacts of these four factors on image quality can help determine the most effective strategy in improving 4D-CBCT for IGRT.

METHODS

Fourteen 4D-CBCT patient projection datasets with various respiratory motion features were reconstructed with the following controllable factors: (i) respiratory signal (real-time position management, projection image intensity analysis, or fiducial marker tracking), (ii) binning method (phase, displacement, or equal-projection-density displacement binning), and (iii) reconstruction algorithm [Feldkamp-Davis-Kress (FDK), McKinnon-Bates (MKB), or adaptive-steepest-descent projection-onto-convex-sets (ASD-POCS)]. The image quality was quantified using signal-to-noise ratio (SNR), contrast-to-noise ratio, and edge-response width in order to assess noise/streaking and blur. The SNR values were also analyzed with respect to the maximum, mean, and root-mean-squared-error (RMSE) projection angular spacing to investigate how projection angular spacing affects image quality.

RESULTS

The choice of respiratory signals was found to have no significant impact on image quality. Displacement-based binning was found to be less prone to motion artifacts compared to phase binning in more than half of the cases, but was shown to suffer from large interbin image quality variation and large projection angular gaps. Both MKB and ASD-POCS resulted in noticeably improved image quality almost 100% of the time relative to FDK. In addition, SNR values were found to increase with decreasing RMSE values of projection angular gaps with strong correlations (r ≈ -0.7) regardless of the reconstruction algorithm used.

CONCLUSIONS

Based on the authors' results, displacement-based binning methods, better reconstruction algorithms, and the acquisition of even projection angular views are the most important factors to consider for improving thoracic 4D-CBCT image quality. In view of the practical issues with displacement-based binning and the fact that projection angular spacing is not currently directly controllable, development of better reconstruction algorithms represents the most effective strategy for improving image quality in thoracic 4D-CBCT for IGRT applications at the current stage.

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.

Actively triggered 4d cone-beam CT acquisition

PURPOSE

4d cone-beam computed tomography (CBCT) scans are usually reconstructed by extracting the motion information from the 2d projections or an external surrogate signal, and binning the individual projections into multiple respiratory phases. In this "after-the-fact" binning approach, however, projections are unevenly distributed over respiratory phases resulting in inefficient utilization of imaging dose. To avoid excess dose in certain respiratory phases, and poor image quality due to a lack of projections in others, the authors have developed a novel 4d CBCT acquisition framework which actively triggers 2d projections based on the forward-predicted position of the tumor.

METHODS

The forward-prediction of the tumor position was independently established using either (i) an electromagnetic (EM) tracking system based on implanted EM-transponders which act as a surrogate for the tumor position, or (ii) an external motion sensor measuring the chest-wall displacement and correlating this external motion to the phase-shifted diaphragm motion derived from the acquired images. In order to avoid EM-induced artifacts in the imaging detector, the authors devised a simple but effective "Faraday" shielding cage. The authors demonstrated the feasibility of their acquisition strategy by scanning an anthropomorphic lung phantom moving on 1d or 2d sinusoidal trajectories.

RESULTS

With both tumor position devices, the authors were able to acquire 4d CBCTs free of motion blurring. For scans based on the EM tracking system, reconstruction artifacts stemming from the presence of the EM-array and the EM-transponders were greatly reduced using newly developed correction algorithms. By tuning the imaging frequency independently for each respiratory phase prior to acquisition, it was possible to harmonize the number of projections over respiratory phases. Depending on the breathing period (3.5 or 5 s) and the gantry rotation time (4 or 5 min), between ∼90 and 145 projections were acquired per respiratory phase resulting in a dose of ∼1.7-2.6 mGy per respiratory phase. Further dose savings and decreases in the scanning time are possible by acquiring only a subset of all respiratory phases, for example, peak-exhale and peak-inhale only scans.

CONCLUSIONS

This study is the first experimental demonstration of a new 4d CBCT acquisition paradigm in which imaging dose is efficiently utilized by actively triggering only those projections that are desired for the reconstruction process.

Evaluation of the effect of respiratory and anatomical variables on a Fourier technique for markerless, self-sorted 4D-CBCT

A novel technique based on Fourier transform theory has been developed that directly extracts respiratory information from projections without the use of external surrogates. While the feasibility has been demonstrated with three patients, a more extensive validation is necessary. Therefore, the purpose of this work is to investigate the effects of a variety of respiratory and anatomical scenarios on the performance of the technique with the 4D digital extended cardiac torso phantom. FT-phase and FT-magnitude methods were each applied to identify peak-inspiration projections and quantitatively compared to the gold standard of visual identification. Both methods proved to be robust across the studied scenarios with average differences in respiratory phase <10% and percentage of projections assigned within 10% of the gold standard >90%, when incorporating minor modifications to region-of-interest (ROI) selection and/or low-frequency location for select cases of DA and lung percentage in the field of view of the projection. Nevertheless, in the instance where one method initially faltered, the other method prevailed and successfully identified peak-inspiration projections. This is promising because it suggests that the two methods provide complementary information to each other. To ensure appropriate clinical adaptation of markerless, self-sorted four-dimensional cone-beam CT (4D-CBCT), perhaps an optimal integration of the two methods can be developed.

Respiratory triggered 4D cone-beam computed tomography: a novel method to reduce imaging dose

PURPOSE

A novel method called respiratory triggered 4D cone-beam computed tomography (RT 4D CBCT) is described whereby imaging dose can be reduced without degrading image quality. RT 4D CBCT utilizes a respiratory signal to trigger projections such that only a single projection is assigned to a given respiratory bin for each breathing cycle. In contrast, commercial 4D CBCT does not actively use the respiratory signal to minimize image dose.

METHODS

To compare RT 4D CBCT with conventional 4D CBCT, 3600 CBCT projections of a thorax phantom were gathered and reconstructed to generate a ground truth CBCT dataset. Simulation pairs of conventional 4D CBCT acquisitions and RT 4D CBCT acquisitions were developed assuming a sinusoidal respiratory signal which governs the selection of projections from the pool of 3600 original projections. The RT 4D CBCT acquisition triggers a single projection when the respiratory signal enters a desired acquisition bin; the conventional acquisition does not use a respiratory trigger and projections are acquired at a constant frequency. Acquisition parameters studied were breathing period, acquisition time, and imager frequency. The performance of RT 4D CBCT using phase based and displacement based sorting was also studied. Image quality was quantified by calculating difference images of the test dataset from the ground truth dataset. Imaging dose was calculated by counting projections.

RESULTS

Using phase based sorting RT 4D CBCT results in 47% less imaging dose on average compared to conventional 4D CBCT. Image quality differences were less than 4% at worst. Using displacement based sorting RT 4D CBCT results in 57% less imaging dose on average, than conventional 4D CBCT methods; however, image quality was 26% worse with RT 4D CBCT.

CONCLUSIONS

Simulation studies have shown that RT 4D CBCT reduces imaging dose while maintaining comparable image quality for phase based 4D CBCT; image quality is degraded for displacement based RT 4D CBCT in its current implementation.

Extracting respiratory signals from thoracic cone beam CT projections

The patient respiratory signal associated with the cone beam CT (CBCT) projections is important for lung cancer radiotherapy. In contrast to monitoring an external surrogate of respiration, such a signal can be extracted directly from the CBCT projections. In this paper, we propose a novel local principal component analysis (LPCA) method to extract the respiratory signal by distinguishing the respiration motion-induced content change from the gantry rotation-induced content change in the CBCT projections. The LPCA method is evaluated by comparing with three state-of-the-art projection-based methods, namely the Amsterdam Shroud method, the intensity analysis method and the Fourier-transform-based phase analysis method. The clinical CBCT projection data of eight patients, acquired under various clinical scenarios, were used to investigate the performance of each method. We found that the proposed LPCA method has demonstrated the best overall performance for cases tested and thus is a promising technique for extracting a respiratory signal. We also identified the applicability of each existing method.

Target-specific optimization of four-dimensional cone beam computed tomography

PURPOSE

Under-sampling artifacts are a major problem in four-dimensional cone-beam computed tomography (4D-CBCT) and may compromise evaluation of target motion. Several studies have addressed scan parameter selection for minimizing under-sampling artifacts; however, the role of the target characteristics in scan parameter selection has not been investigated. In this work, the authors evaluated 4D-CBCT performance by assessing the accuracy of target motion measurements for various target sizes and motions. The results may serve as patient-specific guidelines for scan parameters selection in 4D-CBCT.

METHODS

The authors acquired 4D-CBCT scans of a moving phantom consisting of six water-filled sphere targets of sizes 10-37 mm, using various scan times ranging from 30 s to 3 min., setting the motion to 3-s and 6-s periods. The authors used automatic image registration to extract the target motion trajectories and evaluated these measurements for various target sizes and motions over various combinations of scan parameters including scan time, detector configuration, number of respiration phases, and reconstruction filters.

RESULTS

The most important object parameter to 4D-CBCT performance was the period of motion. Measurements for the 6-s motion were always systematically less accurate than measurements for the 3-s motion for 34 of 36 objects of various sizes and periods of motion. The 6-s motion required a greater scan time than the 3-s motion for equivalent measurement accuracy. The second most influential parameter was the target size. For the 3-s period of motion, objects larger than 13 mm were tracked with sub-millimeter accuracy with a 1-min scan time. For the 6-s period of motion, objects larger than 22-mm were tracked with sub-millimeter accuracy with a 1.5-min scan time. For all sizes and motions, temporal blurring was observed when the number of phases was fewer than 8. Offset detector configuration provided the same performance as centered detector except for small targets (20 mm.) at short scan times (≤1 min). Finally, reconstruction filtering and number of respiratory phases did not affect performance.

CONCLUSIONS

Scan time should be set according to target size and motion. The authors have provided figures that provide the minimum scan time needed to achieve the particular motion measurement accuracy for a particular size and motion period.