Comparison of automated beam hardening correction (ABHC) algorithms for myocardial perfusion imaging using computed tomography.

purposeMyocardial perfusion imaging using computed tomography (MPI-CT) and coronary CT angiography (CTA) have the potential to make CT an ideal noninvasive imaging gatekeeper exam for invasive coronary angiography. However, beam hardening can prevent accurate blood flow estimation in dynamic MPI-CT and can create artifacts that resemble flow deficits in single-shot MPI-CT. In this work, we compare four automatic beam hardening correction algorithms (ABHCs) applied to CT images, for their ability to produce accurate single images of contrast and accurate MPI flow maps using images from conventional CT systems, without energy sensitivity.

methodsPreviously, we reported a method, herein called ABHC-1, where we iteratively optimized a cost function sensitive to beam hardening artifacts in MPI-CT images and used a low order polynomial correction on projections of segmentation-processed CT images. Here, we report results from two new algorithms with higher order polynomial corrections, ABHC-2 and ABHC-3 (with three and seven free parameters, respectively), having potentially better correction but likely reduced estimability. Additionally, we compared results to an algorithm reported by others in the literature (ABHC-NH). Comparisons were made on a digital static phantom with simulated water, bone, and iodine regions; on a digital dynamic anthropomorphic phantom, with simulated blood flow; and on preclinical porcine experiments. We obtained CT images on a prototype spectral detector CT (Philips Healthcare) scanner that provided both conventional and virtual keV images, allowing us to quantitatively compare corrected CT images to virtual keV images. To test these methods' parameter optimization sensitivity to noise, we evaluated results on images obtained using different mAs.

resultsIn images of the static phantom, ABHC-2 reduced beam hardening artifacts better than our previous ABHC-1 algorithm, giving artifacts smaller than 1.8 HU, even in the presence of high noise which should affect parameter optimization. Taken together, the quality of static phantom results ordered ABHC-2> ABHC-3> ABHC-1>> ABHC-NH. In an anthropomorphic MPI-CT simulator with homogeneous myocardial blood flow of 100 ml⋅min

conclusionThe automated algorithm can be used with different parametric BH correction models. ABHC-2 improved MPI-CT blood flow estimation as compared to other approaches and was robust to noisy images. In simulation and preclinical experiments, ABHC-2 gave results approaching gold standard 70 keV measurements.