A strategy to decrease partial scan reconstruction artifacts in myocardial perfusion CT: phantom and in vivo evaluation

Myocardial late enhancement and extracellular volume with single-energy, dual-energy, and photon-counting computed tomography

Computed tomography late enhancement (CT-LE) is emerging as a non-invasive technique for cardiac diagnosis with wider accessibility compared to MRI, despite its typically lower contrast-to-noise ratio. Optimizing CT-LE image quality necessitates a thorough methodology addressing contrast administration, timing, and radiation dose, alongside a robust understanding of extracellular volume (ECV) quantification methods. This review summarizes CT-LE protocols, clinical utility, and advances in ECV measurement through both single-energy and dual-energy CT. It also highlights photon-counting detector CT technology as an innovative means to potentially improve image quality and reduce radiation exposure.

Dynamic Computed Tomography Angiography for capturing vessel wall motion: A phantom study for optimal image reconstruction

BACKGROUND

Reliably capturing sub-millimeter vessel wall motion over time, using dynamic Computed Tomography Angiography (4D CTA), might provide insight in biomechanical properties of these vessels. This may improve diagnosis, prognosis, and treatment decision making in vascular pathologies.

PURPOSE

The aim of this study is to determine the most suitable image reconstruction method for 4D CTA to accurately assess harmonic diameter changes of vessels.

METHODS

An elastic tube (inner diameter 6 mm, wall thickness 2 mm) was exposed to sinusoidal pressure waves with a frequency of 70 beats-per-minute. Five flow amplitudes were set, resulting in increasing sinusoidal diameter changes of the elastic tube, measured during three simulated pulsation cycles, using ECG-gated 4D CTA on a 320-detector row CT system. Tomographic images were reconstructed using one of the following three reconstruction methods: hybrid iterative (Hybrid-IR), model-based iterative (MBIR) and deep-learning based (DLR) reconstruction. The three reconstruction methods where based on 180 degrees (half reconstruction mode) and 360 degrees (full reconstruction mode) raw data. The diameter change, captured by 4D CTA, was computed based on image registration. As a reference metric for diameter change measurement, a 9 MHz linear ultrasound transducer was used. The sum of relative absolute differences (SRAD) between the ultrasound and 4D CTA measurements was calculated for each reconstruction method. The standard deviation was computed across the three pulsation cycles.

RESULTS

MBIR and DLR resulted in a decreased SRAD and standard deviation compared to Hybrid-IR. Full reconstruction mode resulted in a decreased SRAD and standard deviations, compared to half reconstruction mode.

CONCLUSIONS

4D CTA can capture a diameter change pattern comparable to the pattern captured by US. DLR and MBIR algorithms show more accurate results than Hybrid-IR. Reconstruction with DLR is >3 times faster, compared to reconstruction with MBIR. Full reconstruction mode is more accurate than half reconstruction mode.

Effect of image quality on myocardial extracellular volume quantification using cardiac computed tomography: a phantom study

BACKGROUND

The image quality directly affects the accuracy of computed tomography (CT) extracellular volume (ECV) quantification.

PURPOSE

To investigate the effects of image quality and acquisition protocol on the accuracy of ECV quantification.

MATERIAL AND METHODS

One-volume scans were performed on a 320-row multidetector CT volume scanner using a multi-energy CT phantom. To simulate the blood pool and myocardium, solid rods representing blood and soft tissue were used in precontrast CT. Moreover, the solid rods including different iodine concentrations were used in postcontrast CT. The tube voltage was set at 120 kVp, and the tube current was changed from 750 mA (100% dose) to 190 mA (25% dose). All images underwent full- and half-scan reconstructions based on model-based iterative reconstruction. The ECV was calculated from the CT numbers between pre- and postcontrast.

RESULTS

The mean ECV with full- and half-scan reconstructions at the central portion was 0.275 at 100% scan dose to 0.271 at 25% scan dose and 0.276 at 100% scan dose to 0.269 at 25% scan dose. Compared with that in the 100% scan dose, the variation in each ECV increased with decreasing radiation dose. The ECV at the center of the image along the z-axis had lower variation than that at outer portion of the images. On the reconstruction algorithm, there was no statistical difference in ECVs with full- and half-scan reconstructions.

CONCLUSION

For stable ECV quantifications, excessive radiation dose reduction may be inappropriate, and it is better to consider the variations in ECV values depending on the slice location.

Evaluation of Four-Dimensional Computed Tomography as a Technique for Quantifying Carpal Motion

Delayed diagnosis of dynamic carpal instability often occurs because early changes in bone alignment and movement are difficult to detect and manifest mainly during a dynamic/functional task. Current diagnostic tools are only able to examine the carpal bones under static or sequential-static conditions. Four-dimensional (three dimensions + time) computed tomography (4DCT) enables quantification of carpal mechanics through 3D volume sequences of the wrist in motion. A comprehensive understanding of carpal mechanics is needed to define normal function and structure and provide targets for treatment of carpal injuries. In this study, measurements of scaphoid translation and joint congruency were taken by creating models from the CT scans of the carpals in extreme frames of motion, registering those models to the neutral position, transforming the models into a local coordinate system, and using software to calculate the joint surface areas (JSA). Results indicated that the centroid of the scaphoid translated 6.4 ± 1.3 mm and extended from extreme radial to extreme ulnar deviation. Results are consistent with the literature. An additional study was performed to measure the responsiveness of the 4DCT technique presented. Bone models from each frame of motion for radio ulnar deviation (RUD) and flexion extension (FE) were created and distinct differences between their JSA were measured qualitatively and quantitatively. The results show that there was statistically significantly different JSA within carpal joints between RUD and FE. These studies provide the first step in developing the methodology when using 4DCT scanning to measure subtle abnormalities in the wrist.

Evaluation of a Quantitative Method for Carpal Motion Analysis Using Clinical 3-D and 4-D CT Protocols

For wrist complaints related to motion, a 2-D radiograph or CT scan of the static wrist may not always be considered diagnostic. 3-D motion imaging, i.e., multiple 3DCT scans in time (4DCT), enables quantifying carpal motion and comparing motion patterns of the affected wrist with those of the healthy contralateral side. The accuracy and precision of the method, however, is limited by noise and motion artifacts. Although, the technique is considered promising in existing literature, the accuracy and precision of carpal motion analysis has never been investigated systematically. In this paper, we introduce and evaluate a semi-automatic segmentation- and registration-based method for 3-D carpal motion analysis. We investigate the accuracy and precision of the method, and its dependency on motion and scan parameters (angular velocity, dose, gantry revolution angle for image reconstruction, and scanner type) using a wrist phantom. During standstill the positioning error was ≤ 0.23 mm and ≤ 0.78°. A partial gantry revolution for 3-D reconstruction introduced image deformation, contributing to a positioning error of approx. 0.8 mm. This error increased with reduced dose, and with increasing angular velocity of the wrist phantom. In cases where the phantom was rotating about an axis parallel to the rotation axis of the gantry, and in a direction opposite to the gantry, the positioning error increased, probably because of the apparent increase in angular velocity with respect to the gantry. Slow carpal motion 4DCT analysis is feasible using a regular CT scanner. A partial gantry revolution angle for 3-D reconstruction may introduce image deformation, which decreases the accuracy of carpal motion analysis. Knowing the positioning error in 4DCT imaging with the proposed method is considered valuable when investigating wrist injury since it enables discrimination of actual motion from apparent motion caused by methodological error.

Low-dose dynamic myocardial perfusion CT imaging using a motion adaptive sparsity prior

PURPOSE

Dynamic myocardial perfusion computed tomography (DM-PCT) imaging offers benefits over quantitative assessment of myocardial blood flow (MBF) for diagnosis and risk stratification of coronary artery disease. However, one major drawback of DM-PCT imaging is that a high radiation level is imparted by repeated scanning. To address this issue, in this work, we developed a statistical iterative reconstruction algorithm based on the penalized weighted least-squares (PWLS) scheme by incorporating a motion adaptive sparsity prior (MASP) model to achieve high-quality DM-PCT imaging with low tube current dynamic data acquisition. For simplicity, we refer to the proposed algorithm as "PWLS-MASP''.

METHODS

The MASP models both the spatial and temporal structured sparsity of DM-PCT sequence images with the assumption that the differences between adjacent frames after motion correction are sparse in the gradient image domain. To validate and evaluate the effectiveness of the present PWLS-MASP algorithm thoroughly, a modified XCAT phantom and preclinical porcine DM-PCT dataset were used in the study.

RESULTS

The present PWLS-MASP algorithm can obtain high-quality DM-PCT images in both phantom and porcine cases, and outperforms the existing filtered back-projection algorithm and PWLS-based algorithms with total variation regularization (PWLS-TV) and robust principal component analysis regularization (PWLS-RPCA) in terms of noise reduction, streak artifacts mitigation, and time density curve estimation. Moreover, the PWLS-MASP algorithm can yield more accurate diagnostic hemodynamic parametric maps than the PWLS-TV and PWLS-RPCA algorithms.

CONCLUSIONS

The study indicates that there is a substantial advantage in using the present PWLS-MASP algorithm for low-dose DM-PCT, and potentially in other dynamic tomography areas.

Utility of Dual-Energy CT-based Monochromatic Imaging in the Assessment of Myocardial Delayed Enhancement in Patients with Cardiomyopathy

Purpose To investigate the diagnostic utility of dual-energy computed tomography (CT)-based monochromatic imaging for myocardial delayed enhancement (MDE) assessment in patients with cardiomyopathy. Materials and Methods The institutional review board approved this prospective study, and informed consent was obtained from all participants who were enrolled in the study. Forty patients (27 men and 13 women; mean age, 56 years ± 15 [standard deviation]; age range, 22-81 years) with cardiomyopathy underwent cardiac magnetic resonance (MR) imaging and dual-energy CT. Conventional (120-kV) and monochromatic (60-, 70-, and 80-keV) images were reconstructed from the dual-energy CT acquisition. Subjective quality score, contrast-to-noise ratio (CNR), and beam-hardening artifacts were compared pairwise with the Friedman test at post hoc analysis. With cardiac MR imaging as the reference standard, diagnostic performance of dual-energy CT in MDE detection and its predictive ability for pattern classification were compared pairwise by using logistic regression analysis with the generalized estimating equation in a per-segment analysis. The Bland-Altman method was used to find agreement between cardiac MR imaging and CT in MDE quantification. Results Among the monochromatic images, 70-keV CT images resulted in higher subjective quality (mean score, 3.38 ± 0.54 vs 3.15 ± 0.43; P = .0067), higher CNR (mean, 4.26 ± 1.38 vs 3.93 ± 1.33; P = .0047), and a lower value for beam-hardening artifacts (mean, 3.47 ± 1.56 vs 4.15 ± 1.67; P < .0001) when compared with conventional CT. When compared with conventional CT, 70-keV CT showed improved diagnostic performance for MDE detection (sensitivity, 94.6% vs 90.4% [P = .0032]; specificity, 96.0% vs 94.0% [P = .0031]; and accuracy, 95.6% vs 92.7% [P < .0001]) and improved predictive ability for pattern classification (subendocardial, 91.5% vs 84.3% [P = .0111]; epicardial, 94.3% vs 73.5% [P = .0001]; transmural, 93.0% vs 77.7% [P = .0018]; mesocardial, 85.4% vs 69.2% [P = .0047]; and patchy. 84.4% vs 78.4% [P = .1514]). For MDE quantification, 70-keV CT showed a small bias 0.1534% (95% limits of agreement: -4.7013, 5.0080). Conclusion Dual-energy CT-based 70-keV monochromatic images improve MDE assessment in patients with cardiomyopathy via improved image quality and CNR and reduced beam-hardening artifacts when compared with conventional CT images. ^© RSNA, 2017 Online supplemental material is available for this article.

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.

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.

A new redundancy weighting scheme for nonstationary data for computed tomography

PURPOSE

The same projection data (or line integrals) are often measured multiple times, e.g., twice from opposite directions during one gantry rotation. The redundant data must be normalized by applying redundancy weighting such as the halfscan algorithm, which assumes that the noise of the data is uniform. This assumption, however, is not correct when a tube current modulation technique is employed. The variance of line integrals, which is inversely related to the tube current, could vary significantly. The purpose of this work is to improve how the projection data are used during analytical reconstruction when the tube current is modulated during the scan.

METHODS

The authors developed a new redundancy weighting scheme. It not only takes into account the data statistics but also can control how much to weigh the statistics from 100% (αs = 1.0) to 0% (αs = 0.0) by a parameter αs. The proposed weighting scheme reduces to the conventional redundancy weighting scheme when αs = 0.0. The authors evaluated the performance of the proposed scheme using computer simulations targeting at myocardial perfusion CT imaging. The image quality was evaluated in terms of the image noise and halfscan artifacts, and perfusion defect detection performance was evaluated by the positive predictive value (PPV) and the area-under-the-receiver operating characteristic-curve (AUC) value.

RESULTS

Results showed a tradeoff between the image noise and halfscan artifacts. The normalized noise standard deviation was 1.00 with halfscan, 0.89 with αs = 1.0, 0.97 with αs = 0.5, and 1.20 with αs = 0.0 when projections over one rotation (75% of projections are acquired with full dose, 25% with 1/10 of the full dose) are used. The halfscan artifacts were 13.4 Hounsfield unit (HU) with halfscan, 8.2 HU with αs = 1.0, 4.5 HU with αs = 0.5, and 3.1 HU with αs = 0.0. Both the PPVs and AUCs were improved from the halfscan method: PPV, 69.0%-70.6% vs 58.0%, P < 0.003; AUC, 0.935-0.938 vs 0.908, P < 0.003.

CONCLUSIONS

The new redundancy weight allows for decreasing the image noise and controlling the tradeoff between the image noise and artifacts.

Low dose dynamic CT myocardial perfusion imaging using a statistical iterative reconstruction method

PURPOSE

Dynamic CT myocardial perfusion imaging has the potential to provide both functional and anatomical information regarding coronary artery stenosis. However, radiation dose can be potentially high due to repeated scanning of the same region. The purpose of this study is to investigate the use of statistical iterative reconstruction to improve parametric maps of myocardial perfusion derived from a low tube current dynamic CT acquisition.

METHODS

Four pigs underwent high (500 mA) and low (25 mA) dose dynamic CT myocardial perfusion scans with and without coronary occlusion. To delineate the affected myocardial territory, an N-13 ammonia PET perfusion scan was performed for each animal in each occlusion state. Filtered backprojection (FBP) reconstruction was first applied to all CT data sets. Then, a statistical iterative reconstruction (SIR) method was applied to data sets acquired at low dose. Image voxel noise was matched between the low dose SIR and high dose FBP reconstructions. CT perfusion maps were compared among the low dose FBP, low dose SIR and high dose FBP reconstructions. Numerical simulations of a dynamic CT scan at high and low dose (20:1 ratio) were performed to quantitatively evaluate SIR and FBP performance in terms of flow map accuracy, precision, dose efficiency, and spatial resolution.

RESULTS

Forin vivo studies, the 500 mA FBP maps gave -88.4%, -96.0%, -76.7%, and -65.8% flow change in the occluded anterior region compared to the open-coronary scans (four animals). The percent changes in the 25 mA SIR maps were in good agreement, measuring -94.7%, -81.6%, -84.0%, and -72.2%. The 25 mA FBP maps gave unreliable flow measurements due to streaks caused by photon starvation (percent changes of +137.4%, +71.0%, -11.8%, and -3.5%). Agreement between 25 mA SIR and 500 mA FBP global flow was -9.7%, 8.8%, -3.1%, and 26.4%. The average variability of flow measurements in a nonoccluded region was 16.3%, 24.1%, and 937.9% for the 500 mA FBP, 25 mA SIR, and 25 mA FBP, respectively. In numerical simulations, SIR mitigated streak artifacts in the low dose data and yielded flow maps with mean error <7% and standard deviation <9% of mean, for 30 × 30 pixel ROIs (12.9 × 12.9 mm(2)). In comparison, low dose FBP flow errors were -38% to +258%, and standard deviation was 6%-93%. Additionally, low dose SIR achieved 4.6 times improvement in flow map CNR(2) per unit input dose compared to low dose FBP.

CONCLUSIONS

SIR reconstruction can reduce image noise and mitigate streaking artifacts caused by photon starvation in dynamic CT myocardial perfusion data sets acquired at low dose (low tube current), and improve perfusion map quality in comparison to FBP reconstruction at the same dose.

Early changes in myocardial microcirculation in asymptomatic hypercholesterolemic subjects: as detected by perfusion CT

Intramyocardial microvessels show functional changes in early stages of atherosclerosis prior to epicardial coronary artery stenosis. However, clinical CT does not have adequate spatial resolution to resolve the microvessels. To clinically detect changes in the function of the intramyocardial microcirculation, the spatial heterogeneity of the distribution of myocardial perfusion (F) and intramyocardial microcirculatory blood volume (Bv) was determined by perfusion CT. Two human subject groups were studied: (i) a "Control" group (24) with no risk factors nor evidence of coronary artery disease (CAD), and (ii) an "At-Risk" group (24) with hypercholesterolemia, but no evidence of CAD. In the perfusion CT image, a region of interest (ROI) covering the left ventricular myocardium was subdivided into multiple nested ROI (nROI) of equal size and used to compute F and Bv for each nROI. No significant differences between the groups were demonstrable in overall myocardial F, or Bv. The nROI data showed significantly increased spatial heterogeneity in the "At Risk" group when compared to "Control" subjects. This study demonstrates that subresolution distribution at the microcirculatory level can be quantified with myocardial perfusion CT and significant changes in these parameters occur in hypercholesterolemic subjects before they have developed significant changes in conventional perfusion parameters.

Myocardial perfusion imaging with cardiac computed tomography: state of the art

Cardiac computed tomography (CCT) has become an important tool for the anatomic assessment of patients with suspected coronary disease. Its diagnostic accuracy for detecting the presence of underlying coronary artery disease and ability to risk stratify patients are well documented. However, the role of CCT for the physiologic assessment of myocardial perfusion during resting and stress conditions is only now emerging. With the addition of myocardial perfusion imaging to coronary imaging, CCT has the potential to assess both coronary anatomy and its functional significance with a single non-invasive test. In this review, we discuss the current state of CCT myocardial perfusion imaging for the detection of myocardial ischemia and myocardial infarction and examine its complementary role to CCT coronary imaging.