Perhaps the most important aspect of the method presented here is its computational simplicity in producing measures of motion. The basic mathematical operations here include: the inverse Fourier transform to produce a complex image; computing the HARP images; and computing the derivatives of the HARP images. These computations are fast and in few minutes the resulting strain maps are produced. For example, for the canine heart with 20 time frames, it took only three min to compute the strain. Bearing in mind that the code was implemented using MATLAB, there is clearly potential for even faster computations