Restoration of the velocity field of the heart from two-dimensional echocardiograms

A method to quantify the motion of the heart from digitized sequences of two-dimensional echocardiograms (2-D) echos was recently proposed. This method computes on every point of the 2-D echoes, the 2-D apparent velocity vector (or optical flow) which characterizes its interframe motion. However, further analysis is required to determine what part of this motion is due to translation, rotation, contraction, and deformation of the myocardium. A method to locally obtain this information is presented. The proposed method assumes that the interframe velocity field U(xy), V(x,y) can be locally described by linear equations in the form U(x,y)=a+Ax+By; V(x,y)=b+Cx+Dy. The additional constraint was introduced in the computation of the local velocity field by the method of projections onto convex sets. Since this constraint is only valid locally, the myocardium must be first divided into sectors and the velocity fields computed independently for each sector.