Tagged magnetic resonance imaging is a specially developed technique to noninvasively assess contractile function of the heart. Several methods have been developed to estimate myocardial deformation from tagged image data. Most of these methods do not explicitly impose a continuity constraint through time although myocardial motion is a continuous physical phenomenon. In this paper, we propose to model the spatio-temporal myocardial displacement field by a cosine series model fitted to the entire tagged dataset. The method has been implemented in two dimensions (2D)+time. Its accuracy was successively evaluated on actual tagged data and on a simulated two-dimensional (2D) moving heart model. The simulations show that an overall theoretical mean accuracy of 0.1 mm can be attained with adequate model orders. The influence of the tagging pattern was evaluated and computing time is provided as a function of the model complexity and data size. This method provides an analytical and hierarchical model of the 2D+time deformation inside the myocardium. It was applied to actual tagged data from a healthy subject and from a patient with ischemia. The results demonstrate the adequacy of the proposed model for this evaluation.