The analysis of blood flow patterns and the interaction between salient topological flow features and cardiovascular structure plays an important role in the study of cardiovascular function. Flow velocity images acquired by Magnetic Resonance (MR) velocity imaging are generally subject to noise that are intrinsic to system hardware setup and those specific to patient movement in relation to imaging sequence designs. To improve the accuracy of the quantitative analysis of the evolution of topological flow features, it is essential to restore the original flow fields so that the associated critical points can be more reliably detected. In this study, we propose a total variation based variational method for the restoration of flow vector fields. The method is formulated as a constrained optimisation problem by minimizing the total variation energy of the normalized velocity field subject to a constraint that depends on the noise level. The effectiveness of this restoration method great...