Improving motion estimation by accounting for local image distortion

Cardiac elastography is a useful diagnostic technique for detection of heart function abnormalities, based on analysis of echocardiograms. The analysis of the regional heart motion allows assessing the extent of myocardial ischemia and infarction. In this paper, a new two-stage algorithm for cardiac motion estimation is proposed, where the data is taken from a sequence of 2D echocardiograms. The method combines the advantages of block-matching and optical flow techniques. The first stage employs a standard block-matching algorithm (sum of absolute differences) to provide a displacement estimate with accuracy of up to one pixel. At the second stage, this estimate is corrected by estimating the parameters of a local image transform within a test window. The parameters of the image transform are estimated in the least-square sense. In order to account for typical heart motions, like contraction/expansion, translation and rotation, a local affine model is assumed within the test window. The accuracy of the new algorithm is evaluated using a sequence of 500 grayscale B-mode images, which are generated as distorted, but known copies of an original ROI, taken from a real echocardiogram. The accuracy of the motion estimation is expressed in terms of errors: maximum absolute error, root-mean-square error, average error and standard deviation. The errors of the proposed algorithm are compared with these of the known block-matching technique with cross-correlation and interpolation in the sub-pixel space. Statistical analysis of the errors shows that the proposed algorithm provides more accurate estimates of the heart motion than the cross-correlation technique with interpolation in the sub-pixel space.