Dictionary-Driven Ischemia Detection From Cardiac Phase-Resolved Myocardial BOLD MRI at Rest

• Published in: IEEE transactions on medical imaging Vol. 35; no. 1; pp. 282 - 293
• Main Authors: Bevilacqua, Marco, Dharmakumar, Rohan, Tsaftaris, Sotirios A.
• Format: Journal Article
• Language: English
• Published: United States IEEE 01.01.2016; The Institute of Electrical and Electronics Engineers, Inc. (IEEE)
• Subjects: Acute coronary syndromes; Algorithms; Blood-oxygen-level dependent; cardiac MRI; Data models; Dictionaries; dictionary learning; Heart attacks; Humans; Image Processing, Computer-Assisted - methods; Image segmentation; Ischemia; Magnetic resonance imaging; Magnetic Resonance Imaging - methods; Myocardial Ischemia - diagnosis; Myocardium; Oxygen - blood; shift-invariance; sparse representations; Support Vector Machine; Support vector machines; Time series analysis; Unsupervised Machine Learning
• ISSN: 0278-0062; 1558-254X; 1558-254X
• DOI: 10.1109/TMI.2015.2470075

Cardiac Phase-resolved Blood-Oxygen-Level Dependent (CP-BOLD) MRI provides a unique opportunity to image an ongoing ischemia at rest. However, it requires post-processing to evaluate the extent of ischemia. To address this, here we propose an unsupervised ischemia detection (UID) method which relies on the inherent spatio-temporal correlation between oxygenation and wall motion to formalize a joint learning and detection problem based on dictionary decomposition. Considering input data of a single subject, it treats ischemia as an anomaly and iteratively learns dictionaries to represent only normal observations (corresponding to myocardial territories remote to ischemia). Anomaly detection is based on a modified version of One-class Support Vector Machines (OCSVM) to regulate directly the margins by incorporating the dictionary-based representation errors. A measure of ischemic extent (IE) is estimated, reflecting the relative portion of the myocardium affected by ischemia. For visualization purposes an ischemia likelihood map is created by estimating posterior probabilities from the OCSVM outputs, thus obtaining how likely the classification is correct. UID is evaluated on synthetic data and in a 2D CP-BOLD data set from a canine experimental model emulating acute coronary syndromes. Comparing early ischemic territories identified with UID against infarct territories (after several hours of ischemia), we find that IE, as measured by UID, is highly correlated (Pearson's r=0.84) with respect to infarct size. When advances in automated registration and segmentation of CP-BOLD images and full coverage 3D acquisitions become available, we hope that this method can enable pixel-level assessment of ischemia with this truly non-invasive imaging technique.