Myocardial velocities carry important diagnostic information in a range of cardiac diseases, and play an important role in diagnosing and grading left ventricular diastolic dysfunction. Tissue Phase Mapping (TPM) Magnetic Resonance Imaging (MRI) enables discrete sampling of the myocardium’s underlying smooth and continuous velocity field. This paper presents a post-processing framework for constructing a spatially and temporally smooth and continuous representation of the myocardium’s velocity field from TPM data. In the proposed scheme, the velocity field is represented through either linear or cubic B-spline basis functions. The framework facilitates both interpolation and noise reducing approximation. As a proof-of-concept, the framework was evaluated using artificially noisy (i.e., synthetic) velocity fields created by adding different levels of noise to an original TPM data. The framework’s ability to restore the original velocity field was investigated using Bland-Altman statistics. Moreover, we calculated myocardial material point trajectories through temporal integration of the original and synthetic fields. The effect of noise reduction on the calculated trajectories was investigated by assessing the distance between the start and end position of material points after one complete cardiac cycle (end point error). We found that the Bland-Altman limits of agreement between the original and the synthetic velocity fields were reduced after application of the framework. Furthermore, the integrated trajectories exhibited consistently lower end point error. These results suggest that the proposed method generates a realistic continuous representation of myocardial velocity fields from noisy and discrete TPM data. Linear B-splines resulted in narrower limits of agreement between the original and synthetic fields, compared to Cubic B-splines. The end point errors were also consistently lower for Linear B-splines than for cubic. Linear B-splines therefore appear to be more suitable for TPM data.