Reduction of half-scan shading artifact based on full-scan correction

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.

Noise spatial nonuniformity and the impact of statistical image reconstruction in CT myocardial perfusion imaging

PURPOSE

To achieve high temporal resolution in CT myocardial perfusion imaging (MPI), images are often reconstructed using filtered backprojection (FBP) algorithms from data acquired within a short-scan angular range. However, the variation in the central angle from one time frame to the next in gated short scans has been shown to create detrimental partial scan artifacts when performing quantitative MPI measurements. This study has two main purposes. (1) To demonstrate the existence of a distinct detrimental effect in short-scan FBP, i.e., the introduction of a nonuniform spatial image noise distribution; this nonuniformity can lead to unexpectedly high image noise and streaking artifacts, which may affect CT MPI quantification. (2) To demonstrate that statistical image reconstruction (SIR) algorithms can be a potential solution to address the nonuniform spatial noise distribution problem and can also lead to radiation dose reduction in the context of CT MPI.

METHODS

Projection datasets from a numerically simulated perfusion phantom and an in vivo animal myocardial perfusion CT scan were used in this study. In the numerical phantom, multiple realizations of Poisson noise were added to projection data at each time frame to investigate the spatial distribution of noise. Images from all datasets were reconstructed using both FBP and SIR reconstruction algorithms. To quantify the spatial distribution of noise, the mean and standard deviation were measured in several regions of interest (ROIs) and analyzed across time frames. In the in vivo study, two low-dose scans at tube currents of 25 and 50 mA were reconstructed using FBP and SIR. Quantitative perfusion metrics, namely, the normalized upslope (NUS), myocardial blood volume (MBV), and first moment transit time (FMT), were measured for two ROIs and compared to reference values obtained from a high-dose scan performed at 500 mA.

RESULTS

Images reconstructed using FBP showed a highly nonuniform spatial distribution of noise. This spatial nonuniformity led to large fluctuations in the temporal direction. In the numerical phantom study, the level of noise was shown to vary by as much as 87% within a given image, and as much as 110% between different time frames for a ROI far from isocenter. The spatially nonuniform noise pattern was shown to correlate with the source trajectory and the object structure. In contrast, images reconstructed using SIR showed a highly uniform spatial distribution of noise, leading to smaller unexpected noise fluctuations in the temporal direction when a short scan angular range was used. In the numerical phantom study, the noise varied by less than 37% within a given image, and by less than 20% between different time frames. Also, the noise standard deviation in SIR images was on average half of that of FBP images. In the in vivo studies, the deviation observed between quantitative perfusion metrics measured from low-dose scans and high-dose scans was mitigated when SIR was used instead of FBP to reconstruct images.

CONCLUSIONS

(1) Images reconstructed using FBP suffered from nonuniform spatial noise levels. This nonuniformity is another manifestation of the detrimental effects caused by short-scan reconstruction in CT MPI. (2) Images reconstructed using SIR had a much lower and more uniform noise level and thus can be used as a potential solution to address the FBP nonuniformity. (3) Given the improvement in the accuracy of the perfusion metrics when using SIR, it may be desirable to use a statistical reconstruction framework to perform low-dose dynamic CT MPI.

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

PURPOSE

Partial scan reconstruction (PSR) artifacts are present in myocardial perfusion imaging using dynamic multidetector computed tomography (MDCT). PSR artifacts appear as temporal CT number variations due to inconsistencies in the angular data range used to reconstruct images and compromise the quantitative value of myocardial perfusion when using MDCT. The purpose of this work is to present and evaluate a technique termed targeted spatial frequency filtration (TSFF) to reduce CT number variations due to PSR when applied to myocardial perfusion imaging using MDCT.

METHODS

The TSFF algorithm requires acquiring enough X-ray projections to reconstruct both partial (π + fan angle α) and full (2π) scans. Then, using spatial linear filters, the TSFF-corrected image data are created by superimposing the low spatial frequency content of the full scan reconstruction (containing no PSR artifacts, but having low spatial resolution and worse temporal resolution) with the high spatial frequency content of the partial scan reconstruction (containing high spatial frequencies and better temporal resolution). The TSFF method was evaluated both in a static anthropomorphic thoracic phantom and using an in vivo porcine model and compared with a previously validated reference standard technique that avoids PSR artifacts by pacing the animal heart in synchrony with the gantry rotation. CT number variations were quantified by measuring the range and standard deviation of CT numbers in selected regions of interest (ROIs) over time. Myocardial perfusion parameters such as blood volume (BV), mean transit time (MTT), and blood flow (BF) were quantified and compared in the in vivo study.

RESULTS

Phantom experiments demonstrated that TSFF reduced PSR artifacts by as much as tenfold, depending on the location of the ROI. For the in vivo experiments, the TSFF-corrected data showed two- to threefold decrease in CT number variations. Also, after TSFF, the perfusion parameters had an average difference of 13.1% (range 4.5%-25.6%) relative to the reference method, in contrast to an average difference of 31.8% (range 0.3%-54.0%) between the non-TSFF processed data with the reference method.

CONCLUSIONS

TSFF demonstrated consistent reduction in CT number variations due to PSR using controlled phantom and in vivo experiments. TSFF-corrected data provided quantitative measures of perfusion (BV, MTT, and BF) with better agreement to a reference method compared to noncorrected data. Practical implementation of TSFF is expected to incur in an additional radiation exposure of 14%, when tube current is modulated to 20% of its maximum, to complete the needed full scan reconstruction.

Quantitative myocardial CT perfusion: a pictorial review and the current state of technology development

Coronary artery disease (CAD) is one of the leading causes of morbidity and mortality, and is associated with substantial and increasing resource burden. A combined physiologic and anatomic assessment may improve identification of patients with CAD who would benefit from revascularization and reduce unnecessary diagnostic and interventional procedures. Cardiovascular computed tomography (CT) has the potential to provide a comprehensive evaluation of CAD in a single setting. Although coronary CT angiography has been widely implemented for clinical use, the application of myocardial CT perfusion (CTP) has been relatively restricted because of a few limitations, such as beam hardening and the high radiation dose delivered. In this article, we first review the fundamental basis of the qualitative, semiquantitative, and quantitative techniques for myocardial CTP and discussed the strength and weakness of each approach. Beam-hardening correction for myocardial CTP with image-based method and dual-energy CT are then discussed with example cases demonstrating the effectiveness of each method. Initial experiences suggest both techniques can reduce beam-hardening artifact to a satisfactory extent. An overview on dose reduction technologies, such as prospective ECG triggering and iterative reconstruction for myocardial CTP, is also provided. Preclinical studies suggest it is feasible to perform low-dose quantitative myocardial CTP without affecting perfusion measurement. Finally, the impact of scan length on myocardial CTP is addressed.