A parametric level-set approach to simultaneous object identification and background reconstruction for dual-energy computed tomography

Virtual monochromatic imaging with projection-based material decomposition algorithm for metal artifacts reduction in photon-counting detector computed tomography

Metal artifacts present a major challenge to computed tomography (CT) because they reduce the image quality in medical diagnosis and treatment. Several metal artifact reduction (MAR) methods have been proposed to address this issue in previous studies. This study aimed to synthesize a virtual monochromatic image for MAR in CT images using projection-based material decomposition (MD) algorithms. We developed a spectral micro-CT prototype system equipped with a photon-counting detector (PCD) and PCD-CT imaging simulator to assess the performances of different MAR methods. Two projection-based MD algorithms were implemented and evaluated for their MAR performances in CT images and compared with conventional sinogram inpainting MAR methods. Different parts of digital 4D-extended cardiac torso (XCAT) phantoms with metal implants were designed to simulate various real scenarios. A homemade metal artifact evaluation (MAE) phantom was used to evaluate the MAR performance in experiments. The simulated results of the XCAT phantom indicated that the projection-based virtual monochromatic CT (VMCT) images provided better image quality than the conventional MAR images without blurring the normal tissues at the position of the metal artifacts. Various quantitative indicators support this conclusion. Additionally, the experimental results of the MAE phantom reveal that projection-based VMCT images can avoid image distortion caused by metal artifacts, unlike conventional MAR methods. In regards to the projection-based VMCT images, the simulated and experimental results demonstrated that using the linear maximum likelihood estimators with an error correction look-up table algorithm yielded better MAR performance compared to that obtained using a polynomial algorithm. Furthermore, projection-based VMCT images can not only reduce metal artifacts effectively but also simultaneously prevents object blurring at the metal artifact position and image distortion of the metal implants. Hence, the CT image quality can be further improved to increase the abilities for both preoperative and postoperative assessment of metal implants.

Significance of the spectral correction of photon counting detector response in material classification from spectral x-ray CT

Purpose: Photon counting imaging detectors (PCD) has paved the way for spectral x-ray computed tomography (spectral CT), which simultaneously measures a sample's linear attenuation coefficient (LAC) at multiple energies. However, cadmium telluride (CdTe)-based PCDs working under high flux suffer from detector effects, such as charge sharing and photon pileup. These effects result in the severe spectral distortions of the measured spectra and significant deviation of the extracted LACs from the reference attenuation curve. We analyze the influence of the spectral distortion correction on material classification performance. Approach: We employ a spectral correction algorithm to reduce the primary spectral distortions. We use a method for material classification that measures system-independent material properties, such as electron density, ρ e , and effective atomic number, Z eff . These parameters are extracted from the LACs using attenuation decomposition and are independent of the scanner specification. The classification performance with the raw and corrected data is tested on different numbers of energy bins and projections and different radiation dose levels. We use experimental data with a broad range of materials in the range of 6 ≤ Z eff ≤ 15 , acquired with a custom laboratory instrument for spectral CT. Results: We show that using the spectral correction leads to an accuracy increase of 1.6 and 3.8 times in estimating ρ e and Z eff , respectively, when the image reconstruction is performed from only 12 projections and the 15 energy bins approach is used. Conclusions: The correction algorithm accurately reconstructs the measured attenuation curve and thus gives better classification performance.

On the relevance of denoising and artefact reduction in 3D segmentation and classification within complex computed tomography imagery

We evaluate the impact of denoising and Metal Artefact Reduction (MAR) on 3D object segmentation and classification in low-resolution, cluttered dual-energy Computed Tomography (CT). To this end, we present a novel 3D materials-based segmentation technique based on the Dual-Energy Index (DEI) to automatically generate subvolumes for classification. Subvolume classification is performed using an extension of Extremely Randomised Clustering (ERC) forest codebooks, constructed using dense feature-point sampling and multiscale Density Histogram (DH) descriptors. Within this experimental framework, we evaluate the impact on classification accuracy and computational expense of pre-processing by intensity thresholding, Non-Local Means (NLM) filtering, Linear Interpolation-based MAR (LIMar) and Distance-Driven MAR (DDMar) in the domain of 3D baggage security screening. We demonstrate that basic NLM filtering, although removing fewer artefacts, produces state-of-the-art classification results comparable to the more complex DDMar but at a significant reduction in computational cost - bringing into question the importance (in terms of automated CT analysis) of computationally expensive artefact reduction techniques. Overall, it was found that the use of MAR pre-processing approaches produced only a marginal improvement in classification performance (< 1%) at considerable additional computational cost (> 10×) when compared to NLM pre-processing.

A Novel Region Reconstruction Method for Fluorescence Molecular Tomography

Fluorescence molecular tomography (FMT) could exploit the distribution of fluorescent biomarkers that target tumors accurately and effectively, which enables noninvasive real-time 3-D visualization as well as quantitative analysis of small tumors in small animal studies in vivo. Due to the difficulties of reconstruction, continuous efforts are being made to find more practical and efficient approaches to accurately obtain the characteristics of fluorescent regions inside biological tissues. In this paper, we propose a region reconstruction method for FMT, which is defined as an L1-norm regularization piecewise constant level set approach. The proposed approach adopts a priori information including the sparsity of the fluorescent sources and the fluorescent contrast between the target and background. When the contrast of different fluorescent sources is low to a certain degree, our approach can simultaneously solve the detection and characterization problems for the reconstruction of FMT. To evaluate the performance of the region reconstruction method, numerical phantom experiments and in vivo bead-implanted mouse experiments were performed. The results suggested that the proposed region reconstruction method was able to reconstruct the features of the fluorescent regions accurately and effectively, and the proposed method was able to be feasibly adopted in in vivo application.

Tensor-based formulation and nuclear norm regularization for multienergy computed tomography

The development of energy selective, photon counting X-ray detectors allows for a wide range of new possibilities in the area of computed tomographic image formation. Under the assumption of perfect energy resolution, here we propose a tensor-based iterative algorithm that simultaneously reconstructs the X-ray attenuation distribution for each energy. We use a multilinear image model rather than a more standard stacked vector representation in order to develop novel tensor-based regularizers. In particular, we model the multispectral unknown as a three-way tensor where the first two dimensions are space and the third dimension is energy. This approach allows for the design of tensor nuclear norm regularizers, which like its 2D counterpart, is a convex function of the multispectral unknown. The solution to the resulting convex optimization problem is obtained using an alternating direction method of multipliers approach. Simulation results show that the generalized tensor nuclear norm can be used as a standalone regularization technique for the energy selective (spectral) computed tomography problem and when combined with total variation regularization it enhances the regularization capabilities especially at low energy images where the effects of noise are most prominent.

Model-Based Iterative Reconstruction for Dual-Energy X-Ray CT Using a Joint Quadratic Likelihood Model

Dual-energy X-ray CT (DECT) has the potential to improve contrast and reduce artifacts as compared to traditional CT. Moreover, by applying model-based iterative reconstruction (MBIR) to dual-energy data, one might also expect to reduce noise and improve resolution. However, the direct implementation of dual-energy MBIR requires the use of a nonlinear forward model, which increases both complexity and computation. Alternatively, simplified forward models have been used which treat the material-decomposed channels separately, but these approaches do not fully account for the statistical dependencies in the channels. In this paper, we present a method for joint dual-energy MBIR (JDE-MBIR), which simplifies the forward model while still accounting for the complete statistical dependency in the material-decomposed sinogram components. The JDE-MBIR approach works by using a quadratic approximation to the polychromatic log-likelihood and a simple but exact nonnegativity constraint in the image domain. We demonstrate that our method is particularly effective when the DECT system uses fast kVp switching, since in this case the model accounts for the inaccuracy of interpolated sinogram entries. Both phantom and clinical results show that the proposed model produces images that compare favorably in quality to previous decomposition-based methods, including FBP and other statistical iterative approaches.

Modelling the physics in the iterative reconstruction for transmission computed tomography

There is an increasing interest in iterative reconstruction (IR) as a key tool to improve quality and increase applicability of x-ray CT imaging. IR has the ability to significantly reduce patient dose; it provides the flexibility to reconstruct images from arbitrary x-ray system geometries and allows one to include detailed models of photon transport and detection physics to accurately correct for a wide variety of image degrading effects. This paper reviews discretization issues and modelling of finite spatial resolution, Compton scatter in the scanned object, data noise and the energy spectrum. The widespread implementation of IR with a highly accurate model-based correction, however, still requires significant effort. In addition, new hardware will provide new opportunities and challenges to improve CT with new modelling.

Modelling the physics in the iterative reconstruction for transmission computed tomography

There is an increasing interest in iterative reconstruction (IR) as a key tool to improve quality and increase applicability of x-ray CT imaging. IR has the ability to significantly reduce patient dose; it provides the flexibility to reconstruct images from arbitrary x-ray system geometries and allows one to include detailed models of photon transport and detection physics to accurately correct for a wide variety of image degrading effects. This paper reviews discretization issues and modelling of finite spatial resolution, Compton scatter in the scanned object, data noise and the energy spectrum. The widespread implementation of IR with a highly accurate model-based correction, however, still requires significant effort. In addition, new hardware will provide new opportunities and challenges to improve CT with new modelling.

Parametric estimation of 3D tubular structures for diffuse optical tomography

We explore the use of diffuse optical tomography (DOT) for the recovery of 3D tubular shapes representing vascular structures in breast tissue. Using a parametric level set method (PaLS) our method incorporates the connectedness of vascular structures in breast tissue to reconstruct shape and absorption values from severely limited data sets. The approach is based on a decomposition of the unknown structure into a series of two dimensional slices. Using a simplified physical model that ignores 3D effects of the complete structure, we develop a novel inter-slice regularization strategy to obtain global regularity. We report on simulated and experimental reconstructions using realistic optical contrasts where our method provides a more accurate estimate compared to an unregularized approach and a pixel based reconstruction.