Efficient characterization of inhomogeneity in contraction strain pattern

Cardiac dyssynchrony often accompanies patients with heart failure (HF) and can lead to an increase in mortality rate. Cardiac resynchronization therapy (CRT) has been shown to provide substantial benefits to the HF population with ventricular dyssynchrony; however, there still exists a group of patients who do not respond to this treatment. In order to better understand patient response to CRT, it is necessary to quantitatively characterize both electrical and mechanical dyssynchrony. The quantification of mechanical dyssynchrony via characterization of contraction strain field inhomogeneity is the focus of this modeling investigation. Raw data from a 3D finite element (FE) model were received from Roy Kerckhoffs et al. and analyzed in MATLAB. The FE model consisted of canine left and right ventricles coupled to a closed circulation with the effects of the pericardium acting as a pressure on the epicardial surface. For each of three simulations (normal synchronous, SYNC, right ventricular apical pacing, RVA, and left ventricular free wall pacing, LVFW) the Gauss point locations and values were used to generate lookup tables (LUTs) with each entry representing a location in the heart. In essence, we employed piecewise cubic interpolation to generate a fine point cloud (LUTs) from a course point cloud (Gauss points). Strain was calculated in the fiber direction and was then displayed in multiple ways to better characterize strain inhomogeneity. By plotting average strain and standard deviation over time, the point of maximum contraction and the point of maximal inhomogeneity were found for each simulation. Strain values were organized into seven strain bins to show operative strain ranges and extent of inhomogeneity throughout the heart wall. In order to visualize strain propagation, magnitude, and inhomogeneity over time, we created 2D area maps displaying strain over the entire cardiac cycle. To visualize spatial strain distribution at the time point of maximum inhomogeneity, a 3D point cloud was created for each simulation, and a CURE index was calculated. We found that both the RVA and LFVW simulations took longer to reach maximum contraction than the SYNC simulation, while also exhibiting larger disparities in strain values during contraction. Strain in the hoop direction was also analyzed and was found to be similar to the fiber strain results. It was found that our method of analyzing contraction strain pattern yielded more detailed spacial and temporal information about fiber strain in the heart over the cardiac cycle than the more conventional CURE index method. We also observed that our method of strain binning aids in visualization of the strain fields, and in particular, the separation of the mass points into separate images associated with each strain bin allows the strain pattern to be explicitly compartmentalized.