Methods and systems for accelerated Phase-contrast magnetic resonance imaging (PC-MRI). The technique is based on Bayesian inference and provides for fast computation via an approximate message passing algorithm. The Bayesian formulation allows modeling and exploitation of the statistical relationships across space, time, and encodings in order to achieve reproducible estimation of flow from highly undersampled data.