Radiation Dose Reduction in CT Myocardial Perfusion Imaging Using SMART-RECON

DIOR: Deep Iterative Optimization-Based Residual-Learning for Limited-Angle CT Reconstruction

Limited-angle CT is a challenging problem in real applications. Incomplete projection data will lead to severe artifacts and distortions in reconstruction images. To tackle this problem, we propose a novel reconstruction framework termed Deep Iterative Optimization-based Residual-learning (DIOR) for limited-angle CT. Instead of directly deploying the regularization term on image space, the DIOR combines iterative optimization and deep learning based on the residual domain, significantly improving the convergence property and generalization ability. Specifically, the asymmetric convolutional modules are adopted to strengthen the feature extraction capacity in smooth regions for deep priors. Besides, in our DIOR method, the information contained in low-frequency and high-frequency components is also evaluated by perceptual loss to improve the performance in tissue preservation. Both simulated and clinical datasets are performed to validate the performance of DIOR. Compared with existing competitive algorithms, quantitative and qualitative results show that the proposed method brings a promising improvement in artifact removal, detail restoration and edge preservation.

CLEAR: Comprehensive Learning Enabled Adversarial Reconstruction for Subtle Structure Enhanced Low-Dose CT Imaging

X-ray computed tomography (CT) is of great clinical significance in medical practice because it can provide anatomical information about the human body without invasion, while its radiation risk has continued to attract public concerns. Reducing the radiation dose may induce noise and artifacts to the reconstructed images, which will interfere with the judgments of radiologists. Previous studies have confirmed that deep learning (DL) is promising for improving low-dose CT imaging. However, almost all the DL-based methods suffer from subtle structure degeneration and blurring effect after aggressive denoising, which has become the general challenging issue. This paper develops the Comprehensive Learning Enabled Adversarial Reconstruction (CLEAR) method to tackle the above problems. CLEAR achieves subtle structure enhanced low-dose CT imaging through a progressive improvement strategy. First, the generator established on the comprehensive domain can extract more features than the one built on degraded CT images and directly map raw projections to high-quality CT images, which is significantly different from the routine GAN practice. Second, a multi-level loss is assigned to the generator to push all the network components to be updated towards high-quality reconstruction, preserving the consistency between generated images and gold-standard images. Finally, following the WGAN-GP modality, CLEAR can migrate the real statistical properties to the generated images to alleviate over-smoothing. Qualitative and quantitative analyses have demonstrated the competitive performance of CLEAR in terms of noise suppression, structural fidelity and visual perception improvement.

Motion compensation combining with local low rank regularization for low dose dynamic CT myocardial perfusion reconstruction

Dynamic CT myocardial perfusion imaging (DCT-MPI) is a reliable examination tool for the assessment of myocardium and vascular, while its special scan protocol may result in excessive radiation exposure to patients and inevitable inter-frame motion. Lowering the tube current is a simple way to reduce radiation exposure. However, low mAs will certainly cause severe image noise, thus may further impact the accuracy of functional hemodynamic parameters, which are used for the assessment of blood supply. In this work, we present a novel scheme applying motion compensation and local low rank regularization (MC-LLR) for obtaining high quality motion compensated DCT-MPI images. Specifically, motion compensation by using robust data decomposition registration (RDDR) was introduced. Robust principal component analysis coupled with optical flow-based registration algorithm were used in RDDR. Then, the local low rank constraint on the motion compensated time series images was applied for the DCT-MPI reconstruction. One healthy mini pig and two patient datasets were used to evaluate the proposed MC-LLR algorithm. Results show that the present method achieved satisfactory image quality with higher CNRs, smaller rRMSEs, and more accurate hemodynamic parameter maps.

Temporal feature prior-aided separated reconstruction method for low-dose dynamic myocardial perfusion computed tomography

Dynamic myocardial perfusion computed tomography (DMP-CT) is an effective medical imaging technique for coronary artery disease diagnosis and therapy guidance. However, the radiation dose received by the patient during repeated CT scans is a widespread concern of radiologists because of the increased risk of cancer. The sparse few-view CT scanning protocol can be a feasible approach to reduce the radiation dose of DMP-CT imaging; however, an advanced reconstruction algorithm is needed. In this paper, a temporal feature prior-aided separated reconstruction method (TFP-SR) for low-dose DMP-CT images reconstruction from sparse few-view sinograms is proposed. To implement the proposed method, the objective perfusion image is divided into the baseline fraction and the enhancement fraction introduced by the arrival of the contrast agent. The core of the proposed TFP-SR method is the utilization of the temporal evolution information that naturally exists in the DMP-CT image sequence to aid the enhancement image reconstruction from limited data. The temporal feature vector of an image pixel is defined by the intensities of this pixel in the pre-reconstructed enhancement sequence, and the connection between two related features is calculated via a zero-mean Gaussian function. A prior matrix is constructed based on the connections between the extracted temporal features and used in the iterative reconstruction of the enhancement images. To evaluate the proposed method, the conventional filtered back-projection algorithm, the total variation regularized PWLS (PWLS-TV) and the prior image constrained compressed sensing are compared in this paper based on studies on a digital extended cardiac-torso (XCAT) thoracic phantom and a preclinical porcine DMP-CT data set that take image misregistration into account. The experimental results demonstrate that the proposed TFP-SR method has superior performance in sparse DMP-CT images reconstruction in terms of image quality and the analyses of the time attenuation curve and hemodynamic parameters.

Myocardial perfusion assessment in the infarct core and penumbra zones in an in-vivo porcine model of the acute, sub-acute, and chronic infarction

OBJECTIVES

To assess the longitudinal changes of microvascular function in different myocardial regions after myocardial infarction (MI) using myocardial blood flow derived by dynamic CT perfusion (CTP-MBF), and compare CTP-MBF with the results of cardiac magnetic resonance (CMR) and histopathology.

METHODS

The CTP scanning was performed in a MI porcine model 1 day (n = 15), 7 days (n = 10), and 3 months (n = 5) following induction surgery. CTP-MBF was measured in the infarcted myocardium, penumbra, and remote myocardium, respectively. CMR perfusion and histopathology were performed for validation.

RESULTS

From baseline to follow-up scans, CTP-MBF presented a stepwise increase in the infarcted myocardium (68.51 ± 11.04 vs. 86.73 ± 13.32 vs. 109.53 ± 26.64 ml/100 ml/min, p = 0.001) and the penumbra (104.92 ± 29.29 vs. 120.32 ± 24.74 vs. 183.01 ± 57.98 ml/100 ml/min, p = 0.008), but not in the remote myocardium (150.05 ± 35.70 vs. 166.66 ± 38.17 vs. 195.36 ± 49.64 ml/100 ml/min, p = 0.120). The CTP-MBF correlated with max slope (r = 0.584, p < 0.001), max signal intensity (r = 0.357, p < 0.001), and time to max (r = - 0.378, p < 0.001) by CMR perfusion. Moreover, CTP-MBF defined the infarcted myocardium on triphenyl tetrazolium chloride staining (AUC: 0.810, p < 0.001) and correlated with microvascular density on CD31 staining (r = 0.561, p = 0.002).

CONCLUSION

CTP-MBF could quantify the longitudinal changes of microvascular function in different regions of the post-MI myocardium, which demonstrates good agreement with contemporary CMR and histopathological findings.

KEY POINTS

• The CT perfusion-based myocardial blood flow (CTP-MBF) could quantify the microvascular impairment in different myocardial regions after myocardial infarction (MI) and track its recovery over time. • The assessment of CTP-MBF is in good agreement with contemporary cardiac MRI and histopathological findings, which potentially facilitates a rapid approach for pathophysiological insights following MI.

Contrast-Medium Anisotropy-Aware Tensor Total Variation Model for Robust Cerebral Perfusion CT Reconstruction with Low-Dose Scans

Perfusion computed tomography (PCT) is critical in detecting cerebral ischemic lesions. PCT examination with low-dose scans can effectively reduce radiation exposure to patients at the cost of degraded images with severe noise and artifacts. Tensor total variation (TTV) models are powerful tools that can encode the regional continuous structures underlying a PCT object. In a TTV model, the sparsity structures of the contrast-medium concentration (CMC) across PCT frames are assumed to be isotropic with identical and independent distribution. However, this assumption is inconsistent with practical PCT tasks wherein the sparsity has evident variations and correlations. Such modeling deviation hampers the performance of TTV-based PCT reconstructions. To address this issue, we developed a novel contrast-medium anisotropy-aware tensor total variation (CMAA-TTV) model to describe the intrinsic anisotropy sparsity of the CMC in PCT imaging tasks. Instead of directly on the difference matrices, the CMAA-TTV model characterizes sparsity on a low-rank subspace of the difference matrices which are calculated from the input data adaptively, thus naturally encoding the intrinsic variant and correlated anisotropy sparsity structures of the CMC. We further proposed a robust and efficient PCT reconstruction algorithm to improve low-dose PCT reconstruction performance using the CMAA-TTV model. Experimental studies using a digital brain perfusion phantom, patient data with low-dose simulation and clinical patient data were performed to validate the effectiveness of the presented algorithm. The results demonstrate that the CMAA-TTV algorithm can achieve noticeable improvements over state-of-the-art methods in low-dose PCT reconstruction tasks.

Constructing a tissue-specific texture prior by machine learning from previous full-dose scan for Bayesian reconstruction of current ultralow-dose CT images

Purpose: Bayesian theory provides a sound framework for ultralow-dose computed tomography (ULdCT) image reconstruction with two terms for modeling the data statistical property and incorporating a priori knowledge for the image that is to be reconstructed. We investigate the feasibility of using a machine learning (ML) strategy, particularly the convolutional neural network (CNN), to construct a tissue-specific texture prior from previous full-dose computed tomography. Approach: Our study constructs four tissue-specific texture priors, corresponding with lung, bone, fat, and muscle, and integrates the prior with the prelog shift Poisson (SP) data property for Bayesian reconstruction of ULdCT images. The Bayesian reconstruction was implemented by an algorithm called SP-CNN-T and compared with our previous Markov random field (MRF)-based tissue-specific texture prior algorithm called SP-MRF-T. Results: In addition to conventional quantitative measures, mean squared error and peak signal-to-noise ratio, structure similarity index, feature similarity, and texture Haralick features were used to measure the performance difference between SP-CNN-T and SP-MRF-T algorithms in terms of the structure and tissue texture preservation, demonstrating the feasibility and the potential of the investigated ML approach. Conclusions: Both training performance and image reconstruction results showed the feasibility of constructing CNN texture prior model and the potential of improving the structure preservation of the nodule comparing to our previous regional tissue-specific MRF texture prior model.