Although modern ultrasound acquisition systems allow recording of 3D echocardiographic images, tracking anatomical structures from them is still challenging. In addition, since these images are typically created from information obtained across several cardiac cycles, it is not yet possible to acquire high-quality 3D images from patients presenting varying heart rhythms. In this paper, we propose a method to estimate the motion field from multi-plane echocardiographic images of the left ventricle, which are acquired simultaneously during a single cardiac cycle. The method integrates tri-plane B-mode and tissue Doppler images acquired at different rotation angles around the long axis of the left ventricle. It uses a diffeomorphic continuous spatio-temporal transformation model with a spherical data representation for a better interpolation in the circumferential direction. This framework allows exploiting the spatial relation among the acquired planes. In addition, higher temporal resolution of the transformation in the beam direction is achieved by uncoupling the estimation of the different components of the velocity field. The method was validated using a realistic synthetic dataset including healthy and ischemic cases, obtaining errors of 0.14 ± 0.09 mm for displacements, 0.96 ± 1.03% for longitudinal strain and 3.94 ± 4.38% for radial strain estimation. In addition, the method was also demonstrated on a healthy volunteer and two patients with ischemia.