1. A method of image reconstruction for a magnetic resonance imaging (MRI) system comprising a plurality of coils, the method comprising: obtaining k-space scan data captured by the MRI system, the k-space scan data being representative of an undersampled region over time; determining, with a processor, a respective coil sensitivity profile for the region for each coil of the plurality of coils; and iteratively reconstructing, with the processor, dynamic images for the region from the k-space scan data via an optimization of a minimization problem; wherein the minimization problem is a function of the determined coil sensitivity profiles and redundant Haar wavelet transforms of the dynamic images.