New automated Markov–Gibbs random field based framework for myocardial wall viability quantification on agent enhanced cardiac magnetic resonance images

A novel automated framework for detecting and quantifying viability from agent enhanced cardiac magnetic resonance images is proposed. The framework identifies the pathological tissues based on a joint Markov–Gibbs random field (MGRF) model that accounts for the 1st-order visual appearance of the myocardial wall (in terms of the pixel-wise intensities) and the 2nd-order spatial interactions between pixels. The pathological tissue is quantified based on two metrics: the percentage area in each segment with respect to the total area of the segment, and the trans-wall extent of the pathological tissue. This transmural extent is estimated using point-to-point correspondences based on a Laplace partial differential equation. Transmural extent was validated using a simulated phantom. We tested the proposed framework on 14 datasets (168 images) and validated against manual expert delineation of the pathological tissue by two observers. Mean Dice similarity coefficients (DSC) of 0.90 and 0.88 were obtained for the observers, approaching the ideal value, 1. The Bland–Altman statistic of infarct volumes estimated by manual versus the MGRF estimation revealed little bias difference, and most values fell within the 95% confidence interval, suggesting very good agreement. Using the DSC measure we documented statistically significant superior segmentation performance for our MGRF method versus established intensity-based methods (greater DSC, and smaller standard deviation). Our Laplace method showed good operating characteristics across the full range of extent of transmural infarct, outperforming conventional methods. Phantom validation and experiments on patient data confirmed the robustness and accuracy of the proposed framework.