High-fidelity artifact correction for cone-beam CT imaging of the brain

Feasibility of bone marrow edema detection using dual-energy cone-beam computed tomography

BACKGROUND

Dual-energy (DE) detection of bone marrow edema (BME) would be a valuable new diagnostic capability for the emerging orthopedic cone-beam computed tomography (CBCT) systems. However, this imaging task is inherently challenging because of the narrow energy separation between water (edematous fluid) and fat (health yellow marrow), requiring precise artifact correction and dedicated material decomposition approaches.

PURPOSE

We investigate the feasibility of BME assessment using kV-switching DE CBCT with a comprehensive CBCT artifact correction framework and a two-stage projection- and image-domain three-material decomposition algorithm.

METHODS

DE CBCT projections of quantitative BME phantoms (water containers 100-165 mm in size with inserts presenting various degrees of edema) and an animal cadaver model of BME were acquired on a CBCT test bench emulating the standard wrist imaging configuration of a Multitom Rax twin robotic x-ray system. The slow kV-switching scan protocol involved a 60 kV low energy (LE) beam and a 120 kV high energy (HE) beam switched every 0.5° over a 200° angular span. The DE CBCT data preprocessing and artifact correction framework consisted of (i) projection interpolation onto matched LE and HE projections views, (ii) lag and glare deconvolutions, and (iii) efficient Monte Carlo (MC)-based scatter correction. Virtual non-calcium (VNCa) images for BME detection were then generated by projection-domain decomposition into an Aluminium (Al) and polyethylene basis set (to remove beam hardening) followed by three-material image-domain decomposition into water, Ca, and fat. Feasibility of BME detection was quantified in terms of VNCa image contrast and receiver operating characteristic (ROC) curves. Robustness to object size, position in the field of view (FOV) and beam collimation (varied 20-160 mm) was investigated.

RESULTS

The MC-based scatter correction delivered > 69% reduction of cupping artifacts for moderate to wide collimations (> 80 mm beam width), which was essential to achieve accurate DE material decomposition. In a forearm-sized object, a 20% increase in water concentration (edema) of a trabecular bone-mimicking mixture presented as ∼15 HU VNCa contrast using 80-160 mm beam collimations. The variability with respect to object position in the FOV was modest (< 15% coefficient of variation). The areas under the ROC curve were > 0.9. A femur-sized object presented a somewhat more challenging task, resulting in increased sensitivity to object positioning at 160 mm collimation. In animal cadaver specimens, areas of VNCa enhancement consistent with BME were observed in DE CBCT images in regions of MRI-confirmed edema.

CONCLUSION

Our results indicate that the proposed artifact correction and material decomposition pipeline can overcome the challenges of scatter and limited spectral separation to achieve relatively accurate and sensitive BME detection in DE CBCT. This study provides an important baseline for clinical translation of musculoskeletal DE CBCT to quantitative, point-of-care bone health assessment.

Super-resolution dual-layer CBCT imaging with model-guided deep learning

Objective.This study aims at investigating a novel super resolution CBCT imaging approach with a dual-layer flat panel detector (DL-FPD).Approach.With DL-FPD, the low-energy and high-energy projections acquired from the top and bottom detector layers contain over-sampled spatial information, from which super-resolution CT images can be reconstructed. A simple mathematical model is proposed to explain the signal formation procedure in DL-FPD, and a dedicated recurrent neural network, named suRi-Net, is developed based upon the above imaging model to nonlinearly retrieve the high-resolution dual-energy information. Physical benchtop experiments are conducted to validate the performance of this newly developed super-resolution CBCT imaging method.Main Results.The results demonstrate that the proposed suRi-Net can accurately retrieve high spatial resolution information from the low-energy and high-energy projections of low spatial resolution. Quantitatively, the spatial resolution of the reconstructed CBCT images from the top and bottom detector layers is increased by about 45% and 54%, respectively.Significance.In the future, suRi-Net will provide a new approach to perform high spatial resolution dual-energy imaging in DL-FPD-based CBCT systems.

A quantitative CBCT pipeline based on 2D antiscatter grid and grid-based scatter sampling for image-guided radiation therapy

BACKGROUND

Quantitative accuracy is critical for expanding the role of cone beam CT (CBCT) imaging from target localization to quantitative treatment monitoring and plan adaptations in radiation therapy. Despite advances in CBCT image quality improvement methods, quantitative accuracy gap between CBCT and multi-detector CT (MDCT) remains.

PURPOSE

In this work, a physics-driven approach was investigated that combined robust scatter rejection, raw data correction and iterative image reconstruction to further improve CBCT image quality and quantitative accuracy, referred to as quantitative CBCT (qCBCT).

METHODS

QCBCT approach includes tungsten 2D antiscatter grid hardware, residual scatter correction with grid-based scatter sampling, image lag, and beam hardening correction for offset detector geometry linac-mounted CBCT. Images were reconstructed with iterative image reconstruction to reduce image noise. qCBCT was evaluated using a variety of phantoms to investigate the effect of object size and its composition on image quality, and image quality was benchmarked against clinical CBCT and gold standard MDCT images used for treatment planning.

RESULTS

QCBCT provided statistically significant improvement in CT number accuracy and reduced image artifacts when compared to clinical CBCT images. When compared to gold standard MDCT, mean HU errors in qCBCT and clinical CBCT were 17 ± 9 and 38 ± 29 HU, respectively. Magnitude of phantom size dependent HU variations were comparable between MDCT and qCBCT images. With iterative reconstruction, contrast-to-noise ratio improved by 25% when compared to clinical CBCT protocols.

CONCLUSIONS

Combination of novel scatter suppression techniques and other data correction methods in qCBCT provided CT number accuracy comparable to gold standard MDCT used for treatment planning. This approach may potentially improve CBCT's promise in fulfilling the tasks that demand high quantitative accuracy, such as online dose calculations and treatment response assessment, in image guided radiation therapy.

Dual-energy head cone-beam CT using a dual-layer flat-panel detector: Hybrid material decomposition and a feasibility study

BACKGROUND

Flat panel detector (FPD) based cone-beam computed tomography (CT) has made tremendous progress in the last two decades, with many new and advanced medical and industrial applications keeping emerging from diagnostic imaging and image guidance for radiotherapy and interventional surgery. The current cone-beam CT (CBCT), however, is still suboptimal for head CT scan which requires a high standard of image quality. While the dual-layer FPD technology is under extensive development and is promising to further advance CBCT from qualitative anatomic imaging to quantitative dual-energy CT, its potential of enabling head CBCT applications has not yet been fully investigated.

PURPOSE

The relatively moderate energy separation from the dual-layer FPD and the overall low signal level especially at the bottom-layer detector, could raise significant challenges in performing high-quality dual-energy material decomposition (MD). In this work, we propose a hybrid, physics and model guided, MD algorithm that attempts to fully use the detected x-ray signals and prior-knowledge behind head CBCT using dual-layer FPD.

METHODS

Firstly, a regular projection-domain MD is performed as initial results of our approach and for comparison as conventional method. Secondly, based on the combined projection, a dual-layer multi-material spectral correction (dMMSC) is applied to generate beam hardening free images. Thirdly, the dMMSC corrected projections are adopted as a physics-model based guidance to generate a hybrid MD. A set of physics experiments including fan-beam scan and cone-beam scan using a head phantom and a Gammex Multi-Energy CT phantom are conducted to validate our proposed approach.

RESULTS

The combined reconstruction could reduce noise by about 10% with no visible resolution degradation. The fan-beam studies on the Gammex phantom demonstrated an improved MD performance, with the averaged iodine quantification error for the 5-15 mg/ml iodine inserts reduced from about 5.6% to 3.0% by the hybrid method. On fan-beam scan of the head phantom, our proposed hybrid MD could significantly reduce the streak artifacts, with CT number nonuniformity (NU) in the selected regions of interest (ROIs) reduced from 23 Hounsfield Units (HU) to 4.2 HU, and the corresponding noise suppressed from 31 to 6.5 HU. For cone-beam scan, after scatter correction (SC) and cone-beam artifact reduction (CBAR), our approach can also significantly improve image quality, with CT number NU in the selected ROI reduced from 24.2 to 6.6 HU and the noise level suppressed from 22.1 to 8.2 HU.

CONCLUSIONS

Our proposed physics and model guided hybrid MD for dual-layer FPD based head CBCT can significantly improve the robustness of MD and suppress the low-signal artifact. This preliminary feasibility study also demonstrated that the dual-layer FPD is promising to enable head CBCT spectral imaging.

An unsupervised dual contrastive learning framework for scatter correction in cone-beam CT image

PURPOSE

Cone-beam computed tomography (CBCT) is widely utilized in modern radiotherapy; however, CBCT images exhibit increased scatter artifacts compared to planning CT (pCT), compromising image quality and limiting further applications. Scatter correction is thus crucial for improving CBCT image quality.

METHODS

In this study, we proposed an unsupervised contrastive learning method for CBCT scatter correction. Initially, we transformed low-quality CBCT into high-quality synthetic pCT (spCT) and generated forward projections of CBCT and spCT. By computing the difference between these projections, we obtained a residual image containing image details and scatter artifacts. Image details primarily comprise high-frequency signals, while scatter artifacts consist mainly of low-frequency signals. We extracted the scatter projection signal by applying a low-pass filter to remove image details. The corrected CBCT (cCBCT) projection signal was obtained by subtracting the scatter artifacts projection signal from the original CBCT projection. Finally, we employed the FDK reconstruction algorithm to generate the cCBCT image.

RESULTS

To evaluate cCBCT image quality, we aligned the CBCT and pCT of six patients. In comparison to CBCT, cCBCT maintains anatomical consistency and significantly enhances CT number, spatial homogeneity, and artifact suppression. The mean absolute error (MAE) of the test data decreased from 88.0623 ± 26.6700 HU to 17.5086 ± 3.1785 HU. The MAE of fat regions of interest (ROIs) declined from 370.2980 ± 64.9730 HU to 8.5149 ± 1.8265 HU, and the error between their maximum and minimum CT numbers decreased from 572.7528 HU to 132.4648 HU. The MAE of muscle ROIs reduced from 354.7689 ± 25.0139 HU to 16.4475 ± 3.6812 HU. We also compared our proposed method with several conventional unsupervised synthetic image generation techniques, demonstrating superior performance.

CONCLUSIONS

Our approach effectively enhances CBCT image quality and shows promising potential for future clinical adoption.

An Unsupervised Learning-Based Regional Deformable Model for Automated Multi-Organ Contour Propagation

The aim of this study is to evaluate a regional deformable model based on a deep unsupervised learning model for automatic contour propagation in breast cone-beam computed tomography-guided adaptive radiation therapy. A deep unsupervised learning model was introduced to map breast's tumor bed, clinical target volume, heart, left lung, right lung, and spinal cord from planning computed tomography to cone-beam CT. To improve the traditional image registration method's performance, we used a regional deformable framework based on the narrow-band mapping, which can mitigate the effect of the image artifacts on the cone-beam CT. We retrospectively selected 373 anonymized cone-beam CT volumes from 111 patients with breast cancer. The cone-beam CTs are divided into three sets. 311 / 20 / 42 cone-beam CT images were used for training, validating, and testing. The manual contour was used as reference for the testing set. We compared the results between the reference and the model prediction for evaluating the performance. The mean Dice between manual reference segmentations and the model predicted segmentations for breast tumor bed, clinical target volume, heart, left lung, right lung, and spinal cord were 0.78 ± 0.09, 0.90 ± 0.03, 0.88 ± 0.04, 0.94 ± 0.03, 0.95 ± 0.02, and 0.77 ± 0.07, respectively. The results demonstrated a good agreement between the reference and the proposed contours. The proposed deep learning-based regional deformable model technique can automatically propagate contours for breast cancer adaptive radiotherapy. Deep learning in contour propagation was promising, but further investigation was warranted.

A novel beam stopper-based approach for scatter correction in digital planar radiography

X-ray scatter in planar radiography degrades the contrast resolution of the image, thus reducing its diagnostic utility. Antiscatter grids partially block scattered photons at the cost of increasing the dose delivered by two- to four-fold and posing geometrical restrictions that hinder their use for other acquisition settings, such as portable radiography. The few software-based approaches investigated for planar radiography mainly estimate the scatter map from a low-frequency version of the image. We present a novel method for scatter correction in planar imaging based on direct patient measurements. Samples from the shadowed regions of an additional partially obstructed projection acquired with a beam stopper placed between the X-ray source and the patient are used to estimate the scatter map. Evaluation with simulated and real data showed an increase in contrast resolution for both lung and spine and recovery of ground truth values superior to those of three recently proposed methods. Our method avoids the biases of post-processing methods and yields results similar to those for an antiscatter grid while removing geometrical restrictions at around half the radiation dose. It can be used in unconventional imaging techniques, such as portable radiography, where training datasets needed for deep-learning approaches would be very difficult to obtain.

The butterfly effect: improving brain cone-beam CT image artifacts for stroke assessment using a novel dual-axis trajectory

BACKGROUND

Cone-beam computed tomography (CBCT) imaging of the brain can be performed in the angiography suite to support various neurovascular procedures. Relying on CBCT brain imaging solely, however, still lacks full diagnostic confidence due to the inferior image quality compared with CT and various imaging artifacts that persist even with modern CBCT.

OBJECTIVE

To perform a detailed evaluation of image artifact improvement using a new CBCT protocol which implements a novel dual-axis 'butterfly' trajectory.

METHODS

Our study included 94 scans from 47 patients who received CBCT imaging for assessment of either ischemia or hemorrhage during a neurovascular procedure. Both a traditional uni-axis 'circular' and novel dual-axis 'butterfly' protocol were performed on each patient (same-patient control). Each brain scan was divided into six regions and scored out of 3 based on six artifacts originating from various physics-based and patient-based sources.

RESULTS

The dual-axis trajectory produces CBCT images with significantly fewer image artifacts than the traditional circular scan (whole brain average artifact score, AS: 0.20 vs 0.33), with the greatest improvement in bone beam hardening (AS: 0.13 vs 0.78) and cone-beam artifacts (AS: 0.04 vs 0.55).

CONCLUSIONS

Recent developments in CBCT imaging protocols have significantly improved image artifacts, which has improved diagnostic confidence for stroke and supports a direct-to-angiography suite transfer approach for patients with acute ischemic stroke.

A correlated sampling-based Monte Carlo simulation for fast CBCT iterative scatter correction

BACKGROUND

In recent years, cone-beam computed tomography (CBCT) has played an important role in medical imaging. However, the applications of CBCT are limited due to the severe scatter contamination. Conventional Monte Carlo (MC) simulation can provide accurate scatter estimation for scatter correction, but the expensive computational cost has always been the bottleneck of MC method in clinical application.

PURPOSE

In this work, an MC simulation method combined with a variance reduction technique called correlated sampling is proposed for fast iterative scatter correction.

METHODS

Correlated sampling exploits correlation between similar simulation systems to reduce the variance of interest quantities. Specifically, conventional MC simulation is first performed on the scatter-contaminated CBCT to generate the initial scatter signal. In the subsequent correction iterations, scatter estimation is then updated by applying correlated MC sampling to the latest corrected CBCT images by reusing the random number sequences of the task-related photons in conventional MC. Afterward, the corrected projections obtained by subtracting the scatter estimation from raw projections are utilized for FDK reconstruction. These steps are repeated until an adequate scatter correction is obtained. The performance of the proposed framework is evaluated by the accuracy of the scatter estimation, the quality of corrected CBCT images and efficiency.

RESULTS

Overall, the difference in mean absolute percentage error between scatter estimation with and without correlated sampling is 0.25% for full-fan case and 0.34% for half-fan case, respectively. In simulation studies, scatter artifacts are substantially eliminated, where the mean absolute error value is reduced from 15 to 2 HU in full-fan case and from 53 to 13 HU in half-fan case. Scatter-to-primary ratio is reduced to 0.02 for full-fan and 0.04 for half-fan, respectively. In phantom study, the contrast-to-noise ratio (CNR) is increased by a factor of 1.63, and the contrast is increased by a factor of 1.77. As for clinical studies, the CNR is improved by 11% and 14% for half-fan and full-fan, respectively. The contrast after correction is increased by 19% for half-fan and 44% for full-fan. Furthermore, root mean square error is also effectively reduced, especially from 78 to 4 HU for full-fan. Experimental results demonstrate that the figure of merit is improved between 23 and 43 folds when using correlated sampling. The proposed method takes less than 25 s for the whole iterative scatter correction process.

CONCLUSIONS

The proposed correlated sampling-based MC simulation method can achieve fast and accurate scatter correction for CBCT, making it suitable for real-time clinical use.

Dual-Energy Cone-Beam CT with Three-Material Decomposition for Bone Marrow Edema Imaging

We investigate the feasibility of bone marrow edema (BME) detection using a kV-switching Dual-Energy (DE) Cone-Beam CT (CBCT) protocol. This task is challenging due to unmatched x-ray paths in the low-energy (LE) and high-energy (HE) spectral channels, CBCT non-idealities such as x-ray scatter, and narrow spectral separation between fat (bone marrow) and water (BME). We propose a comprehensive DE decomposition framework consisting of projection interpolation onto matching LE and HE view angles, fast Monte Carlo scatter correction with low number of tracked photons and Gaussian denoising, and two-stage three-material decompositions involving two-material (fat-Aluminium) Projection-Domain Decomposition (PDD) followed by image-domain three-material (fat-water-bone) base-change. Performance in BME detection was evaluated in simulations and experiments emulating a kV-switching CBCT wrist imaging protocol on a robotic x-ray system with 60 kV LE beam, 120 kV HE beam, and 0.5° angular shift between the LE and HE views. Cubic B-spline interpolation was found to be adequate to resample HE and LE projections of a wrist onto common view angles required by PDD. The DE decomposition maintained acceptable BME detection specificity (<0.2 mL erroneously detected BME volume compared to 0.85 mL true BME volume) over +/-10% range of scatter magnitude errors, as long as the scatter shape was estimated without major distortions. Physical test bench experiments demonstrated successful discrimination of ~20% change in fat concentrations in trabecular bone-mimicking solutions of varying water and fat content.

Feasibility of Dual-Energy Cone-Beam CT of Bone Marrow Edema Using Dual-Layer Flat Panel Detectors

Purpose

We investigated the feasibility of detection and quantification of bone marrow edema (BME) using dual-energy (DE) Cone-Beam CT (CBCT) with a dual-layer flat panel detector (FPD) and three-material decomposition.

Methods

A realistic CBCT system simulator was applied to study the impact of detector quantization, scatter, and spectral calibration errors on the accuracy of fat-water-bone decompositions of dual-layer projections. The CBCT system featured 975 mm source-axis distance, 1,362 mm source-detector distance and a 430 × 430 mm2 dual-layer FPD (top layer: 0.20 mm CsI:Tl, bottom layer: 0.55 mm CsI:Tl; a 1 mm Cu filter between the layers to improve spectral separation). Tube settings were 120 kV (+2 mm Al, +0.2 mm Cu) and 10 mAs per exposure. The digital phantom consisted of a 160 mm water cylinder with inserts containing mixtures of water (volume fraction ranging 0.18 to 0.46) - fat (0.5 to 0.7) - Ca (0.04 to 0.12); decreasing fractions of fat indicated increasing degrees of BME. A two-stage three-material DE decomposition was applied to DE CBCT projections: first, projection-domain decomposition (PDD) into fat-aluminum basis, followed by CBCT reconstruction of intermediate base images, followed by image-domain change of basis into fat, water and bone. Sensitivity to scatter was evaluated by i) adjusting source collimation (12 to 400 mm width) and ii) subtracting various fractions of the true scatter from the projections at 400 mm collimation. The impact of spectral calibration was studied by shifting the effective beam energy (± 2 keV) when creating the PDD lookup table. We further simulated a realistic BME imaging framework, where the scatter was estimated using a fast Monte Carlo (MC) simulation from a preliminary decomposition of the object; the object was a realistic wrist phantom with an 0.85 mL BME stimulus in the radius.

Results

The decomposition is sensitive to scatter: approx. <20 mm collimation width or <10% error of scatter correction in a full field-of-view setting is needed to resolve BME. A mismatch in PDD decomposition calibration of ± 1 keV results in ~25% error in fat fraction estimates. In the wrist phantom study with MC scatter corrections, we were able to achieve ~0.79 mL true positive and ~0.06 mL false positive BME detection (compared to 0.85 mL true BME volume).

Conclusions

Detection of BME using DE CBCT with dual-layer FPD is feasible, but requires scatter mitigation, accurate scatter estimation, and robust spectral calibration.

Flat-panel conebeam CT in the clinic: history and current state

Research into conebeam CT concepts began as soon as the first clinical single-slice CT scanner was conceived. Early implementations of conebeam CT in the 1980s focused on high-contrast applications where concurrent high resolution ( < 200    μ m ), for visualization of small contrast-filled vessels, bones, or teeth, was an imaging requirement that could not be met by the contemporaneous CT scanners. However, the use of nonlinear imagers, e.g., x-ray image intensifiers, limited the clinical utility of the earliest diagnostic conebeam CT systems. The development of consumer-electronics large-area displays provided a technical foundation that was leveraged in the 1990s to first produce large-area digital x-ray detectors for use in radiography and then compact flat panels suitable for high-resolution and high-frame-rate conebeam CT. In this review, we show the concurrent evolution of digital flat panel (DFP) technology and clinical conebeam CT. We give a brief summary of conebeam CT reconstruction, followed by a brief review of the correction approaches for DFP-specific artifacts. The historical development and current status of flat-panel conebeam CT in four clinical areas-breast, fixed C-arm, image-guided radiation therapy, and extremity/head-is presented. Advances in DFP technology over the past two decades have led to improved visualization of high-contrast, high-resolution clinical tasks, and image quality now approaches the soft-tissue contrast resolution that is the standard in clinical CT. Future technical developments in DFPs will enable an even broader range of clinical applications; research in the arena of flat-panel CT shows no signs of slowing down.

Slot-scan dual-energy bone densitometry using motorized X-ray systems

PURPOSE

We investigate the feasibility of slot-scan dual-energy (DE) bone densitometry on motorized radiographic equipment. This approach will enable fast quantitative measurements of areal bone mineral density (aBMD) for opportunistic evaluation of osteoporosis.

METHODS

We investigated DE slot-scan protocols to obtain aBMD measurements at the lumbar spine (L-spine) and hip using a motorized x-ray platform capable of synchronized translation of the x-ray source and flat-panel detector (FPD). The slot dimension was 5 × 20 cm² . The DE slot views were processed as follows: (1) convolution kernel-based scatter correction, (2) unfiltered backprojection to tile the slots into long-length radiographs, and (3) projection-domain DE decomposition, consisting of an initial adipose-water decomposition in a bone-free region followed by water-CaHA decomposition with adjustment for adipose content. The accuracy and reproducibility of slot-scan aBMD measurements were investigated using a high-fidelity simulator of a robotic x-ray system (Siemens Multitom Rax) in a total of 48 body phantom realizations: four average bone density settings (cortical bone mass fraction: 10-40%), four body sizes (waist circumference, WC = 70-106 cm), and three lateral shifts of the body within the slot field of view (FOV) (centered and ±1 cm off-center). Experimental validations included: (1) x-ray test-bench feasibility study of adipose-water decomposition and (2) initial demonstration of slot-scan DE bone densitometry on the robotic x-ray system using the European Spine Phantom (ESP) with added attenuation (polymethyl methacrylate [PMMA] slabs) ranging 2 to 6 cm thick.

RESULTS

For the L-spine, the mean aBMD error across all WC settings ranged from 0.08 g/cm² for phantoms with average cortical bone fraction w_cortical  = 10% to ∼0.01 g/cm² for phantoms with w_cortical  = 40%. The L-spine aBMD measurements were fairly robust to changes in body size and positioning, e.g., coefficient of variation (CV) for L1 with w_cortical  = 30% was ∼0.034 for various WC and ∼0.02 for an obese patient (WC = 106 cm) changing lateral shift. For the hip, the mean aBMD error across all phantom configurations was about 0.07 g/cm² for a centered patient. The reproducibility of hip aBMD was slightly worse than in the L-spine (e.g., in the femoral neck, the CV with respect to changing WC was ∼0.13 for phantom realizations with w_cortical  = 30%) due to more challenging scatter estimation in the presence of an air-tissue interface within the slot FOV. The aBMD of the hip was therefore sensitive to lateral positioning of the patient, especially for obese patients: e.g., the CV with respect to patient lateral shift for femoral neck with WC = 106 cm and w_cortical  = 30% was 0.14. Empirical evaluations confirmed substantial reduction in aBMD errors with the proposed adipose estimation procedure and demonstrated robust aBMD measurements on the robotic x-ray system, with aBMD errors of ∼0.1 g/cm² across all three simulated ESP vertebrae and all added PMMA attenuator settings.

CONCLUSIONS

We demonstrated that accurate aBMD measurements can be obtained on a motorized FPD-based x-ray system using DE slot-scans with kernel-based scatter correction, backprojection-based slot view tiling, and DE decomposition with adipose correction.

Model-based dual-energy tomographic image reconstruction of objects containing known metal components

Dual-energy (DE) decomposition has been adopted in orthopedic imaging to measure bone composition and visualize intraarticular contrast enhancement. One of the potential applications involves monitoring of callus mineralization for longitudinal assessment of fracture healing. However, fracture repair usually involves internal fixation hardware that can generate significant artifacts in reconstructed images. To address this challenge, we develop a novel algorithm that combines simultaneous reconstruction-decomposition using a previously reported method for model-based material decomposition (MBMD) augmented by the known-component (KC) reconstruction framework to mitigate metal artifacts. We apply the proposed algorithm to simulated DE data representative of a dedicated extremity cone-beam CT (CBCT) employing an x-ray unit with three vertically arranged sources. The scanner generates DE data with non-coinciding high- and low-energy projection rays when the central source is operated at high tube potential and the peripheral sources at low potential. The proposed algorithm was validated using a digital extremity phantom containing varying concentrations of Ca-water mixtures and Ti implants. Decomposition accuracy was compared to MBMD without the KC model. The proposed method suppressed metal artifacts and yielded estimated Ca concentrations that approached the reconstructions of an implant-free phantom for most mixture regions. In the vicinity of simple components, the errors of Ca density estimates obtained by incorporating KC in MBMD were ∼1.5-5× lower than the errors of conventional MBMD; for cases with complex implants, the errors were ∼3-5× lower. In conclusion, the proposed method can achieve accurate bone mineral density measurements in the presence of metal implants using non-coinciding DE projections acquired on a multisource CBCT system.

Contextual loss based artifact removal method on CBCT image

Purpose

Cone beam computed tomography (CBCT) offers advantages such as high ray utilization rate, the same spatial resolution within and between slices, and high precision. It is one of the most actively studied topics in international computed tomography (CT) research. However, its application is hindered owing to scatter artifacts. This paper proposes a novel scatter artifact removal algorithm that is based on a convolutional neural network (CNN), where contextual loss is employed as the loss function.

Methods

In the proposed method, contextual loss is added to a simple CNN network to correct the CBCT artifacts in the pelvic region. The algorithm aims to learn the mapping from CBCT images to planning CT images. The 627 CBCT‐CT pairs of 11 patients were used to train the network, and the proposed algorithm was evaluated in terms of the mean absolute error (MAE), average peak signal‐to‐noise ratio (PSNR) and so on. The proposed method was compared with other methods to illustrate its effectiveness.

Results

The proposed method can remove artifacts (including streaking, shadowing, and cupping) in the CBCT image. Furthermore, key details such as the internal contours and texture information of the pelvic region are well preserved. Analysis of the average CT number, average MAE, and average PSNR indicated that the proposed method improved the image quality. The test results obtained with the chest data also indicated that the proposed method could be applied to other anatomies.

Conclusions

Although the CBCT‐CT image pairs are not completely matched at the pixel level, the method proposed in this paper can effectively correct the artifacts in the CBCT slices and improve the image quality. The average CT number of the regions of interest (including bones, skin) also exhibited a significant improvement. Furthermore, the proposed method can be applied to enhance the performance on such applications as dose estimation and segmentation.

Scatter correction based on adaptive photon path-based Monte Carlo simulation method in Multi-GPU platform

Monte Carlo (MC)-based simulation is the most precise method in scatter correction for Cone-beam CT (CBCT). Nonetheless, the existing MC methods cannot be fully applied in clinical due to its low efficiency. The traditional MC simulations perform calculations via a particle-by-particle scheme, which leads to high computation costs because abundant photons do not reach the X-ray detector in transport. The conventional approaches cannot control where the particle ends. Hence, it unavoidably waste lots of time in transporting numerous photons that have no contribution to the signal at the detector, yielding a low computational efficiency. To solve the problem, an innovative GPU-based Metropolis MC (gMMC) method was proposed. Compared with the traditional ones, the Metropolis based algorithm utilizes a path-by-path sampling method. The method can automatically control each particle path and eventually accelerate the convergence. In this paper, we firstly take planning CT image as prior information because of its precise CT value, and utilize gMMC to estimate scatter signal. Then the scatter signals are removed from the raw CBCT projections. Afterwards, FDK reconstruction is performed to obtain the corrected image,some accelerating strategies including reducing photon history number, pixels sampling, projection angles sampling and reconstructed image down-sampling achieve adaptive fast CBCT image reconstruction. For having high computational efficiency, we implemented the whole workflow on a 4-GPU workstation. In order to verify the feasibility of the the method, the experiment of several cases are conducted including simulation, phantom, and real patient cases. Results indicate that the image contrast becomes better, the scatter artifacts are eliminated. The maximum error (e_max), the minimum error (e_min), the 95th percentile error (e_95%), average error (¯e) are reduced from 264, 56, 14 and 21 HU to 28, 10, 3 and 7 HU in full-fan case, and from 387, 5, 19 and 95 HU to 39, 2, 2 and 6 HU in the half-fan case. In terms of computation time, the MC simulation time of all cases is within 2.5 seconds, and the total time is within 15 seconds.

Cone-beam imaging with tilted rotation axis: Method and performance evaluation

PURPOSE

The recently introduced robotic x-ray systems provide the flexibility to acquire cone-beam computed tomography (CBCT) data using customized, application-specific source-detector trajectories. We exploit this capability to mitigate the effects of x-ray scatter and noise in CBCT imaging of weight-bearing foot and cervical spine (C-spine) using scan orbits with a tilted rotation axis.

METHODS

We used an advanced CBCT simulator implementing accurate models of x-ray scatter, primary attenuation, and noise to investigate the effects of the orbital tilt angle in upright foot and C-spine imaging. The system model was parameterized using a laboratory version of a three-dimensional (3D) robotic x-ray system (Multitom RAX, Siemens Healthineers). We considered a generalized tilted axis scan configuration, where the detector remained parallel to patient's long body axis during the acquisition, but the elevation of source and detector was changing. A modified Feldkamp-Davis-Kress (FDK) algorithm was developed for reconstruction in this configuration, which departs from the FDK assumption of a detector that is perpendicular to the scan plane. The simulated foot scans involved source-detector distance (SDD) of 1386 mm, orbital tilt angles ranging 10° to 40°, and 400 views at 1 mAs/view and 0.5° increment; the C-spine scans involved -25° to -45° tilt angles, SDD of 1090 mm, and 202 views at 1.3 mAs and 1° increment The imaging performance was assessed by projection-domain measurements of the scatter-to-primary ratio (SPR) and by reconstruction-domain measurements of contrast, noise and generalized contrast-to-noise ratio (gCNR, accounting for both image noise and background nonuniformity) of the metatarsals (foot imaging) and cervical vertebrae (spine imaging). The effects of scatter correction were also compared for horizontal and tilted scans using an ideal Monte Carlo (MC)-based scatter correction and a frame-by-frame mean scatter correction.

RESULTS

The proposed modified FDK, involving projection resampling, mitigated streak artifacts caused by the misalignment between the filtering direction and the detector rows. For foot imaging (no grids), an optimized 20° tilted orbit reduced the maximum SPR from ~1.5 in a horizontal scan to <0.5. The gCNR of the second metatarsal was enhanced twofold compared to a horizontal orbit. For the C-spine (with vertical grids), imaging with a tilted orbit avoided highly attenuating x-ray paths through the lower cervical vertebrae and shoulders. A -35° tilted orbit yielded improved image quality and visualization of the lower cervical spine: the SPR of lower cervical vertebrae was reduced from ~10 (horizontal orbit) to <6 (tilted orbit), and the gCNR for C5-C7 increased by a factor of 2. Furthermore, tilted orbits showed potential benefits over horizontal orbits by enabling scatter correction with a simple frame-by-frame mean correction without substantial increase in noise-induced artifacts after the correction.

CONCLUSIONS

Tilted scan trajectories, enabled by the emerging robotic x-ray system technology, were optimized for CBCT imaging of foot and cervical spine using an advanced simulation framework. The results demonstrated the potential advantages of tilted axis orbits in mitigation of scatter artifacts and improving contrast-to-noise ratio in CBCT reconstructions.

Artificial intelligence in radiotherapy: a technological review

Radiation therapy (RT) is widely used to treat cancer. Technological advances in RT have occurred in the past 30 years. These advances, such as three-dimensional image guidance, intensity modulation, and robotics, created challenges and opportunities for the next breakthrough, in which artificial intelligence (AI) will possibly play important roles. AI will replace certain repetitive and labor-intensive tasks and improve the accuracy and consistency of others, particularly those with increased complexity because of technological advances. The improvement in efficiency and consistency is important to manage the increasing cancer patient burden to the society. Furthermore, AI may provide new functionalities that facilitate satisfactory RT. The functionalities include superior images for real-time intervention and adaptive and personalized RT. AI may effectively synthesize and analyze big data for such purposes. This review describes the RT workflow and identifies areas, including imaging, treatment planning, quality assurance, and outcome prediction, that benefit from AI. This review primarily focuses on deep-learning techniques, although conventional machine-learning techniques are also mentioned.

Cone-beam CT for imaging of the head/brain: Development and assessment of scanner prototype and reconstruction algorithms

PURPOSE

Our aim was to develop a high-quality, mobile cone-beam computed tomography (CBCT) scanner for point-of-care detection and monitoring of low-contrast, soft-tissue abnormalities in the head/brain, such as acute intracranial hemorrhage (ICH). This work presents an integrated framework of hardware and algorithmic advances for improving soft-tissue contrast resolution and evaluation of its technical performance with human subjects.

METHODS

Four configurations of a CBCT scanner prototype were designed and implemented to investigate key aspects of hardware (including system geometry, antiscatter grid, bowtie filter) and technique protocols. An integrated software pipeline (c.f., a serial cascade of algorithms) was developed for artifact correction (image lag, glare, beam hardening and x-ray scatter), motion compensation, and three-dimensional image (3D) reconstruction [penalized weighted least squares (PWLS), with a hardware-specific statistical noise model]. The PWLS method was extended in this work to accommodate multiple, independently moving regions with different resolution (to address both motion compensation and image truncation). Imaging performance was evaluated quantitatively and qualitatively with 41 human subjects in the neurosciences critical care unit (NCCU) at our institution.

RESULTS

The progression of four scanner configurations exhibited systematic improvement in the quality of raw data by variations in system geometry (source-detector distance), antiscatter grid, and bowtie filter. Quantitative assessment of CBCT images in 41 subjects demonstrated: ~70% reduction in image nonuniformity with artifact correction methods (lag, glare, beam hardening, and scatter); ~40% reduction in motion-induced streak artifacts via the multi-motion compensation method; and ~15% improvement in soft-tissue contrast-to-noise ratio (CNR) for PWLS compared to filtered backprojection (FBP) at matched resolution. Each of these components was important to improve contrast resolution for point-of-care cranial imaging.

CONCLUSIONS

This work presents the first application of a high-quality, point-of-care CBCT system for imaging of the head/ brain in a neurological critical care setting. Hardware configuration iterations and an integrated software pipeline for artifacts correction and PWLS reconstruction mitigated artifacts and noise to achieve image quality that could be valuable for point-of-care detection and monitoring of a variety of intracranial abnormalities, including ICH and hydrocephalus.

Quantitative Assessment of Weight-Bearing Fracture Biomechanics Using Extremity Cone-Beam CT

PURPOSE

We investigate an application of multisource extremity Cone-Beam CT (CBCT) with capability of weight-bearing tomographic imaging to obtain quantitative measurements of load-induced deformation of metal internal fixation hardware (e.g. tibial plate). Such measurements are desirable to improve the detection of delayed fusion or non-union of fractures, potentially facilitating earlier return to weight-bearing activities.

METHODS

To measure the deformation, we perform a deformable 3D-2D registration of a prior model of the implant to its CBCT projections under load-bearing. This Known-Component Registration (KC-Reg) framework avoids potential errors that emerge when the deformation is estimated directly from 3D reconstructions with metal artifacts. The 3D-2D registration involves a free-form deformable (FFD) point cloud model of the implant and a 3D cubic B-spline representation of the deformation. Gradient correlation is used as the optimization metric for the registration. The proposed approach was tested in experimental studies on the extremity CBCT system. A custom jig was designed to apply controlled axial loads to a fracture model, emulating weight-bearing imaging scenarios. Performance evaluation involved a Sawbone tibia phantom with an ~4 mm fracture gap. The model was fixed with a locking plate and imaged under five loading conditions. To investigate performance in the presence of confounding background gradients, additional experiments were performed with a pre-deformed femoral plate placed in a water bath with Ca bone mineral density inserts. Errors were measured using eight reference BBs for the tibial plate, and surface point distances for the femoral plate, where a prior model of deformed implant was available for comparison.

RESULTS

Both in the loaded tibial plate case and for the femoral plate with confounding background gradients, the proposed KC-Reg framework estimated implant deformations with errors of <0.2 mm for the majority of the investigated deformation magnitudes (error range 0.14 - 0.25 mm). The accuracy was comparable between 3D-2D registrations performed from 12 x-ray views and registrations obtained from as few as 3 views. This was likely enabled by the unique three-source x-ray unit on the extremity CBCT scanner, which implements two off-central-plane focal spots that provided oblique views of the field-of-view to aid implant pose estimation.

CONCLUSION

Accurate measurements of fracture hardware deformations under physiological weight-bearing are feasible using an extremity CBCT scanner and FFD 3D-2D registration. The resulting deformed implant models can be incorporated into tomographic reconstructions to reduce metal artifacts and improve quantification of the mineral content of fracture callus in CBCT volumes.

Technical Principles of Dual-Energy Cone Beam Computed Tomography and Clinical Applications for Radiation Therapy

PURPOSE

Medical imaging is an indispensable tool in radiotherapy for dose planning, image guidance and treatment monitoring. Cone beam CT (CBCT) is a low dose imaging technique with high spatial resolution capability as a direct by-product of using flat-panel detectors. However, certain issues such as x-ray scatter, beam hardening and other artifacts limit its utility to the verification of patient positioning using image-guided radiotherapy.

METHODS AND MATERIALS

Dual-energy (DE)-CBCT has recently demonstrated promise as an improved tool for tumor visualization in benchtop applications. It has the potential to improve soft-tissue contrast and reduce artifacts caused by beam hardening and metal. In this review, the practical aspects of developing a DE-CBCT based clinical and technical workflow are presented based on existing DE-CBCT literature and concepts adapted from the well-established library of work in DE-CT. Furthermore, the potential applications of DE-CBCT on its future role in radiotherapy are discussed.

RESULTS AND CONCLUSIONS

Based on current literature and an investigation of future applications, there is a clear potential for DE-CBCT technologies to be incorporated into radiotherapy. The applications of DE-CBCT include (but are not limited to): adaptive radiotherapy, brachytherapy, proton therapy, radiomics and theranostics.

Motion compensation in extremity cone-beam computed tomography

OBJECTIVES

To evaluate the improvement in extremity cone-beam computed tomography (CBCT) image quality in datasets with motion artifact using a motion compensation method based on maximizing image sharpness.

METHODS

Following IRB approval, retrospective analysis of 308 CBCT scans of lower extremities was performed by a fellowship-trained musculoskeletal radiologist to identify images with moderate to severe motion artifact. Twenty-four scans of 22 patients (18 male, four female; mean, 32 years old, range, 21-74 years old) were chosen for inclusion. Sharp (bone) and smooth (soft tissue) reconstructions were processed using the motion compensation algorithm. Two experts rated visualization of trabecular bone, cortical bone, joint spaces, and tendon on a nine-level Likert scale with and without motion compensation (a total of 96 datasets). Visual grading characteristics (VGC) was used to quantitatively determine the difference in image quality following motion compensation. Intra-class correlation coefficient (ICC) was obtained to assess inter-observer agreement.

RESULTS

Motion-compensated images exhibited appreciable reduction in artifacts. The observer study demonstrated the associated improvement in diagnostic quality. The fraction of cases receiving scores better than "Fair" increased from less than 10% without compensation to 40-70% following compensation, depending on the task. The area under the VGC curve was 0.75 (tendon) to 0.85 (cortical bone), confirming preference for motion compensated images. ICC values showed excellent agreement between readers before (ICC range, 0.8-0.91) and after motion compensation (ICC range, 0.92-0.97).

CONCLUSIONS

The motion compensation algorithm significantly improved the visualization of bone and soft tissue structures in extremity CBCT for cases exhibiting patient motion.

Quantitative Cone-Beam CT of Bone Mineral Density Using Model-Based Reconstruction

PURPOSE

We develop and validate a model-based framework for artifact correction and image reconstruction to enable application of Cone-Beam CT (CBCT) in quantitative assessment of bone mineral density (BMD). Compared to conventional quantitative CT, this approach does not require a BMD calibration phantom in the field-of-view during an object scan.

METHODS

The quantitative CBCT (qCBCT) imaging framework combined fast Monte Carlo (MC) scatter estimation, accurate models of detector response, and polyenergetic Poisson likelihood (PolyPL, Elbakri et al 2003). The underlying object model assumed that the tissues were ideal mixtures of water and calcium carbonate (CaCO3). Accuracy and reproducibility of qCBCT was evaluated in benchtop test-retest studies emulating a compact extremity CBCT system (axis-detector distance=56 cm, 90 kVp x-ray beam, ~16 mGy central dose). Various arrangements of Ca inserts (50-500 mg/mL) were placed in water cylinders of ~11 cm to ~15 cm diameter and scanned at multiple positions inside the field-of-view for a total of 20 configurations. In addition, a cadaveric ankle was imaged in five configurations (with and without Ca inserts and water bath). Coefficient of variation (CV) of BMD values across different experimental configurations was used to assess reproducibility under varying imaging conditions. The performance of the model-based qCBCT framework (MC + PolyPL) was compared to FDK with water beam hardening correction and MC scatter correction.

RESULTS

The PolyPL framework achieved accuracy of 20 mg/mL or better across all insert densities and experimental configurations. By comparison, the accuracy of the FDK-based BMD estimates deteriorated with higher mineralization, resulting in ~120 mg/mL error for a 500 mg/mL Ca insert. Additionally, the model-based approach mitigated residual streaks that were present in FDK reconstructions. The CV of both methods was ~15% at 50 mg/mL Ca and less than ~8% for higher density inserts, where the PolyPL framework achieved 20-25% lower CV than the FDK-based approach.

CONCLUSION

Accurate and reproducible BMD measurements can be achieved in extremity CBCT, supporting clinical applications in quantitative monitoring of fracture risk, osteoporosis treatment, and early osteoarthritis.

Personalization of X-Ray Tube Motion in Digital Breast Tomosynthesis Using Virtual Defrise Phantoms

In digital breast tomosynthesis (DBT), projection images are acquired as the x-ray tube rotates in the plane of the chest wall. We constructed a prototype next-generation tomosynthesis (NGT) system that has an additional component of tube motion in the perpendicular direction (i.e., posteroanterior motion). Our previous work demonstrated the advantages of the NGT system using the Defrise phantom. The reconstruction shows higher contrast and fewer blurring artifacts. To expand upon that work, this paper analyzes how image quality can be further improved by customizing the motion path of the x-ray tube based on the object being imaged. In simulations, phantoms are created with realistic 3D breast outlines based on an established model of the breast under compression. The phantoms are given an internal structure similar to a Defrise phantom. Two tissue types (fibroglandular and adipose) are arranged in a square-wave pattern. The reconstruction is analyzed as a binary classification task using thresholding to segment the two tissue types. At various thresholds, the classification of each voxel in the reconstruction is compared against the phantom, and receiver operating characteristic (ROC) curves are calculated. It is shown that the area under the ROC curve (AUC) is dependent on the x-ray tube trajectory. The trajectory that maximizes AUC differs between phantoms. In conclusion, this paper demonstrates that the acquisition geometry in DBT should be personalized to the object being imaged in order to optimize the image quality.

Quantitative cone-beam CT reconstruction with polyenergetic scatter model fusion

Scatter can account for large errors in cone-beam CT (CBCT) due to its wide field of view, and its complicated nature makes its compensation difficult. Iterative polyenergetic reconstruction algorithms offer the potential to provide quantitative imaging in CT, but they are usually incompatible with scatter contaminated measurements. In this work, we introduce a polyenergetic convolutional scatter model that is directly fused into the reconstruction process, and exploits information readily available at each iteration for a fraction of additional computational cost. We evaluate this method with numerical and real CBCT measurements, and show significantly enhanced electron density estimation and artifact mitigation over pre-calculated fast adaptive scatter kernel superposition (fASKS). We demonstrate our approach has two levels of benefit: reducing the bias introduced by estimating scatter prior to reconstruction; and adapting to the spectral and spatial properties of the specimen.

Cone-Beam CT image contrast and attenuation-map linearity improvement (CALI) for brain stereotactic radiosurgery procedures

A Contrast and Attenuation‐map Linearity Improvement (CALI) framework is proposed for cone‐beam CT (CBCT) images used for brain stereotactic radiosurgery (SRS). The proposed framework is tailored to improve soft tissue contrast of a new point‐of‐care image‐guided SRS system that employs a challenging half cone beam geometry, but can be readily reproduced on any CBCT platform. CALI includes a pre‐ and post‐processing step. In pre‐processing we apply a shading and beam hardening artifact correction to the projections, and in post‐processing step we correct the dome/capping artifact on reconstructed images caused by the spatial variations in X‐ray energy generated by the bowtie‐filter. The shading reduction together with the beam hardening and dome artifact correction algorithms aim to improve the linearity and accuracy of the CT‐numbers (CT#). The CALI framework was evaluated using CatPhan to quantify linearity, contrast‐to‐noise (CNR), and CT# accuracy, as well as subjectively on patient images acquired on a clinical system. Linearity of the reconstructed attenuation‐map was improved from 0.80 to 0.95. The CT# mean absolute measurement error was reduced from 76.1 to 26.9 HU. The CNR of the acrylic insert in the sensitometry module was improved from 1.8 to 7.8. The resulting clinical brain images showed substantial improvements in soft tissue contrast visibility, revealing structures such as ventricles which were otherwise undetectable in the original clinical images obtained from the system. The proposed reconstruction framework also improved CT# accuracy compared to the original images acquired on the system. For frameless image‐guided SRS, improving soft tissue visibility can facilitate evaluation of MR to CBCT co‐registration. Moreover, more accurate CT# may enable the use of CBCT for daily dose delivery measurements.

Image quality and dose characteristics for an O-arm intraoperative imaging system with model-based image reconstruction

PURPOSE

To assess the imaging performance and radiation dose characteristics of the O-arm CBCT imaging system (Medtronic Inc., Littleton MA) and demonstrate the potential for improved image quality and reduced dose via model-based image reconstruction (MBIR).

METHODS

Two main studies were performed to investigate previously unreported characteristics of the O-arm system. First is an investigation of dose and 3D image quality achieved with filtered back-projection (FBP) - including enhancements in geometric calibration, handling of lateral truncation and detector saturation, and incorporation of an isotropic apodization filter. Second is implementation of an MBIR algorithm based on Huber-penalized likelihood estimation (PLH) and investigation of image quality improvement at reduced dose. Each study involved measurements in quantitative phantoms as a basis for analysis of contrast-to-noise ratio and spatial resolution as well as imaging of a human cadaver to test the findings under realistic imaging conditions.

RESULTS

View-dependent calibration of system geometry improved the accuracy of reconstruction as quantified by the full-width at half maximum of the point-spread function - from 0.80 to 0.65 mm - and yielded subtle but perceptible improvement in high-contrast detail of bone (e.g., temporal bone). Standard technique protocols for the head and body imparted absorbed dose of 16 and 18 mGy, respectively. For low-to-medium contrast (<100 HU) imaging at fixed spatial resolution (1.3 mm edge-spread function) and fixed dose (6.7 mGy), PLH improved CNR over FBP by +48% in the head and +35% in the body. Evaluation at different dose levels demonstrated 30% increase in CNR at 62% of the dose in the head and 90% increase in CNR at 50% dose in the body.

CONCLUSIONS

A variety of improvements in FBP implementation (geometric calibration, truncation and saturation effects, and isotropic apodization) offer the potential for improved image quality and reduced radiation dose on the O-arm system. Further gains are possible with MBIR, including improved soft-tissue visualization, low-dose imaging protocols, and extension to methods that naturally incorporate prior information of patient anatomy and/or surgical instrumentation.

Multienergy element-resolved cone beam CT (MEER-CBCT) realized on a conventional CBCT platform

PURPOSE

Cone beam CT (CBCT) has been widely used in radiation therapy. However, its main application is still to acquire anatomical information for patient positioning. This study proposes a multienergy element-resolved (MEER) CBCT framework that employs energy-resolved data acquisition on a conventional CBCT platform and then simultaneously reconstructs images of x-ray attenuation coefficients, electron density relative to water (rED), and elemental composition (EC) to support advanced applications.

METHODS

The MEER-CBCT framework is realized on a Varian TrueBeam CBCT platform using a kVp-switching scanning scheme. A simultaneous image reconstruction and elemental decomposition model is formulated as an optimization problem. The objective function uses a least square term to enforce fidelity between x-ray attenuation coefficients and projection measurements. Spatial regularization is introduced via sparsity under a tight wavelet-frame transform. Consistency is imposed among rED, EC, and attenuation coefficients and inherently serves as a regularization term along the energy direction. The EC is further constrained by a sparse combination of ECs in a dictionary containing tissues commonly existing in humans. The optimization problem is solved by a novel alternating-direction minimization scheme. The MEER-CBCT framework was tested in a simulation study using an NCAT phantom and an experimental study using a Gammex phantom.

RESULTS

MEER-CBCT framework was successfully realized on a clinical Varian TrueBeam onboard CBCT platform with three energy channels of 80, 100, and 120 kVp. In the simulation study, the attenuation coefficient image achieved a structural similarity index of 0.98, compared to 0.61 for the image reconstructed by the conventional conjugate gradient least square (CGLS) algorithm, primarily because of reduction in artifacts. In the experimental study, the attenuation image obtained a contrast-to-noise ratio ≥60, much higher than that of CGLS results (~16) because of noise reduction. The median errors in rED and EC were 0.5% and 1.4% in the simulation study and 1.4% and 2.3% in the experimental study.

CONCLUSION

We proposed a novel MEER-CBCT framework realized on a clinical CBCT platform. Simulation and experimental studies demonstrated its capability to simultaneously reconstruct x-ray attenuation coefficient, rED, and EC images accurately.

Effect of detruncation on the accuracy of Monte Carlo-based scatter estimation in truncated CBCT

PURPOSE

The purpose of this study is to investigate the necessity of detruncation for scatter estimation of truncated cone-beam CT (CBCT) data and to evaluate different detruncation algorithms. Scattered radiation results in some of the most severe artifacts in CT and depends strongly on the size and the shape of the scanned object. Especially in CBCT systems the large cone-angle and the small detector-to-isocenter distance lead to a large amount of scatter detected, resulting in cupping artifacts, streak artifacts, and inaccurate CT-values. If a small field of measurement (FOM) is used, as it is often the case in CBCT systems, data are truncated in longitudinal and lateral direction. Since only truncated data are available as input for the scatter estimation, the already challenging correction of scatter artifacts becomes even more difficult.

METHODS

The following detruncation methods are compared and evaluated with respect to scatter estimation: constant detruncation, cosine detruncation, adaptive detruncation, and prior-based detruncation using anatomical data from a similar phantom or patient, also compared to the case where no detruncation was performed. Each of the resulting, detruncated reconstructions serve as input volume for a Monte Carlo (MC) scatter estimation and subsequent scatter correction. An evaluation is performed on a head simulation, measurements of a head phantom and a patient using a dental CBCT geometry with a FOM diameter of 11 cm. Additionally, a thorax phantom is measured to assess performance in a C-Arm geometry with a FOM of up to 20 cm.

RESULTS

If scatter estimation is based on simple detruncation algorithms like a constant or a cosine detruncation scatter is estimated inaccurately, resulting in incorrect CT-values as well as streak artifacts in the corrected volume. For the dental CBCT phantom measurement CT-values for soft tissue were corrected from -204 HU (no scatter correction) to -87 HU (no detruncation), -218 HU (constant detruncation), -141 HU (cosine detruncation), -91 HU (adaptive detruncation), -34 HU (prior-based detruncation using a different prior) and -24 HU (prior-based detruncation using the identical prior) for a reference value of -26 HU measured in slit scan mode. In all cases the prior-based detruncation results in the best scatter correction, followed by the adaptive detruncation, as these algorithms provide a rather accurate model of high-density structures outside the FOM, compared to a simple constant or a cosine detruncation.

CONCLUSIONS

Our contribution is twofold: first we give a comprehensive comparison of various detruncation methods for the purpose of scatter estimation. We find that the choice of the detruncation method has a significant influence on the quality of MC-based scatter correction. Simple or no detruncation is often insufficient for artifact removal and results in inaccurate CT-values. On the contrary, prior-based detruncation can achieve a high CT-value accuracy and nearly artifact-free volumes from truncated CBCT data when combined with other state-of-the-art artifact corrections. Secondly, we show that prior-based detruncation is effective even with data from a different patient or phantom. The fact that data completion does not require data from the same patient dramatically increases the applicability and usability of this scatter estimation.

Acuros CTS: A fast, linear Boltzmann transport equation solver for computed tomography scatter - Part I: Core algorithms and validation

PURPOSE

To describe Acuros^® CTS, a new software tool for rapidly and accurately estimating scatter in x-ray projection images by deterministically solving the linear Boltzmann transport equation (LBTE).

METHODS

The LBTE describes the behavior of particles as they interact with an object across spatial, energy, and directional (propagation) domains. Acuros CTS deterministically solves the LBTE by modeling photon transport associated with an x-ray projection in three main steps: (a) Ray tracing photons from the x-ray source into the object where they experience their first scattering event and form scattering sources. (b) Propagating photons from their first scattering sources across the object in all directions to form second scattering sources, then repeating this process until all high-order scattering sources are computed using the source iteration method. (c) Ray-tracing photons from scattering sources within the object to the detector, accounting for the detector's energy and anti-scatter grid responses. To make this process computationally tractable, a combination of analytical and discrete methods is applied. The three domains are discretized using the Linear Discontinuous Finite Elements, Multigroup, and Discrete Ordinates methods, respectively, which confer the ability to maintain the accuracy of a continuous solution. Furthermore, through the implementation in CUDA, we sought to exploit the parallel computing capabilities of graphics processing units (GPUs) to achieve the speeds required for clinical utilization. Acuros CTS was validated against Geant4 Monte Carlo simulations using two digital phantoms: (a) a water phantom containing lung, air, and bone inserts (WLAB phantom) and (b) a pelvis phantom derived from a clinical CT dataset. For these studies, we modeled the TrueBeam^® (Varian Medical Systems, Palo Alto, CA) kV imaging system with a source energy of 125 kVp. The imager comprised a 600 μm-thick Cesium Iodide (CsI) scintillator and a 10:1 one-dimensional anti-scatter grid. For the WLAB studies, the full-fan geometry without a bowtie filter was used (with and without the anti-scatter grid). For the pelvis phantom studies, a half-fan geometry with bowtie was used (with the anti-scatter grid). Scattered and primary photon fluences and energies deposited in the detector were recorded.

RESULTS

The Acuros CTS and Monte Carlo results demonstrated excellent agreement. For the WLAB studies, the average percent difference between the Monte Carlo- and Acuros-generated scattered photon fluences at the face of the detector was -0.7%. After including the detector response, the average percent differences between the Monte Carlo- and Acuros-generated scatter fractions (SF) were -0.1% without the grid and 0.6% with the grid. For the digital pelvis simulation, the Monte Carlo- and Acuros-generated SFs agreed to within 0.1% on average, despite the scatter-to-primary ratios (SPRs) being as high as 5.5. The Acuros CTS computation time for each scatter image was ~1 s using a single GPU.

CONCLUSIONS

Acuros CTS enables a fast and accurate calculation of scatter images by deterministically solving the LBTE thus offering a computationally attractive alternative to Monte Carlo methods. Part II describes the application of Acuros CTS to scatter correction of CBCT scans on the TrueBeam system.

Low-Dose CT Perfusion of the Liver using Reconstruction of Difference

Liver CT perfusion (CTP) is used in the detection, staging, and treatment response analysis of hepatic diseases. Unfortunately, CTP radiation exposures is significant, limiting more widespread use. Traditional CTP data processing reconstructs individual temporal samples, ignoring a large amount of shared anatomical information between temporal samples, suggesting opportunities for improved data processing. We adopt a prior-image-based reconstruction approach called Reconstruction of Difference (RoD) to enable low-exposure CTP acquisition. RoD differs from many algorithms by directly estimating the attenuation changes between the current patient state and a prior CT volume. We propose to use a high-fidelity unenhanced baseline CT image to integrate prior anatomical knowledge into subsequent data reconstructions. Using simulation studies based on a 4D digital anthropomorphic phantom with realistic time-attenuation curves, we compare RoD with conventional filtered-backprojection, penalized-likelihood estimation, and prior image penalized-likelihood estimation. We evaluate each method in comparisons of reconstructions at individual time points, accuracy of estimated time-attenuation curves, and in an analysis of common perfusion metric maps including hepatic arterial perfusion, hepatic portal perfusion, perfusion index, and time-to-peak. Results suggest that RoD enables significant exposure reductions, outperforming standard and more sophisticated model-based reconstruction, making RoD a potentially important tool to enable low-dose liver CTP.

Modeling and evaluation of a high-resolution CMOS detector for cone-beam CT of the extremities

PURPOSE

Quantitative assessment of trabecular bone microarchitecture in extremity cone-beam CT (CBCT) would benefit from the high spatial resolution, low electronic noise, and fast scan time provided by complementary metal-oxide semiconductor (CMOS) x-ray detectors. We investigate the performance of CMOS sensors in extremity CBCT, in particular with respect to potential advantages of thin (<0.7 mm) scintillators offering higher spatial resolution.

METHODS

A cascaded systems model of a CMOS x-ray detector incorporating the effects of CsI:Tl scintillator thickness was developed. Simulation studies were performed using nominal extremity CBCT acquisition protocols (90 kVp, 0.126 mAs/projection). A range of scintillator thickness (0.35-0.75 mm), pixel size (0.05-0.4 mm), focal spot size (0.05-0.7 mm), magnification (1.1-2.1), and dose (15-40 mGy) was considered. The detectability index was evaluated for both CMOS and a-Si:H flat-panel detector (FPD) configurations for a range of imaging tasks emphasizing spatial frequencies associated with feature size aobj. Experimental validation was performed on a CBCT test bench in the geometry of a compact orthopedic CBCT system (SAD = 43.1 cm, SDD = 56.0 cm, matching that of the Carestream OnSight 3D system). The test-bench studies involved a 0.3 mm focal spot x-ray source and two CMOS detectors (Dalsa Xineos-3030HR, 0.099 mm pixel pitch) - one with the standard CsI:Tl thickness of 0.7 mm (C700) and one with a custom 0.4 mm thick scintillator (C400). Measurements of modulation transfer function (MTF), detective quantum efficiency (DQE), and CBCT scans of a cadaveric knee (15 mGy) were obtained for each detector.

RESULTS

Optimal detectability for high-frequency tasks (feature size of ~0.06 mm, consistent with the size of trabeculae) was ~4× for the C700 CMOS detector compared to the a-Si:H FPD at nominal system geometry of extremity CBCT. This is due to ~5× lower electronic noise of a CMOS sensor, which enables input quantum-limited imaging at smaller pixel size. Optimal pixel size for high-frequency tasks was <0.1 mm for a CMOS, compared to ~0.14 mm for an a-Si:H FPD. For this fine pixel pitch, detectability of fine features could be improved by using a thinner scintillator to reduce light spread blur. A 22% increase in detectability of 0.06 mm features was found for the C400 configuration compared to C700. An improvement in the frequency at 50% modulation (f₅₀ ) of MTF was measured, increasing from 1.8 lp/mm for C700 to 2.5 lp/mm for C400. The C400 configuration also achieved equivalent or better DQE as C700 for frequencies above ~2 mm^-1 . Images of cadaver specimens confirmed improved visualization of trabeculae with the C400 sensor.

CONCLUSIONS

The small pixel size of CMOS detectors yields improved performance in high-resolution extremity CBCT compared to a-Si:H FPDs, particularly when coupled with a custom 0.4 mm thick scintillator. The results indicate that adoption of a CMOS detector in extremity CBCT can benefit applications in quantitative imaging of trabecular microstructure in humans.

Task-based statistical image reconstruction for high-quality cone-beam CT

Task-based analysis of medical imaging performance underlies many ongoing efforts in the development of new imaging systems. In statistical image reconstruction, regularization is often formulated in terms to encourage smoothness and/or sharpness (e.g. a linear, quadratic, or Huber penalty) but without explicit formulation of the task. We propose an alternative regularization approach in which a spatially varying penalty is determined that maximizes task-based imaging performance at every location in a 3D image. We apply the method to model-based image reconstruction (MBIR-viz., penalized weighted least-squares, PWLS) in cone-beam CT (CBCT) of the head, focusing on the task of detecting a small, low-contrast intracranial hemorrhage (ICH), and we test the performance of the algorithm in the context of a recently developed CBCT prototype for point-of-care imaging of brain injury. Theoretical predictions of local spatial resolution and noise are computed via an optimization by which regularization (specifically, the quadratic penalty strength) is allowed to vary throughout the image to maximize local task-based detectability index ([Formula: see text]). Simulation studies and test-bench experiments were performed using an anthropomorphic head phantom. Three PWLS implementations were tested: conventional (constant) penalty; a certainty-based penalty derived to enforce constant point-spread function, PSF; and the task-based penalty derived to maximize local detectability at each location. Conventional (constant) regularization exhibited a fairly strong degree of spatial variation in [Formula: see text], and the certainty-based method achieved uniform PSF, but each exhibited a reduction in detectability compared to the task-based method, which improved detectability up to ~15%. The improvement was strongest in areas of high attenuation (skull base), where the conventional and certainty-based methods tended to over-smooth the data. The task-driven reconstruction method presents a promising regularization method in MBIR by explicitly incorporating task-based imaging performance as the objective. The results demonstrate improved ICH conspicuity and support the development of high-quality CBCT systems.

A review of GPU-based medical image reconstruction

Tomographic image reconstruction is a computationally demanding task, even more so when advanced models are used to describe a more complete and accurate picture of the image formation process. Such advanced modeling and reconstruction algorithms can lead to better images, often with less dose, but at the price of long calculation times that are hardly compatible with clinical workflows. Fortunately, reconstruction tasks can often be executed advantageously on Graphics Processing Units (GPUs), which are exploited as massively parallel computational engines. This review paper focuses on recent developments made in GPU-based medical image reconstruction, from a CT, PET, SPECT, MRI and US perspective. Strategies and approaches to get the most out of GPUs in image reconstruction are presented as well as innovative applications arising from an increased computing capacity. The future of GPU-based image reconstruction is also envisioned, based on current trends in high-performance computing.

Technical assessment of a prototype cone-beam CT system for imaging of acute intracranial hemorrhage

PURPOSE

A cone-beam CT scanner has been developed for detection and monitoring of traumatic brain injury and acute intracranial hemorrhage (ICH) at the point of care. This work presents a technical assessment of imaging performance and dose for the scanner in phantom and cadaver studies as a prerequisite to clinical translation.

METHODS

The scanner incorporates a compact, rotating-anode x-ray source and a flat-panel detector (43 × 43 cm²) on a mobile U-arm gantry with source-axis distance = 550 mm and source-detector distance = 1000 mm. Central and peripheral doses were measured in 16 cm diameter CTDI phantoms using a 0.6 cm³ Farmer ionization chamber for various scan techniques and as a function of longitudinal position, including out of field. Spatial resolution, contrast, noise, and image uniformity were assessed in quantitative and anthropomorphic head phantoms. Two reconstruction protocols were evaluated, including filtered backprojection (FBP) for high-resolution bone imaging and penalized weighted least squares (PWLS) reconstruction for low-contrast soft tissue (ICH) visualization. A fresh cadaver was imaged with and without simulated ICH using the scanner as well as a diagnostic multidetector CT (MDCT) scanner using a standard head protocol. Images were interpreted by a fellowship-trained neuroradiologist for imaging tasks of ICH detection, gray-white-CSF differentiation, detection of midline shift, and fracture detection.

RESULTS

The nominal scan protocol involved 720 projections acquired over a 360° orbit at 100 kV and 216 mAs, giving a dose (weighted CTDI) of 22.8 mGy (∼1.2 mSv effective dose). Out-of-field dose decreased to <10% within 6 cm of the field edge (approximate to the thyroid position). Image uniformity demonstrated <1% variation between the edge of the field (near the cranium) and center of the image. The high-resolution FBP reconstruction protocol showed ∼0.9 mm point spread function (PSF) full-width at half-maximum (FWHM). The smooth PWLS reconstruction protocol yielded ∼1.2 mm PSF FWHM and contrast-to-noise ratio exceeding 5.7 in ∼50 HU spherical ICH, resulting in conspicuous depiction of ICH down to ∼2 mm (the smallest diameter investigated). Cadaver images demonstrated good differentiation of brain and CSF (sufficient, but inferior to MDCT, recognizing that the CBCT dose was one-third that of MDCT), excellent visualization of cranial sutures and fracture (potentially superior to MDCT), clear detection of midline shift, and conspicuous detection of ICH.

CONCLUSIONS

Technical assessment of the prototype demonstrates dose characteristics and imaging performance consistent with point-of-care detection and monitoring of head injury-most notably, conspicuous detection of ICH-and supports translation of the system to clinical studies.

Multi-resolution statistical image reconstruction for mitigation of truncation effects: application to cone-beam CT of the head

A prototype cone-beam CT (CBCT) head scanner featuring model-based iterative reconstruction (MBIR) has been recently developed and demonstrated the potential for reliable detection of acute intracranial hemorrhage (ICH), which is vital to diagnosis of traumatic brain injury and hemorrhagic stroke. However, data truncation (e.g. due to the head holder) can result in artifacts that reduce image uniformity and challenge ICH detection. We propose a multi-resolution MBIR method with an extended reconstruction field of view (RFOV) to mitigate truncation effects in CBCT of the head. The image volume includes a fine voxel size in the (inner) nontruncated region and a coarse voxel size in the (outer) truncated region. This multi-resolution scheme allows extension of the RFOV to mitigate truncation effects while introducing minimal increase in computational complexity. The multi-resolution method was incorporated in a penalized weighted least-squares (PWLS) reconstruction framework previously developed for CBCT of the head. Experiments involving an anthropomorphic head phantom with truncation due to a carbon-fiber holder were shown to result in severe artifacts in conventional single-resolution PWLS, whereas extending the RFOV within the multi-resolution framework strongly reduced truncation artifacts. For the same extended RFOV, the multi-resolution approach reduced computation time compared to the single-resolution approach (viz. time reduced by 40.7%, 83.0%, and over 95% for an image volume of 600³, 800³, 1000³ voxels). Algorithm parameters (e.g. regularization strength, the ratio of the fine and coarse voxel size, and RFOV size) were investigated to guide reliable parameter selection. The findings provide a promising method for truncation artifact reduction in CBCT and may be useful for other MBIR methods and applications for which truncation is a challenge.

Evaluation of detector readout gain mode and bowtie filters for cone-beam CT imaging of the head

The effects of detector readout gain mode and bowtie filters on cone-beam CT (CBCT) image quality and dose were characterized for a new CBCT system developed for point-of-care imaging of the head, with potential application to diagnosis of traumatic brain injury, intracranial hemorrhage (ICH), and stroke. A detector performance model was extended to include the effects of detector readout gain on electronic digitization noise. The noise performance for high-gain (HG), low-gain (LG), and dual-gain (DG) detector readout was evaluated, and the benefit associated with HG mode in regions free from detector saturation was quantified. Such benefit could be realized (without detector saturation) either via DG mode or by incorporation of a bowtie filter. Therefore, three bowtie filters were investigated that varied in thickness and curvature. A polyenergetic gain correction method was developed to equalize the detector response between the flood-field and projection data in the presence of a bowtie. The effect of bowtie filters on dose, scatter-to-primary ratio, contrast, and noise was quantified in phantom studies, and results were compared to a high-speed Monte Carlo (MC) simulation to characterize x-ray scatter and dose distributions in the head. Imaging in DG mode improved the contrast-to-noise ratio (CNR) by ~15% compared to LG mode at a dose (D 0, measured at the center of a 16 cm CTDI phantom) of 19 mGy. MC dose calculations agreed with CTDI measurements and showed that bowtie filters reduce peripheral dose by as much as 50% at the same central dose. Bowtie filters were found to increase the CNR per unit square-root dose near the center of the image by ~5-20% depending on bowtie thickness, but reduced CNR in the periphery by ~10-40%. Images acquired at equal CTDIw with and without a bowtie demonstrated a 24% increase in CNR at the center of an anthropomorphic head phantom. Combining a thick bowtie filter with a short arc (180°  +  fan angle) scan centered on the posterior of the head reduced dose to the eye lens by up to 90%. Acquisition in DG mode (without a bowtie filter) was beneficial to the detection of small, low contrast lesions (e.g. subtle ICH) in CBCT. While bowtie filters were found to reduce dose, mitigate sensor saturation at the periphery in HG mode, and improve CNR at the center of the image, the image quality at the periphery was slightly reduced compared to DG mode, and the use of a bowtie required careful implementation of the polyenergetic flood-field correction to avoid artifacts.

X-Ray Scatter Correction on Soft Tissue Images for Portable Cone Beam CT

Soft tissue images from portable cone beam computed tomography (CBCT) scanners can be used for diagnosis and detection of tumor, cancer, intracerebral hemorrhage, and so forth. Due to large field of view, X-ray scattering which is the main cause of artifacts degrades image quality, such as cupping artifacts, CT number inaccuracy, and low contrast, especially on soft tissue images. In this work, we propose the X-ray scatter correction method for improving soft tissue images. The X-ray scatter correction scheme to estimate X-ray scatter signals is based on the deconvolution technique using the maximum likelihood estimation maximization (MLEM) method. The scatter kernels are obtained by simulating the PMMA sheet on the Monte Carlo simulation (MCS) software. In the experiment, we used the QRM phantom to quantitatively compare with fan-beam CT (FBCT) data in terms of CT number values, contrast to noise ratio, cupping artifacts, and low contrast detectability. Moreover, the PH3 angiography phantom was also used to mimic human soft tissues in the brain. The reconstructed images with our proposed scatter correction show significant improvement on image quality. Thus the proposed scatter correction technique has high potential to detect soft tissues in the brain.

Model-based iterative reconstruction for flat-panel cone-beam CT with focal spot blur, detector blur, and correlated noise

While model-based reconstruction methods have been successfully applied to flat-panel cone-beam CT (FP-CBCT) systems, typical implementations ignore both spatial correlations in the projection data as well as system blurs due to the detector and focal spot in the x-ray source. In this work, we develop a forward model for flat-panel-based systems that includes blur and noise correlation associated with finite focal spot size and an indirect detector (e.g. scintillator). This forward model is used to develop a staged reconstruction framework where projection data are deconvolved and log-transformed, followed by a generalized least-squares reconstruction that utilizes a non-diagonal statistical weighting to account for the correlation that arises from the acquisition and data processing chain. We investigate the performance of this novel reconstruction approach in both simulated data and in CBCT test-bench data. In comparison to traditional filtered backprojection and model-based methods that ignore noise correlation, the proposed approach yields a superior noise-resolution tradeoff. For example, for a system with 0.34 mm FWHM scintillator blur and 0.70 FWHM focal spot blur, using the correlated noise model instead of an uncorrelated noise model increased resolution by 42% (with variance matched at 6.9  ×  10(-8) mm(-2)). While this advantage holds across a wide range of systems with differing blur characteristics, the improvements are greatest for systems where source blur is larger than detector blur.

A scatter correction method for contrast-enhanced dual-energy digital breast tomosynthesis

Contrast-enhanced dual energy digital breast tomosynthesis (CE-DE-DBT) is designed to image iodinated masses while suppressing breast anatomical background. Scatter is a problem, especially for high energy acquisition, in that it causes severe cupping artifact and iodine quantitation errors. We propose a patient specific scatter correction (SC) algorithm for CE-DE-DBT. The empirical algorithm works by interpolating scatter data outside the breast shadow into an estimate within the breast shadow. The interpolated estimate is further improved by operations that use an easily obtainable (from phantoms) table of scatter-to-primary-ratios (SPR)--a single SPR value for each breast thickness and acquisition angle. We validated our SC algorithm for two breast emulating phantoms by comparing SPR from our SC algorithm to that measured using a beam-passing pinhole array plate. The error in our SC computed SPR, averaged over acquisition angle and image location, was about 5%, with slightly worse errors for thicker phantoms. The SC projection data, reconstructed using OS-SART, showed a large degree of decupping. We also observed that SC removed the dependence of iodine quantitation on phantom thickness. We applied the SC algorithm to a CE-DE-mammographic patient image with a biopsy confirmed tumor at the breast periphery. In the image without SC, the contrast enhanced tumor was masked by the cupping artifact. With our SC, the tumor was easily visible. An interpolation-based SC was proposed by (Siewerdsen et al 2006 Med. Phys. 33 187-97) for cone-beam CT (CBCT), but our algorithm and application differ in several respects. Other relevant SC techniques include Monte-Carlo and convolution-based methods for CBCT, storage of a precomputed library of scatter maps for DBT, and patient acquisition with a beam-passing pinhole array for breast CT. Our SC algorithm can be accomplished in clinically acceptable times, requires no additional imaging hardware or extra patient dose and is easily transportable between sites.

Statistical reconstruction for cone-beam CT with a post-artifact-correction noise model: application to high-quality head imaging

Non-contrast CT reliably detects fresh blood in the brain and is the current front-line imaging modality for intracranial hemorrhage such as that occurring in acute traumatic brain injury (contrast ~40-80 HU, size  >  1 mm). We are developing flat-panel detector (FPD) cone-beam CT (CBCT) to facilitate such diagnosis in a low-cost, mobile platform suitable for point-of-care deployment. Such a system may offer benefits in the ICU, urgent care/concussion clinic, ambulance, and sports and military theatres. However, current FPD-CBCT systems face significant challenges that confound low-contrast, soft-tissue imaging. Artifact correction can overcome major sources of bias in FPD-CBCT but imparts noise amplification in filtered backprojection (FBP). Model-based reconstruction improves soft-tissue image quality compared to FBP by leveraging a high-fidelity forward model and image regularization. In this work, we develop a novel penalized weighted least-squares (PWLS) image reconstruction method with a noise model that includes accurate modeling of the noise characteristics associated with the two dominant artifact corrections (scatter and beam-hardening) in CBCT and utilizes modified weights to compensate for noise amplification imparted by each correction. Experiments included real data acquired on a FPD-CBCT test-bench and an anthropomorphic head phantom emulating intra-parenchymal hemorrhage. The proposed PWLS method demonstrated superior noise-resolution tradeoffs in comparison to FBP and PWLS with conventional weights (viz. at matched 0.50 mm spatial resolution, CNR = 11.9 compared to CNR = 5.6 and CNR = 9.9, respectively) and substantially reduced image noise especially in challenging regions such as skull base. The results support the hypothesis that with high-fidelity artifact correction and statistical reconstruction using an accurate post-artifact-correction noise model, FPD-CBCT can achieve image quality allowing reliable detection of intracranial hemorrhage.

A practical cone-beam CT scatter correction method with optimized Monte Carlo simulations for image-guided radiation therapy

Cone-beam CT (CBCT) has become the standard image guidance tool for patient setup in image-guided radiation therapy. However, due to its large illumination field, scattered photons severely degrade its image quality. While kernel-based scatter correction methods have been used routinely in the clinic, it is still desirable to develop Monte Carlo (MC) simulation-based methods due to their accuracy. However, the high computational burden of the MC method has prevented routine clinical application. This paper reports our recent development of a practical method of MC-based scatter estimation and removal for CBCT. In contrast with conventional MC approaches that estimate scatter signals using a scatter-contaminated CBCT image, our method used a planning CT image for MC simulation, which has the advantages of accurate image intensity and absence of image truncation. In our method, the planning CT was first rigidly registered with the CBCT. Scatter signals were then estimated via MC simulation. After scatter signals were removed from the raw CBCT projections, a corrected CBCT image was reconstructed. The entire workflow was implemented on a GPU platform for high computational efficiency. Strategies such as projection denoising, CT image downsampling, and interpolation along the angular direction were employed to further enhance the calculation speed. We studied the impact of key parameters in the workflow on the resulting accuracy and efficiency, based on which the optimal parameter values were determined. Our method was evaluated in numerical simulation, phantom, and real patient cases. In the simulation cases, our method reduced mean HU errors from 44 to 3 HU and from 78 to 9 HU in the full-fan and the half-fan cases, respectively. In both the phantom and the patient cases, image artifacts caused by scatter, such as ring artifacts around the bowtie area, were reduced. With all the techniques employed, we achieved computation time of less than 30 s including the time for both the scatter estimation and CBCT reconstruction steps. The efficacy of our method and its high computational efficiency make our method attractive for clinical use.

Cone-Beam CT of Traumatic Brain Injury Using Statistical Reconstruction with a Post-Artifact-Correction Noise Model

Traumatic brain injury (TBI) is a major cause of death and disability. The current front-line imaging modality for TBI detection is CT, which reliably detects intracranial hemorrhage (fresh blood contrast 30-50 HU, size down to 1 mm) in non-contrast-enhanced exams. Compared to CT, flat-panel detector (FPD) cone-beam CT (CBCT) systems offer lower cost, greater portability, and smaller footprint suitable for point-of-care deployment. We are developing FPD-CBCT to facilitate TBI detection at the point-of-care such as in emergent, ambulance, sports, and military applications. However, current FPD-CBCT systems generally face challenges in low-contrast, soft-tissue imaging. Model-based reconstruction can improve image quality in soft-tissue imaging compared to conventional filtered backprojection (FBP) by leveraging high-fidelity forward model and sophisticated regularization. In FPD-CBCT TBI imaging, measurement noise characteristics undergo substantial change following artifact correction, resulting in non-negligible noise amplification. In this work, we extend the penalized weighted least-squares (PWLS) image reconstruction to include the two dominant artifact corrections (scatter and beam hardening) in FPD-CBCT TBI imaging by correctly modeling the variance change following each correction. Experiments were performed on a CBCT test-bench using an anthropomorphic phantom emulating intra-parenchymal hemorrhage in acute TBI, and the proposed method demonstrated an improvement in blood-brain contrast-to-noise ratio (CNR = 14.2) compared to FBP (CNR = 9.6) and PWLS using conventional weights (CNR = 11.6) at fixed spatial resolution (1 mm edge-spread width at the target contrast). The results support the hypothesis that FPD-CBCT can fulfill the image quality requirements for reliable TBI detection, using high-fidelity artifact correction and statistical reconstruction with accurate post-artifact-correction noise models.