Reduction of noise-induced streak artifacts in X-ray computed tomography through spline-based penalized-likelihood sinogram smoothing

Deep learning-based motion compensation for four-dimensional cone-beam computed tomography (4D-CBCT) reconstruction

BACKGROUND

Motion-compensated (MoCo) reconstruction shows great promise in improving four-dimensional cone-beam computed tomography (4D-CBCT) image quality. MoCo reconstruction for a 4D-CBCT could be more accurate using motion information at the CBCT imaging time than that obtained from previous 4D-CT scans. However, such data-driven approaches are hampered by the quality of initial 4D-CBCT images used for motion modeling.

PURPOSE

This study aims to develop a deep-learning method to generate high-quality motion models for MoCo reconstruction to improve the quality of final 4D-CBCT images.

METHODS

A 3D artifact-reduction convolutional neural network (CNN) was proposed to improve conventional phase-correlated Feldkamp-Davis-Kress (PCF) reconstructions by reducing undersampling-induced streaking artifacts while maintaining motion information. The CNN-generated artifact-mitigated 4D-CBCT images (CNN enhanced) were then used to build a motion model which was used by MoCo reconstruction (CNN+MoCo). The proposed procedure was evaluated using in-vivo patient datasets, an extended cardiac-torso (XCAT) phantom, and the public SPARE challenge datasets. The quality of reconstructed images for XCAT phantom and SPARE datasets was quantitatively assessed using root-mean-square-error (RMSE) and normalized cross-correlation (NCC).

RESULTS

The trained CNN effectively reduced the streaking artifacts of PCF CBCT images for all datasets. More detailed structures can be recovered using the proposed CNN+MoCo reconstruction procedure. XCAT phantom experiments showed that the accuracy of estimated motion model using CNN enhanced images was greatly improved over PCF. CNN+MoCo showed lower RMSE and higher NCC compared to PCF, CNN enhanced and conventional MoCo. For the SPARE datasets, the average (± standard deviation) RMSE in mm^-1 for body region of PCF, CNN enhanced, conventional MoCo and CNN+MoCo were 0.0040 ± 0.0009, 0.0029 ± 0.0002, 0.0024 ± 0.0003 and 0.0021 ± 0.0003. Corresponding NCC were 0.84 ± 0.05, 0.91 ± 0.05, 0.91 ± 0.05 and 0.93 ± 0.04.

CONCLUSIONS

CNN-based artifact reduction can substantially reduce the artifacts in the initial 4D-CBCT images. The improved images could be used to enhance the motion modeling and ultimately improve the quality of the final 4D-CBCT images reconstructed using MoCo.

STEDNet: Swin transformer-based encoder-decoder network for noise reduction in low-dose CT

BACKGROUND

Low-dose computed tomography (LDCT) can reduce the dose of X-ray radiation, making it increasingly significant for routine clinical diagnosis and treatment planning. However, the noise introduced by low-dose X-ray exposure degrades the quality of CT images, affecting the accuracy of clinical diagnosis. Purpose The noises, artifacts, and high-frequency components are similarly distributed in LDCT images. Transformer can capture global context information in an attentional manner to create distant dependencies on targets and extract more powerful features. In this paper, we reduce the impact of image errors on the ability to retain detailed information and improve the noise suppression performance by fully mining the distribution characteristics of image information.

METHODS

This paper proposed an LDCT noise and artifact suppressing network based on Swin Transformer. The network includes a noise extraction sub-network and a noise removal sub-network. The noise extraction and removal capability are improved using a coarse extraction network of high-frequency features based on full convolution. The noise removal sub-network improves the network's ability to extract relevant image features by using a Swin Transformer with a shift window as an encoder-decoder and skip connections for global feature fusion. Also, the perceptual field is extended by extracting multi-scale features of the images to recover the spatial resolution of the feature maps. The network uses a loss constraint with a combination of L1 and MS-SSIM to improve and ensure the stability and denoising effect of the network.

RESULTS

The denoising ability and clinical applicability of the methods were tested using clinical datasets. Compared with DnCNN, RED-CNN, CBDNet and TSCN, the STEDNet method shows a better denoising effect on RMSE and PSNR. The STEDNet method effectively removes image noise and preserves the image structure to the maximum extent, making the reconstructed image closest to the NDCT image. The subjective and objective analysis of several sets of experiments shows that the method in this paper can effectively maintain the structure, edges, and textures of the denoised images while having good noise suppression performance. In the real data evaluation, the RMSE of this method is reduced by 18.82%, 15.15%, 2.25%, and 1.10% on average compared with DnCNN, RED-CNN, CBDNet, and TSCNN, respectively. The average improvement of PSNR is 9.53%, 7.33%, 2.65%, and 3.69%, respectively.

CONCLUSIONS

This paper proposed a LDCT image denoising algorithm based on end-to-end training. The method in this paper can effectively improve the diagnostic performance of CT images by constraining the details of the images and restoring the LDCT image structure. The problem of increased noise and artifacts in CT images can be solved while maintaining the integrity of CT image tissue structure and pathological information. Compared with other algorithms, this method has better denoising effects both quantitatively and qualitatively.

A fidelity-embedded learning for metal artifact reduction in dental CBCT

PURPOSE

Dental cone-beam computed tomography (CBCT) has been increasingly used for dental and maxillofacial imaging. However, the presence of metallic inserts, such as implants, crowns, and dental braces, violates the CT model assumption, which leads to severe metal artifacts in the reconstructed CBCT image, resulting in the degradation of diagnostic performance. In this study, we used deep learning to reduce metal artifacts.

METHODS

The metal artifacts, appearing as streaks and shadows, are non-local and highly associated with various factors, including the geometry of metallic inserts, energy-dependent attenuation, and energy spectrum of the incident X-ray beam, making it difficult to learn their complicated structures directly. To provide a step-by-step environment in which deep learning can be trained, we propose an iterative learning approach in which the network at each iteration step learns the correction error caused by the previous network, while enforcing the data fidelity in the projection domain. To generate a realistic paired training dataset, metal-free CBCT scans were collected from patients without metallic inserts, and then simulated metal projection data were added to generate the corresponding metal-corrupted projection data.

RESULTS

The feasibility of the proposed method was investigated in clinical metal-affected CBCT scans, as well as simulated metal-affected CBCT scans. The results show that the proposed method significantly reduces metal artifacts while preserving the morphological structures near metallic objects and outperforms direct image domain learning.

CONCLUSION

The proposed fidelity-embedded learning can effectively reduce metal artifacts in dental CBCT compared with direct image domain learning. This article is protected by copyright. All rights reserved.

PWLS-PR: low-dose computed tomography image reconstruction using a patch-based regularization method based on the penalized weighted least squares total variation approach

BACKGROUND

Radiation exposure computed tomography (CT) scans and the associated risk of cancer in patients have been major clinical concerns. Existing research can achieve low-dose CT imaging by reducing the X-ray current and the number of projections per rotation of the human body. However, this method may produce excessive noise and fringe artifacts in the traditional filtered back projection (FBP)-reconstructed image.

METHODS

To solve this problem, iterative image reconstruction is a promising option to obtain high-quality images from low-dose scans. This paper proposes a patch-based regularization method based on penalized weighted least squares total variation (PWLS-PR) for iterative image reconstruction. This method uses neighborhood patches instead of single pixels to calculate the nonquadratic penalty. The proposed regularization method is more robust than the conventional regularization method in identifying random fluctuations caused by sharp edges and noise. Each iteration of the proposed algorithm can be described in the following three steps: image updating via the total variation based on penalized weighted least squares (PWLS-TV), image smoothing, and pixel-by-pixel image fusion.

RESULTS

Simulation and real-world projection experiments show that the proposed PWLS-PR algorithm achieves a higher image reconstruction performance than similar algorithms. Through the qualitative and quantitative evaluation of simulation experiments, the effectiveness of the method is also verified.

CONCLUSIONS

Furthermore, this study shows that the PWLS-PR method reduces the amount of projection data required for repeated CT scans and has the useful potential to reduce the radiation dose in clinical medical applications.

An iterative reconstruction method for sparse-projection data for low-dose CT

Reducing X-ray radiation is beneficial for reducing the risk of cancer in patients. There are two main approaches for achieving this goal namely, one is to reduce the X-ray current, and another is to apply sparse-view protocols to do image scanning and projections. However, these techniques usually lead to degradation of the reconstructed image quality, resulting in excessive noise and severe edge artifacts, which seriously affect the diagnosis result. In order to overcome such limitation, this study proposes and tests an algorithm based on guided kernel filtering. The algorithm combines the characteristics of anisotropic edges between adjacent image voxels, expresses the relevant weights with an exponential function, and adjusts the weights adaptively through local gray gradients to better preserve the image structure while suppressing noise information. Experiments show that the proposed method can effectively suppress noise and preserve the image structure. Comparing with similar algorithms, the proposed algorithm greatly improves the peak signal-to-noise ratio (PSNR), structural similarity (SSIM), and root mean square error (RMSE) of the reconstructed image. The proposed algorithm has the best effect in quantitative analysis, which verifies the effectiveness of the proposed method and good image reconstruction performance. Overall, this study demonstrates that the proposed method can reduce the number of projections required for repeated CT scans and has potential for medical applications in reducing radiation doses.

DeepOrganNet: On-the-Fly Reconstruction and Visualization of 3D / 4D Lung Models from Single-View Projections by Deep Deformation Network

This paper introduces a deep neural network based method, i.e., DeepOrganNet, to generate and visualize fully high-fidelity 3D / 4D organ geometric models from single-view medical images with complicated background in real time. Traditional 3D / 4D medical image reconstruction requires near hundreds of projections, which cost insufferable computational time and deliver undesirable high imaging / radiation dose to human subjects. Moreover, it always needs further notorious processes to segment or extract the accurate 3D organ models subsequently. The computational time and imaging dose can be reduced by decreasing the number of projections, but the reconstructed image quality is degraded accordingly. To our knowledge, there is no method directly and explicitly reconstructing multiple 3D organ meshes from a single 2D medical grayscale image on the fly. Given single-view 2D medical images, e.g., 3D / 4D-CT projections or X-ray images, our end-to-end DeepOrganNet framework can efficiently and effectively reconstruct 3D / 4D lung models with a variety of geometric shapes by learning the smooth deformation fields from multiple templates based on a trivariate tensor-product deformation technique, leveraging an informative latent descriptor extracted from input 2D images. The proposed method can guarantee to generate high-quality and high-fidelity manifold meshes for 3D / 4D lung models; while, all current deep learning based approaches on the shape reconstruction from a single image cannot. The major contributions of this work are to accurately reconstruct the 3D organ shapes from 2D single-view projection, significantly improve the procedure time to allow on-the-fly visualization, and dramatically reduce the imaging dose for human subjects. Experimental results are evaluated and compared with the traditional reconstruction method and the state-of-the-art in deep learning, by using extensive 3D and 4D examples, including both synthetic phantom and real patient datasets. The efficiency of the proposed method shows that it only needs several milliseconds to generate organ meshes with 10K vertices, which has great potential to be used in real-time image guided radiation therapy (IGRT).

Known-component 3D image reconstruction for improved intraoperative imaging in spine surgery: A clinical pilot study

PURPOSE

Intraoperative imaging plays an increased role in support of surgical guidance and quality assurance for interventional approaches. However, image quality sufficient to detect complications and provide quantitative assessment of the surgical product is often confounded by image noise and artifacts. In this work, we translated a three-dimensional model-based image reconstruction (referred to as "Known-Component Reconstruction," KC-Recon) for the first time to clinical studies with the aim of resolving both limitations.

METHODS

KC-Recon builds upon a penalized weighted least-squares (PWLS) method by incorporating models of surgical instrumentation ("known components") within a joint image registration-reconstruction process to improve image quality. Under IRB approval, a clinical pilot study was conducted with 17 spine surgery patients imaged under informed consent using the O-arm cone-beam CT system (Medtronic, Littleton MA) before and after spinal instrumentation. Volumetric images were generated for each patient using KC-Recon in comparison to conventional filtered backprojection (FBP). Imaging performance prior to instrumentation ("preinstrumentation") was evaluated in terms of soft-tissue contrast-to-noise ratio (CNR) and spatial resolution. The quality of images obtained after the instrumentation ("postinstrumentation") was assessed by quantifying the magnitude of metal artifacts (blooming and streaks) arising from pedicle screws. The potential low-dose advantages of the algorithm were tested by simulating low-dose data (down to one-tenth of the dose of standard protocols) from images acquired at normal dose.

RESULTS

Preinstrumentation images (at normal clinical dose and matched resolution) exhibited an average 24.0% increase in soft-tissue CNR with KC-Recon compared to FBP (N = 16, P = 0.02), improving visualization of paraspinal muscles, major vessels, and other soft-tissues about the spine and abdomen. For a total of 72 screws in postinstrumentation images, KC-Recon yielded a significant reduction in metal artifacts: 66.3% reduction in overestimation of screw shaft width due to blooming (P < 0.0001) and reduction in streaks at the screw tip (65.8% increase in attenuation accuracy, P < 0.0001), enabling clearer depiction of the screw within the pedicle and vertebral body for an assessment of breach. Depending on the imaging task, dose reduction up to an order of magnitude appeared feasible while maintaining soft-tissue visibility and metal artifact reduction.

CONCLUSIONS

KC-Recon offers a promising means to improve visualization in the presence of surgical instrumentation and reduce patient dose in image-guided procedures. The improved soft-tissue visibility could facilitate the use of cone-beam CT to soft-tissue surgeries, and the ability to precisely quantify and visualize instrument placement could provide a valuable check against complications in the operating room (cf., postoperative CT).

CT sinogram-consistency learning for metal-induced beam hardening correction

PURPOSE

This paper proposes a sinogram-consistency learning method to deal with beam hardening-related artifacts in polychromatic computerized tomography (CT). The presence of highly attenuating materials in the scan field causes an inconsistent sinogram that does not match the range space of the Radon transform. When the mismatched data are entered into the range space during CT reconstruction, streaking and shading artifacts are generated owing to the inherent nature of the inverse Radon transform METHODS: The proposed learning method aims to repair inconsistent sinogram by removing the primary metal-induced beam hardening factors along the metal trace in the sinogram. Taking account of the fundamental difficulty in obtaining sufficient training data in a medical environment, the learning method is designed to use simulated training data and a patient's implant type-specific learning model is used to simplify the learning process.

RESULTS

The feasibility of the proposed method is investigated using a dataset, consisting of real CT scans of pelvises containing simulated hip prostheses. The anatomical areas in training and test data are different, in order to demonstrate that the proposed method extracts the beam hardening features, selectively. The results show that our method successfully corrects sinogram inconsistency by extracting beam hardening sources by means of deep learning.

CONCLUSION

This paper proposed a deep learning method of sinogram correction for beam hardening reduction in CT for the first time. Conventional methods for beam hardening reduction are based on regularizations, and have the fundamental drawback of being not easily able to use manifold CT images, while a deep learning approach has the potential to do so.

Quantitative accuracy of CT numbers: Theoretical analyses and experimental studies

PURPOSE

The CT number accuracy, that is, CT number bias, plays an important role in clinical diagnosis. When strategies to reduce radiation dose are discussed, it is important to make sure that the CT number bias is controlled within an acceptable range. The purpose of this paper was to investigate the dependence of CT number bias on radiation dose level and on image contrast (i.e., the difference in CT number between the ROI and the background) in Computed Tomography (CT).

METHODS

A lesion-background model was introduced to theoretically study how the CT number bias changes with radiation exposure level and with CT number contrast when a simple linear reconstruction algorithm such as filtered backprojection (FBP) is used. The theoretical results were validated with experimental studies using a benchtop CT system equipped with a photon-counting detector (XC-HYDRA FX50, XCounter AB, Sweden) and a clinical diagnostic MDCT scanner (Discovery CT750 HD, GE Healthcare, Waukesha, WI, USA) equipped with an energy-integrating detector. The Catphan phantom (Catphan 600, the Phantom Laboratory, Salem, NY, USA) was scanned at different mAs levels and 50 scans were performed for each mAs. The bias of CT number was evaluated for each combination of mAs and ROIs with different contrast levels. An anthropomorphic phantom (ATOM 10-year-old phantom, Model 706, CIRS Inc. Norfolk, VA, USA) with much more heterogeneous object content was used to test the applicability of the theory to the more general image object cases.

RESULTS

Both theoretical and experimental studies showed that the CT number bias is inversely proportional to the radiation exposure level yet linearly dependent on the CT number contrast between the lesion and the background, that is, Bias ( μ ^ 1 FBP ) = α mAs ( 1 + β Δ H U ) .

CONCLUSIONS

The quantitative accuracy of CT numbers can be problematic and thus needs some extra attention when radiation dose is reduced. In this work, we showed that the bias of the FBP reconstruction increases as mAs is reduced; both positive and negative bias can be observed depending on the contrast difference between a targeted ROI and its surrounding background tissues.

Low-dose cone-beam CT via raw counts domain low-signal correction schemes: Performance assessment and task-based parameter optimization (Part II. Task-based parameter optimization)

PURPOSE

Different low-signal correction (LSC) methods have been shown to efficiently reduce noise streaks and noise level in CT to provide acceptable images at low-radiation dose levels. These methods usually result in CT images with highly shift-variant and anisotropic spatial resolution and noise, which makes the parameter optimization process highly nontrivial. The purpose of this work was to develop a local task-based parameter optimization framework for LSC methods.

METHODS

Two well-known LSC methods, the adaptive trimmed mean (ATM) filter and the anisotropic diffusion (AD) filter, were used as examples to demonstrate how to use the task-based framework to optimize filter parameter selection. Two parameters, denoted by the set P, for each LSC method were included in the optimization problem. For the ATM filter, these parameters are the low- and high-signal threshold levels p_l and p_h ; for the AD filter, the parameters are the exponents δ and γ in the brightness gradient function. The detectability index d' under the non-prewhitening (NPW) mathematical observer model was selected as the metric for parameter optimization. The optimization problem was formulated as an unconstrained optimization problem that consisted of maximizing an objective function d'(P), where i and j correspond to the i-th imaging task and j-th spatial location, respectively. Since there is no explicit mathematical function to describe the dependence of d' on the set of parameters P for each LSC method, the optimization problem was solved via an experimentally measured d' map over a densely sampled parameter space. In this work, three high-contrast-high-frequency discrimination imaging tasks were defined to explore the parameter space of each of the LSC methods: a vertical bar pattern (task I), a horizontal bar pattern (task II), and a multidirectional feature (task III). Two spatial locations were considered for the analysis, a posterior region-of-interest (ROI) located within the noise streaks region and an anterior ROI, located further from the noise streaks region. Optimal results derived from the task-based detectability index metric were compared to other operating points in the parameter space with different noise and spatial resolution trade-offs.

RESULTS

The optimal operating points determined through the d' metric depended on the interplay between the major spatial frequency components of each imaging task and the highly shift-variant and anisotropic noise and spatial resolution properties associated with each operating point in the LSC parameter space. This interplay influenced imaging performance the most when the major spatial frequency component of a given imaging task coincided with the direction of spatial resolution loss or with the dominant noise spatial frequency component; this was the case of imaging task II. The performance of imaging tasks I and III was influenced by this interplay in a smaller scale than imaging task II, since the major frequency component of task I was perpendicular to imaging task II, and because imaging task III did not have strong directional dependence. For both LSC methods, there was a strong dependence of the overall d' magnitude and shape of the contours on the spatial location within the phantom, particularly for imaging tasks II and III. The d' value obtained at the optimal operating point for each spatial location and imaging task was similar when comparing the LSC methods studied in this work.

CONCLUSIONS

A local task-based detectability framework to optimize the selection of parameters for LSC methods was developed. The framework takes into account the potential shift-variant and anisotropic spatial resolution and noise properties to maximize the imaging performance of the CT system. Optimal parameters for a given LSC method depend strongly on the spatial location within the image object.

Discriminative Feature Representation to Improve Projection Data Inconsistency for Low Dose CT Imaging

In low dose computed tomography (LDCT) imaging, the data inconsistency of measured noisy projections can significantly deteriorate reconstruction images. To deal with this problem, we propose here a new sinogram restoration approach, the sinogram- discriminative feature representation (S-DFR) method. Different from other sinogram restoration methods, the proposed method works through a 3-D representation-based feature decomposition of the projected attenuation component and the noise component using a well-designed composite dictionary containing atoms with discriminative features. This method can be easily implemented with good robustness in parameter setting. Its comparison to other competing methods through experiments on simulated and real data demonstrated that the S-DFR method offers a sound alternative in LDCT.

Low-dose 4D cone-beam CT via joint spatiotemporal regularization of tensor framelet and nonlocal total variation

Excessive radiation exposure is still a major concern in 4D cone-beam computed tomography (4D-CBCT) due to its prolonged scanning duration. Radiation dose can be effectively reduced by either under-sampling the x-ray projections or reducing the x-ray flux. However, 4D-CBCT reconstruction under such low-dose protocols is prone to image artifacts and noise. In this work, we propose a novel joint regularization-based iterative reconstruction method for low-dose 4D-CBCT. To tackle the under-sampling problem, we employ spatiotemporal tensor framelet (STF) regularization to take advantage of the spatiotemporal coherence of the patient anatomy in 4D images. To simultaneously suppress the image noise caused by photon starvation, we also incorporate spatiotemporal nonlocal total variation (SNTV) regularization to make use of the nonlocal self-recursiveness of anatomical structures in the spatial and temporal domains. Under the joint STF-SNTV regularization, the proposed iterative reconstruction approach is evaluated first using two digital phantoms and then using physical experiment data in the low-dose context of both under-sampled and noisy projections. Compared with existing approaches via either STF or SNTV regularization alone, the presented hybrid approach achieves improved image quality, and is particularly effective for the reconstruction of low-dose 4D-CBCT data that are not only sparse but noisy.

A new CT reconstruction technique using adaptive deformation recovery and intensity correction (ADRIC)

PURPOSE

Sequential same-patient CT images may involve deformation-induced and non-deformation-induced voxel intensity changes. An adaptive deformation recovery and intensity correction (ADRIC) technique was developed to improve the CT reconstruction accuracy, and to separate deformation from non-deformation-induced voxel intensity changes between sequential CT images.

MATERIALS AND METHODS

ADRIC views the new CT volume as a deformation of a prior high-quality CT volume, but with additional non-deformation-induced voxel intensity changes. ADRIC first applies the 2D-3D deformation technique to recover the deformation field between the prior CT volume and the new, to-be-reconstructed CT volume. Using the deformation-recovered new CT volume, ADRIC further corrects the non-deformation-induced voxel intensity changes with an updated algebraic reconstruction technique ("ART-dTV"). The resulting intensity-corrected new CT volume is subsequently fed back into the 2D-3D deformation process to further correct the residual deformation errors, which forms an iterative loop. By ADRIC, the deformation field and the non-deformation voxel intensity corrections are optimized separately and alternately to reconstruct the final CT. CT myocardial perfusion imaging scenarios were employed to evaluate the efficacy of ADRIC, using both simulated data of the extended-cardiac-torso (XCAT) digital phantom and experimentally acquired porcine data. The reconstruction accuracy of the ADRIC technique was compared to the technique using ART-dTV alone, and to the technique using 2D-3D deformation alone. The relative error metric and the universal quality index metric are calculated between the images for quantitative analysis. The relative error is defined as the square root of the sum of squared voxel intensity differences between the reconstructed volume and the "ground-truth" volume, normalized by the square root of the sum of squared "ground-truth" voxel intensities. In addition to the XCAT and porcine studies, a physical lung phantom measurement study was also conducted. Water-filled balloons with various shapes/volumes and concentrations of iodinated contrasts were put inside the phantom to simulate both deformations and non-deformation-induced intensity changes for ADRIC reconstruction. The ADRIC-solved deformations and intensity changes from limited-view projections were compared to those of the "gold-standard" volumes reconstructed from fully sampled projections.

RESULTS

For the XCAT simulation study, the relative errors of the reconstructed CT volume by the 2D-3D deformation technique, the ART-dTV technique, and the ADRIC technique were 14.64%, 19.21%, and 11.90% respectively, by using 20 projections for reconstruction. Using 60 projections for reconstruction reduced the relative errors to 12.33%, 11.04%, and 7.92% for the three techniques, respectively. For the porcine study, the corresponding results were 13.61%, 8.78%, and 6.80% by using 20 projections; and 12.14%, 6.91%, and 5.29% by using 60 projections. The ADRIC technique also demonstrated robustness to varying projection exposure levels. For the physical phantom study, the average DICE coefficient between the initial prior balloon volume and the new "gold-standard" balloon volumes was 0.460. ADRIC reconstruction by 21 projections increased the average DICE coefficient to 0.954.

CONCLUSION

The ADRIC technique outperformed both the 2D-3D deformation technique and the ART-dTV technique in reconstruction accuracy. The alternately solved deformation field and non-deformation voxel intensity corrections can benefit multiple clinical applications, including tumor tracking, radiotherapy dose accumulation, and treatment outcome analysis.

Iterative reconstruction for sparse-view x-ray CT using alpha-divergence constrained total generalized variation minimization

BCKGROUND

Accurate statistical model of the measured projection data is essential for computed tomography (CT) image reconstruction. The transmission data can be described by a compound Poisson distribution upon an electronic noise background. However, such a statistical distribution is numerically intractable for image reconstruction.

OBJECTIVE

Although the sinogram data is easily manipulated, it lacks a statistical description for image reconstruction. To address this problem, we present an alpha-divergence constrained total generalized variation (AD-TGV) method for sparse-view x-ray CT image reconstruction.

METHODS

The AD-TGV method is formulated as an optimization problem, which balances the alpha-divergence (AD) fidelity and total generalized variation (TGV) regularization in one framework. The alpha-divergence is used to measure the discrepancy between the measured and estimated projection data. The TGV regularization can effectively eliminate the staircase and patchy artifacts which is often observed in total variation (TV) regularization. A modified proximal forward-backward splitting algorithm was proposed to minimize the associated objective function.

RESULTS

Qualitative and quantitative evaluations were carried out on both phantom and patient data. Compared with the original TV-based method, the evaluations clearly demonstrate that the AD-TGV method achieves higher accuracy and lower noise, while preserving structural details.

CONCLUSIONS

The experimental results show that the presented AD-TGV method can achieve more gains over the AD-TV method in preserving structural details and suppressing image noise and undesired patchy artifacts. The authors can draw the conclusion that the presented AD-TGV method is potential for radiation dose reduction by lowering the milliampere-seconds (mAs) and/or reducing the number of projection views.

Comparison Between Pre-Log and Post-Log Statistical Models in Ultra-Low-Dose CT Reconstruction

X-ray detectors in clinical computed tomography (CT) usually operate in current-integrating mode. Their complicated signal statistics often lead to intractable likelihood functions for practical use in model-based image reconstruction (MBIR). It is therefore desirable to design simplified statistical models without losing the essential factors. Depending on whether the CT transmission data are logarithmically transformed, pre-log and post-log models are two major categories of choices in CT MBIR. Both being approximations, it remains an open question whether one model can notably improve image quality over the other on real scanners. In this study, we develop and compare several pre-log and post-log MBIR algorithms under a unified framework. Their reconstruction accuracy based on simulation and clinical datasets are evaluated. The results show that pre-log MBIR can achieve notably better quantitative accuracy than post-log MBIR in ultra-low-dose CT, although in less extreme cases, post-log MBIR with handcrafted pre-processing remains a competitive alternative. Pre-log MBIR could play a growing role in emerging ultra-low-dose CT applications.

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

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

Cerebral perfusion computed tomography deconvolution via structure tensor total variation regularization

PURPOSE

Cerebral perfusion computed tomography (PCT) imaging as an accurate and fast acute ischemic stroke examination has been widely used in clinic. Meanwhile, a major drawback of PCT imaging is the high radiation dose due to its dynamic scan protocol. The purpose of this work is to develop a robust perfusion deconvolution approach via structure tensor total variation (STV) regularization (PD-STV) for estimating an accurate residue function in PCT imaging with the low-milliampere-seconds (low-mAs) data acquisition.

METHODS

Besides modeling the spatio-temporal structure information of PCT data, the STV regularization of the present PD-STV approach can utilize the higher order derivatives of the residue function to enhance denoising performance. To minimize the objective function, the authors propose an effective iterative algorithm with a shrinkage/thresholding scheme. A simulation study on a digital brain perfusion phantom and a clinical study on an old infarction patient were conducted to validate and evaluate the performance of the present PD-STV approach.

RESULTS

In the digital phantom study, visual inspection and quantitative metrics (i.e., the normalized mean square error, the peak signal-to-noise ratio, and the universal quality index) assessments demonstrated that the PD-STV approach outperformed other existing approaches in terms of the performance of noise-induced artifacts reduction and accurate perfusion hemodynamic maps (PHM) estimation. In the patient data study, the present PD-STV approach could yield accurate PHM estimation with several noticeable gains over other existing approaches in terms of visual inspection and correlation analysis.

CONCLUSIONS

This study demonstrated the feasibility and efficacy of the present PD-STV approach in utilizing STV regularization to improve the accuracy of residue function estimation of cerebral PCT imaging in the case of low-mAs.

Sinogram restoration in computed tomography with an edge-preserving penalty

PURPOSE

With the goal of producing a less computationally intensive alternative to fully iterative penalized-likelihood image reconstruction, our group has explored the use of penalized-likelihood sinogram restoration for transmission tomography. Previously, we have exclusively used a quadratic penalty in our restoration objective function. However, a quadratic penalty does not excel at preserving edges while reducing noise. Here, we derive a restoration update equation for nonquadratic penalties. Additionally, we perform a feasibility study to extend our sinogram restoration method to a helical cone-beam geometry and clinical data.

METHODS

A restoration update equation for nonquadratic penalties is derived using separable parabolic surrogates (SPS). A method for calculating sinogram degradation coefficients for a helical cone-beam geometry is proposed. Using simulated data, sinogram restorations are performed using both a quadratic penalty and the edge-preserving Huber penalty. After sinogram restoration, Fourier-based analytical methods are used to obtain reconstructions, and resolution-noise trade-offs are investigated. For the fan-beam geometry, a comparison is made to image-domain SPS reconstruction using the Huber penalty. The effects of varying object size and contrast are also investigated. For the helical cone-beam geometry, we investigate the effect of helical pitch (axial movement/rotation). Huber-penalty sinogram restoration is performed on 3D clinical data, and the reconstructed images are compared to those generated with no restoration.

RESULTS

We find that by applying the edge-preserving Huber penalty to our sinogram restoration methods, the reconstructed image has a better resolution-noise relationship than an image produced using a quadratic penalty in the sinogram restoration. However, we find that this relatively straightforward approach to edge preservation in the sinogram domain is affected by the physical size of imaged objects in addition to the contrast across the edge. This presents some disadvantages of this method relative to image-domain edge-preserving methods, although the computational burden of the sinogram-domain approach is much lower. For a helical cone-beam geometry, we found applying sinogram restoration in 3D was reasonable and that pitch did not make a significant difference in the general effect of sinogram restoration. The application of Huber-penalty sinogram restoration to clinical data resulted in a reconstruction with less noise while retaining resolution.

CONCLUSIONS

Sinogram restoration with the Huber penalty is able to provide better resolution-noise performance than restoration with a quadratic penalty. Additionally, sinogram restoration with the Huber penalty is feasible for helical cone-beam CT and can be applied to clinical data.

3D-2D Deformable Image Registration Using Feature-Based Nonuniform Meshes

By using prior information of planning CT images and feature-based nonuniform meshes, this paper demonstrates that volumetric images can be efficiently registered with a very small portion of 2D projection images of a Cone-Beam Computed Tomography (CBCT) scan. After a density field is computed based on the extracted feature edges from planning CT images, nonuniform tetrahedral meshes will be automatically generated to better characterize the image features according to the density field; that is, finer meshes are generated for features. The displacement vector fields (DVFs) are specified at the mesh vertices to drive the deformation of original CT images. Digitally reconstructed radiographs (DRRs) of the deformed anatomy are generated and compared with corresponding 2D projections. DVFs are optimized to minimize the objective function including differences between DRRs and projections and the regularity. To further accelerate the above 3D-2D registration, a procedure to obtain good initial deformations by deforming the volume surface to match 2D body boundary on projections has been developed. This complete method is evaluated quantitatively by using several digital phantoms and data from head and neck cancer patients. The feature-based nonuniform meshing method leads to better results than either uniform orthogonal grid or uniform tetrahedral meshes.

Iterative image reconstruction for ultra-low-dose CT with a combined low-mAs and sparse-view protocol

Ultra-low-dose x-ray computed tomography (CT) imaging is needed in CT fields. Through a scan protocol by lowering the milliampere-seconds (mAs) and reducing the number of projections per rotation around the body, we can realize low-dose CT imaging. However, the resulting noisy and insufficient measurements will unavoidably cause the degradation of desired-image. To solve this problem, iterative image reconstruction is a promising choice for achieving high-quality image with a low-dose scan. In this study, we are focusing on ultra-low-dose CT image reconstruction by using penalized weighted least-square (PWLS) criteria with a combined low-mAs and sparse-view protocol. Specifically, the sinogram data acquired with a combined low-mAs and sparse-view protocol is first restored by using a PWLS based sinogram restoration method. Then, the restored sinogram data is hereafter used to reconstruct image by using a PWLS based total variation (PWLS-TV) method. Qualitative and quantitative evaluations by simulations were carried out to validate the present method.

Statistical image reconstruction for low-dose CT using nonlocal means-based regularization. Part II: An adaptive approach

To reduce radiation dose in X-ray computed tomography (CT) imaging, one common strategy is to lower the tube current and exposure time settings during projection data acquisition. However, this strategy would inevitably increase the projection data noise, and the resulting image by the conventional filtered back-projection (FBP) method may suffer from excessive noise and streak artifacts. The well-known edge-preserving nonlocal means (NLM) filtering can reduce the noise-induced artifacts in the FBP reconstructed image, but it sometimes cannot completely eliminate the artifacts, especially under the very low-dose circumstance when the image is severely degraded. Instead of taking NLM filtering, we proposed a NLM-regularized statistical image reconstruction scheme, which can effectively suppress the noise-induced artifacts and significantly improve the reconstructed image quality. From our previous investigation on NLM-based strategy, we noted that using a spatially invariant filtering parameter in the regularization was rarely optimal for the entire field of view (FOV). Therefore, in this study we developed a novel strategy for designing spatially variant filtering parameters which are adaptive to the local characteristics of the image to be reconstructed. This adaptive NLM-regularized statistical image reconstruction method was evaluated with low-contrast phantoms and clinical patient data to show (1) the necessity in introducing the spatial adaptivity and (2) the efficacy of the adaptivity in achieving superiority in reconstructing CT images from low-dose acquisitions.

Deriving adaptive MRF coefficients from previous normal-dose CT scan for low-dose image reconstruction via penalized weighted least-squares minimization

PURPOSE

Repeated computed tomography (CT) scans are required for some clinical applications such as image-guided interventions. To optimize radiation dose utility, a normal-dose scan is often first performed to set up reference, followed by a series of low-dose scans for intervention. One common strategy to achieve the low-dose scan is to lower the x-ray tube current and exposure time (mAs) or tube voltage (kVp) setting in the scanning protocol, but the resulted image quality by the conventional filtered back-projection (FBP) method may be severely degraded due to the excessive noise. Penalized weighted least-squares (PWLS) image reconstruction has shown the potential to significantly improve the image quality from low-mAs acquisitions, where the penalty plays an important role. In this work, the authors' explore an adaptive Markov random field (MRF)-based penalty term by utilizing previous normal-dose scan to improve the subsequent low-dose scans image reconstruction.

METHODS

In this work, the authors employ the widely-used quadratic-form MRF as the penalty model and explore a novel idea of using the previous normal-dose scan to obtain the MRF coefficients for adaptive reconstruction of the low-dose images. In the coefficients determination, the authors further explore another novel idea of using the normal-dose scan to obtain a scale map, which describes an optimal neighborhood for the coefficients determination such that a local uniform region has a small spread of frequency spectrum and, therefore, a small MRF window, and vice versa. The proposed penalty term is incorporated into the PWLS image reconstruction framework, and the low-dose images are reconstructed via the PWLS minimization.

RESULTS

The presented adaptive MRF based PWLS algorithm was validated by physical phantom and patient data. The experimental results demonstrated that the presented algorithm is superior to the PWLS reconstruction using the conventional Gaussian MRF penalty or the edge-preserving Huber penalty and the conventional FBP method, in terms of image noise reduction and edge/detail/contrast preservation.

CONCLUSIONS

This study demonstrated the feasibility and efficacy of the proposed scheme in utilizing previous normal-dose CT scan to improve the subsequent low-dose scans.

Statistical image reconstruction for low-dose CT using nonlocal means-based regularization

Low-dose computed tomography (CT) imaging without sacrifice of clinical tasks is desirable due to the growing concerns about excessive radiation exposure to the patients. One common strategy to achieve low-dose CT imaging is to lower the milliampere-second (mAs) setting in data scanning protocol. However, the reconstructed CT images by the conventional filtered back-projection (FBP) method from the low-mAs acquisitions may be severely degraded due to the excessive noise. Statistical image reconstruction (SIR) methods have shown potentials to significantly improve the reconstructed image quality from the low-mAs acquisitions, wherein the regularization plays a critical role and an established family of regularizations is based on the Markov random field (MRF) model. Inspired by the success of nonlocal means (NLM) in image processing applications, in this work, we propose to explore the NLM-based regularization for SIR to reconstruct low-dose CT images from low-mAs acquisitions. Experimental results with both digital and physical phantoms consistently demonstrated that SIR with the NLM-based regularization can achieve more gains than SIR with the well-known Gaussian MRF regularization or the generalized Gaussian MRF regularization and the conventional FBP method, in terms of image noise reduction and resolution preservation.

Deformable image registration with local rigidity constraints for cone-beam CT-guided spine surgery

Image-guided spine surgery (IGSS) is associated with reduced co-morbidity and improved surgical outcome. However, precise localization of target anatomy and adjacent nerves and vessels relative to planning information (e.g., device trajectories) can be challenged by anatomical deformation. Rigid registration alone fails to account for deformation associated with changes in spine curvature, and conventional deformable registration fails to account for rigidity of the vertebrae, causing unrealistic distortions in the registered image that can confound high-precision surgery. We developed and evaluated a deformable registration method capable of preserving rigidity of bones while resolving the deformation of surrounding soft tissue. The method aligns preoperative CT to intraoperative cone-beam CT (CBCT) using free-form deformation (FFD) with constraints on rigid body motion imposed according to a simple intensity threshold of bone intensities. The constraints enforced three properties of a rigid transformation-namely, constraints on affinity (AC), orthogonality (OC), and properness (PC). The method also incorporated an injectivity constraint (IC) to preserve topology. Physical experiments involving phantoms, an ovine spine, and a human cadaver as well as digital simulations were performed to evaluate the sensitivity to registration parameters, preservation of rigid body morphology, and overall registration accuracy of constrained FFD in comparison to conventional unconstrained FFD (uFFD) and Demons registration. FFD with orthogonality and injectivity constraints (denoted FFD+OC+IC) demonstrated improved performance compared to uFFD and Demons. Affinity and properness constraints offered little or no additional improvement. The FFD+OC+IC method preserved rigid body morphology at near-ideal values of zero dilatation (D = 0.05, compared to 0.39 and 0.56 for uFFD and Demons, respectively) and shear (S = 0.08, compared to 0.36 and 0.44 for uFFD and Demons, respectively). Target registration error (TRE) was similarly improved for FFD+OC+IC (0.7 mm), compared to 1.4 and 1.8 mm for uFFD and Demons. Results were validated in human cadaver studies using CT and CBCT images, with FFD+OC+IC providing excellent preservation of rigid morphology and equivalent or improved TRE. The approach therefore overcomes distortions intrinsic to uFFD and could better facilitate high-precision IGSS.

Soft-tissue imaging with C-arm cone-beam CT using statistical reconstruction

The potential for statistical image reconstruction methods such as penalized-likelihood (PL) to improve C-arm cone-beam CT (CBCT) soft-tissue visualization for intraoperative imaging over conventional filtered backprojection (FBP) is assessed in this work by making a fair comparison in relation to soft-tissue performance. A prototype mobile C-arm was used to scan anthropomorphic head and abdomen phantoms as well as a cadaveric torso at doses substantially lower than typical values in diagnostic CT, and the effects of dose reduction via tube current reduction and sparse sampling were also compared. Matched spatial resolution between PL and FBP was determined by the edge spread function of low-contrast (∼ 40-80 HU) spheres in the phantoms, which were representative of soft-tissue imaging tasks. PL using the non-quadratic Huber penalty was found to substantially reduce noise relative to FBP, especially at lower spatial resolution where PL provides a contrast-to-noise ratio increase up to 1.4-2.2 × over FBP at 50% dose reduction across all objects. Comparison of sampling strategies indicates that soft-tissue imaging benefits from fully sampled acquisitions at dose above ∼ 1.7 mGy and benefits from 50% sparsity at dose below ∼ 1.0 mGy. Therefore, an appropriate sampling strategy along with the improved low-contrast visualization offered by statistical reconstruction demonstrates the potential for extending intraoperative C-arm CBCT to applications in soft-tissue interventions in neurosurgery as well as thoracic and abdominal surgeries by overcoming conventional tradeoffs in noise, spatial resolution, and dose.

Iterative image reconstruction for sparse-view CT using normal-dose image induced total variation prior

X-ray computed tomography (CT) iterative image reconstruction from sparse-view projection data has been an important research topic for radiation reduction in clinic. In this paper, to relieve the requirement of misalignment reduction operation of the prior image constrained compressed sensing (PICCS) approach introduced by Chen et al, we present an iterative image reconstruction approach for sparse-view CT using a normal-dose image induced total variation (ndiTV) prior. The associative objective function of the present approach is constructed under the penalized weighed least-square (PWLS) criteria, which contains two terms, i.e., the weighted least-square (WLS) fidelity and the ndiTV prior, and is referred to as "PWLS-ndiTV". Specifically, the WLS fidelity term is built based on an accurate relationship between the variance and mean of projection data in the presence of electronic background noise. The ndiTV prior term is designed to reduce the influence of the misalignment between the desired- and prior- image by using a normal-dose image induced non-local means (ndiNLM) filter. Subsequently, a modified steepest descent algorithm is adopted to minimize the associative objective function. Experimental results on two different digital phantoms and an anthropomorphic torso phantom show that the present PWLS-ndiTV approach for sparse-view CT image reconstruction can achieve noticeable gains over the existing similar approaches in terms of noise reduction, resolution-noise tradeoff, and low-contrast object detection.

SR-NLM: a sinogram restoration induced non-local means image filtering for low-dose computed tomography

Radiation dose has raised significant concerns to patients and operators in modern X-ray computed tomography (CT) examinations. A simple and cost-effective means to perform a low-dose CT scan is to lower the milliampere-seconds (mAs) as low as reasonably achievable in data acquisition. However, the associated image quality with lower-mAs scans (or low-dose scans) will be unavoidably degraded due to the excessive data noise, if no adequate noise control is applied during image reconstruction. For image reconstruction with low-dose scans, sinogram restoration algorithms based on modeling the noise properties of measurement can produce an image with noise-induced artifact suppression, but they often suffer noticeable resolution loss. As an alternative technique, the noise-reduction algorithms via edge-preserving image filtering can yield an image without noticeable resolution loss, but they often do not completely eliminate the noise-induced artifacts. With above observations, in this paper, we present a sinogram restoration induced non-local means (SR-NLM) image filtering algorithm to retain the CT image quality by fully considering the advantages of the sinogram restoration and image filtering algorithms in low-dose image reconstruction. Extensive experimental results show that the present SR-NLM algorithm outperforms the existing methods in terms of cross profile, noise reduction, contrast-to-ratio measure, noise-resolution tradeoff and receiver operating characteristic (ROC) curves.

Iterative image reconstruction for cerebral perfusion CT using a pre-contrast scan induced edge-preserving prior

Cerebral perfusion x-ray computed tomography (PCT) imaging, which detects and characterizes the ischemic penumbra, and assesses blood-brain barrier permeability with acute stroke or chronic cerebrovascular diseases, has been developed extensively over the past decades. However, due to its sequential scan protocol, the associated radiation dose has raised significant concerns to patients. Therefore, in this study we developed an iterative image reconstruction algorithm based on the maximum a posterior (MAP) principle to yield a clinically acceptable cerebral PCT image with lower milliampere-seconds (mA s). To preserve the edges of the reconstructed image, an edge-preserving prior was designed using a normal-dose pre-contrast unenhanced scan. For simplicity, the present algorithm was termed as 'MAP-ndiNLM'. Evaluations with the digital phantom and the simulated low-dose clinical brain PCT datasets clearly demonstrate that the MAP-ndiNLM method can achieve more significant gains than the existing FBP and MAP-Huber algorithms with better image noise reduction, low-contrast object detection and resolution preservation. More importantly, the MAP-ndiNLM method can yield more accurate kinetic enhanced details and diagnostic hemodynamic parameter maps than the MAP-Huber method.

Variance analysis of x-ray CT sinograms in the presence of electronic noise background

PURPOSE

Low-dose x-ray computed tomography (CT) is clinically desired. Accurate noise modeling is a fundamental issue for low-dose CT image reconstruction via statistics-based sinogram restoration or statistical iterative image reconstruction. In this paper, the authors analyzed the statistical moments of low-dose CT data in the presence of electronic noise background.

METHODS

The authors first studied the statistical moment properties of detected signals in CT transmission domain, where the noise of detected signals is considered as quanta fluctuation upon electronic noise background. Then the authors derived, via the Taylor expansion, a new formula for the mean-variance relationship of the detected signals in CT sinogram domain, wherein the image formation becomes a linear operation between the sinogram data and the unknown image, rather than a nonlinear operation in the CT transmission domain. To get insight into the derived new formula by experiments, an anthropomorphic torso phantom was scanned repeatedly by a commercial CT scanner at five different mAs levels from 100 down to 17.

RESULTS

The results demonstrated that the electronic noise background is significant when low-mAs (or low-dose) scan is performed.

CONCLUSIONS

The influence of the electronic noise background should be considered in low-dose CT imaging.

Noise reduction with low dose CT data based on a modified ROF model

In order to reduce the radiation exposure caused by Computed Tomography (CT) scanning, low dose CT has gained much interest in research as well as in industry. One fundamental difficulty for low dose CT lies in its heavy noise pollution in the raw data which leads to quality deterioration for reconstructed images. In this paper, we propose a modified ROF model to denoise low dose CT measurement data in light of Poisson noise model. Experimental results indicate that the reconstructed CT images based on measurement data processed by our model are in better quality, compared to the original ROF model or bilateral filtering.

Comparison of sinogram- and image-domain penalized-likelihood image reconstruction estimators

PURPOSE

In recent years, the authors and others have been exploring the use of penalized-likelihood sinogram-domain smoothing and restoration approaches for emission and transmission tomography. The motivation for this strategy was initially pragmatic: to provide a more computationally feasible alternative to fully iterative penalized-likelihood image reconstruction involving expensive backprojections and reprojections, while still obtaining some of the benefits of the statistical modeling employed in penalized-likelihood approaches. In this work, the authors seek to compare the two approaches in greater detail.

METHODS

The sinogram-domain strategy entails estimating the "ideal" line integrals needed for reconstruction of an activity or attenuation distribution from the set of noisy, potentially degraded tomographic measurements by maximizing a penalized-likelihood objective function. The objective function models the data statistics as well as any degradation that can be represented in the sinogram domain. The estimated line integrals can then be input to analytic reconstruction algorithms such as filtered backprojection (FBP). The authors compare this to fully iterative approaches maximizing similar objective functions.

RESULTS

The authors present mathematical analyses based on so-called equivalent optimization problems that establish that the approaches can be made precisely equivalent under certain restrictive conditions. More significantly, by use of resolution-variance tradeoff studies, the authors show that they can yield very similar performance under more relaxed, realistic conditions.

CONCLUSIONS

The sinogram- and image-domain approaches are equivalent under certain restrictive conditions and can perform very similarly under more relaxed conditions. The match is particularly good for fully sampled, high-resolution CT geometries. One limitation of the sinogram-domain approach relative to the image-domain approach is the difficulty of imposing additional constraints, such as image non-negativity.

Computer-assisted scan protocol and reconstruction (CASPAR)-reduction of image noise and patient dose

X-ray computed tomography is a powerful medical imaging device. It allows high-resolution 3-D visualization of the human body. However, one drawback is the health risk associated with ionizing radiation. Simply downscaling the radiation intensities over the entire scan results in increased quantum noise. This paper proposes the concept of computer-assisted scan protocol and reconstruction. More specifically, we propose a method to compute patient and task-specific intensity profiles that achieve an optimal tradeoff between radiation dose and image quality. Therefore, reasonable image variance and dose metrics are derived. Conventional third-generation systems as well as inverted geometry concepts are considered. Two dose/noise minimization problems are formulated and solved by an efficient algorithm providing optimized milliampere (mA)-profiles. Thorax phantom simulations demonstrate the promising advantage of this technique: in this particular example, the dose is reduced by 53% for third-generation systems and by 86% for an inverted geometry in comparison to a sinusoidal mA-profile at a constant upper noise limit.

Electronic noise modeling in statistical iterative reconstruction

We consider electronic noise modeling in tomographic image reconstruction when the measured signal is the sum of a Gaussian distributed electronic noise component and another random variable whose log-likelihood function satisfies a certain linearity condition. Examples of such likelihood functions include the Poisson distribution and an exponential dispersion (ED) model that can approximate the signal statistics in integration mode X-ray detectors. We formulate the image reconstruction problem as a maximum-likelihood estimation problem. Using an expectation-maximization approach, we demonstrate that a reconstruction algorithm can be obtained following a simple substitution rule from the one previously derived without electronic noise considerations. To illustrate the applicability of the substitution rule, we present examples of a fully iterative reconstruction algorithm and a sinogram smoothing algorithm both in transmission CT reconstruction when the measured signal contains additive electronic noise. Our simulation studies show the potential usefulness of accurate electronic noise modeling in low-dose CT applications.

Dose reduction for kilovotage cone-beam computed tomography in radiation therapy

Kilovotage cone-beam computed tomography (kV-CBCT) has shown potentials to improve the accuracy of a patient setup in radiotherapy. However, daily and repeated use of CBCT will deliver high extra radiation doses to patients. One way to reduce the patient dose is to lower mAs when acquiring projection data. This, however, degrades the quality of low mAs CBCT images dramatically due to excessive noises. In this work, we aim to improve the CBCT image quality from low mAs scans. Based on the measured noise properties of the sinogram, a penalized weighted least-squares (PWLS) objective function was constructed, and the ideal sinogram was then estimated by minimizing the PWLS objection function. To preserve edge information in the projection data, an anisotropic penalty term was designed using the intensity difference between neighboring pixels. The effectiveness of the presented algorithm was demonstrated by two experimental phantom studies. Noise in the reconstructed CBCT image acquired with a low mAs protocol was greatly suppressed after the proposed sinogram domain image processing, without noticeable sacrifice of the spatial resolution.

Multiscale penalized weighted least-squares sinogram restoration for low-dose x-ray computed tomography

In this paper, we propose a novel multiscale penalized weighted least-squares (PWLS) method for restoration of low-dose computed tomography (CT) sinogram. The method utilizes wavelet transform for the multiscale or multiresolution analysis on the sinogram. Specifically, the Mallat-Zhong's wavelet transform is applied to decompose the sinogram to different resolution levels. At each decomposed resolution level, a PWLS criterion is applied to restore the noise-contaminated wavelet coefficients, where the penalty is adaptive to each resolution scale and the weight is updated by an exponential relationship between the data variance and mean at each scale and location. The proposed PWLS method is based on the observations that 1) noise in the CT sinogram after logarithm transform and calibration can be modeled as signal-dependent variables and the sample variance depends on the sample mean by an exponential relationship; and 2) noise reduction can be more effective when it is adaptive to different resolution levels. The effectiveness of the proposed multiscale PWLS method is validated by both computer simulations and experimental studies. The gain by multiscale approach over single scale means is quantified by noise-resolution tradeoff measures.

Multiscale penalized weighted least-squares sinogram restoration for low-dose X-ray computed tomography

We propose a novel multiscale penalized weighted least-squares (PWLS) method for restoration of low-dose computed tomography (CT) sinogram. The method utilizes wavelet transform for the multiscale or multi-resolution analysis on the sinogram. Specifically the Mallat-Zhong's wavelet transform is applied to decompose the sinogram to different resolution levels. At each decomposed resolution level, a PWLS criterion is applied to restore the noise-contaminated wavelet coefficients, where the penalty is adaptive to each resolution scale and the weight is adaptive to each scale and each location. The proposed PWLS method is based on the observation that (1) the noisy sinogram of low-dose CT after logarithm transform can be modeled as signal-dependent Gaussian variables and the sample variance depends on the sample mean; and (2) the noise restoration can be more effective when it is adaptive to different resolution levels. The effectiveness of the proposed multiscale PWLS method is validated by an experimental study. The gain by multiscale approach over single-scale means is quantified by noise-resolution tradeoff measures.

Sparseness prior based iterative image reconstruction for retrospectively gated cardiac micro-CT

Recent advances in murine cardiac studies with three-dimensional (3D) cone beam micro-CT used a retrospective gating technique. However, this sampling technique results in a limited number of projections with an irregular angular distribution due to the temporal resolution requirements and radiation dose restrictions. Both angular irregularity and undersampling complicate the reconstruction process, since they cause significant streaking artifacts. This work provides an iterative reconstruction solution to address this particular challenge. A sparseness prior regularized weighted l2 norm optimization is proposed to mitigate streaking artifacts based on the fact that most medical images are compressible. Total variation is implemented in this work as the regularizer for its simplicity. Comparison studies are conducted on a 3D cardiac mouse phantom generated with experimental data. After optimization, the method is applied to in vivo cardiac micro-CT data.

Noise Reduction for Low-Dose Single-Slice Helical CT Sinograms

Helical computed tomography (HCT) has several advantages over conventional step-and-shoot CT for imaging a relatively large object, especially for dynamic studies. However, HCT may increase X-ray exposure significantly. This work aims to reduce the radiation by lowering X-ray tube current (mA) and filtering low-mA (or dose) sinogram noise of HCT. The noise reduction method is based on three observations on HCT: (1) the axial sampling of HCT projections is nearly continuous as detection system rotates; (2) the noise distribution in sinogram space is nearly a Gaussian after system calibration (including logarithmic transform); and (3) the relationship between the calibrated data mean and variance can be expressed as an exponential functional across the field-of-view. Based on the second and third observations, a penalized weighted least-squares (PWLS) solution is an optimal choice, where the weight is given by the mean-variance relationship. The first observation encourages the use of Karhunen-Loève (KL) transform along the axial direction because of the associated correlation. In the KL domain, the eigenvalue of each principal component and the derived data variance provide the signal-to-noise ratio (SNR) information, resulting in a SNR-adaptive noise reduction. The KL-PWLS noise-reduction method was implemented analytically for efficient restoration of large volume HCT sinograms. Simulation studies showed a noticeable improvement, in terms of image quality and defect detectability, of the proposed noise-reduction method over the Ordered-Subsets Expectation-Maximization reconstruction and the conventional low-pass noise filtering with optimal cutoff frequency and/or other filter parameters.

Local noise weighted filtering for emphysema scoring of low-dose CT images

Computed tomography (CT) has become the new reference standard for quantification of emphysema. The most popular measure of emphysema derived from CT is the pixel index (PI), which expresses the fraction of the lung volume with abnormally low intensity values. As PI is calculated from a single, fixed threshold on intensity, this measure is strongly influenced by noise. This effect shows up clearly when comparing the PI score of a high-dose scan to the PI score of a low-dose (i.e., noisy) scan of the same subject. In this paper, the noise variance (NOVA) filter is presented: a general framework for (iterative) nonlinear filtering, which uses an estimate of the spatially dependent noise variance in an image. The NOVA filter iteratively estimates the local image noise and filters the image. For the specific purpose of emphysema quantification of low-dose CT images, a dedicated, noniterative NOVA filter is constructed by using prior knowledge of the data to obtain a good estimate of the spatially dependent noise in an image. The performance of the NOVA filter is assessed by comparing characteristics of pairs of high-dose and low-dose scans. The compared characteristics are the PI scores for different thresholds and the size distributions of emphysema bullae. After filtering, the PI scores of high-dose and low-dose images agree to within 2%-3% points. The reproducibility of the high-dose bullae size distribution is also strongly improved. NOVA filtering of a CT image of typically 400 x 512 x 512 voxels takes only a couple of minutes which makes it suitable for routine use in clinical practice.

Comparison of quadratic- and median-based roughness penalties for penalized-likelihood sinogram restoration in computed tomography

We have compared the performance of two different penalty choices for a penalized-likelihood sinogram-restoration strategy we have been developing. One is a quadratic penalty we have employed previously and the other is a new median-based penalty. We compared the approaches to a noniterative adaptive filter that loosely but not explicitly models data statistics. We found that the two approaches produced similar resolution-variance tradeoffs to each other and that they outperformed the adaptive filter in the low-dose regime, which suggests that the particular choice of penalty in our approach may be less important than the fact that we are explicitly modeling data statistics at all. Since the quadratic penalty allows for derivation of an algorithm that is guaranteed to monotonically increase the penalized-likelihood objective function, we find it to be preferable to the median-based penalty.

Analytical reconstruction of deconvolved Fourier rebinned PET sinograms

Fully 3D PET data are often rebinned into 2D data sets in order to avoid computationally intensive fully 3D reconstruction. Then, conventional 2D reconstruction techniques are employed to obtain images from the rebinned data. In a common scenario, 2D filtered back projection (FBP) is applied to Fourier rebinned (FORE) data. This approach is suboptimal because FBP is based on an idealized mathematical model of the data and cannot account for the statistical structure of data and noise. FORE data contain some blur in all three dimensions in comparison to conventional 2D PET data. In this work, we propose methods for approximating this blur in the sinogram domain due to FORE through its point spread function (PSF). We also explore simple methods for deconvolving the rebinned data with this PSF to restore it to a more ideal state prior to FBP. Our results show that deconvolution of the approximate transaxial PSF yields no improvement. When low image noise levels are required for detection tasks, the deconvolution of the axial PSF does not provide adequate resolution or quantitative benefits to justify its application. When accurate quantitation is required and higher noise levels are acceptable, the deconvolution of the axial PSF leads to considerable gains (30%) in accuracy over conventional FORE+FBP at matched noise levels.