Why do commercial CT scanners still employ traditional, filtered back-projection for image reconstruction?

Prior-image-based low-dose CT reconstruction for adaptive radiation therapy

Objective. The study aims to reduce the imaging radiation dose in Adaptive Radiotherapy (ART) while maintaining high-quality CT images, critical for effective treatment planning and monitoring.Approach. We developed the Prior-aware Learned Primal-Dual Network (pLPD-UNet), which uses prior CT images to enhance reconstructions from low-dose scans. The network was separately trained on thorax and abdomen datasets to accommodate the unique imaging requirements of each anatomical region.Main results. The pLPD-UNet demonstrated improved reconstruction accuracy and robustness in handling sparse data compared to traditional methods. It effectively maintained image quality essential for precise organ delineation and dose calculation, while achieving a significant reduction in radiation exposure.Significance. This method offers a significant advancement in the practice of ART by integrating prior imaging data, potentially setting a new standard for balancing radiation safety with the need for high-resolution imaging in cancer treatment planning.

Enhancing global sensitivity and uncertainty quantification in medical image reconstruction with Monte Carlo arbitrary-masked mamba

Deep learning has been extensively applied in medical image reconstruction, where Convolutional Neural Networks (CNNs) and Vision Transformers (ViTs) represent the predominant paradigms, each possessing distinct advantages and inherent limitations: CNNs exhibit linear complexity with local sensitivity, whereas ViTs demonstrate quadratic complexity with global sensitivity. The emerging Mamba has shown superiority in learning visual representation, which combines the advantages of linear scalability and global sensitivity. In this study, we introduce MambaMIR, an Arbitrary-Masked Mamba-based model with wavelet decomposition for joint medical image reconstruction and uncertainty estimation. A novel Arbitrary Scan Masking (ASM) mechanism "masks out" redundant information to introduce randomness for further uncertainty estimation. Compared to the commonly used Monte Carlo (MC) dropout, our proposed MC-ASM provides an uncertainty map without the need for hyperparameter tuning and mitigates the performance drop typically observed when applying dropout to low-level tasks. For further texture preservation and better perceptual quality, we employ the wavelet transformation into MambaMIR and explore its variant based on the Generative Adversarial Network, namely MambaMIR-GAN. Comprehensive experiments have been conducted for multiple representative medical image reconstruction tasks, demonstrating that the proposed MambaMIR and MambaMIR-GAN outperform other baseline and state-of-the-art methods in different reconstruction tasks, where MambaMIR achieves the best reconstruction fidelity and MambaMIR-GAN has the best perceptual quality. In addition, our MC-ASM provides uncertainty maps as an additional tool for clinicians, while mitigating the typical performance drop caused by the commonly used dropout.

Compressed-sensing-inspired reconstruction algorithms in low-dose computed tomography: A review

BACKGROUND

Optimization of the dose the patient receives during scanning is an important problem in modern medical X-ray computed tomography (CT). One of the basic ways to its solution is to reduce the number of views. Compressed sensing theory helped promote the development of a new class of effective reconstruction algorithms for limited data CT. These compressed-sensing-inspired (CSI) algorithms optimize the Lp (0 ≤ p ≤ 1) norm of images and can accurately reconstruct CT tomograms from a very few views. The paper presents a review of the CSI algorithms and discusses prospects for their further use in commercial low-dose CT.

METHODS

Many literature references with the CSI algorithms have been were searched. To structure the material collected the author gives a classification framework within which he describes Lp regularization methods, the basic CSI algorithms that are used most often in few-view CT, and some of their derivatives. Lots of examples are provided to illustrate the use of the CSI algorithms in few-view and low-dose CT.

RESULTS

A list of the CSI algorithms is compiled from the literature search. For better demonstrativeness they are summarized in a table. The inference is done that already today some of the algorithms are capable of reconstruction from 20 to 30 views with acceptable quality and dose reduction by a factor of 10.

DISCUSSION

In conclusion the author discusses how soon the CSI reconstruction algorithms can be introduced in the practice of medical diagnosis and used in commercial CT scanners.

Imperfect refractive index matching in scanning laser optical tomography and a method for digital correction

Significance

Scanning laser optical tomography (SLOT) is a volumetric multi-modal imaging technique that is comparable to optical projection tomography and computer tomography. Image quality is crucially dependent on matching the refractive indexes (RIs) of the sample and surrounding medium, but RI matching often requires some effort and is never perfect.

Aim

Reducing the burden of RI matching between the immersion medium and sample in biomedical imaging is a challenging and interesting task. We aim at implementing a post processing strategy for correcting SLOT measurements that have errors caused by RI mismatch.

Approach

To better understand the problems with poorly matched Ris, simulated SLOT measurements with imperfect RI matching of the sample and medium are performed and presented here. A method to correct distorted measurements was developed and is presented and evaluated. This method is then applied to a sample containing fluorescent polystyrene beads and a sample made of olydimethylsiloxane with embedded fluorescent nanoparticles.

Results

From the simulations, it is evident that measurements with an RI mismatch larger than 0.02 and no correction yield considerably worse results compared to perfectly matched measurements. RI mismatches larger than 0.05 make it almost impossible to resolve finer details and structures. By contrast, the simulations imply that a measurement with an RI mismatch of up to 0.1 can still yield reasonable results if the presented correction method is applied. The experiments validate the simulated results for an RI mismatch of about 0.09.

Conclusions

The method significantly improves the SLOT image quality for samples with imperfectly matched Ris. Although the absolutely best imaging quality will be achieved with perfect RI matching, these results pave the way for imaging in SLOT with RI mismatches while maintaining high image quality.

Spectral CT image reconstruction using a constrained optimization approach-An algorithm for AAPM 2022 spectral CT grand challenge and beyond

BACKGROUND

CT reconstruction is of essential importance in medical imaging. In 2022, the American Association of Physicists in Medicine (AAPM) sponsored a Grand Challenge to investigate the challenging inverse problem of spectral CT reconstruction, with the aim of achieving the most accurate reconstruction results. The authors of this paper participated in the challenge and won as a runner-up team.

PURPOSE

This paper reports details of our PROSPECT algorithm (Prior-based Restricted-variable Optimization for SPEctral CT) and follow-up studies regarding the algorithm's accuracy and enhancement of its convergence speed.

METHODS

We formulated the reconstruction task as an optimization problem. PROSPECT employed a one-step backward iterative scheme to solve this optimization problem by allowing estimation of and correction for the difference between the actual polychromatic projection model and the monochromatic model used in the optimization problem. PROSPECT incorporated various forms of prior information derived by analyzing training data provided by the Grand Challenge to reduce the number of unknown variables. We investigated the impact of projection data precision on the resulting solution accuracy and improved convergence speed of the PROSPECT algorithm by incorporating a beam-hardening correction (BHC) step in the iterative process. We also studied the algorithm's performance under noisy projection data.

RESULTS

Prior knowledge allowed a reduction of the number of unknown variables by85.9 % $85.9\%$ . PROSPECT algorithm achieved the average root of mean square error (RMSE) of3.3 × 10 - 6 $3.3\,\times \,10^{-6}$ in the test data set provided by the Grand Challenge. Performing the reconstruction with the same algorithm but using double-precision projection data reduced RMSE to1.2 × 10 - 11 $1.2\,\times \,10^{-11}$ . Including the BHC step in the PROSPECT algorithm accelerated the iteration process with a 40% reduction in computation time.

CONCLUSIONS

PROSPECT algorithm achieved a high degree of accuracy and computational efficiency.

Hierarchical decomposed dual-domain deep learning for sparse-view CT reconstruction

Objective. X-ray computed tomography employing sparse projection views has emerged as a contemporary technique to mitigate radiation dose. However, due to the inadequate number of projection views, an analytic reconstruction method utilizing filtered backprojection results in severe streaking artifacts. Recently, deep learning (DL) strategies employing image-domain networks have demonstrated remarkable performance in eliminating the streaking artifact caused by analytic reconstruction methods with sparse projection views. Nevertheless, it is difficult to clarify the theoretical justification for applying DL to sparse view computed tomography (CT) reconstruction, and it has been understood as restoration by removing image artifacts, not reconstruction.Approach. By leveraging the theory of deep convolutional framelets (DCF) and the hierarchical decomposition of measurement, this research reveals the constraints of conventional image and projection-domain DL methodologies, subsequently, the research proposes a novel dual-domain DL framework utilizing hierarchical decomposed measurements. Specifically, the research elucidates how the performance of the projection-domain network can be enhanced through a low-rank property of DCF and a bowtie support of hierarchical decomposed measurement in the Fourier domain.Main results. This study demonstrated performance improvement of the proposed framework based on the low-rank property, resulting in superior reconstruction performance compared to conventional analytic and DL methods.Significance. By providing a theoretically justified DL approach for sparse-view CT reconstruction, this study not only offers a superior alternative to existing methods but also opens new avenues for research in medical imaging. It highlights the potential of dual-domain DL frameworks to achieve high-quality reconstructions with lower radiation doses, thereby advancing the field towards safer and more efficient diagnostic techniques. The code is available athttps://github.com/hanyoseob/HDD-DL-for-SVCT.

Dual-Domain Reconstruction Network Incorporating Multi-Level Wavelet Transform and Recurrent Convolution for Sparse View Computed Tomography Imaging

Sparse view computed tomography (SVCT) aims to reduce the number of X-ray projection views required for reconstructing the cross-sectional image of an object. While SVCT significantly reduces X-ray radiation dose and speeds up scanning, insufficient projection data give rise to issues such as severe streak artifacts and blurring in reconstructed images, thereby impacting the diagnostic accuracy of CT detection. To address this challenge, a dual-domain reconstruction network incorporating multi-level wavelet transform and recurrent convolution is proposed in this paper. The dual-domain network is composed of a sinogram domain network (SDN) and an image domain network (IDN). Multi-level wavelet transform is employed in both IDN and SDN to decompose sinograms and CT images into distinct frequency components, which are then processed through separate network branches to recover detailed information within their respective frequency bands. To capture global textures, artifacts, and shallow features in sinograms and CT images, a recurrent convolution unit (RCU) based on convolutional long and short-term memory (Conv-LSTM) is designed, which can model their long-range dependencies through recurrent calculation. Additionally, a self-attention-based multi-level frequency feature normalization fusion (MFNF) block is proposed to assist in recovering high-frequency components by aggregating low-frequency components. Finally, an edge loss function based on the Laplacian of Gaussian (LoG) is designed as the regularization term for enhancing the recovery of high-frequency edge structures. The experimental results demonstrate the effectiveness of our approach in reducing artifacts and enhancing the reconstruction of intricate structural details across various sparse views and noise levels. Our method excels in both performance and robustness, as evidenced by its superior outcomes in numerous qualitative and quantitative assessments, surpassing contemporary state-of-the-art CNNs or Transformer-based reconstruction methods.

Advances in medical image analysis with vision Transformers: A comprehensive review

The remarkable performance of the Transformer architecture in natural language processing has recently also triggered broad interest in Computer Vision. Among other merits, Transformers are witnessed as capable of learning long-range dependencies and spatial correlations, which is a clear advantage over convolutional neural networks (CNNs), which have been the de facto standard in Computer Vision problems so far. Thus, Transformers have become an integral part of modern medical image analysis. In this review, we provide an encyclopedic review of the applications of Transformers in medical imaging. Specifically, we present a systematic and thorough review of relevant recent Transformer literature for different medical image analysis tasks, including classification, segmentation, detection, registration, synthesis, and clinical report generation. For each of these applications, we investigate the novelty, strengths and weaknesses of the different proposed strategies and develop taxonomies highlighting key properties and contributions. Further, if applicable, we outline current benchmarks on different datasets. Finally, we summarize key challenges and discuss different future research directions. In addition, we have provided cited papers with their corresponding implementations in https://github.com/mindflow-institue/Awesome-Transformer.

Exploring the Limitations of Hybrid Adiabatic Quantum Computing for Emission Tomography Reconstruction

Our study explores the feasibility of quantum computing in emission tomography reconstruction, addressing a noisy ill-conditioned inverse problem. In current clinical practice, this is typically solved by iterative methods minimizing a L2 norm. After reviewing quantum computing principles, we propose the use of a commercially available quantum annealer and employ corresponding hybrid solvers, which combine quantum and classical computing to handle more significant problems. We demonstrate how to frame image reconstruction as a combinatorial optimization problem suited for these quantum annealers and hybrid systems. Using a toy problem, we analyze reconstructions of binary and integer-valued images with respect to their image size and compare them to conventional methods. Additionally, we test our method's performance under noise and data underdetermination. In summary, our method demonstrates competitive performance with traditional algorithms for binary images up to an image size of 32×32 on the toy problem, even under noisy and underdetermined conditions. However, scalability challenges emerge as image size and pixel bit range increase, restricting hybrid quantum computing as a practical tool for emission tomography reconstruction until significant advancements are made to address this issue.

Iterative image reconstruction with polar coordinate discretized system matrix for optical CT radiochromic gel dosimetry

BACKGROUND

Gel dosimeters are a potential tool for measuring the complex dose distributions that characterize modern radiotherapy. A prototype tabletop solid-tank fan-beam optical CT scanner for readout of gel dosimeters was recently developed. This scanner does not have a straight raypath from source to detector, thus images cannot be reconstructed using filtered backprojection (FBP) and iterative techniques are required. Iterative image reconstruction requires a system matrix that describes the geometry of the imaging system. Stored system matrices can become immensely large, making them impractical for storage on a typical desktop computer.

PURPOSE

Here we develop a method to reduce the storage size of optical CT system matrices through use of polar coordinate discretization while accounting for the refraction in optical CT systems.

METHODS

A ray tracing simulator was developed to track the path of light rays as they traverse the different mediums of the optical CT scanner. Cartesian coordinate discretized system matrices (CCDSMs) and polar coordinate discretized system matrices (PCDSMs) were generated by discretizing the reconstruction area of the optical CT scanner into a Cartesian pixel grid and a polar coordinate pixel grid, respectively. The length of each ray through each pixel was calculated and used to populate the system matrices. To ensure equal weighting during iterative reconstruction, the radial rings of PCDSMs were asymmetrically spaced such that the area of each polar pixel was constant. Two clinical phantoms and several synthetic phantoms were produced and used to evaluate the reconstruction techniques under known conditions. Reconstructed images were analyzed in terms of spatial resolution, signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), signal nonuniformity (SNU), and Gamma map pass percentage.

RESULTS

A storage size reduction of 99.72% was found when comparing a PCDSM to a CCDSM with the same total number of pixels. Images reconstructed with a PCDSM were found to have superior SNR, CNR, SNU, and Gamma (1 mm, 1%) pass percentage compared to those reconstructed with a CCDSM. Increasing spatial resolution in the radial direction with increasing radial distance was found in both PCDSM and CCDSM reconstructions due to the outer regions refracting light more severely. Images reconstructed with a PCDSM showed a decrease in spatial resolution in the azimuthal directions as radial distance increases, due to the widening of the polar pixels. However, this can be mitigated with only a slight increase in storage size by increasing the number of projections. A loss of spatial resolution in the radial direction within 5 mm radially from center was found when reconstructing with a PCDSM, due to the large innermost pixels. However, this was remedied by increasing the number of radial rings within the PCDSM, yielding radial spatial resolution on par with images reconstructed with a CCDSM and a storage size reduction of 99.26%.

CONCLUSIONS

Discretizing the image pixel elements in polar coordinates achieved a system matrix storage size reduction of 99.26% with only minimal reduction in the image quality.

A novel reconstruction method combining multi-detector SPECT with an elliptical orbit and computer tomography for cardiac imaging

The myocardial single photon emission computed tomography (SPECT) is a good study due to its clinical significance in the diagnosis of myocardial disease and the requirement for improving image quality. However, SPECT imaging faces challenges related to low spatial resolution and significant statistical noise, which concerns patient radiation safety. In this paper, a novel reconstruction system combining multi-detector elliptical SPECT (ME-SPECT) and computer tomography (CT) is proposed to enhance spatial resolution and sensitivity. The hybrid imaging system utilizes a slit-slat collimator and elliptical orbit to improve sensitivity and signal-to-noise ratio (SNR), obtains accurate attenuation mapping matrices, and requires prior information from integrated CT. Collimator parameters are corrected based on CT reconstruction results. The SPECT imaging system employs an iterative reconstruction algorithm that utilizes prior knowledge. An iterative reconstruction algorithm based on prior knowledge is applied to the SPECT imaging system, and a method for prioritizing the reconstruction of regions of interest (ROI) is introduced to deal with severely truncated data from ME-SPECT. Simulation results show that the proposed method can significantly improve the system's spatial resolution, SNR, and image fidelity. The proposed method can effectively suppress distortion and artifacts with the higher spatial resolution ordered subsets expectation maximization (OSEM); slit-slat collimation.

Generalized-equiangular geometry CT: Concept and shift-invariant FBP algorithms

BACKGROUND

With advanced x-ray source and detector technologies being continuously developed, non-traditional CT geometries have been widely explored. Generalized-Equiangular Geometry CT (GEGCT) architecture, in which an x-ray source might be positioned radially far away from the focus of arced detector array that is equiangularly spaced, is of importance in many novel CT systems and designs.

PURPOSE

GEGCT, unfortunately, has no theoretically exact and shift-invariant analytical image reconstruction algorithm in general. In this study, to obtain fast and accurate reconstruction from GEGCT and to promote its system design and optimization, an in-depth investigation on a group of approximate Filtered Back-Projection (FBP) algorithms with a variety of weighting strategies has been conducted.

METHODS

The architecture of GEGCT is first presented and characterized by using a normalized-radial-offset distance (NROD). Next, shift-invariant weighted FBP-type algorithms are derived in a unified framework, with pre-filtering, filtering, and post-filtering weights, for both fixed and dynamic NROD configurations. Three viable weighting strategies are then presented including a classic one developed by Besson in the literature and two new ones generated from a curvature fitting and from an empirical formula, where all of the three weights can be expressed as certain functions of NROD. After that, an analysis of reconstruction accuracy is conducted with a wide range of NROD. Finally, the weighted FBP algorithm for GEGCT is extended to a three-dimensional form in the case of cone-beam scan with a cylindrical detector array.

RESULTS

Theoretical analysis and numerical study show that weights in the shift-invariant FBP algorithms can guarantee highly accurate reconstruction for GEGCT. A simulation of Shepp-Logan phantom and a GEGCT scan of lung mimicked by using a clinical lung CT dataset both demonstrate that FBP reconstructions with Besson and polynomial weights can achieve excellent image quality, with Peak Signal to Noise Ratio and Structural Similarity being at the same level as that from the standard equiangular fan-beam CT scan. Reconstruction of a cylinder object with multiple contrasts from simulated GEGCT scan with dynamic NROD is also highly consistent with fixed ones when using the Besson and polynomial weights, with root mean square error less than 7 hounsfield units, demonstrating the robustness and flexibility of the presented FBP algorithms. In terms of resolution, the direct FBP methods for GEGCT could achieve 1.35 lp/mm of spatial resolution at 10% modulation transfer functions point, higher than that of the rebinning method which can only reach 1.14 lp/mm. Moreover, 3D reconstructions of a disc phantom reveal that a greater value of NROD for GEGCT will bring less cone beam artifacts as expected.

CONCLUSIONS

We propose the concept of GEGCT and investigate the feasibility of using shift-invariant weighted FBP-type algorithms for reconstruction from GEGCT data without rebinning. A comprehensive analysis and phantom studies have been conducted to validate the effectiveness of proposed weighting strategies in a wide range of NROD for GEGCT with fixed and dynamic NROD.

Dual-scale similarity-guided cycle generative adversarial network for unsupervised low-dose CT denoising

Removing the noise in low-dose CT (LDCT) is crucial to improving the diagnostic quality. Previously, many supervised or unsupervised deep learning-based LDCT denoising algorithms have been proposed. Unsupervised LDCT denoising algorithms are more practical than supervised ones since they do not need paired samples. However, unsupervised LDCT denoising algorithms are rarely used clinically due to their unsatisfactory denoising ability. In unsupervised LDCT denoising, the lack of paired samples makes the direction of gradient descent full of uncertainty. On the contrary, paired samples used in supervised denoising allow the parameters of networks to have a clear direction of gradient descent. To bridge the gap in performance between unsupervised and supervised LDCT denoising, we propose dual-scale similarity-guided cycle generative adversarial network (DSC-GAN). DSC-GAN uses similarity-based pseudo-pairing to better accomplish unsupervised LDCT denoising. We design a Vision Transformer-based global similarity descriptor and a residual neural network-based local similarity descriptor for DSC-GAN to effectively describe the similarity between two samples. During training, pseudo-pairs, i.e., similar LDCT samples and normal-dose CT (NDCT) samples, dominate parameter updates. Thus, the training can achieve equivalent effect as training with paired samples. Experiments on two datasets demonstrate that DSC-GAN beats the state-of-the-art unsupervised algorithms and reaches a level close to supervised LDCT denoising algorithms.

Tunable neural networks for CT image formation

Optimization of CT image quality typically involves balancing variance and bias. In traditional filtered back-projection, this trade-off is controlled by the filter cutoff frequency. In model-based iterative reconstruction, the regularization strength parameter often serves the same function. Deep neural networks (DNNs) typically do not provide this tunable control over output image properties. Models are often trained to minimize the expected mean squared error, which penalizes both variance and bias in image outputs but does not offer any control over the trade-off between the two. We propose a method for controlling the output image properties of neural networks with a new loss function called weighted covariance and bias (WCB). Our proposed method uses multiple noise realizations of the input images during training to allow for separate weighting matrices for the variance and bias penalty terms. Moreover, we show that tuning these weights enables targeted penalization of specific image features with spatial frequency domain penalties. To evaluate our method, we present a simulation study using digital anthropomorphic phantoms, physical simulation of CT measurements, and image formation with various algorithms. We show that the WCB loss function offers a greater degree of control over trade-offs between variance and bias, whereas mean-squared error provides only one specific image quality configuration. We also show that WCB can be used to control specific image properties including variance, bias, spatial resolution, and the noise correlation of neural network outputs. Finally, we present a method to optimize the proposed weights for a spiculated lung nodule shape discrimination task. Our results demonstrate this new image quality can control the image properties of DNN outputs and optimize image quality for task-specific performance.

A Simple but Universal Fully Linearized ADMM Algorithm for Optimization Based Image Reconstruction

Background and Objective

Optimization based image reconstruction algorithm is an advanced algorithm in medical imaging. However, the corresponding solving algorithm is challenging because the optimization model is usually large-scale and non-smooth. This work aims to devise a simple but universal solver for optimization models.

Methods

The alternating direction method of multipliers (ADMM) algorithm is a simple and effective solver of the optimization models. However, there always exists a sub-problem that has not closed-form solution. One may use gradient descent algorithm to solve this sub-problem, but the step-size selection via line search is time-consuming. Or, one may use fast Fourier transform (FFT) to get a closed-form solution if the system matrix and the sparse transform matrix are both of special structure. In this work, we propose a simple but universal fully linearized ADMM (FL-ADMM) algorithm that avoids line search to determine step-size and applies to system matrix and sparse transform of any structures.

Results

We derive the FL-ADMM algorithm instances for three total variation (TV) models in 2D computed tomography (CT). Further, we validate and evaluate one FL-ADMM algorithm and explore how the two important factors impact convergence rate. Also, we compare this algorithm with the Chambolle-Pock algorithm via real CT phantom reconstructions. These studies show that the FL-ADMM algorithm may accurately solve optimization models in image reconstruction.

Conclusion

The FL-ADMM algorithm is a simple, effective, convergent and universal solver of optimization models in image reconstruction. Compared to the existing ADMM algorithms, the new algorithm does not need time-consuming step-size line-search or special demand to system matrix and sparse transform. It is a rapid prototyping tool for optimization based image reconstruction.

Improved image quality with deep learning reconstruction - a study on a semi-anthropomorphic upper-abdomen phantom

Purpose

To assess image quality of a deep learning reconstruction (DLR) algorithm across dose levels using a semi-anthropomorphic upper-abdominal phantom, and compare with filtered back projection (FBP) and hybrid iterative reconstruction (IR).

Material and methods

CT scans obtained at five dose levels (CTDIvol 5, 10, 15, 20 and 25 mGy) were reconstructed with FBP, hybrid IR (IR50, IR70 and IR90) and DLR of low (DLL), medium (DLM) and high strength (DLH) in 0.625 mm and 2.5 mm slices. CT number, homogeneity, noise, contrast, contrast-to-noise ratio (CNR), noise texture deviation (NTD; a measure of IR-specific artifacts), noise power spectrum (NPS) and task-based transfer function (TTF) were compared between reconstruction algorithms.

Results

CT numbers were highly consistent across reconstruction algorithms. Image noise was significantly reduced with higher levels of DLR. Noise texture (NPS and NTD) was with DLR maintained at comparable levels to FBP, contrary to increasing levels of hybrid IR. Images reconstructed with DLR of low and high strength in 0.625 mm slices showed similar noise characteristics to 2.5 mm slice FBP and IR50, respectively. Dose-reduction potential based on image noise with IR50 as reference was estimated to 35% for DLM and 74% for DLH.

Conclusions

The novel DLR algorithm demonstrates robust noise reduction with maintained noise texture characteristics despite higher algorithm strength, and may have overcome important limitations of IR. There may be potential for dose reduction and additional benefit from thin-slice reconstruction.

Absolute Grüneisen parameter measurement in deep tissue based on X-ray-induced acoustic computed tomography

The Grüneisen parameter is a primary parameter of the initial sound pressure signal in the photoacoustic effect, which can provide unique biological information and is related to the temperature change information of an object. The accurate measurement of this parameter is of great significance in biomedical research. Combining X-ray-induced acoustic tomography and conventional X-ray computed tomography, we proposed a method to obtain the absolute Grüneisen parameter. The theory development, numerical simulation, and biomedical application scenarios are discussed. The results reveal that our method not only can determine the Grüneisen parameter but can also obtain the body internal temperature distribution, presenting its potential in the diagnosis of a broad range of diseases.

Dual-domain sparse-view CT reconstruction with Transformers

PURPOSE

Computed Tomography (CT) has been widely used in the medical field. Sparse-view CT is an effective and feasible method to reduce the radiation dose. However, the conventional filtered back projection (FBP) algorithm will suffer from severe artifacts in sparse-view CT. Iterative reconstruction algorithms have been adopted to remove artifacts, but they are time-consuming due to repeated projection and back projection and may cause blocky effects. To overcome the difficulty in sparse-view CT, we proposed a dual-domain sparse-view CT algorithm CT Transformer (CTTR) and paid attention to sinogram information.

METHODS

CTTR treats sinograms as sentences and enhances reconstructed images with sinogram's characteristics. We qualitatively evaluate the CTTR, an iterative method TVM-POCS, a convolutional neural network based method FBPConvNet in terms of a reduction in artifacts and a preservation of details. Besides, we also quantitatively evaluate these methods in terms of RMSE, PSNR and SSIM.

RESULTS

We evaluate our method on the Lung Image Database Consortium image collection with different numbers of projection views and noise levels. Experiment studies show that, compared with other methods, CTTR can reduce more artifacts and preserve more details on various scenarios. Specifically, CTTR improves the FBPConvNet performance of PSNR by 0.76dB with 30 projections.

CONCLUSIONS

The performance of our proposed CTTR is better than the method based on CNN in the case of extremely sparse views both on visual results and quantitative evaluation. Our proposed method provides a new idea for the application of Transformers to CT image processing.

Reflective Tomography Lidar Image Reconstruction for Long Distance Non-Cooperative Target

In the long-distance space target detection, the technique of laser reflection tomography (LRT) shows great power and attracts more attention for further study and real use. However, space targets are often non-cooperative, and normally a 360° complete view of reflection projections cannot be obtained. Therefore, this article firstly introduces an improved LRT system design with more advanced laser equipment for long-distance non-cooperative detection to ensure the high quality of the lidar beam and the lidar projection data. Then, the LRT image reconstruction is proposed and focused on the laser image reconstruction method utilizing the total variation (TV) minimization approach based on the sparse algebraic reconstruction technique (ART) model, in order to reconstruct the laser image in a sparse or incomplete view of projections. At last, comparative experiments with the system are performed to validate the advantages of this method with the LRT system. In both near and far field experiments, the effectiveness and superiority of the proposed method are verified for different types of projection data through comparison to typical methods.

End-to-end deep learning for interior tomography with low-dose x-ray CT

Objective. There are several x-ray computed tomography (CT) scanning strategies used to reduce radiation dose, such as (1) sparse-view CT, (2) low-dose CT and (3) region-of-interest (ROI) CT (called interior tomography). To further reduce the dose, sparse-view and/or low-dose CT settings can be applied together with interior tomography. Interior tomography has various advantages in terms of reducing the number of detectors and decreasing the x-ray radiation dose. However, a large patient or a small field-of-view (FOV) detector can cause truncated projections, and then the reconstructed images suffer from severe cupping artifacts. In addition, although low-dose CT can reduce the radiation exposure dose, analytic reconstruction algorithms produce image noise. Recently, many researchers have utilized image-domain deep learning (DL) approaches to remove each artifact and demonstrated impressive performances, and the theory of deep convolutional framelets supports the reason for the performance improvement. Approach. In this paper, we found that it is difficult to solve coupled artifacts using an image-domain convolutional neural network (CNN) based on deep convolutional framelets. Significance. To address the coupled problem, we decouple it into two sub-problems: (i) image-domain noise reduction inside the truncated projection to solve low-dose CT problem and (ii) extrapolation of the projection outside the truncated projection to solve the ROI CT problem. The decoupled sub-problems are solved directly with a novel proposed end-to-end learning method using dual-domain CNNs. Main results. We demonstrate that the proposed method outperforms the conventional image-domain DL methods, and a projection-domain CNN shows better performance than the image-domain CNNs commonly used by many researchers.

NeuRec: Incorporating Interpatient prior to Sparse-View Image Reconstruction for Neurorehabilitation

Medical imaging technologies such as computed tomography (CT) and magnetic resonance imaging (MRI) imaging are indispensable for contemporary neurorehabilitation diagnostics, intervention, and monitoring. It would be desirable to reconstruct images from sparse measurements to reduce the ionizing radiation and motion artifacts. Although recent coordinate-based representation methods have shown promise advances for sparse-view reconstruction, they overfit a single MLP on a single patient. In this work, we generalize it across many patients by incorporating an interpatient prior into the ill-posed inverse/reconstruction problem, which is the missing ingredient in the previous works. The experiment demonstrates that our method significantly improves image quality over the state-of-the-art both qualitatively and quantitatively. Thus, our method provides a powerful and principled means to deal with the measurement-scarce problem.

Combining convolutional sparse coding with total variation for sparse-view CT reconstruction

Conventional dictionary-learning-based computed tomography (CT) reconstruction methods extract patches from an original image to train, ignoring the consistency of pixels in overlapping patches. To address the problem, this paper proposes a method combining convolutional sparse coding (CSC) with total variation (TV) for sparse-view CT reconstruction. The proposed method inherits the advantages of CSC by directly processing the whole image without dividing it into overlapping patches, which preserves more details and reduces artifacts caused by patch aggregation. By introducing a TV regularization term to enhance the constraint of the image domain, the noise can be effectively further suppressed. The alternating direction method of multipliers algorithm is employed to solve the objective function. Numerous experiments are conducted to validate the performance of the proposed method in different views. Qualitative and quantitative results show the superiority of the proposed method in terms of noise suppression, artifact reduction, and image details recovery.

Control of Variance and Bias in CT Image Processing with Variational Training of Deep Neural Networks

Optimization of CT image quality typically involves balancing noise and bias. In filtered back-projection, this trade-off is controlled by the particular filter and cutoff frequency. In penalized-likelihood iterative reconstruction, the penalty weight serves the same function. Deep neural networks typically do not provide this tuneable control over output image properties. Models are often trained to minimize mean squared error which penalizes both variance and bias in image outputs but does not offer any control over the trade-off between the two. In this work, we propose a method for controlling the output image properties of neural networks with a new loss function called weighted covariance and bias (WCB). Our proposed method includes separate weighting parameters to control the relative importance of noise or bias reduction. Moreover, we show that tuning these weights enables targeted penalization of specific image features (e.g. spatial frequencies). To evaluate our method, we present a simulation study using digital anthropormorphic phantoms, physical simulation of non-ideal CT data, and image formation with various algorithms. We show that WCB offers a greater degree of control over trade-offs between variance and bias whereas MSE has only one configuration. We also show that WCB can be used to control specific image properties including variance, bias, spatial resolution, and the noise correlation of neural network outputs. Finally, we present a method to optimize the proposed weights for stimulus detectability. Our results demonstrate the potential for this new capability to control the image properties of DNN outputs and optimize image quality for the task-specific applications.

Performance evaluation of segmentation methods for assessing the lens of the frog Thoropa miliaris from synchrotron-based phase-contrast micro-CT images

PURPOSE

In the context of synchrotron microtomography using propagation-based phase-contrast imaging (XSPCT), we evaluated the performance of semiautomatic and automatic image segmentation of soft biological structures by means of Dice Similarity Coefficient (DSC) and volume quantification.

METHODS

We took advantage of the phase-contrast effects of XSPCT to provide enhanced object boundaries and improved visualization of the lenses of the frog Thoropa miliaris. Then, we applied semiautomatic segmentation methods 1 and 2 (Interpolation and Watershed, respectively) and method 3, an automatic segmentation algorithm using the U-Net architecture, to the reconstructed images. DSC and volume quantification of the lenses were used to quantify the performance of image segmentation methods.

RESULTS

Comparing the lenses segmented by the three methods, the most pronounced difference in volume quantification was between methods 1 and 3: a reduction of 4.24%. Method 1, 2 and 3 obtained the global average DSC of 97.02%, 95.41% and 89.29%, respectively. Although it obtained the lowest DSC, method 3 performed the segmentation in a matter of seconds, while the semiautomatic methods had the average time to segment the lenses around 1 h and 30 min.

CONCLUSIONS

Our results suggest that the performance of U-Net was impaired due to the irregularities of the ROI edges mainly in its lower and upper regions, but it still showed high accuracy (DSC = 89.29%) with significantly reduced segmentation time compared to the semiautomatic methods. Besides, with the present work we have established a baseline for future assessments of Deep Neural Networks applied to XSPCT volumes.

Low-contrast detectability and potential for radiation dose reduction using deep learning image reconstruction—A 20-reader study on a semi-anthropomorphic liver phantom

Background

A novel deep learning image reconstruction (DLIR) algorithm for CT has recently been clinically approved.

Purpose

To assess low-contrast detectability and dose reduction potential for CT images reconstructed with the DLIR algorithm and compare with filtered back projection (FBP) and hybrid iterative reconstruction (IR).

Material and methods

A customized upper-abdomen phantom containing four cylindrical liver inserts with low-contrast lesions was scanned at CT dose indexes of 5, 10, 15, 20 and 25 mGy. Images were reconstructed with FBP, 50% hybrid IR (IR50), and DLIR of low strength (DLL), medium strength (DLM) and high strength (DLH). Detectability was assessed by 20 independent readers using a two-alternative forced choice approach. Dose reduction potential was estimated separately for each strength of DLIR using a fitted model, with the detectability performance of FBP and IR50 as reference.

Results

For the investigated dose levels of 5 and 10 mGy, DLM improved detectability compared to FBP by 5.8 and 6.9 percentage points (p.p.), and DLH improved detectability by 9.6 and 12.3 p.p., respectively (all p < .007). With IR50 as reference, DLH improved detectability by 5.2 and 9.8 p.p. for the 5 and 10 mGy dose level, respectively (p < .03). With respect to this low-contrast detectability task, average dose reduction potential relative to FBP was estimated to 39% for DLM and 55% for DLH. Relative to IR50, average dose reduction potential was estimated to 21% for DLM and 42% for DLH.

Conclusions:

Low-contrast detectability performance is improved when applying a DLIR algorithm, with potential for radiation dose reduction.

•

Deep learning image reconstruction improves low-contrast detectability in CT.

•

Performance improved with increasing strength of deep learning image reconstruction.

•

Results suggest potential for CT radiation dose reduction.

Foam-like phantoms for comparing tomography algorithms

Tomographic algorithms are often compared by evaluating them on certain benchmark datasets. For fair comparison, these datasets should ideally (i) be challenging to reconstruct, (ii) be representative of typical tomographic experiments, (iii) be flexible to allow for different acquisition modes, and (iv) include enough samples to allow for comparison of data-driven algorithms. Current approaches often satisfy only some of these requirements, but not all. For example, real-world datasets are typically challenging and representative of a category of experimental examples, but are restricted to the acquisition mode that was used in the experiment and are often limited in the number of samples. Mathematical phantoms are often flexible and can sometimes produce enough samples for data-driven approaches, but can be relatively easy to reconstruct and are often not representative of typical scanned objects. In this paper, we present a family of foam-like mathematical phantoms that aims to satisfy all four requirements simultaneously. The phantoms consist of foam-like structures with more than 100000 features, making them challenging to reconstruct and representative of common tomography samples. Because the phantoms are computer-generated, varying acquisition modes and experimental conditions can be simulated. An effectively unlimited number of random variations of the phantoms can be generated, making them suitable for data-driven approaches. We give a formal mathematical definition of the foam-like phantoms, and explain how they can be generated and used in virtual tomographic experiments in a computationally efficient way. In addition, several 4D extensions of the 3D phantoms are given, enabling comparisons of algorithms for dynamic tomography. Finally, example phantoms and tomographic datasets are given, showing that the phantoms can be effectively used to make fair and informative comparisons between tomography algorithms.

A Review of Sensing Technologies for Non-Destructive Evaluation of Structural Composite Materials

The growing demand and diversity in the application of industrial composites and the current inability of present non-destructive evaluation (NDE) methods to perform detailed inspection of these composites has motivated this comprehensive review of sensing technologies. NDE has the potential to be a versatile tool for maintaining composite structures deployed in hazardous and inaccessible areas, such as offshore wind farms and nuclear power plants. Therefore, the future composite solutions need to take into consideration the niche requirements of these high-value/critical applications. Composite materials are intrinsically complex due to their anisotropic and non-homogeneous characteristics. This presents a significant challenge for evaluation and the associated data analysis for NDEs. For example, the quality assurance, certification of composite structures, and early detection of the failure is complex due to the variability and tolerances involved in the composite manufacturing. Adapting existing NDE methods to detect and locate the defects at multiple length scales in the complex materials represents a significant challenge, resulting in a delayed and incorrect diagnosis of the structural health. This paper presents a comprehensive review of the NDE techniques, that includes a detailed discussion of their working principles, setup, advantages, limitations, and usage level for the structural composites. A comparison between these techniques is also presented, providing an insight into the future trends for composites’ prognostic and health management (PHM). Current research trends show the emergence of the non-contact-type NDE (including digital image correlation, infrared tomography, as well as disruptive frequency-modulated continuous wave techniques) for structural composites, and the reasons for their choice over the most popular contact-type (ultrasonic, acoustic, and piezoelectric testing) NDE methods is also discussed. The analysis of this new sensing modality for composites’ is presented within the context of the state-of-the-art and projected future requirements.

Iterative reconstruction improves image quality and reduces radiation dose in trauma protocols; A human cadaver study

Background

Radiation-related cancer risk is an object of concern in CT of trauma patients, as these represent a young population. Different radiation reducing methods, including iterative reconstruction (IR), and spilt bolus techniques have been introduced in the recent years in different large scale trauma centers.

Purpose

To compare image quality in human cadaver exposed to thoracoabdominal computed tomography using IR and standard filtered back-projection (FBP) at different dose levels.

Material and methods

Ten cadavers were scanned at full dose and a dose reduction in CTDIvol of 5 mGy (low dose 1) and 7.5 mGy (low dose 2) on a Siemens Definition Flash 128-slice computed tomography scanner. Low dose images were reconstructed with FBP and Sinogram affirmed iterative reconstruction (SAFIRE) level 2 and 4. Quantitative image quality was analyzed by comparison of contrast-to-noise ratio (CNR) and signal-to-noise ratio (SNR). Qualitative image quality was evaluated by use of visual grading regression (VGR) by four radiologists.

Results

Readers preferred SAFIRE reconstructed images over FBP at a dose reduction of 40% (low dose 1) and 56% (low dose 2), with significant difference in overall impression of image quality. CNR and SNR showed significant improvement for images reconstructed with SAFIRE 2 and 4 compared to FBP at both low dose levels.

Conclusions

Iterative image reconstruction, SAFIRE 2 and 4, resulted in equal or improved image quality at a dose reduction of up to 56% compared to full dose FBP and may be used a strong radiation reduction tool in the young trauma population.

Superiorization versus regularization: A comparison of algorithms for solving image reconstruction problems with applications in computed tomography

PURPOSE

A system matrix can be built in order to account for the refractions in an optical computed tomography (CT) system. In order to utilize this system matrix, iterative methods are employed to solve the image reconstruction problem. The purpose of this study is to compare potential iterative algorithms to solve this image reconstruction problem. Comparisons examine both solution time and the quality of the reconstructed image. While our work is motivated by optical CT, the results can be extended more generally to CT.

METHODS

A collection of 21 algorithms for solving the image reconstruction problem were evaluated. Specifically, algorithms using (i) superiorization techniques and (ii) regularization to avoid overfitting were compared. Multiple test problems are investigated using 18 different image phantoms, parallel-beam and fan-beam system matrices, and varying noise levels. Comparison of the algorithms is done using performance profiles on three different performance measures.

RESULTS

The results for both the synthetic and clinical test problems show that there is not one single algorithm outperforming all others, but instead a set of top algorithms that give the best values on the performance profiles. When qualitative analyses such as reliance on stopping conditions, number of input parameters, and run time are also considered, FISTA-TV shows slight advantages over the other top algorithms.

CONCLUSIONS

There is a set of top algorithms that all show good results in the performance profiles with a mix of superiorized and regularized model algorithms. As to which of these top algorithms outperforms the rest is undetermined and further research needs to be investigated.

Tomosipo: fast, flexible, and convenient 3D tomography for complex scanning geometries in Python

Tomography is a powerful tool for reconstructing the interior of an object from a series of projection images. Typically, the source and detector traverse a standard path (e.g., circular, helical). Recently, various techniques have emerged that use more complex acquisition geometries. Current software packages require significant handwork, or lack the flexibility to handle such geometries. Therefore, software is needed that can concisely represent, visualize, and compute reconstructions of complex acquisition geometries. We present tomosipo, a Python package that provides these capabilities in a concise and intuitive way. Case studies demonstrate the power and flexibility of tomosipo.

From EMI to AI: a brief history of commercial CT reconstruction algorithms

Computed tomography was one of the first imaging modalities to require a computerized solution of an inverse problem to produce a useful image from the data acquired by the sensor hardware. The computerized solutions, which are known as image reconstruction algorithms, have thus been a critical component of every CT scanner ever sold. We review the history of commercially deployed CT reconstruction algorithms and consider the forces that led, at various points, both to innovation and to convergence around certain broadly useful algorithms. The forces include the emergence of new hardware capabilities, competitive pressures, the availability of computational power, and regulatory considerations. We consider four major historical periods and turning points. The original EMI scanner was developed with an iterative reconstruction algorithm, but an explosion of innovation coupled with rediscovery of an older literature led to the development of alternative algorithms throughout the early 1970s. Most CT vendors quickly converged on the use of the filtered back-projection (FBP) algorithm, albeit layered with a variety of proprietary corrections in both projection data and image domains to improve image quality. Innovations such as helical scanning and multi-row detectors were both enabled by and drove the development of additional applications of FBP in the 1990s and 2000s. Finally, the last two decades have seen a return of iterative reconstruction and the introduction of artificial intelligence approaches that benefit from increased computational power to reduce radiation dose and improve image quality.

Ultra-Low-Dose Spectral CT Based on a Multi-level Wavelet Convolutional Neural Network

Spectral computed tomography (CT) based on a photon-counting detector (PCD) is a promising technique with the potential to improve lesion detection, tissue characterization, and material decomposition. PCD-based scanners have several technical issues including operation in the step-and-scan mode and long data acquisition time. One straightforward solution to these issues is to reduce the number of projection views. However, if the projection data are under-sampled or noisy, it would be challenging to produce a correct solution without precise prior information. Recently, deep-learning approaches have demonstrated impressive performance for under-sampled CT reconstruction. In this work, the authors present a multilevel wavelet convolutional neural network (MWCNN) to address the limitations of PCD-based scanners. Data properties of the proposed method in under-sampled spectral CT are analyzed with respect to the proposed deep-running-network-based image reconstruction using two measures: sampling density and data incoherence. This work presents the proposed method and four different methods to restore sparse sampling. We investigate and compare these methods through a simulation and real experiments. In addition, data properties are quantitatively analyzed and compared for the effect of sparse sampling on the image quality. Our results indicate that both sampling density and data incoherence affect the image quality in the studied methods. Among the different methods, the proposed MWCNN shows promising results. Our method shows the highest performance in terms of various evaluation parameters such as the structural similarity, root mean square error, and resolution. Based on the results of imaging and quantitative evaluation, this study confirms that the proposed deep-running network structure shows excellent image reconstruction in sparse-view PCD-based CT. These results demonstrate the feasibility of sparse-view PCD-based CT using the MWCNN. The advantage of sparse view CT is that it can significantly reduce the radiation dose and obtain images with several energy bands by fusing PCDs. These results indicate that the MWCNN possesses great potential for sparse-view PCD-based CT.

Autoregression and Structured Low-Rank Modeling of Sinogram Neighborhoods

Sinograms are commonly used to represent the raw data from tomographic imaging experiments. Although it is already well-known that sinograms posess some amount of redundancy, in this work, we present novel theory suggesting that sinograms will often possess substantial additional redundancies that have not been explicitly exploited by previous methods. Specifically, we derive that sinograms will often satisfy multiple simple data-dependent autoregression relationships. This kind of autoregressive structure enables missing/degraded sinogram samples to be linearly predicted using a simple shift-invariant linear combination of neighboring samples. Our theory also further implies that if sinogram samples are assembled into a structured Hankel/Toeplitz matrix, then the matrix will be expected to have low-rank characteristics. As a result, sinogram restoration problems can be formulated as structured low-rank matrix recovery problems. Illustrations of this approach are provided using several different (real and simulated) X-ray imaging datasets, including comparisons against a state-of-the-art deep learning approach. Results suggest that structured low-rank matrix methods for sinogram recovery can have comparable performance to state-of-the-art approaches. Although our evaluation focuses on competitive comparisons against other approaches, we believe that autoregressive constraints are actually complementary to existing approaches with strong potential synergies.

Reducing residual-motion artifacts in iterative 3D CBCT reconstruction in image-guided radiation therapy

PURPOSE

Recent evaluations of a 3D iterative cone-beam computed tomography (iCBCT) reconstruction method available on Varian radiation treatment devices demonstrated that iCBCT provides superior image quality when compared to analytical Feldkamp-Davis-Kress (FDK) method. However, iCBCT employs statistical penalized likelihood (PL) that is known to be highly sensitive to inconsistencies due to physiological motion occurring during the acquisition. We propose a computationally inexpensive extension of iCBCT addressing this deficiency.

METHODS

During the iterative process, the gradients of PL are modified to avoid the generation of motion-related artifacts. To assess the impact of this modification, we propose a motion simulation generating CBCT projections of a moving anatomy together with artifact-free images used as ground truth. Contrast-to-noise ratio and power spectra of difference images are computed to quantify the impact of the motion on reconstructed CBCT volumes as well as the effect of the proposed modification.

RESULTS

Using both simulated and clinical data, it is shown that the motion of patient's abdominal wall during the acquisition results in artifacts that can be quantified as low-frequency components in volumes reconstructed with iCBCT. Further, a quantitative evaluation demonstrates that the proposed modification of PL reduces these low-frequency components. While preserving the advantages of PL, it effectively suppresses the propagation of motion-related artifacts into clinically important regions, thus increasing the motion resiliency of iCBCT.

CONCLUSIONS

The proposed modified iterative reconstruction method significantly improves the quality of CBCT images of anatomies suffering from residual motion.

Synchrotron X-Ray Tomography for Rechargeable Battery Research: Fundamentals, Setups and Applications

Understanding the complicated interplay of the continuously evolving electrode materials in their inherent 3D states during the battery operating condition is of great importance for advancing rechargeable battery research. In this regard, the synchrotron X-ray tomography technique, which enables non-destructive, multi-scale, and 3D imaging of a variety of electrode components before/during/after battery operation, becomes an essential tool to deepen this understanding. The past few years have witnessed an increasingly growing interest in applying this technique in battery research. Hence, it is time to not only summarize the already obtained battery-related knowledge by using this technique, but also to present a fundamental elucidation of this technique to boost future studies in battery research. To this end, this review firstly introduces the fundamental principles and experimental setups of the synchrotron X-ray tomography technique. After that, a user guide to its application in battery research and examples of its applications in research of various types of batteries are presented. The current review ends with a discussion of the future opportunities of this technique for next-generation rechargeable batteries research. It is expected that this review can enhance the reader's understanding of the synchrotron X-ray tomography technique and stimulate new ideas and opportunities in battery research.

Core Imaging Library - Part I: a versatile Python framework for tomographic imaging

We present the Core Imaging Library (CIL), an open-source Python framework for tomographic imaging with particular emphasis on reconstruction of challenging datasets. Conventional filtered back-projection reconstruction tends to be insufficient for highly noisy, incomplete, non-standard or multi-channel data arising for example in dynamic, spectral and in situ tomography. CIL provides an extensive modular optimization framework for prototyping reconstruction methods including sparsity and total variation regularization, as well as tools for loading, preprocessing and visualizing tomographic data. The capabilities of CIL are demonstrated on a synchrotron example dataset and three challenging cases spanning golden-ratio neutron tomography, cone-beam X-ray laminography and positron emission tomography. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 2'.

Core Imaging Library - Part I: a versatile Python framework for tomographic imaging

We present the Core Imaging Library (CIL), an open-source Python framework for tomographic imaging with particular emphasis on reconstruction of challenging datasets. Conventional filtered back-projection reconstruction tends to be insufficient for highly noisy, incomplete, non-standard or multi-channel data arising for example in dynamic, spectral and in situ tomography. CIL provides an extensive modular optimization framework for prototyping reconstruction methods including sparsity and total variation regularization, as well as tools for loading, preprocessing and visualizing tomographic data. The capabilities of CIL are demonstrated on a synchrotron example dataset and three challenging cases spanning golden-ratio neutron tomography, cone-beam X-ray laminography and positron emission tomography. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 2'.

Manual kidney stone size measurements in computed tomography are most accurate using multiplanar image reformatations and bone window settings

Computed tomography in suspected urolithiasis provides information about the presence, location and size of stones. Particularly stone size is a key parameter in treatment decision; however, data on impact of reformatation and measurement strategies is sparse. This study aimed to investigate the influence of different image reformatations, slice thicknesses and window settings on stone size measurements. Reference stone sizes of 47 kidney stones representative for clinically encountered compositions were measured manually using a digital caliper (Man-M). Afterwards stones were placed in a 3D-printed, semi-anthropomorphic phantom, and scanned using a low dose protocol (CTDIvol 2 mGy). Images were reconstructed using hybrid-iterative and model-based iterative reconstruction algorithms (HIR, MBIR) with different slice thicknesses. Two independent readers measured largest stone diameter on axial (2 mm and 5 mm) and multiplanar reformatations (based upon 0.67 mm reconstructions) using different window settings (soft-tissue and bone). Statistics were conducted using ANOVA ± correction for multiple comparisons. Overall stone size in CT was underestimated compared to Man-M (8.8 ± 2.9 vs. 7.7 ± 2.7 mm, p < 0.05), yet closely correlated (r = 0.70). Reconstruction algorithm and slice thickness did not significantly impact measurements (p > 0.05), while image reformatations and window settings did (p < 0.05). CT measurements using multiplanar reformatation with a bone window setting showed closest agreement with Man-M (8.7 ± 3.1 vs. 8.8 ± 2.9 mm, p < 0.05, r = 0.83). Manual CT-based stone size measurements are most accurate using multiplanar image reformatation with a bone window setting, while measurements on axial planes with different slice thicknesses underestimate true stone size. Therefore, this procedure is recommended when impacting treatment decision.

A geometry-guided deep learning technique for CBCT reconstruction

Purpose.Although deep learning (DL) technique has been successfully used for computed tomography (CT) reconstruction, its implementation on cone-beam CT (CBCT) reconstruction is extremely challenging due to memory limitations. In this study, a novel DL technique is developed to resolve the memory issue, and its feasibility is demonstrated for CBCT reconstruction from sparsely sampled projection data.Methods.The novel geometry-guided deep learning (GDL) technique is composed of a GDL reconstruction module and a post-processing module. The GDL reconstruction module learns and performs projection-to-image domain transformation by replacing the traditional single fully connected layer with an array of small fully connected layers in the network architecture based on the projection geometry. The DL post-processing module further improves image quality after reconstruction. We demonstrated the feasibility and advantage of the model by comparing ground truth CBCT with CBCT images reconstructed using (1) GDL reconstruction module only, (2) GDL reconstruction module with DL post-processing module, (3) Feldkamp, Davis, and Kress (FDK) only, (4) FDK with DL post-processing module, (5) ray-tracing only, and (6) ray-tracing with DL post-processing module. The differences are quantified by peak-signal-to-noise ratio (PSNR), structural similarity index measure (SSIM), and root-mean-square error (RMSE).Results.CBCT images reconstructed with GDL show improvements in quantitative scores of PSNR, SSIM, and RMSE. Reconstruction time per image for all reconstruction methods are comparable. Compared to current DL methods using large fully connected layers, the estimated memory requirement using GDL is four orders of magnitude less, making DL CBCT reconstruction feasible.Conclusion.With much lower memory requirement compared to other existing networks, the GDL technique is demonstrated to be the first DL technique that can rapidly and accurately reconstruct CBCT images from sparsely sampled data.

A Novel Iterative MLEM Image Reconstruction Algorithm Based on Beltrami Filter: Application to ECT Images

The implementation of emission-computed tomography (ECT), including positron emission tomography and single-photon emission-computed tomography, has been an important research topic in recent years and is of significant and practical importance. However, the slow rate of convergence and the computational complexity have severely impeded the efficient implementation of iterative reconstruction. By combining the maximum-likelihood expectation maximization (MLEM) iteratively along with the Beltrami filter, this paper proposes a new approach to reformulate the MLEM algorithm. Beltrami filtering is applied to an image obtained using the MLEM algorithm for each iteration. The role of Beltrami filtering is to remove mainly out-of-focus slice blurs, which are artifacts present in most existing images. To improve the quality of an image reconstructed using MLEM, the Beltrami filter employs similar structures, which in turn reduce the number of errors in the reconstructed image. Numerical image reconstruction tomography experiments have demonstrated the performance capability of the proposed algorithm in terms of an increase in signal-to-noise ratio (SNR) and the recovery of fine details that can be hidden in the data. The SNR and visual inspections of the reconstructed images are significantly improved compared to those of a standard MLEM. We conclude that the proposed algorithm provides an edge-preserving image reconstruction and substantially suppress noise and edge artifacts.

Photoacoustic computed tomography for functional human brain imaging [Invited]

The successes of magnetic resonance imaging and modern optical imaging of human brain function have stimulated the development of complementary modalities that offer molecular specificity, fine spatiotemporal resolution, and sufficient penetration simultaneously. By virtue of its rich optical contrast, acoustic resolution, and imaging depth far beyond the optical transport mean free path (∼1 mm in biological tissues), photoacoustic computed tomography (PACT) offers a promising complementary modality. In this article, PACT for functional human brain imaging is reviewed in its hardware, reconstruction algorithms, in vivo demonstration, and potential roadmap.

StaticCodeCT: single coded aperture tensorial X-ray CT

Coded aperture X-ray CT (CAXCT) is a new low-dose imaging technology that promises far-reaching benefits in industrial and clinical applications. It places various coded apertures (CA) at a time in front of the X-ray source to partially block the radiation. The ill-posed inverse reconstruction problem is then solved using l1-norm-based iterative reconstruction methods. Unfortunately, to attain high-quality reconstructions, the CA patterns must change in concert with the view-angles making the implementation impractical. This paper proposes a simple yet radically different approach to CAXCT, which is coined StaticCodeCT, that uses a single-static CA in the CT gantry, thus making the imaging system amenable for practical implementations. Rather than using conventional compressed sensing algorithms for recovery, we introduce a new reconstruction framework for StaticCodeCT. Namely, we synthesize the missing measurements using low-rank tensor completion principles that exploit the multi-dimensional data correlation and low-rank nature of a 3-way tensor formed by stacking the 2D coded CT projections. Then, we use the FDK algorithm to recover the 3D object. Computational experiments using experimental projection measurements exhibit up to 10% gains in the normalized root mean square distance of the reconstruction using the proposed method compared with those attained by alternative low-dose systems.

Improved image quality in abdominal computed tomography reconstructed with a novel Deep Learning Image Reconstruction technique - initial clinical experience

BACKGROUND

A novel Deep Learning Image Reconstruction (DLIR) technique for computed tomography has recently received clinical approval.

PURPOSE

To assess image quality in abdominal computed tomography reconstructed with DLIR, and compare with standardly applied iterative reconstruction.

MATERIAL AND METHODS

Ten abdominal computed tomography scans were reconstructed with iterative reconstruction and DLIR of medium and high strength, with 0.625 mm and 2.5 mm slice thickness. Image quality was assessed using eight visual grading criteria in a side-by-side comparative setting. All series were presented twice to evaluate intraobserver agreement. Reader scores were compared using univariate logistic regression. Image noise and contrast-to-noise ratio were calculated for quantitative analyses.

RESULTS

For 2.5 mm slice thickness, DLIR images were more frequently perceived as equal or better than iterative reconstruction across all visual grading criteria (for both DLIR of medium and high strength, p < 0.001). Correspondingly, DLIR images were more frequently perceived as better (as opposed to equal or in favor of iterative reconstruction) for visual reproduction of liver parenchyma, intrahepatic vascular structures as well as overall impression of image noise and texture (p < 0.001). This improved image quality was also observed for 0.625 mm slice images reconstructed with DLIR of high strength when directly comparing to traditional iterative reconstruction in 2.5 mm slices. Image noise was significantly lower and contrast-to-noise ratio measurements significantly higher for images reconstructed with DLIR compared to iterative reconstruction (p < 0.01).

CONCLUSIONS

Abdominal computed tomography images reconstructed using a DLIR technique shows improved image quality when compared to standardly applied iterative reconstruction across a variety of clinical image quality criteria.

Limited-view binary tomography reconstruction assisted by shape centroid

In this paper, the binary tomographic reconstruction problem for very limited projection data availability is considered. Being this inverse problem highly ill-posed, we propose a new reconstruction model that uses a shape centroid-based regularization term, i.e., we assume that the center of gravity of the object of interest is known, at least approximately, in advance. Motivation for this regularization is found in the close connection between the projection data and the object centroid, as we will show. Experimental evaluation underpins that reasonable results can be obtained from practically minimal amount of projection data, gathered from just one projection direction.

A feasibility study of realizing low-dose abdominal CT using deep learning image reconstruction algorithm

OBJECTIVES

To explore the feasibility of achieving diagnostic images in low-dose abdominal CT using a Deep Learning Image Reconstruction (DLIR) algorithm.

METHODS

Prospectively enrolled 47 patients requiring contrast-enhanced abdominal CT scans. The late-arterial phase scan was added and acquired using lower-dose mode (tube current range, 175-545 mA; 80 kVp for patients with BMI ≤24 kg/m2 and 100 kVp for patients with BMI > 24 kg/m2) and reconstructed with DLIR at medium setting (DLIR-M) and high setting (DLIR-H), ASIR-V at 0% (FBP), 40% and 80% strength. Both the quantitative measurement and qualitative analysis of the five types of reconstruction methods were compared. In addition, radiation dose and image quality between the early-arterial phase ASIR-V images using standard-dose and the late-arterial phase DLIR images using low-dose were compared.

RESULTS

For the late-arterial phase, all five reconstructions had similar CT value (P > 0.05). DLIR-H, DLIR-M and ASIR-V80% images significantly reduced the image noise and improved the image contrast noise ratio, compared with the standard ASIR-V40% images (P < 0.05). ASIR-V80% images had undesirable image characteristics with obvious "waxy" artifacts, while DLIR-H images maintained high spatial resolution and had the highest subjective image quality. Compared with the early-arterial scans, the late-arterial phase scans significantly reduced the radiation dose (P < 0.05), while the DLIR-H images exhibited lower image noise and good display of the specific image details of lesions.

CONCLUSIONS

DLIR algorithm improves image quality under low-dose scan condition and may be used to reduce the radiation dose without adversely affecting the image quality.

Region-Based Segmentation and Wiener Pilot-Based Novel Amoeba Denoising Scheme for CT Imaging

Computed tomography (CT) is one of the most common and beneficial medical imaging schemes, but the associated high radiation dose injurious to the patient is always a concern. Therefore, postprocessing-based enhancement of a CT reconstructed image acquired using a reduced dose is an active research area. Amoeba- (or spatially variant kernel-) based filtering is a strong candidate scheme for postprocessing of the CT image, which adapts its shape according to the image contents. In the reported research work, the amoeba filtering is customized for postprocessing of CT images acquired at a reduced X-ray dose. The proposed scheme modifies both the pilot image formation and amoeba shaping mechanism of the conventional amoeba implementation. The proposed scheme uses a Wiener filter-based pilot image, while region-based segmentation is used for amoeba shaping instead of the conventional amoeba distance-based approach. The merits of the proposed scheme include being more suitable for CT images because of the similar region-based and symmetric nature of the human body anatomy, image smoothing without compromising on the edge details, and being adaptive in nature and more robust to noise. The performance of the proposed amoeba scheme is compared to the traditional amoeba kernel in the image denoising application for CT images using filtered back projection (FBP) on sparse-view projections. The scheme is supported by computer simulations using fan-beam projections of clinically reconstructed and simulated head CT phantoms. The scheme is tested using multiple image quality matrices, in the presence of additive projection noise. The scheme implementation significantly improves the image quality visually and statistically, providing better contrast and image smoothing without compromising on edge details. Promising results indicate the efficacy of the proposed scheme.