Applications of nonlocal means algorithm in low-dose X-ray CT image processing and reconstruction: A review

4D-Precise: Learning-based 3D motion estimation and high temporal resolution 4DCT reconstruction from treatment 2D+t X-ray projections

BACKGROUND AND OBJECTIVE

In radiotherapy treatment planning, respiration-induced motion introduces uncertainty that, if not appropriately considered, could result in dose delivery problems. 4D cone-beam computed tomography (4D-CBCT) has been developed to provide imaging guidance by reconstructing a pseudo-motion sequence of CBCT volumes through binning projection data into breathing phases. However, it suffers from artefacts and erroneously characterizes the averaged breathing motion. Furthermore, conventional 4D-CBCT can only be generated post-hoc using the full sequence of kV projections after the treatment is complete, limiting its utility. Hence, our purpose is to develop a deep-learning motion model for estimating 3D+t CT images from treatment kV projection series.

METHODS

We propose an end-to-end learning-based 3D motion modelling and 4DCT reconstruction model named 4D-Precise, abbreviated from Probabilistic reconstruction of image sequences from CBCT kV projections. The model estimates voxel-wise motion fields and simultaneously reconstructs a 3DCT volume at any arbitrary time point of the input projections by transforming a reference CT volume. Developing a Torch-DRR module, it enables end-to-end training by computing Digitally Reconstructed Radiographs (DRRs) in PyTorch. During training, DRRs with matching projection angles to the input kVs are automatically extracted from reconstructed volumes and their structural dissimilarity to inputs is penalised. We introduced a novel loss function to regulate spatio-temporal motion field variations across the CT scan, leveraging planning 4DCT for prior motion distribution estimation.

RESULTS

The model is trained patient-specifically using three kV scan series, each including over 1200 angular/temporal projections, and tested on three other scan series. Imaging data from five patients are analysed here. Also, the model is validated on a simulated paired 4DCT-DRR dataset created using the Surrogate Parametrised Respiratory Motion Modelling (SuPReMo). The results demonstrate that the reconstructed volumes by 4D-Precise closely resemble the ground-truth volumes in terms of Dice, volume similarity, mean contour distance, and Hausdorff distance, whereas 4D-Precise achieves smoother deformations and fewer negative Jacobian determinants compared to SuPReMo.

CONCLUSIONS

Unlike conventional 4DCT reconstruction techniques that ignore breath inter-cycle motion variations, the proposed model computes both intra-cycle and inter-cycle motions. It represents motion over an extended timeframe, covering several minutes of kV scan series.

Iterative reconstruction for limited-angle CT using implicit neural representation

Objective.Limited-angle computed tomography (CT) presents a challenge due to its ill-posed nature. In such scenarios, analytical reconstruction methods often exhibit severe artifacts. To tackle this inverse problem, several supervised deep learning-based approaches have been proposed. However, they are constrained by limitations such as generalization issue and the difficulty of acquiring a large amount of paired CT images.Approach.In this work, we propose an iterative neural reconstruction framework designed for limited-angle CT. By leveraging a coordinate-based neural representation, we formulate tomographic reconstruction as a convex optimization problem involving a deep neural network. We then employ differentiable projection layer to optimize this network by minimizing the discrepancy between the predicted and measured projection data. In addition, we introduce a prior-based weight initialization method to ensure the network starts optimization with an informed initial guess. This strategic initialization significantly improves the quality of iterative reconstruction by stabilizing the divergent behavior in ill-posed neural fields. Our method operates in a self-supervised manner, thereby eliminating the need for extensive data.Main results.The proposed method outperforms other iterative and learning-based methods. Experimental results on XCAT and Mayo Clinic datasets demonstrate the effectiveness of our approach in restoring anatomical features as well as structures. This finding was substantiated by visual inspections and quantitative evaluations using NRMSE, PSNR, and SSIM. Moreover, we conduct a comprehensive investigation into the divergent behavior of iterative neural reconstruction, thus revealing its suboptimal convergence when starting from scratch. In contrast, our method consistently produced accurate images by incorporating an initial estimate as informed initialization.Significance.This work showcases the feasibility to reconstruct high-fidelity CT images from limited-angle x-ray projections. The proposed methodology introduces a novel data-free approach to enhance medical imaging, holding promise across various clinical applications.

CT image denoising methods for image quality improvement and radiation dose reduction

With the ever-increasing use of computed tomography (CT), concerns about its radiation dose have become a significant public issue. To address the need for radiation dose reduction, CT denoising methods have been widely investigated and applied in low-dose CT images. Numerous noise reduction algorithms have emerged, such as iterative reconstruction and most recently, deep learning (DL)-based approaches. Given the rapid advancements in Artificial Intelligence techniques, we recognize the need for a comprehensive review that emphasizes the most recently developed methods. Hence, we have performed a thorough analysis of existing literature to provide such a review. Beyond directly comparing the performance, we focus on pivotal aspects, including model training, validation, testing, generalizability, vulnerability, and evaluation methods. This review is expected to raise awareness of the various facets involved in CT image denoising and the specific challenges in developing DL-based models.

Optimization Method to Predict Optimal Noise Reduction Parameters for the Non-Local Means Algorithm Based on the Scintillator Thickness in Radiography

The resulting image obtained from an X-ray imaging system depends significantly on the characteristics of the detector. In particular, when an X-ray image is acquired by thinning the detector, a relatively large amount of noise inevitably occurs. In addition, when a thick detector is used to reduce noise in X-ray images, blurring increases and the ability to distinguish target areas deteriorates. In this study, we aimed to derive the optimal X-ray image quality by deriving the optimal noise reduction parameters based on the non-local means (NLM) algorithm. The detectors used were of two thicknesses (96 and 140 μm), and images were acquired based on the IEC 62220-1-1:2015 RQA-5 protocol. The optimal parameters were derived by calculating the edge preservation index and signal-to-noise ratio according to the sigma value of the NLM algorithm. As a result, a sigma value of the optimized NLM algorithm (0.01) was derived, and this algorithm was applied to a relatively thin X-ray detector system to obtain appropriate noise level and spatial resolution data. The no-reference-based blind/referenceless image spatial quality evaluator value, which analyzes the overall image quality, was best when using the proposed method. In conclusion, we propose an optimized NLM algorithm based on a new method that can overcome the noise amplification problem in thin X-ray detector systems and is expected to be applied in various photon imaging fields in the future.

Affine Iterative Closest Point Algorithm Based on Color Information and Correntropy for Precise Point Set Registration

In this paper, we propose a novel affine iterative closest point algorithm based on color information and correntropy, which can effectively deal with the registration problems with a large number of noise and outliers and small deformations in RGB-D datasets. Firstly, to alleviate the problem of low registration accuracy for data with weak geometric structures, we consider introducing color features into traditional affine algorithms to establish more accurate and reliable correspondences. Secondly, we introduce the correntropy measurement to overcome the influence of a large amount of noise and outliers in the RGB-D datasets, thereby further improving the registration accuracy. Experimental results demonstrate that the proposed registration algorithm has higher registration accuracy, with error reduction of almost 10 times, and achieves more stable robustness than other advanced algorithms.

An optimized EBRSA-Bi LSTM model for highly undersampled rapid CT image reconstruction

COVID-19 has spread all over the world, causing serious panic around the globe. Chest computed tomography (CT) images are integral in confirming COVID positive patients. Several investigations were conducted to improve or maintain the image reconstruction quality for the sample image reconstruction. Deep learning (DL) methods have recently been proposed to achieve fast reconstruction, but many have focused on a single domain, such as the image domain of k-space. In this research, the highly under-sampled enhanced battle royale self-attention based bi-directional long short-term (EBRSA-bi LSTM) CT image reconstruction model is proposed to reconstruct the image from the under-sampled data. The research is adapted with two phases, namely, pre-processing and reconstruction. The extended cascaded filter (ECF) is proposed for image pre-processing and tends to suppress the noise and enhance the reconstruction accuracy. In the reconstruction model, the battle royale optimization (BrO) is intended to diminish the loss function of the reconstruction network model and weight updation. The proposed model is tested with two datasets, COVID-CT- and SARS-CoV-2 CT. The reconstruction accuracy of the proposed model with two datasets is 93.5 % and 97.7 %, respectively. Also, the image quality assessment parameters such as Peak-Signal to Noise Ratio (PSNR), Root Mean Square Error (RMSE) and Structural Similarity Index metric (SSIM) are evaluated, and it yields an outcome of (45 and 46 dB), (0.0026 and 0.0022) and (0.992, 0.996) with two datasets.

An efficient deep neural network to segment lung nodule using optimized HDCCARUNet model

Lung cancer is a severe disease that may lead to death if left undiagnosed and untreated. Lung cancer recognition and segmentation is a difficult task in medical image processing. The study of Computed Tomography (CT) is an important phase for detecting abnormal tissues in the lung. The size of a nodule as well as the fine details of nodule can be varied for various images. Radiologists face a difficult task in diagnosing nodules from multiple images. Deep learning approaches outperform traditional learning algorithms when the data amount is large. One of the most common deep learning architectures is convolutional neural networks. Convolutional Neural Networks use pre-trained models like LeNet, AlexNet, GoogleNet, VGG16, VGG19, Resnet50, and others for learning features. This study proposes an optimized HDCCARUNet (Hybrid Dilated Convolutional Channel Attention Res-UNet) architecture, which combines an improved U-Net with a modified channel attention (MCA) block, and a HDAC (hybrid dilated attention convolutional) layer to accurately and effectively do medical image segmentation for various tasks. The attention mechanism aids in focusing on the desired outcome. The ability to dynamically allot input weights to neurons allows it to focus only on the most important information. In order to gather key details about different object features and infer a finer channel-wise attention, the proposed system uses a modified channel attention (MCA) block. The experiment is conducted on LIDC-IDRI dataset. The noises present in the dataset images are denoised by enhanced DWT filter and the performance is analysed at various noise levels. The proposed method achieves an accuracy rate of 99.58 % . Performance measures like accuracy, sensitivity, specificity, and ROC curves are evaluated and the system significantly outperforms other state-of-the-art systems.

Synthetic Data Generation for Automatic Segmentation of X-ray Computed Tomography Reconstructions of Complex Microstructures

The greatest challenge when using deep convolutional neural networks (DCNNs) for automatic segmentation of microstructural X-ray computed tomography (XCT) data is the acquisition of sufficient and relevant data to train the working network. Traditionally, these have been attained by manually annotating a few slices for 2D DCNNs. However, complex multiphase microstructures would presumably be better segmented with 3D networks. However, manual segmentation labeling for 3D problems is prohibitive. In this work, we introduce a method for generating synthetic XCT data for a challenging six-phase Al-Si alloy composite reinforced with ceramic fibers and particles. Moreover, we propose certain data augmentations (brightness, contrast, noise, and blur), a special in-house designed deep convolutional neural network (Triple UNet), and a multi-view forwarding strategy to promote generalized learning from synthetic data and therefore achieve successful segmentations. We obtain an overall Dice score of 0.77. Lastly, we prove the detrimental effects of artifacts in the XCT data on achieving accurate segmentations when synthetic data are employed for training the DCNNs. The methods presented in this work are applicable to other materials and imaging techniques as well. Successful segmentation coupled with neural networks trained with synthetic data will accelerate scientific output.

Superiorization-inspired unrolled SART algorithm with U-Net generated perturbations for sparse-view and limited-angle CT reconstruction

Objective.Unrolled algorithms are a promising approach for reconstruction of CT images in challenging scenarios, such as low-dose, sparse-view and limited-angle imaging. In an unrolled algorithm, a fixed number of iterations of a reconstruction method are unrolled into multiple layers of a neural network, and interspersed with trainable layers. The entire network is then trained end-to-end in a supervised fashion, to learn an appropriate regularizer from training data. In this paper we propose a novel unrolled algorithm, and compare its performance with several other approaches on sparse-view and limited-angle CT.Approach.The proposed algorithm is inspired by the superiorization methodology, an optimization heuristic in which iterates of a feasibility-seeking method are perturbed between iterations, typically using descent directions of a model-based penalty function. Our algorithm instead uses a modified U-net architecture to introduce the perturbations, allowing a network to learn beneficial perturbations to the image at various stages of the reconstruction, based on the training data.Main Results.In several numerical experiments modeling sparse-view and limited angle CT scenarios, the algorithm provides excellent results. In particular, it outperforms several competing unrolled methods in limited-angle scenarios, while providing comparable or better performance on sparse-view scenarios.Significance.This work represents a first step towards exploiting the power of deep learning within the superiorization methodology. Additionally, it studies the effect of network architecture on the performance of unrolled methods, as well as the effectiveness of the unrolled approach on both limited-angle CT, where previous studies have primarily focused on the sparse-view and low-dose cases.

Automated development of the contrast-detail curve based on statistical low-contrast detectability in CT images

Purpose

We have developed a software to automatically find the contrast–detail (C–D) curve based on the statistical low‐contrast detectability (LCD) in images of computed tomography (CT) phantoms at multiple cell sizes and to generate minimum detectable contrast (MDC) characteristics.

Methods

A simple graphical user interface was developed to set the initial parameters needed to create multiple grid region of interest of various cell sizes with a 2‐pixel increment. For each cell in the grid, the average CT number was calculated to obtain the standard deviation (SD). Detectability was then calculated by multiplying the SD of the mean CT numbers by 3.29. This process was automatically repeated as many times as the cell size was set at initialization. Based on the obtained LCD, the C–D curve was obtained and the target size at an MDC of 0.6% (i.e., 6‐HU difference) was determined. We subsequently investigated the consistency of the target sizes for a 0.6% MDC at four locations within the homogeneous image. We applied the software to images with six noise levels, images of two modules of the American College of Radiology CT phantom, images of four different phantoms, and images of four different CT scanners. We compared the target sizes at a 0.6% MDC based on the statistical LCD and the results from a human observer.

Results

The developed system was able to measure C–D curves from different phantoms and scanners. We found that the C–D curves follow a power‐law fit. We found that higher noise levels resulted in a higher MDC for a target of the same size. The low‐contrast module image had a slightly higher MDC than the distance module image. The minimum size of an object detected by visual observation was slightly larger than the size using statistical LCD.

Conclusions

The statistical LCD measurement method can generate a C–D curve automatically, quickly, and objectively.

Improved Ultrafast Power Doppler Imaging by Using Spatiotemporal Non-Local Means Filtering

The change of microvasculature is associated with the occurrence and development of many diseases. Ultrafast power Doppler imaging (uPDI) is an emerging technology for the visualization of microvessels due to the development of ultrafast plane wave (PW) imaging and advanced clutter filters. However, the low signal-to-noise ratio (SNR) caused by unfocused transmit of PW imaging deteriorates the subsequent imaging of microvasculature. Nonlocal means (NLM) filtering has been demonstrated to be effective in the denoising of both natural and medical images, including ultrasound power Doppler images. However, the feasibility and performance of applying an NLM filter on the ultrasound radio frequency (RF) data have not been investigated so far. In this study, we propose to apply an NLM filter on the spatiotemporal domain of clutter filtered blood flow RF data (St-NLM) to improve the quality of uPDI. Experiments were conducted to compare the proposed method with three different methods (under various similarity window sizes), including conventional uPDI without NLM filtering (Non-NLM), NLM filtering on the obtained power Doppler images (PD-NLM), and NLM filtering on the spatial domain of clutter filtered blood flow RF data (S-NLM). Phantom experiments, in vivo contrast-enhanced human spinal cord tumor experiments, and in vivo contrast-free human liver experiments were performed to demonstrate the superiority of the proposed St-NLM method over the other three methods. Qualitative and quantitative results show that the proposed St-NLM method can effectively suppress the background noise, improve the contrast between vessels and background, and preserve the details of small vessels at the same time. In the human liver study, the proposed St-NLM method achieves 31.05-, 24.49-, and 11.15-dB higher contrast-to-noise ratios (CNRs) and 36.86-, 36.86-, and 15.22-dB lower noise powers than Non-NLM, PD-NLM, and S-NLM, respectively. In the human spinal cord tumor, the full-width at half-maximums (FWHMs) of vessel cross Section are 76, 201, and [Formula: see text] for St-NLM, Non-NLM, and S-NLM, respectively. The proposed St-NLM method can enhance the microvascular visualization in uPDI and has the potential for the diagnosis of many microvessel-change-related diseases.

Low-Dose CT Image Post-Processing Based on Learn-Type Sparse Transform

As a detection method, X-ray Computed Tomography (CT) technology has the advantages of clear imaging, short detection time, and low detection cost. This makes it more widely used in clinical disease screening, detection, and disease tracking. This study exploits the ability of sparse representation to learn sparse transformations of information and combines it with image decomposition theory. The structural information of low-dose CT images is separated from noise and artifact information, and the sparse expression of sparse transformation is used to improve the imaging effect. In this paper, two different learned sparse transformations are used. The first covers more organizational information about the scanned object. The other can cover more noise artifacts. Both methods can improve the ability to learn sparse transformations to express various image information. Experimental results show that the algorithm is effective.

Artificial Intelligence in Radiation Therapy

Artificial intelligence (AI) has great potential to transform the clinical workflow of radiotherapy. Since the introduction of deep neural networks, many AI-based methods have been proposed to address challenges in different aspects of radiotherapy. Commercial vendors have started to release AI-based tools that can be readily integrated to the established clinical workflow. To show the recent progress in AI-aided radiotherapy, we have reviewed AI-based studies in five major aspects of radiotherapy including image reconstruction, image registration, image segmentation, image synthesis, and automatic treatment planning. In each section, we summarized and categorized the recently published methods, followed by a discussion of the challenges, concerns, and future development. Given the rapid development of AI-aided radiotherapy, the efficiency and effectiveness of radiotherapy in the future could be substantially improved through intelligent automation of various aspects of radiotherapy.

Intelligent Algorithm-Based Multislice Spiral Computed Tomography to Diagnose Coronary Heart Disease

In this study, dictionary learning and expectation maximization reconstruction (DLEM) was combined to denoise 64-slice spiral CT images, and results of coronary angiography (CAG) were used as standard to evaluate its clinical value in diagnosing coronary artery diseases. 120 patients with coronary heart disease (CHD) confirmed by CAG examination were retrospectively selected as the research subjects. According to the random number table method, the patients were divided into two groups: the control group was diagnosed by conventional 64-slice spiral CT images, and the observation group was diagnosed by 64-slice spiral CT images based on the DLEM algorithm, with 60 cases in both groups. With CAG examination results as the standard, the diagnostic effects of the two CT examination methods were compared. The results showed that when the number of iterations of maximum likelihood expectation maximization (MLEM) algorithm reached 50, the root mean square error (RMSE) and peak signal to noise ratio (PSNR) values were similar to the results obtained by the DLEM algorithm under a number of iterations of 10 when the RMSE and PSNR values were 18.9121 dB and 74.9911 dB, respectively. In the observation group, 28.33% (17/60) images were of grade 4 or above before processing; after processing, it was 70% (42/60), significantly higher than the proportion of high image quality before processing. The overall diagnostic consistency, sensitivity, specificity, and accuracy (88.33%, 86.67%, 80%, and 85%) of the observation group were better than those in the control group (60.46%, 62.5%, 58.33%, and 61.66%). In conclusion, the DLEM algorithm has good denoising effect on 64-slice spiral CT images, which significantly improves the accuracy in the diagnosis of coronary artery stenosis and has good clinical diagnostic value and is worth promoting.

Total variation combining nonlocal means filtration for image reconstruction in X-ray computed tomography

BACKGROUND

Image reconstruction for realistic medical images under incomplete observation is still one of the core tasks for computed tomography (CT). However, the stair-case artifacts of Total variation (TV) based ones have restricted the usage of the reconstructed images.

OBJECTIVE

This work aims to propose and test an accurate and efficient algorithm to improve reconstruction quality under the idea of synergy between local and nonlocal regularizations.

METHODS

The total variation combining the nonlocal means filtration is proposed and the alternating direction method of multipliers is utilized to develop an efficient algorithm. The first order approximation of linear expansion at intermediate point is applied to overcome the computation of the huge CT system matrix.

RESULTS

The proposed method improves root mean squared error by 25.6% compared to the recent block-matching sparsity regularization (BMSR) on simulation dataset of 19 views. The structure similarities of image of the new method is higher than 0.95, while that of BMSR is about 0.92. Moreover, on real rabbit dataset of 20 views, the peak signal-to-noise ratio (PSNR) of the new method is 36.84, while using other methods PSNR are lower than 35.81.

CONCLUSIONS

The proposed method shows advantages on noise suppression and detail preservations over the competing algorithms used in CT image reconstruction.

MAGIC: Manifold and Graph Integrative Convolutional Network for Low-Dose CT Reconstruction

Low-dose computed tomography (LDCT) scans, which can effectively alleviate the radiation problem, will degrade the imaging quality. In this paper, we propose a novel LDCT reconstruction network that unrolls the iterative scheme and performs in both image and manifold spaces. Because patch manifolds of medical images have low-dimensional structures, we can build graphs from the manifolds. Then, we simultaneously leverage the spatial convolution to extract the local pixel-level features from the images and incorporate the graph convolution to analyze the nonlocal topological features in manifold space. The experiments show that our proposed method outperforms both the quantitative and qualitative aspects of state-of-the-art methods. In addition, aided by a projection loss component, our proposed method also demonstrates superior performance for semi-supervised learning. The network can remove most noise while maintaining the details of only 10% (40 slices) of the training data labeled.

Noise Conscious Training of Non Local Neural Network Powered by Self Attentive Spectral Normalized Markovian Patch GAN for Low Dose CT Denoising

The explosive rise of the use of Computer tomography (CT) imaging in medical practice has heightened public concern over the patient's associated radiation dose. On the other hand, reducing the radiation dose leads to increased noise and artifacts, which adversely degrades the scan's interpretability. In recent times, the deep learning-based technique has emerged as a promising method for low dose CT(LDCT) denoising. However, some common bottleneck still exists, which hinders deep learning-based techniques from furnishing the best performance. In this study, we attempted to mitigate these problems with three novel accretions. First, we propose a novel convolutional module as the first attempt to utilize neighborhood similarity of CT images for denoising tasks. Our proposed module assisted in boosting the denoising by a significant margin. Next, we moved towards the problem of non-stationarity of CT noise and introduced a new noise aware mean square error loss for LDCT denoising. The loss mentioned above also assisted to alleviate the laborious effort required while training CT denoising network using image patches. Lastly, we propose a novel discriminator function for CT denoising tasks. The conventional vanilla discriminator tends to overlook the fine structural details and focus on the global agreement. Our proposed discriminator leverage self-attention and pixel-wise GANs for restoring the diagnostic quality of LDCT images. Our method validated on a publicly available dataset of the 2016 NIH-AAPM-Mayo Clinic Low Dose CT Grand Challenge performed remarkably better than the existing state of the art method. The corresponding source code is available at: https://github.com/reach2sbera/ldct_nonlocal.

Computed Tomography Image Texture under Feature Extraction Algorithm in the Diagnosis of Effect of Specific Nursing Intervention on Mycoplasma Pneumonia in Children

To evaluate the effect of specific nursing intervention in children with mycoplasma pneumonia (MP), a feature extraction algorithm based on gray level co-occurrence matrix (GLCM) was proposed and combined with computed tomography (CT) image texture features. Then, 98 children with MP were rolled into the observation group with 49 cases (specific nursing) and the control group with 49 cases (routine nursing). CT images based on feature extraction algorithm of optimized GLCM were used to examine the children before and after nursing intervention, and the recovery of the two groups of children was discussed. The results showed that the proportion of lung texture increase, rope shadow, ground glass shadow, atelectasis, and pleural effusion in the observation group (24.11%, 3.86%, 8.53%, 15.03%, and 3.74%) was significantly lower than that in the control group (28.53%, 10.23%, 13.34%, 21.15%, and 8.13%) after nursing (P < 0.05). There were no significant differences in the proportion of small patchy shadows, large patchy consolidation shadows, and bronchiectasis between the observation group and the control group (P > 0.05). In the course of nursing intervention, in the observation group, the disappearance time of cough, normal temperature, disappearance time of lung rales, and absorption time of lung shadow (2.15 ± 0.86 days, 4.81 ± 1.14 days, 3.64 ± 0.55 days, and 5.96 ± 0.62 days) were significantly shorter than those in the control group (2.87 ± 0.95 days, 3.95 ± 1.06 days, 4.51 ± 1.02 days, and 8.14 ± 1.35 days) (P < 0.05). After nursing intervention, the proportion of satisfaction and total satisfaction in the experimental group (67.08% and 28.66%) was significantly higher than that in the control group (40.21% and 47.39%), while the proportion of dissatisfaction (4.26%) was significantly lower than that in the control group (12.4%) (P < 0.05). To sum up, specific nursing intervention was more beneficial to improve the progress of characterization recovery and the overall recovery effect of children with MP relative to conventional nursing. CT image based on feature extraction algorithm of optimized GLCM was of good adoption value in the diagnosis and treatment of MP in children.

Implementing a non-local means method to CTA data of aortic dissection

It is necessary to conserve important information, like edges, details, and textures, in CT aortic dissection images, as this helps the radiologist examine and diagnose the disease. Hence, a less noisy image is required to support medical experts in performing better diagnoses. In this work, the non-local means (NLM) method is conducted to minimize the noise in CT images of aortic dissection patients as a preprocessing step to produce accurate aortic segmentation results. The method is implemented in an existing segmentation system using six different kernel functions, and the evaluation is done by assessing DSC, precision, and recall of segmentation results. Furthermore, the visual quality of denoised images is also taken into account to be determined. Besides, a comparative analysis between NLM and other denoising methods is done in this experiment. The results showed that NLM yields encouraging segmentation results, even though the visualization of denoised images is unacceptable. Applying the NLM algorithm with the flat function provides the highest DSC, precision, and recall values of 0.937101, 0.954835, and 0.920517 consecutively.

General simultaneous motion estimation and image reconstruction (G-SMEIR)

To achieve better performance for 4D multi-frame reconstruction with the parametric motion model (MF-PMM), a general simultaneous motion estimation and image reconstruction (G-SMEIR) method is proposed. In G-SMEIR, projection domain motion estimation and image domain motion estimation are performed alternatively to achieve better 4D reconstruction. This method can mitigate the local optimum trapping problem in either domain. To improve computational efficiency, the image domain motion estimation is accelerated by adapting fast convergent algorithms and graphics processing unit (GPU) computing. The proposed G-SMEIR method is tested using a cone-beam computed tomography (CBCT) simulation study of 4D XCAT phantom at different dose levels and compared with 3D total variation-based reconstruction (3D TV), 4D reconstruction with image domain motion estimation (IM4D), and SMEIR. G-SMEIR shows strong denoising capability and achieves similar performance at regular dose and half dose. The root mean squared error (RMSE) of G-SMEIR is the best among the four methods and improved about 12% over SMEIR for all respiratory phase images at full dose. G-SMEIR also achieved the best structural similarity index (SSIM) values among all methods. More importantly, G-SMEIR leads to more than 40% improvement of the mean deviation from the phantom tumor motion over SMEIR. A preliminary patient CBCT image reconstruction also shows better image quality of G-SMEIR than that of the frame-by-frame reconstruction (3D TV) and MF-PMM either using image domain motion estimation (IM4D) or using projection domain motion estimation (SMEIR) alone. G-SMEIR with a flexible combination of image domain and projection domain motion estimation provides an effective tool for 4D tomographic reconstruction.

Prior-image-based CT reconstruction using attenuation-mismatched priors

Prior-image-based reconstruction (PIBR) methods are powerful tools for reducing radiation doses and improving the image quality of low-dose computed tomography (CT). Apart from anatomical changes, prior and current images can also have different attenuations because they originated from different scanners or from the same scanner but with different x-ray beam qualities (e.g., kVp settings, beam filters) during data acquisition. In such scenarios, with attenuation-mismatched priors, PIBR is challenging. In this work, we investigate a specific PIBR method, called statistical image reconstruction, using normal-dose image-induced nonlocal means regularization (SIR-ndiNLM), to address PIBR with such attenuation-mismatched priors and achieve quantitative low-dose CT imaging. We propose two corrective schemes for the original SIR-ndiNLM method, (1) a global histogram-matching approach and (2) a local attenuation correction approach, to account for the attenuation differences between the prior and current images in PIBR. We validate the efficacy of the proposed schemes using images acquired from dual-energy CT scanners to simulate attenuation mismatches. Meanwhile, we utilize different CT slices to simulate anatomical mismatches or changes between the prior and the current low-dose image. We observe that the original SIR-ndiNLM introduces artifacts to the reconstruction when an attenuation-mismatched prior is used. Furthermore, we find that a larger attenuation mismatch between the prior and current images results in more severe artifacts in the SIR-ndiNLM reconstruction. Our two proposed corrective schemes enable SIR-ndiNLM to effectively handle the attenuation mismatch and anatomical changes between the two images and successfully eliminate the artifacts. We demonstrate that the proposed techniques permit SIR-ndiNLM to leverage the attenuation-mismatched prior and achieve quantitative low-dose CT reconstruction from both low-flux and sparse-view data acquisitions. This work permits robust and reliable PIBR for CT data acquired using different beam settings.

Computed Tomography- (CT-) Based Virtual Surgery Planning for Spinal Intervertebral Foraminal Assisted Clinical Treatment

With the development of minimally invasive spine concepts and the introduction of new minimally invasive instruments, minimally invasive spine technology, represented by foraminoscopy, has flourished, and percutaneous foraminoscopy has become one of the most reliable minimally invasive procedures for the treatment of lumbar disc herniation. Percutaneous foraminoscopy is a safe and effective minimally invasive spinal endoscopic surgical technique. It fully protects the paravertebral muscles and soft tissues as well as the posterior column structure of the spine, provides precise treatment of the target nucleus pulposus tissue, with the advantages of less surgical trauma, fewer postoperative complications, and rapid postoperative recovery, and is widely promoted and used in clinical practice. In this paper, we can view the location, morphology, structure, alignment, and adjacency relationships by performing coronary, CT, and diagonal reconstruction along the attachment of the yellow ligaments and performing 3D reconstruction or processing techniques after performing CT scans. This allows clinicians to observe the laminoplasty and the stenosis of the vertebral canal in a more intuitive and overall manner. It has clinical significance for the display of the sublaminar spine as well as the physician's judgment of the disease and the choice of surgery.

OMEGA-open-source emission tomography software

In this paper we present OMEGA, an open-source software, for efficient and fast image reconstruction in positron emission tomography (PET). OMEGA uses the scripting language of MATLAB and GNU Octave allowing reconstruction of PET data with a MATLAB or GNU Octave interface. The goal of OMEGA is to allow easy and fast reconstruction of any PET data, and to provide a computationally efficient, easy-access platform for development of new PET algorithms with built-in forward and backward projection operations available to the user as a MATLAB/Octave class. OMEGA also includes direct support for GATE simulated data, facilitating easy evaluation of the new algorithms using Monte Carlo simulated PET data. OMEGA supports parallel computing by utilizing OpenMP for CPU implementations and OpenCL for GPU allowing any hardware to be used. OMEGA includes built-in function for the computation of normalization correction and allows several other corrections to be applied such as attenuation, randoms or scatter. OMEGA includes several different maximum-likelihood and maximum a posteriori (MAP) algorithms with several different priors. The user can also input their own priors to the built-in MAP functions. The image reconstruction in OMEGA can be computed either by using an explicitly computed system matrix or with a matrix-free formalism, where the latter can be accelerated with OpenCL. We provide an overview on the software and present some examples utilizing the different features of the software.

Parameter-Transferred Wasserstein Generative Adversarial Network (PT-WGAN) for Low-Dose PET Image Denoising

Due to the widespread use of positron emission tomography (PET) in clinical practice, the potential risk of PET-associated radiation dose to patients needs to be minimized. However, with the reduction in the radiation dose, the resultant images may suffer from noise and artifacts that compromise diagnostic performance. In this paper, we propose a parameter-transferred Wasserstein generative adversarial network (PT-WGAN) for low-dose PET image denoising. The contributions of this paper are twofold: i) a PT-WGAN framework is designed to denoise low-dose PET images without compromising structural details, and ii) a task-specific initialization based on transfer learning is developed to train PT-WGAN using trainable parameters transferred from a pretrained model, which significantly improves the training efficiency of PT-WGAN. The experimental results on clinical data show that the proposed network can suppress image noise more effectively while preserving better image fidelity than recently published state-of-the-art methods. We make our code available at https://github.com/90n9-yu/PT-WGAN.

Anisotropic structure property based image reconstruction method for limited-angle computed tomography

Limited-angle computed tomography (CT) may appear in restricted CT scans. Since the available projection data is incomplete, the images reconstructed by filtered back-projection (FBP) or algebraic reconstruction technique (ART) often encounter shading artifacts. However, using the anisotropy property of the shading artifacts that coincide with the characteristic of limited-angle CT images can reduce the shading artifacts. Considering this concept, we combine the anisotropy property of the shading artifacts with the anisotropic structure property of an image to develop a new algorithm for image reconstruction. Specifically, we propose an image reconstruction method based on adaptive weighted anisotropic total variation (AwATV). This method, termed as AwATV method for short, is designed to preserve image structures and then remove the shading artifacts. It characterizes both of above properties. The anisotropy property of the shading artifacts accounts for reducing artifacts, and the anisotropic structure property of an image accounts for preserving structures. In order to evaluate the performance of AwATV, we use the simulation projection data of FORBILD head phantom and real CT data for image reconstruction. Experimental results show that AwATV can always reconstruct images with higher SSIM and PSNR, and smaller RMSE, which means that AwATV enables to reconstruct images with higher quality in term of artifact reduction and structure preservation.

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

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

Qualitative and Quantitative Assessment of Nonlocal Means Reconstruction Algorithm in a Flexible PET Scanner

OBJECTIVE. Flexible PET (fxPET) was designed to fit existing MRI systems. The newly modified nonlocal means (NLM) algorithm is combined with the 3D dynamic row-action maximum likelihood algorithm (DRAMA). We investigated qualitative and quantitative acceptability of fxPET images reconstructed by modified NLM compared with whole-body (WB) PET/CT images and conventional 3D DRAMA reconstruction alone. MATERIALS AND METHODS. Fifty-nine patients with known or suspected malignancies underwent WB PET/CT scanning approximately 1 hour after the injection of ¹⁸F-FDG, after which they underwent fxPET scanning. Two readers rated the quality of fxPET images by consensus. Detection rate (the proportion of lesions found on PET), maximal standardized uptake value (SUV_max), metabolic tumor volume (MTV), total lesion glycolysis (TLG), tumor-to-normal liver ratio (TNR), and background liver signal-to-noise ratio (SNR) were compared among the three datasets. RESULTS. Higher image quality was obtained by modified NLM reconstruction than by conventional reconstruction without statistical significance. The detection rate was comparable among three datasets. SUV_max was significantly higher, and MTV and TLG were significantly lower in the modified NLM dataset (p < 0.002) than in the other two datasets, with significantly positive correlations (p < 0.001; Spearman rank correlation coefficient, 0.87-0.99). The TNRs in modified NLM images were significantly larger than in the other datasets (p < 0.05). The background SNRs in modified NLM images were comparable with those in WB PET/CT images, and significantly higher than in the conventional fxPET images (p < 0.005). CONCLUSION. The modified NLM algorithm was clinically acceptable, yielding higher TNR and background SNR compared with conventional reconstruction. Image quality and the lesion detection rate were comparable in this population.

Noise Reduction in CT Images Using a Selective Mean Filter

Background

Noise reduction is a method for reducing CT dose; however, it can reduce image quality.

Objective

This study aims to propose a selective mean filter (SMF) and evaluate its effectiveness for noise suppression in CT images.

Material and Methods

This experimental study proposed and implemented the new noise reduction algorithm. The proposed algorithm is based on a mean filter (MF), but the calculation of the mean pixel value using the neighboring pixels in a kernel selectively applied a threshold value based on the noise of the image. The SMF method was evaluated using images of phantoms. The dose reduction was estimated by comparing the image noise acquired with a lower dose after implementing the SMF method and the noise in the original image acquired with a higher dose. For comparison, the images were also filtered with an adaptive mean filter (AMF) and a bilateral filter (BF).

Results

The spatial resolution of the image filtered with the SMF was similar to the original images and the images filtered with the BF. While using the AMF, spatial resolution was significantly corrupted. The noise reduction achieved using the SMF was up to 75%, while it was up to 50% using the BF.

Conclusion

SMF significantly reduces the noise and preserves the spatial resolution of the image. The noise reduction was more pronounced with BF, and less pronounced with AMF.

VVBP-Tensor in the FBP Algorithm: Its Properties and Application in Low-Dose CT Reconstruction

For decades, commercial X-ray computed tomography (CT) scanners have been using the filtered backprojection (FBP) algorithm for image reconstruction. However, the desire for lower radiation doses has pushed the FBP algorithm to its limit. Previous studies have made significant efforts to improve the results of FBP through preprocessing the sinogram, modifying the ramp filter, or postprocessing the reconstructed images. In this paper, we focus on analyzing and processing the stacked view-by-view backprojections (named VVBP-Tensor) in the FBP algorithm. A key challenge for our analysis lies in the radial structures in each backprojection slice. To overcome this difficulty, a sorting operation was introduced to the VVBP-Tensor in its z direction (the direction of the projection views). The results show that, after sorting, the tensor contains structures that are similar to those of the object, and structures in different slices of the tensor are correlated. We then analyzed the properties of the VVBP-Tensor, including structural self-similarity, tensor sparsity, and noise statistics. Considering these properties, we have developed an algorithm using the tensor singular value decomposition (named VVBP-tSVD) to denoise the VVBP-Tensor for low-mAs CT imaging. Experiments were conducted using a physical phantom and clinical patient data with different mAs levels. The results demonstrate that the VVBP-tSVD is superior to all competing methods under different reconstruction schemes, including sinogram preprocessing, image postprocessing, and iterative reconstruction. We conclude that the VVBP-Tensor is a suitable processing target for improving the quality of FBP reconstruction, and the proposed VVBP-tSVD is an effective algorithm for noise reduction in low-mAs CT imaging. This preliminary work might provide a heuristic perspective for reviewing and rethinking the FBP algorithm.

SPULTRA: Low-Dose CT Image Reconstruction With Joint Statistical and Learned Image Models

Low-dose CT image reconstruction has been a popular research topic in recent years. A typical reconstruction method based on post-log measurements is called penalized weighted-least squares (PWLS). Due to the underlying limitations of the post-log statistical model, the PWLS reconstruction quality is often degraded in low-dose scans. This paper investigates a shifted-Poisson (SP) model based likelihood function that uses the pre-log raw measurements that better represents the measurement statistics, together with a data-driven regularizer exploiting a Union of Learned TRAnsforms (SPULTRA). Both the SP induced data-fidelity term and the regularizer in the proposed framework are nonconvex. The proposed SPULTRA algorithm uses quadratic surrogate functions for the SP induced data-fidelity term. Each iteration involves a quadratic subproblem for updating the image, and a sparse coding and clustering subproblem that has a closed-form solution. The SPULTRA algorithm has a similar computational cost per iteration as its recent counterpart PWLS-ULTRA that uses post-log measurements, and it provides better image reconstruction quality than PWLS-ULTRA, especially in low-dose scans.

Advanced 4-dimensional cone-beam computed tomography reconstruction by combining motion estimation, motion-compensated reconstruction, biomechanical modeling and deep learning

4-Dimensional cone-beam computed tomography (4D-CBCT) offers several key advantages over conventional 3D-CBCT in moving target localization/delineation, structure de-blurring, target motion tracking, treatment dose accumulation and adaptive radiation therapy. However, the use of the 4D-CBCT in current radiation therapy practices has been limited, mostly due to its sub-optimal image quality from limited angular sampling of cone-beam projections. In this study, we summarized the recent developments of 4D-CBCT reconstruction techniques for image quality improvement, and introduced our developments of a new 4D-CBCT reconstruction technique which features simultaneous motion estimation and image reconstruction (SMEIR). Based on the original SMEIR scheme, biomechanical modeling-guided SMEIR (SMEIR-Bio) was introduced to further improve the reconstruction accuracy of fine details in lung 4D-CBCTs. To improve the efficiency of reconstruction, we recently developed a U-net-based deformation-vector-field (DVF) optimization technique to leverage a population-based deep learning scheme to improve the accuracy of intra-lung DVFs (SMEIR-Unet), without explicit biomechanical modeling. Details of each of the SMEIR, SMEIR-Bio and SMEIR-Unet techniques were included in this study, along with the corresponding results comparing the reconstruction accuracy in terms of CBCT images and the DVFs. We also discussed the application prospects of the SMEIR-type techniques in image-guided radiation therapy and adaptive radiation therapy, and presented potential schemes on future developments to achieve faster and more accurate 4D-CBCT imaging.

New nonlocal forward model for diffuse optical tomography

The forward model in diffuse optical tomography (DOT) describes how light propagates through a turbid medium. It is often approximated by a diffusion equation (DE) that is numerically discretized by the classical finite element method (FEM). We propose a nonlocal diffusion equation (NDE) as a new forward model for DOT, the discretization of which is carried out with an efficient graph-based numerical method (GNM). To quantitatively evaluate the new forward model, we first conduct experiments on a homogeneous slab, where the numerical accuracy of both NDE and DE is compared against the existing analytical solution. We further evaluate NDE by comparing its image reconstruction performance (inverse problem) to that of DE. Our experiments show that NDE is quantitatively comparable to DE and is up to 64% faster due to the efficient graph-based representation that can be implemented identically for geometries in different dimensions. The proposed discretization method can be easily applied to other imaging techniques like diffuse correlation spectroscopy which are normally modeled by the diffusion equation.

poly-DART: A discrete algebraic reconstruction technique for polychromatic X-ray CT

The discrete algebraic reconstruction technique (DART) is a tomographic method to reconstruct images from X-ray projections in which prior knowledge on the number of object materials is exploited. In monochromatic X-ray CT (e.g., synchrotron), DART has been shown to lead to high-quality reconstructions, even with a low number of projections or a limited scanning view. However, most X-ray sources are polychromatic, leading to beam hardening effects, which significantly degrade the performance of DART. In this work, we propose a new discrete tomography algorithm, poly-DART, that exploits sparsity in the attenuation values using DART and simultaneously accounts for the polychromatic nature of the X-ray source. The results show that poly-DART leads to a vastly improved segmentation on polychromatic data obtained from Monte Carlo simulations as well as on experimental data, compared to DART.

Feasibility of fast non local means filter in pediatric chest x-ray for increasing of pulmonary nodule detectability with 3D printed lung nodule phantom

General x-ray images have a lower probability of nodule detection than other modalities. Especially in children, the probability of nodule detection can likely drop due to poor image quality from using low radiation dose. To demonstrate the effectiveness of fast non-local means (FNLM) filter to increase the probability of nodule detection in pediatric chest x-ray images and reduce radiation dose while maintaining image quality. Quantitative assessment of normalised noise power spectrum (NNPS), coefficient of variation (COV) and contrast to noise ratio (CNR) were performed after applying four filters (median, Wiener, total variation and FNLM) on a 1-year-old child phantom. A 3D-printed patient nodule phantom was inserted into the phantom. Assessment was performed on AP and LAT view images acquired with the tube voltage reduced to 38 and 27%, and tube current reduced to 84 and 61%, respectively. The results showed the lowest NNPS and COV values and the highest CNR value when the FNLM filter applied. Moreover, the AP view results showed 37% decrease in COV and 30% increase in CNR in images with the FNLM filter applied (images exposed with the tube voltage and current reduced to 29% and 50%, respectively). The LAT view results showed 5% decrease in COV and 36% increase in CNR in images with the FNLM filter applied (images exposed with the tube current reduced by 27%). By applying the FNLM filter, the probability of nodule detection could be increased by denoising and contrast enhancement. Moreover, using the FNLM filter could reduce cancer risk in pediatric patients by reducing radiation dose about 30% to 44%.

Learning-based CBCT correction using alternating random forest based on auto-context model

PURPOSE

Quantitative Cone Beam CT (CBCT) imaging is increasing in demand for precise image-guided radiotherapy because it provides a foundation for advanced image-guided techniques, including accurate treatment setup, online tumor delineation, and patient dose calculation. However, CBCT is currently limited only to patient setup in the clinic because of the severe issues in its image quality. In this study, we develop a learning-based approach to improve CBCT's image quality for extended clinical applications.

MATERIALS AND METHODS

An auto-context model is integrated into a machine learning framework to iteratively generate corrected CBCT (CCBCT) with high-image quality. The first step is data preprocessing for the built training dataset, in which uninformative image regions are removed, noise is reduced, and CT and CBCT images are aligned. After a CBCT image is divided into a set of patches, the most informative and salient anatomical features are extracted to train random forests. Within each patch, alternating RF is applied to create a CCBCT patch as the output. Moreover, an iterative refinement strategy is exercised to enhance the image quality of CCBCT. Then, all the CCBCT patches are integrated to reconstruct final CCBCT images.

RESULTS

The learning-based CBCT correction algorithm was evaluated using the leave-one-out cross-validation method applied on a cohort of 12 patients' brain data and 14 patients' pelvis data. The mean absolute error (MAE), peak signal-to-noise ratio (PSNR), normalized cross-correlation (NCC) indexes, and spatial nonuniformity (SNU) in the selected regions of interest (ROIs) were used to quantify the proposed algorithm's correction accuracy and generat the following results: mean MAE = 12.81 ± 2.04 and 19.94 ± 5.44 HU, mean PSNR = 40.22 ± 3.70 and 31.31 ± 2.85 dB, mean NCC = 0.98 ± 0.02 and 0.95 ± 0.01, and SNU = 2.07 ± 3.36% and 2.07 ± 3.36% for brain and pelvis data.

CONCLUSION

Preliminary results demonstrated that the novel learning-based correction method can significantly improve CBCT image quality. Hence, the proposed algorithm is of great potential in improving CBCT's image quality to support its clinical utility in CBCT-guided adaptive radiotherapy.

Statistical CT reconstruction using region-aware texture preserving regularization learning from prior normal-dose CT image

In some clinical applications, prior normal-dose CT (NdCT) images are available, and the valuable textures and structure features in them may be used to promote follow-up low-dose CT (LdCT) reconstruction. This study aims to learn texture information from the NdCT images and leverage it for follow-up LdCT image reconstruction to preserve textures and structure features. Specifically, the proposed reconstruction method first learns the texture information from those patches with similar structures in NdCT image, and the similar patches can be clustered by searching context features efficiently from the surroundings of the current patch. Then it utilizes redundant texture information from the similar patches as a priori knowledge to describe specific regions in the LdCT image. The advanced region-aware texture preserving prior is termed as 'RATP'. The main advantage of the PATP prior is that it can properly learn the texture features from available NdCT images and adaptively characterize the region-specific structures in the LdCT image. The experiments using patient data were performed to evaluate the performance of the proposed method. The proposed RATP method demonstrated superior performance in LdCT imaging compared to the filtered back projection (FBP) and statistical iterative reconstruction (SIR) methods using Gaussian regularization, Huber regularization and the original texture preserving regularization.

Iterative quality enhancement via residual-artifact learning networks for low-dose CT

Radiation exposure and the associated risk of cancer for patients in computed tomography (CT) scans have been major clinical concerns. The radiation exposure can be reduced effectively via lowering the x-ray tube current (mA). However, this strategy may lead to excessive noise and streak artifacts in the conventional filtered back-projection reconstructed images. To address this issue, some deep convolutional neural network (ConvNet) based approaches have been developed for low-dose CT imaging inspired by the recent development of machine learning. Nevertheless, some of the image textures reconstructed by the ConvNet could be corrupted by the severe streaks, especially in ultra-low-dose cases, which could be close to prostheses and hamper diagnosis. Therefore, in this work, we propose an iterative residual-artifact learning ConvNet (IRLNet) approach to improve the reconstruction performance over the ConvNet based approaches. Specifically, the proposed IRLNet estimates the high-frequency details within the noise and then removes them iteratively; after eliminating severe streaks in the low-dose CT images, the residual low-frequency details can be processed through the conventional network. Moreover, the proposed IRLNet scheme can be extended for robust handling of quantitative dual energy CT/cerebral perfusion CT imaging, and statistical iterative reconstruction. Real patient data are used to evaluate the proposed IRLNet, and the experimental results demonstrate that the proposed IRLNet approach outperforms the previous ConvNet based approaches in reducing the image noise and streak artifacts efficiently at the same time as preserving edge details well, suggesting that the proposed IRLNet approach can be used to improve the CT image quality, especially in ultra-low-dose cases.

Regularization strategies in statistical image reconstruction of low-dose x-ray CT: A review

Statistical image reconstruction (SIR) methods have shown potential to substantially improve the image quality of low-dose x-ray computed tomography (CT) as compared to the conventional filtered back-projection (FBP) method. According to the maximum a posteriori (MAP) estimation, the SIR methods are typically formulated by an objective function consisting of two terms: (a) a data-fidelity term that models imaging geometry and physical detection processes in projection data acquisition, and (b) a regularization term that reflects prior knowledge or expectations of the characteristics of the to-be-reconstructed image. SIR desires accurate system modeling of data acquisition, while the regularization term also has a strong influence on the quality of reconstructed images. A variety of regularization strategies have been proposed for SIR in the past decades, based on different assumptions, models, and prior knowledge. In this paper, we review the conceptual and mathematical bases of these regularization strategies and briefly illustrate their efficacies in SIR of low-dose CT.

Assessment of prior image induced nonlocal means regularization for low-dose CT reconstruction: Change in anatomy

PURPOSE

Repeated computed tomography (CT) scans are prescribed for some clinical applications such as lung nodule surveillance. Several studies have demonstrated that incorporating a high-quality prior image into the reconstruction of subsequent low-dose CT (LDCT) acquisitions can either improve image quality or reduce data fidelity requirements. Our proposed previous normal-dose image induced nonlocal means (ndiNLM) regularization method for LDCT is an example of such a method. However, one major concern with prior image based methods is that they might produce false information when the prior image and the current LDCT image show different structures (for example, if a lung nodule emerges, grows, shrinks, or disappears over time). This study aims to assess the performance of the ndiNLM regularization method in situations with change in anatomy.

METHOD

We incorporated the ndiNLM regularization into the statistical image reconstruction (SIR) framework for reconstruction of subsequent LDCT images. Because of its patch-based search mechanism, a rough registration between the prior image and the current LDCT image is adequate for the SIR-ndiNLM method. We assessed the performance of the SIR-ndiNLM method in lung nodule surveillance for two different scenarios: (a) the nodule was not found in a baseline exam but appears in a follow-up LDCT scan; (b) the nodule was present in a baseline exam but disappears in a follow-up LDCT scan. We further investigated the effect of nodule size on the performance of the SIR-ndiNLM method.

RESULTS

We found that a relatively large search-window (e.g., 33 × 33) should be used for the SIR-ndiNLM method to account for misalignment between the prior image and the current LDCT image, and to ensure that enough similar patches can be found in the prior image. With proper selection of other parameters, experimental results with two patient datasets demonstrated that the SIR-ndiNLM method did not miss true nodules nor introduce false nodules in the lung nodule surveillance scenarios described above. We also found that the SIR-ndiNLM reconstruction shows improved image quality when the prior image is similar to the current LDCT image in anatomy. These gains in image quality might appear small upon visual inspection, but they can be detected using quantitative measures. Finally, the SIR-ndiNLM method also performed well in ultra-low-dose conditions and with different nodule sizes.

CONCLUSIONS

This study assessed the performance of the SIR-ndiNLM method in situations in which the prior image and the current LDCT image show substantial anatomical differences, specifically, changes in lung nodules. The experimental results demonstrate that the SIR-ndiNLM method does not introduce false lung nodules nor miss true nodules, which relieves the concern that this method might produce false information. However, there is insufficient evidence that these findings will hold true for all kinds of anatomical changes.

PWLS-ULTRA: An Efficient Clustering and Learning-Based Approach for Low-Dose 3D CT Image Reconstruction

The development of computed tomography (CT) image reconstruction methods that significantly reduce patient radiation exposure, while maintaining high image quality is an important area of research in low-dose CT imaging. We propose a new penalized weighted least squares (PWLS) reconstruction method that exploits regularization based on an efficient Union of Learned TRAnsforms (PWLS-ULTRA). The union of square transforms is pre-learned from numerous image patches extracted from a dataset of CT images or volumes. The proposed PWLS-based cost function is optimized by alternating between a CT image reconstruction step, and a sparse coding and clustering step. The CT image reconstruction step is accelerated by a relaxed linearized augmented Lagrangian method with ordered-subsets that reduces the number of forward and back projections. Simulations with 2-D and 3-D axial CT scans of the extended cardiac-torso phantom and 3-D helical chest and abdomen scans show that for both normal-dose and low-dose levels, the proposed method significantly improves the quality of reconstructed images compared to PWLS reconstruction with a nonadaptive edge-preserving regularizer. PWLS with regularization based on a union of learned transforms leads to better image reconstructions than using a single learned square transform. We also incorporate patch-based weights in PWLS-ULTRA that enhance image quality and help improve image resolution uniformity. The proposed approach achieves comparable or better image quality compared to learned overcomplete synthesis dictionaries, but importantly, is much faster (computationally more efficient).

Low-dose x-ray tomography through a deep convolutional neural network

Synchrotron-based X-ray tomography offers the potential for rapid large-scale reconstructions of the interiors of materials and biological tissue at fine resolution. However, for radiation sensitive samples, there remain fundamental trade-offs between damaging samples during longer acquisition times and reducing signals with shorter acquisition times. We present a deep convolutional neural network (CNN) method that increases the acquired X-ray tomographic signal by at least a factor of 10 during low-dose fast acquisition by improving the quality of recorded projections. Short-exposure-time projections enhanced with CNNs show signal-to-noise ratios similar to long-exposure-time projections. They also show lower noise and more structural information than low-dose short-exposure acquisitions post-processed by other techniques. We evaluated this approach using simulated samples and further validated it with experimental data from radiation sensitive mouse brains acquired in a tomographic setting with transmission X-ray microscopy. We demonstrate that automated algorithms can reliably trace brain structures in low-dose datasets enhanced with CNN. This method can be applied to other tomographic or scanning based X-ray imaging techniques and has great potential for studying faster dynamics in specimens.