Sinogram smoothing techniques for myocardial blood flow estimation from dose-reduced dynamic computed tomography

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

In this paper, a newly developed statistical model-based image reconstruction [referred to as Simultaneous Multiple Artifacts Reduction in Tomographic RECONstruction (SMART-RECON)] is applied to low dose computer tomography (CT) myocardial perfusion imaging (CT-MPI). This method uses the nuclear norm of the spatial-temporal image matrix of the CT-MPI images as a regularizer, rather than a conventional spatial regularizer that incorporates image smoothness, edge preservation, or spatial sparsity into the reconstruction. In addition to providing the needed noise reduction for low-dose CT-MPI, SMART-RECON provides images with spatial resolution and noise power spectrum (NPS) properties, which are independent of contrast and dose levels. Both numerical simulations and in vivo animal studies were performed to validate the proposed method. In these studies, it was found that: 1) quantitative accuracy of perfusion maps in CT-MPI was well maintained for radiation dose level as low as 10 mAs per image frame, compared with the reference standard of 200 mAs for conventional filtered backprojection; 2) flow-occluded myocardium in the porcine heart was well delineated by SMART-RECON at 10 mAs per frame when compared with model-based image reconstruction using spatial total variation (TV) as the regularizer (referred to as TV-SIR) or spatial-temporal TV (ST-TV-SIR); the CT-MPI results were confirmed with positron-emission tomography imaging; 3) image sharpness in SMART-RECON images was nearly independent of image contrast level and radiation dose level, in stark contrast to TV-SIR and ST-TV-SIR, which displayed a strong dependence on both image contrast and radiation dose levels; and 4) the structure of the dose-normalized NPS for the SMART-RECON method did not depend on dose, while the TV-SIR and ST-TV-SIR NPS structure was dose-dependent.

Variable temporal sampling and tube current modulation for myocardial blood flow estimation from dose-reduced dynamic computed tomography

Quantification of myocardial blood flow (MBF) can aid in the diagnosis and treatment of coronary artery disease. However, there are no widely accepted clinical methods for estimating MBF. Dynamic cardiac perfusion computed tomography (CT) holds the promise of providing a quick and easy method to measure MBF quantitatively. However, the need for repeated scans can potentially result in a high patient radiation dose, limiting the clinical acceptance of this approach. In our previous work, we explored techniques to reduce the patient dose by either uniformly reducing the tube current or by uniformly reducing the number of temporal frames in the dynamic CT sequence. These dose reduction techniques result in noisy time-attenuation curves (TACs), which can give rise to significant errors in MBF estimation. We seek to investigate whether nonuniformly varying the tube current and/or sampling intervals can yield more accurate MBF estimates for a given dose. Specifically, we try to minimize the dose and obtain the most accurate MBF estimate by addressing the following questions: when in the TAC should the CT data be collected and at what tube current(s)? We hypothesize that increasing the sampling rate and/or tube current during the time frames when the myocardial CT number is most sensitive to the flow rate, while reducing them elsewhere, can achieve better estimation accuracy for the same dose. We perform simulations of contrast agent kinetics and CT acquisitions to evaluate the relative MBF estimation performance of several clinically viable variable acquisition methods. We find that variable temporal and tube current sequences can be performed that impart an effective dose of 5.5 mSv and allow for reductions in MBF estimation root-mean-square error on the order of 20% compared to uniform acquisition sequences with comparable or higher radiation doses.

Low-dose dynamic myocardial perfusion CT image reconstruction using pre-contrast normal-dose CT scan induced structure tensor total variation regularization

Dynamic myocardial perfusion CT (DMP-CT) imaging provides quantitative functional information for diagnosis and risk stratification of coronary artery disease by calculating myocardial perfusion hemodynamic parameter (MPHP) maps. However, the level of radiation delivered by dynamic sequential scan protocol can be potentially high. The purpose of this work is to develop a pre-contrast normal-dose scan induced structure tensor total variation regularization based on the penalized weighted least-squares (PWLS) criteria to improve the image quality of DMP-CT with a low-mAs CT acquisition. For simplicity, the present approach was termed as 'PWLS-ndiSTV'. Specifically, the ndiSTV regularization takes into account the spatial-temporal structure information of DMP-CT data and further exploits the higher order derivatives of the objective images to enhance denoising performance. Subsequently, an effective optimization algorithm based on the split-Bregman approach was adopted to minimize the associative objective function. Evaluations with modified dynamic XCAT phantom and preclinical porcine datasets have demonstrated that the proposed PWLS-ndiSTV approach can achieve promising gains over other existing approaches in terms of noise-induced artifacts mitigation, edge details preservation, and accurate MPHP maps calculation.

Robust dynamic myocardial perfusion CT deconvolution for accurate residue function estimation via adaptive-weighted tensor total variation regularization: a preclinical study

Dynamic myocardial perfusion computed tomography (MPCT) is a promising technique for quick diagnosis and risk stratification of coronary artery disease. However, one major drawback of dynamic MPCT imaging is the heavy radiation dose to patients due to its dynamic image acquisition protocol. In this work, to address this issue, we present a robust dynamic MPCT deconvolution algorithm via adaptive-weighted tensor total variation (AwTTV) regularization for accurate residue function estimation with low-mA s data acquisitions. For simplicity, the presented method is termed 'MPD-AwTTV'. More specifically, the gains of the AwTTV regularization over the original tensor total variation regularization are from the anisotropic edge property of the sequential MPCT images. To minimize the associative objective function we propose an efficient iterative optimization strategy with fast convergence rate in the framework of an iterative shrinkage/thresholding algorithm. We validate and evaluate the presented algorithm using both digital XCAT phantom and preclinical porcine data. The preliminary experimental results have demonstrated that the presented MPD-AwTTV deconvolution algorithm can achieve remarkable gains in noise-induced artifact suppression, edge detail preservation, and accurate flow-scaled residue function and MPHM estimation as compared with the other existing deconvolution algorithms in digital phantom studies, and similar gains can be obtained in the porcine data experiment.