A method for rapid computation of three-dimensional displacement and Lagrange strain in a high resolution filed of phase data obtained with Displacement Encoding with Stimulated Echoes (DENSE) in magnetic resonance images of the myocardium. The method includes semi-automated segmentation of a region of a heart, phase unwrapping the images in three dimensions, and a custom radial point interpolation method (RPIM). The RPIM is a meshfree numerical analysis method that uses radial basis functions and polynomial functions to calculate the Lagrange strain of DENSE displacement data acquired from the myocardium.