Calibration‐free beam hardening correction for myocardial perfusion imaging using CT

PURPOSE: Computed tomography myocardial perfusion imaging (CT‐MPI) and coronary CTA have the potential to make CT an ideal noninvasive imaging gatekeeper exam for invasive coronary angiography. However, beam hardening (BH) artifacts prevent accurate blood flow calculation in CT‐MPI. BH correction methods require either energy‐sensitive CT, not widely available, or typically, a calibration‐based method in conventional CT. We propose a calibration‐free, automatic BH correction (ABHC) method suitable for CT‐MPI and evaluate its ability to reduce BH artifacts in single “static‐perfusion” images and to create accurate myocardial blood flow (MBF) in dynamic CT‐MPI. METHODS: In the algorithm, we used input CT DICOM images and iteratively optimized parameters in a polynomial BH correction until a BH‐sensitive cost function was minimized on output images. An input image was segmented into a soft tissue image and a highly attenuating material (HAM) image containing bones and regions of high iodine concentrations, using mean HU and temporal enhancement properties. We forward projected HAM, corrected projection values according to a polynomial correction, and reconstructed a correction image to obtain the current iteration's BH corrected image. The cost function was sensitive to BH streak artifacts and cupping. We evaluated the algorithm on simulated CT and physical phantom images, and on preclinical porcine with optional coronary obstruction and clinical CT‐MPI data. Assessments included measures of BH artifact in single images as well as MBF estimates. We obtained CT images on a prototype spectral detector CT (SDCT, Philips Healthcare) scanner that provided both conventional and virtual keV images, allowing us to quantitatively compare corrected CT images to virtual keV images. To stress test the method, we evaluated results on images from a different scanner (iCT, Philips Healthcare) and different kVp values. RESULTS: In a CT‐simulated digital phantom consisting of water with iodine cylinder insets, BH streak artifacts between simulated iodine inserts were reduced from 13 ± 2 to 0 ± 1 HU. In a similar physical phantom having higher iodine concentrations, BH streak artifacts were reduced from 48 ± 6 to 1 ± 5 HU and cupping was reduced by 86%, from 248 to 23 HU. In preclinical CT‐MPI images without coronary obstruction, BH artifact was reduced from 24 ± 6 HU to less than 5 ± 4 HU at peak enhancement. Standard deviation across different regions of interest (ROI) along the myocardium was reduced from 13.26 to 6.86 HU for ABHC, comparing favorably to measurements in the corresponding virtual keV image. Corrections greatly reduced variations in preclinical MBF maps as obtained in normal animals without obstruction (FFR = 1). Coefficients of variations were 22% (conventional CT), 9% (ABHC), and 5% (virtual keV). Moreover, variations in flow tended to be localized after ABHC, giving result which would not be confused with a flow deficit in a coronary vessel territory. CONCLUSION: The automated algorithm can be used to reduce BH artifact in conventional CT and improve CT‐MPI accuracy particularly by removing regions of reduced estimated flow which might be misinterpreted as flow deficits.