 Research
 Open Access
 Published:
Probability mapping of scarred myocardium using texture and intensity features in CMR images
BioMedical Engineering OnLine volume 12, Article number: 91 (2013)
Abstract
Background
The myocardium exhibits heterogeneous nature due to scarring after Myocardial Infarction (MI). In Cardiac Magnetic Resonance (CMR) imaging, Late Gadolinium (LG) contrast agent enhances the intensity of scarred area in the myocardium.
Methods
In this paper, we propose a probability mapping technique using Texture and Intensity features to describe heterogeneous nature of the scarred myocardium in Cardiac Magnetic Resonance (CMR) images after Myocardial Infarction (MI). Scarred tissue and nonscarred tissue are represented with high and low probabilities, respectively. Intermediate values possibly indicate areas where the scarred and healthy tissues are interwoven. The probability map of scarred myocardium is calculated by using a probability function based on Bayes rule. Any set of features can be used in the probability function.
Results
In the present study, we demonstrate the use of two different types of features. One is based on the mean intensity of pixel and the other on underlying texture information of the scarred and nonscarred myocardium. Examples of probability maps computed using the mean intensity of pixel and the underlying texture information are presented. We hypothesize that the probability mapping of myocardium offers alternate visualization, possibly showing the details with physiological significance difficult to detect visually in the original CMR image.
Conclusion
The probability mapping obtained from the two features provides a way to define different cardiac segments which offer a way to identify areas in the myocardium of diagnostic importance (like core and border areas in scarred myocardium).
Introduction
Patients who have suffered but survived myocardial infarction (MI) may subsequently suffer a possibly disabling or fatal cardiac arrhythmia. Late Gadolinium (LG) enhanced Cardiac Magnetic Resonance imaging (CMR) is used for assessing morphology of the myocardium in patients after MI. An example image is depicted in Figure 1 with myocardium and scar delineated by green and red lines, respectively. The myocardial muscle fibers are completely dead at the core area of the scar. Thus, the core area does not react to the electrical signals propagating through the heart muscle telling the heart to contract.
At the border areas, however, the muscle fibers are not all dead. The electrical signals will be disturbed in these areas possibly causing reentries and sometimes arrhythmias. The latter gives us reason to believe that the border areas are a very important part of the scar, and thus good ways of defining and visualizing border areas would be beneficial in many situations. Previous works [1, 2] emphasize the importance of the border area (also referred to as periinfarct area, or gray zone area) being a major determinant for identifying patients at risk of arrhythmia. The heterogeneous nature of cardiac segments is related to various diagnoses and prognosis. For example, Risk of serious arrhythmias and better indication for implantation of implantable cardioverter defibrillator (ICD) [3], heart rate of ventricular tachycardia [4], risk and degree of heart failure, indication of implantation of pacemakers in order to synchronize cardiac contractions [5] etc. Our main focus in this work is to develop methods to visualize, define and analyze the myocardium to identify regions of diagnostic importance in the myocardium.
Earlier reported work
To our knowledge, visualization of the heterogeneous nature of the myocardium is mostly based on thresholding methods described in the literature where thresholds are defined at intensity levels corresponding to some percentage of the maximum intensity level inside the scar area [1–3]. Recently, Shokrollahi et al. was able to divide the scarred myocardium of animals into core and gray zones using fuzzy clustering algorithm and validated the results using invivo electroanatomical voltage mapping [6]. Some previous publications reported automatic segmentation of scarred myocardium in late enhanced magnetic resonance images [7–9]. The segmentation of scar itself is an important application as the scar size contains information useful for finding the patients with high and low risk of arrhythmia [3]. Also, scar size is largely responsible in ventricular remodeling [10]. Dicki et al. use support vector machine (SVM) to classify the scar and the myocardium in delayed enhancement magnetic resonance images in [7] using three features. The three features used in [7] depend on the intensity of scar, myocardium and the whole CMR image.
From our group’s previous work [11], it is shown that patients with high and low risk of getting arrhythmia can be classified using textural measures from graylevel cooccurrence matrices. In our recent work [12], we presented some preliminary results where dictionary learning and sparse representations were used for representing the texture in the myocardial segments followed by a classification.
The proposed probability mapping
All the above discussed previous work have sought to identify methods to segment the hyperenhanced area of the myocardium. However, by only delineating the hyperenhanced myocardium, important information concerning the tissue within the infarct territory is not provided. Hence, our focus is on developing a method that provides probability values as a means to distinguish the areas of the myocardium rather than the crisp definitions offered by a scar vs no scar segmentation.
In this work, we propose a technique for transforming the scarred myocardium in the CMR image to a probability map to indicate varying degrees of its functioning capabilities: where high and low probabilities will indicate the scarred tissue and normal tissue, respectively. Section ‘Calculation of probability maps using Bayes rule’ describes the general principles for calculation of the probability maps, where each pixel is the posterior probability of that pixel being scar compared to healthy myocardium. The posterior probabilities are calculated using Bayes rule with the help of features from the myocardium. Features are calculated pixel by pixel based on a local neighborhood and is explained in Section ‘Extraction of features’. We have based our experiments on two features; a local mean (direct current, DC) that will correspond well to how the doctors manually segment the scar, and a textural feature using sparse representations, that might reveal information regarding gray zone areas.
Section ‘Experiments and results’ shows the experimental framework and examples of probability maps for both features. Also in section ‘Experiments and results’, cardiac segments example is demonstrated to show the diagnostic importance of texture feature probability maps. Finally, the paper is concluded in section ‘Conclusion’.
Calculation of probability maps using Bayes rule
The class specific probability density functions (PDFs) modeling specific feature vector values known to be located within scar or healthy myocardial areas are estimated as p(vscar) and p(vmyo), respectively. According to Bayes rule [13], the probability of a pixel being in the scar area is related to the density function and the prior probability function:
where p(v) = P(scar)p(vscar)+P(myo)p(vmyo), and v is a feature vector. The prior of scar, P(scar), reflects the mean value of the scar size in the training data set, and P(scar) is calculated as the ratio of the number of pixels present in the scar with respect to the myocardium in the entire training data.
The parameters in the class specific PDFs can be estimated using parametric or nonparametric methods, and here they are calculated using maximum likelihood (ML) estimation [14] which is a method of low complexity.
Our primary focus is deriving the probability mapping, and there is a strong link to segmentation as we define the functions so that they should be optimal in the sense of discriminatively expressing the properties of the scar and normal tissue. The minimum error classifier corresponds to classifying as scar when P(scarv)>P(myov). Therefore, we evaluate the features and feature combinations used in various mapping with respect to their segmentation capability. As we will discuss later, this might not be entirely fair as in our case; features reflecting the cardiologist’s crisp decisions are favored.
The extraction of feature vectors for myocardium is illustrated in section ‘Extraction of features’. Feature vector v of the myocardium in general can be represented as:
where d is the dimension of the feature vector v. Finally, the probability maps are obtained by regenerating a new image such that each pixel in the myocardium is given by the value of P(scarv) as calculated for that pixel. The probability values range from 0–1 and they are displayed using a color code shown in Figure 2(b) and (c). The shades of red and yellow in the color code represents high probability of scar and shades of green and blue represents low probability.
Training and testing
The process of generating probability maps is depicted in Figure 3. The whole CMR data set is divided into training and testing images. In the training process, training features are extracted along with some training parameters required in the testing phase. In the next step, the calculated training feature vectors are used for estimating the PDFs of scar and myocardium region using ML. The training process and its outputs are indicated by the dashed lines. The probability maps are generated on test images, different from the training images, and this process is indicated by solid lines. The extraction of features for training and testing phases are explained in section ‘Extraction of features’.
Extraction of features
As discussed in section ‘Introduction’, two features, DC value and texture feature are used to generate probability maps of the scarred myocardium in this work. The detailed process of extracting the two features are discussed in the following sections ‘DC Feature’ and ‘Dictionarybased textural feature’.
DC Feature
Historically, in electronics field, the mean is commonly referred to as DC (direct current) value [15]. Because of the intensity differences between the healthy and scarred myocardium tissues in the LG enhanced CMR images, DC value dc(i,j) is one of the features explored for developing the probability maps. This is probably the feature that would mimic the cardiologists segmentation best since they use the intensity values of scarred myocardium in CMR image to segment the scarred areas manually. We define the feature $\mathit{\text{dc}}(i,j)=\mathit{\text{mean}}\left({I}_{\sqrt{N}\times \sqrt{N}}\right(i,j\left)\right)$ (where ${I}_{\sqrt{N}\times \sqrt{N}}$ is the neighborhood around the pixel I(i,j)), such that a new image I_{ dc }, with dc(i,j) as values at pixel position (i,j) is made with a sliding window averaging over the image.
Dictionarybased textural feature
A dictionary D, is an ensemble of a finite number of atoms, d^{(k)}, and can be used to represent signals. A linear combination of some of the atoms in the dictionary gives exact or approximate representation of the original signal. Let a column vector x of finite length N represent the original signal. The dictionary atoms are arranged as columns in an N × K matrix D where K > N. The representation or the approximation of the signal, $\stackrel{~}{x}$ and the representation error, r can be expressed as:
where w is sparse coefficient vector. In a sparse representation, only a small number of the coefficients w(k) are allowed to be nonzero. Finding the sparse coefficient vector can be a NPhard (Nondeterministic Polynomialtime hard) problem and formulated as:
The pseudonorm ∥.∥_{0} is the number of nonzero elements and, s denotes the sparseness. A vector selection algorithm called Order Recursive Matching Pursuit (ORMP) algorithm is used for sparse approximation [16].
Dictionary learning is the task of learning or training a dictionary on an available training set such that it adapts well to represent that specific class of signals. The training vectors, L and the sparse coefficients are arranged as columns in the matrices X and W, respectively. The objective in dictionary learning is to give a sparse representation of the training set X in order to minimize the sum of the squared error. The cost function is formulated as:
This is a very hard optimization problem. Recursive Least Squares Dictionary Learning Algorithm (RLSDLA), [17] is used for dictionary learning in all our experiments. A brief explanation of the algorithm is given in Appendix.
Sparse representations and learned dictionaries have been shown to work well for texture classification by Skretting and Husøy in [18] and by Mairal et al.[19]. Texture in a small image patch can be modeled as sparse linear combination of dictionary atoms in Frame Texture Classification Method (FTCM) presented by Skretting et al. in [18]. FTCM is developed by modeling a texture as a tiled floor where all the tiles are identical. The color or the gray level, at a given position in the floor is given by an underlying continuous periodic twodimensional function. It is shown based on this model that a vector from spatial neighborhood is indeed a sparse combination of finite dictionary atoms.
The algorithm for sparse representation in our work proceeds as follows in the training phase: Consider the myocardium in a CMR image, I that contains two texture classes: healthy and scarred myocardium. The training vector, y_{ l } of length N is made from that specific pixel and its neighborhood, $\sqrt{N}\times \sqrt{N}$ for each pixel in the training image. In the training set, each pixel is classified into a specified texture class. Then, the dictionaries, D_{ s } and D_{ m } are trained for the predefined texture classes (scar and myocardium) using RLSDLA. Every training vector, y_{ l } is then represented sparsely using ORMP vector selection algorithm [16] twice, both with D_{ s } and D_{ m }. The residuals (or representation errors), R_{ s } and R_{ m }, are calculated for each pixel in the myocardium of a training image as:
where ${w}_{l}^{s}$ and${w}_{l}^{m}$ are sparse coefficient vectors.
The residuals, R_{ s } and R_{ m }, are combined to form the texture feature used in this work. The texture feature is calculated as the ratio of the residual of D_{ m } to the sum of the residuals of D_{ s } and D_{ m }:
Pixel by pixel R_{ p } value can be interpreted as the scaling of residuals R_{ s } and R_{ m }. Smaller values of R_{ p } means that the pixel is not likely to be scar (i.e. healthy myocardium), and larger values means that it is likely to be scar or border area. The texture feature, R_{ p }(i,j) from the training set is used to estimate the PDFs for the probability maps computation. In the testing phase, texture feature vectors are collected from the test image set, in the same way as training feature vectors. In Figure 3, training and classification parameters are indicated. In fact, these training parameters are the dictionaries, D_{ s } and D_{ m }. Using D_{ s } and D_{ m } learned in the training phase, the residual images, R_{ s } and R_{ m } are calculated for test images. Texture is not pixelbypixel local except in the edge between two textures; thus some sort of smoothing is used on the test residual images before calculating the probability maps [18]. The test residual images R_{ s }, and R_{ m } are smoothed using an a × a pixels separable Gaussian low pass filter with variance σ^{2} before calculating R_{ p }.
Experiments and results
The CMR images used in our experiment were provided by the Department of Cardiology in Stavanger University Hospital and were obtained from 44 patients with myocardial infarction. 24 patients had old MI and ICD was implanted in all patients as primary or secondary prophylaxis. This group was defined as the high risk arrhythmic group. The remaining 20 patients had MI one year prior to CMR imaging and during this observation period, no reported incidents of serious arrhythmias occurred. This group was defined as the low risk arrhythmic group. The study was approved by the Regional Ethics Committee and informed consent was obtained from all patients.
All the CMR images were obtained from 1.5 Tesla Philips Intera machine. Images were acquired with a pixel size of 0.82 × 0.82 mm^{2}, covering the whole ventricle with shortaxis slices of 10 m thickness, without interslice gap. The LG enhanced CMR images were stored according to the digital imaging and communications in medicine (DICOM) format with 512 × 512 pixel resolution. The number of image slices with visible scar in each patient varies depending on the size of the scar and the scar size varies from one slice to the other. Only shortaxis image slices were used in our experiments. Out of the total 444 CMR images belonging to the 44 patients, 11 CMR images of poor quality due to motion and artifacts were not included. An experienced cardiologist delineated the borders of the myocardium and the infarction tissues as shown in Figure 1. The resulting set of segmented images were further split into labeled training and test sets. All our experiments were carried out in MATLAB. The CMR training and test images were cropped from 512×512 original CMR image so that they contain only left ventricle with myocardium and infarction areas to save execution time. Preprocessing of any kind was not used on these cropped images.
In all the CMR images, we took into account only the myocardium segmented by the cardiologists. Patients with high risk of getting arrhythmia have profound scars with possibly large grayzone areas. We therefore anticipated that images from these patients contained the information we wanted to represent in the dictionaries: information from the myocardium, the core area of scars, and the grayzone areas. Hence, the high risk patients group was used in the training phase. From the group of 24 high risk patients, 6 patients were randomly selected to form a training set, whereas the remaining 18 high risk patients were used to form a test set. Two sets of training vectors were generated, one from the scar, and one from the healthy myocardium.
This section is organized as follows: Sections ‘Probability maps obtained using DC feature, dc(i,j)’ and ‘Probability maps obtained using texture feature, R_{ p }(i,j)’ explains the processes of calculating the probability maps using DC and texture feature respectively. Section ‘Examples of probability maps obtained from DC and texture feature’ discusses the examples of probability maps obtained in our experiments. Section ‘A possible application  cardiac segments’ illustrates how probability maps obtained with texture feature can be used to define cardiac segments. Section ‘The discriminatory power of the DC and texture features’ considers the discriminative power of the DC and texture features.
Probability maps obtained using DC feature, dc(i,j)
The neighborhood size 3×3 was used to form training vectors as explained in ‘DC Feature’. The same neighborhood size must be used while training and finding the DC images I_{ dc }. The DC images obtained from the training images were used to form the training feature set. The parameters of the class specific PDFs: the mean and the standard deviation were found using ML estimation using the training feature set. DCvalues were scaled to have zero mean and unit variance before finding the ML estimates. The scaling coefficients from the training were stored to scale the test vectors. The probability maps of the scarred myocardium were calculated using Bayes rule as explained in section ‘Calculation of probability maps using Bayes rule’. An example of a probability map obtained using the DC feature with color code is shown in Figure 2(b). Note the sharp transition from the segmented scar to the myocardium.
Probability maps obtained using texture feature, R_{ p }(i,j)
The training vector and test vectors were generated in the same way as in the DC feature experiment using the same neighborhood size 3×3. Consider the pixels on the border zones, their neighborhood extends into other regions that are not under consideration. If we use training vectors from border regions, the dictionaries might learn the texture properties of other regions along with the texture properties they are intended to learn. So, the training vectors for the pixels whose neighborhood span other regions were not considered in our experiments. This is depicted in Figure 4.
The dictionaries were learned on the training vectors generated from scar and healthy myocardium using RLSDLA. The dictionary size of 9 × 90 atoms was used in our experiments. The initial dictionaries were formed by randomly selecting 90 vectors of length 9 from the training sets. Dictionary size was selected based on previous experience, and it should be K > N < < L. The forgetting factor is initialized to 0.995 and slowly increased towards 1 according to the exponential method described in [17]. If the number of dictionary atoms used to represent the image patch increases, then the residual decreases, and a large number of atoms will provide a good approximation with any full rank dictionary. Therefore, higher sparsity lowers the difference between the residuals of healthy myocardium and scar which in turn decreases the differentiation between the classes. The sparsity s used in our work is two, i.e. we used two vectors from the dictionary to represent the original image patch.
After learning the dictionaries, they were used to obtain the residual images. For each pixel in the myocardium of training images, the scaled residual R_{ p }(i,j) was found from residual images as explained in section ‘Dictionarybased textural feature’. The training feature set was then generated from these scaled residual values. The parameters of the PDFs to be used with the classifier was estimated from the labeled training vector set using ML estimation. The test vector set in this experiment was calculated after smoothing (using low pass Gaussian filter with σ = 5 and 9 × 9 window size) the test residual images. Using the ML estimates and Bayes rule, the probability maps for texture features of the scarred myocardium were obtained. An example of such a probability map with color code is shown in Figure 2(c). The color code used is explained under Figure 2. Here, the dictionaries were learned without the removal of the DC value in each image patch.
Dictionary learning was performed with and without the removal of the DC factor from the image patch of the training set. In order to illustrate some of the texture captured by the dictionary atoms, they were plotted as image patches in Figure 5. The difference between the scar and the healthy myocardium dictionaries was visible when the DC value was not removed, and this is according to the known intensity difference between scar and healthy myocardium. To investigate the variation in the dictionary atoms used as image patches, we had calculated a measure from the first derivative of dictionary atoms in both vertical and horizontal directions. The absolute mean of the first derivative of the dictionary atoms was used to compare the scar and healthy myocardium dictionaries learned with and without DC value. The calculated absolute means are shown under Figure 5. The difference between the scar and healthy dictionaries without DC removal was observed clearly, both from the depicted dictionaries and the absolute means. The difference between the scar and healthy myocardium dictionaries trained by removing the DC value from image patches could not be observed visually. A larger block for DC removal prior to textural capturing is probably more appropriate, more similar to suppression of difference in illuminance conditions in computer vision.
Figure 6 shows probability maps of some example CMR images and all of them were calculated using the texture feature, R_{ p }(i,j). In the middle column, the DC value was not removed from each image patch before training dictionaries and calculating the probabilities. In the last column, the DC value was removed from every image patch both before learning dictionaries and before finding probabilities. The DC value is obviously important, so we do not expect this experiment to provide competitive results, but we see a tendency of higher values in the scar and border areas, but not necessarily in the core of the scar. We expect the scar core to be very bright and could be captured easily by some sort of intensity based measure (for example the DC feature). If the physical nature of the tissue is reflected in the image texture, we can expect the scar core to be quite homogeneous. Thus, we are not surprised that the regions that are the brightest in column II were actually quite dark in column III (core), but the intermediate regions were captured fairly. We have continued with the experiments in this work by letting the dictionaries incorporate the DC value as we did not combine the DC feature, dc and texture feature, R_{ p } into a 2D dimensional feature for computing the probability maps. Unless it is mentioned under figures of probability maps, all the probability maps of texture feature, R_{ p } were computed with dictionaries learned without removing DC from every image patch.
Examples of probability maps obtained from DC and texture feature
The second and the third column of Figure 7, shows more examples of probability maps of scarred myocardium using the DC, dc(i,j), and texture, R_{ p }(i,j), features, respectively. Figure 8 shows examples of probability maps of scarred myocardium for both DC and texture feature for the patients without ICD. This shows that the PDFs estimated for the DC and texture features from the patient group implanted with ICD were also able to produce probability maps for low risk arrhythmia patients (i.e. patients without ICD). Both the features were able to identify the healthy myocardium without scars, and an example of this is shown in Figure 9. In order to compare the probability maps obtained by the two methods, they are plotted with the same color map. Experiments of probability mapping combining the DC feature and texture feature gave almost identical probability maps as the DC feature because of the dominance of the DC feature over the texture feature. Improved ways of combining both the features will be explored in our future work.
A possible application  cardiac segments
A cardiac segment, CS can be defined as a region of interest in the myocardium which might have diagnostic importance. We define a candidate cardiac segment by a lower and upper probability so that the segment consists of all pixels with probability values in the range between these. Different candidate cardiac segments might be evaluated, for example, calculating one or several features from the candidate segment and correlate result to a clinical meaningful hypothesis. We demonstrate the concept by making a comprehensive search for significant cardiac segments by comparing two patient groups and testing for statistical significance. The demonstration experiment is meant to illustrate how the probability mapping can be used to localize cardiac segments containing information discriminating between the high and low risk patient groups.
For each patient, cardiac segments were determined by a lower boundary L and an upper boundary U. These boundaries were varied in the range 0–1 in steps of 0.025 with the restriction that U  L ≥ 0.1. Figure 10(a) illustrates this by visualizing each of the resulting candidate cardiac segment as a black dot with coordinate (L,U). For example, the three yellow dots in coordinates (0.15,0.3), (0.3,0.7), and (0.7,0.825) correspond to the cardiac segments CS_{1} = 0.15  0.3, CS_{2} = 0.3  0.7, and CS_{3} = 0.7  0.825, respectively. Figure 10(b) and (c) show one of the original image slices and the corresponding probability mapping (texture feature) for one of the high risk patients. The probability maps, which do not have 0  1 probability range, was extended using a sigmoid function for easy comparability. Figure 10(d–f) shows the cardiac segments corresponding to CS_{1}, CS_{2} and CS_{3} for the image slice shown in Figure 10(b).
Preliminary experiments indicate that corresponding segments determined without using the probability mappings depends on the signal intensity values and seems to lack the ability to define contiguous areas in the same way as the probability mapping for some of the cardiac segments. For example, the cardiac segment determined directly from thresholding relative to the maximum 3D scar signal intensity is shown in Figure 10(g–i) for boundaries corresponding to CS_{1}CS_{2} and CS_{3}, respectively. These plots are comparable to the corresponding cardiac segments based on the probability mapping in Figure 10(d–f), and show less contiguous cardiac segment, especially CS_{1}. Figure 11 shows similar plots of the cardiac segments corresponding to CS_{1}, CS_{2} and CS_{3} for one of the low risk patients.
For each candidate cardiac segment, its relative size (accumulated for all slices) to the total myocardium was calculated. The computed relative sizes for high and low risk patient groups were compared using a MannWhithney test at significance level 0.05. The LUplots in Figure 12 show the pvalues for the test of all the candidate cardiac segments derived using the texture based (a) and DC based (b) probability map and those derived without the use of probability mapping (c). The blue colored dots indicate cardiac segments for which the high and low risk groups showed significant differences of the relative size.
Interpreting the comparison results for the texture feature probability map based cardiac segments in Figure 12 (a), the significant CS with L = 0 and U > 0.5 correspond to the cardiac segment containing all of the healthy myocardium and gradually including the scar area. More interesting is the localization in the upper right corner which makes more specific localization in the myocardial tissue in the regions close to the scar. Note that CS_{3} is identified as significant while CS_{1} and CS_{2} are not. Figure 12 (b) shows the LUplot from the probability maps of DC feature and, the plot shows that cardiac segments does not give any localization. This is due to the narrow transition between high and low probabilities of the DC based probability maps. The DC probability maps will predominantly have extreme values close to 0 or 1 and few values in between 0–1. This provides little information to explore. These plots (Figure 12 (a) and Figure 12 (b)) of cardiac segments clearly shows that the probability maps obtained from the DC features reflect the way cardiologist perceive the scarred myocardium in CMR images and, the probability maps from texture features gives information about the characteristics of myocardial tissue and might emphasize information that is not easy for inspection with a human eye in the original CMR image. The difference between the narrow and broad transitions for the DC and texturebased probability mappings might explain this. As the narrow transition gives little information besides the crisp core/scar discrimination corresponding to the cardiologist’s delineation based on the visual information. Thus, most of the candidate cardiac segments will be empty as there are virtually no intermediate probability values. This is not the case for the texture based probability mappings, which discloses the information present in this intermediate probability range.
In comparison the results for the cardiac segments derived directly without using the probability plots (Figure 12 (c)), the localization seems vaguer, not specifically identifying any coherent region in the myocardium. It is our belief that probability mapping might be used to identify cardiac segments that might be used in future studies.
The discriminatory power of the DC and texture features
The probability maps calculated for the two features, DC and textural feature, give different information. The discriminatory power of the two features were explored by comparing to the manual segmentation of scar and healthy myocardium to obtain Receiver Operating Characteristic (ROC) curves. The ROC curves calculated for the DC and texture features on the CMR images of the 18 test patients with high risk of getting arrhythmia are shown in Figure 13. The area under the average ROC curves is 0.9052 and 0.8428 for the DC and texture features, respectively. From ROC curves, it is observed that the discriminative power of the DC feature is comparatively higher than that of the texture features. The sensitivity, and specificity reported in [7] are plotted as a single point in Figure 13. This reported result is not calculated on our CMR data. This single point lies below within 95% confidence interval of DC feature. This reported segmentation result performed well compared to our average DC and texture ROC curves as Dikici et al.[7] does not include all the CMRI slices of a patient for training and testing where as in our work we include all the CMRI slices i.e. the scar volume.
Conclusion
In this paper, a new technique has been proposed for enhanced visualization of scarred myocardium in LG enhanced CMR images using probability mapping. Using Bayes rule and ML estimation, we obtained the probability maps of scarred myocardium using DC and texture features from each pixel in the myocardium. The probability values obtained from the DC features can be used for determining the scar size, and they reflect the way cardiologist perceive the scarred myocardium in CMR images. The probability maps from texture features gives information about the characteristics of myocardial tissue and emphasize information that is not easy for inspection with a human eye in the original CMR image. Using textural features in the probability mapping provide us with a tool for defining cardiac segments as well as give us features that can be used to distinguish between different groups of patients (patients high and low risk of getting arrhythmia). Our belief that the probability maps calculated by the texture feature captures the heterogeneous nature of scarred myocardium is justified with the example of the cardiac segments experiment. The two features tested in this paper gave valuable, but different, information. In future work, benefit of combining these two features as well as other features will be explored.
Intellectual Property Rights Patent pending. The Patent Application Number is 1121307.1. All Intellectual Property Rights are reserved by the University of Stavanger and Stavanger University Hospital. Prekubator TTO is responsible for IP management and commercialization.
Appendix
Recursive least squares dictionary learning algorithm
The dictionary learning problem can be solved in a two steps algorithm. Step 1) Find sparse coefficient matrix, W, using vector selection algorithms, keeping dictionary, D, fixed. Step 2) Keeping W fixed, find D. The dictionary update step depends on the dictionary learning method we choose to use. D and C are initialized randomly.
In this paper, an online dictionary learning algorithm, Recursive Least Squares Dictionary Learning Algorithm (RLSDLA) presented in [17] is used for dictionary learning and ORMP vector selection algorithm presented in [16] is used for sparse representation. The RLSDLA updates the dictionary with the arrival of each new training vector. In deriving updating rules for RLSDLA, X_{ t } = [x_{1},x_{2},…,x_{ t }] of size S × t, W_{ t } = [w_{1},w_{2},…,w_{ t }] of size K × t and ${C}_{t}={\left({W}_{t}{W}_{t}^{T}\right)}^{1}$ are defined for the ’time step’ t. At each time step the dictionary D_{t1} and the C matrix are updated so that they obey the least squares solution ${D}_{t}=\left({X}_{t}{W}_{t}^{T}\right){\left({W}_{t}{W}_{t}^{T}\right)}^{1}$. The matrix inversion lemma (Woodbury matrix identity) is applied on C_{ t } to get the following update rules:
where u = C_{t1}w_{ t } and $\alpha =1/(1+{w}_{t}^{T}u)$, r_{ t } = x_{ t }  D_{t1}w_{ t } is the approximation error.
With the inclusion of an adaptive forgetting factor λ_{ t } in RLSDLA, the update equation 8 is changed to:
where $u=\left({\lambda}_{t}^{1}{C}_{t1}\right){w}_{t}$. The remaining equations remain the same. Due to introduction of the forgetting factor, RLSDLA has good converging properties and is less dependent on the initial dictionary.
Abbreviations
 CMR:

Cardiac magnetic resonance
 MI:

Myocardial infarction
 LG:

Late gadolinium
 ICD:

Implantable cardioverter defibrillator
 SVM:

Support vector machine
 PDFs:

Probability density functions
 ML:

Maximum likelihood
 DC:

Direct current
 ORMP:

Order recursive matching pursuit
 RLSDLA:

Recursive least squares dictionary learning algorithm
 NPhard:

Nondeterministic polynomialtime hard
 FTCM:

Frame texture classification method
 VT:

Ventricular tachycardia
 DICOM:

Digital imaging and communications in medicine
 ROC:

Receiver operating characteristic.
References
 1.
Roes SD, et al.: Infarct tissue heterogeneity assessed with contrastenhanced MRI predicts spontaneous ventricular arrhythmia in patients with ischemic cardiomyopathy and implantable cardioverterdefibrillator. Circ Cardiovasc Imaging 2009, 2: 183–190. 10.1161/CIRCIMAGING.108.826529
 2.
Schmidt A, et al.: Infarct tissue heterogeneity by magnetic resonance imaging identifies enhanced cardiac arrhythmia susceptibility in patients with left ventricular dysfunction. Circulation 2007, 115: 2006–2014. 10.1161/CIRCULATIONAHA.106.653568
 3.
Yan A, et al.: Characterization of periinfarct zone by contrastenchanced cardiac magentic resonance imaging is a powerful predictor of postmyocardial infarction mortality. Circulation 2006, 114: 32–39. 10.1161/CIRCULATIONAHA.106.613414
 4.
Woie L, et al.: The heart rate of ventricular tachycardia following an old myocardial infarction is inversely related to the size of scarring. Europace 2011, 13(6):864–868. 10.1093/europace/euq466
 5.
Leyva F, et al.: Cardiac resynchronization therapy guided by late gadoliniumenhancement cardiovascular magnetic resonance. Cardiovasc Magn Reson 2011., 13(29):
 6.
Shokrollahi E, et al.: Invivo MRI and invivo electroanatomical voltage map characteristics of infarct heterogeneity in a swine model. In Pro. EMBC conference. Boston: IEEE; 2011:2792–2795.
 7.
Dikici E, et al.: Quantification of delayed Enchancement MR images. In Pro. MICCAI 2004. Berlin, Heidelberg: SpringerLink; 2004:250–257.
 8.
Tao Q, et al.: Automated segmentation of myocardial scar in late enhancement MRI using combined intensity and spatial information. Magn Reson Med 2010, 64: 586–594.
 9.
Corsi C, et al.: Automatic quantification of cardiac scar extent from late gadolinium enhancement MRI. In Proc. Computing in Cardiology. China: Hangzhou; 2011. http://www.cinc.org/2011/preprints/205.pdf
 10.
Ørn S, et al.: Effect of left ventricular scar size, location, and transmurality on left ventricular remodeling with healed myocardial Infarction. Am J Cardiol 2007, 99(8):1109–1114. 10.1016/j.amjcard.2006.11.059
 11.
Engan K, et al.: Exploratory data analysis of image texture and statiscal features on myocardium and infarction areas in cardiac magnetic resonance images. In Pro. EMBS Conference 2010. Buenos Aires: IEEE; 2010:5728–5731.
 12.
Kotu L, et al.: Texture classification of scarred and nonscarred myocardium in cardiac MRI using learned dictionaries. In Pro. ICIP Conference 2011. Brussels, Belgium: IEEE; 2011:65–68.
 13.
Duda R, et al.: Pattern Classification. New York: John Wiley & Sons; 2001.
 14.
Theodoridis S, Koutroumbas K: Pattern Recognition. London: Academic Press; 2008.
 15.
Sayood K: Data Compression. London: Academic Press; 2000.
 16.
GharaviAlkhansari, et al.: A fast orthogonal matching pursuit algorithm. In Pro. ICASSP 1998. Seattle: IEEE; 1998:1389–1392.
 17.
Skretting K, et al.: Recursive least squares dictionary learning algorithm. In Signal Processing. Piscataway: IEEE; 2010:2121–2130.
 18.
Skretting K, et al.: Texture classification using sparse framebased representations. EURASIP Appl Signal Process 2006, 2006: 1–11. Heidelberg, Germany
 19.
Mairal J, et al.: Discriminative learned dictionaries for local image analysis. In Pro. Computer Vision and Pattern Recognition. Anchorage, AK: IEEE; 2008:1–8.
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
LPK carried out the feature extraction (DC and texture), probability mapping experiments and drafted the manuscript. KE and KS helped with the dictionary learning and sparse representation. TE participated in probability mapping and designed and performed the cardiac segment experiments. FM collaborated in finding the cardiac segments with out using the probability mapping. SØ and LW were responsible for the clinical studies including patient recruitment and acquiring of CMR images. They performed the segmentation of the healthy and scarred myocardium and contributed to the drafting and revision of the clinical and CMR related parts of the manuscript. KE, TE and KS also helped in drafting the manuscript. All the authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Kotu, L.P., Engan, K., Skretting, K. et al. Probability mapping of scarred myocardium using texture and intensity features in CMR images. BioMed Eng OnLine 12, 91 (2013). https://doi.org/10.1186/1475925X1291
Received:
Accepted:
Published:
Keywords
 Cardiac Magnetic Resonance
 Texture Feature
 Implantable Cardioverter Defibrillator
 Sparse Representation
 Cardiac Magnetic Resonance Image