Ultrasound speckle tracking provides robust estimates of fine tissue displacements due to the analytic nature of echo data. Extensions to 2D have been proposed. Due to the relatively coarse lateral sampling, several solutions relied on lateral interpolation in order to achieve subsample accuracy. A multi-dimensional speckle tracking method (MDST) with subsample accuracy in all dimensions is used with two algorithms. The first algorithm is based on the gradient of the magnitude and the zero-phase crossing of 2D normalized complex correlation of the generalized analytic signal while the second algorithm is based on the polynomial fit. The generalization method utilizes the Riesz transform which is the vector extension of the Hilbert transform.