Assessing Reliability of Myocardial Blood Flow After Motion Correction With Dynamic PET Using a Bayesian Framework

Assessment of motion correction in dynamic rubidium-82 cardiac PET with and without frame-by-frame adjustment of attenuation maps for calculation of myocardial blood flow

BACKGROUND

Patient motion during pharmacological stressing can have substantial impact on myocardial blood flow (MBF) estimated from dynamic PET. This work evaluated a motion correction algorithm with and without adjustment of the PET attenuation map.

METHODS

Frame-by-frame motion correction was performed by three users on 30 rubidium-82 studies. Data were divided equally into three groups of motion severity [mild (M1), moderate (M2) and severe (M3)]. MBF data were compared for non-motion corrected (NC), motion-corrected-only (MC) and with adjustment of the attenuation map (MCAC). Percentage differences of MBF were calculated in the coronary territories and 17-segment polar plots. Polar plots of spill-over were also generated from the data.

RESULTS

Median differences of 23% were seen in the RCA and 18% for the LAD in the M3 category for MC vs NC images. Differences for MCAC vs MC images were considerably smaller and typically < 10%. Spill-over plots for MC and MCAC were notably more uniform compared with NC images.

CONCLUSION

Motion correction for dynamic rubidium data is desirable for future MBF software updates. Adjustment of the PET attenuation map results in only marginal differences and therefore is unlikely to be an essential requirement. Assessing the uniformity of spill-over plots is a useful visual aid for verifying motion correction techniques.

Lung 4D CT Image Registration Based on High-Order Markov Random Field

To solve the problem that traditional image registration methods based on continuous optimization for large motion lung 4D CT image sequences are easy to fall into local optimal solutions and lead to serious misregistration, a novel image registration method based on high-order Markov Random Field (MRF) is proposed. By analyzing the effect of the deformation field constraint of the potential functions with different order cliques in MRF model, energy functions with high-order cliques form are designed separately for 2D and 3D images to preserve topology of the deformation field. In order to preserve the topology of the deformation field more effectively, it is necessary to apply a smooth term and a topology preservation term simultaneously in the energy function and use logarithmic function to impose a penalty on the Jacobian matrix with high-order cliques in the topology preservation term. For the complexity of the designed energy function with high-order cliques form, Markov Chain Monte Carlo (MCMC) algorithm is used to solve the optimization problem of the designed energy function. To address the high computational requirements in lung 4D CT image registration, a multi-level processing strategy is adopted to reduce the space complexity of the proposed registration method and promotes the computational efficiency. In the DIR-lab dataset with 4D CT images and the COPD (Chronic Obstructive Pulmonary Disease) dataset with 3D CT images, the average target registration error (TRE) of our proposed method can reach 0.95 mm respectively.