Data sampling and constrained reconstruction for high-dimensional MRI