MR imaging can leverage a wide variety of intrinsic contrast mechanisms to provide detailed information regarding the anatomy, function, physiology, and metabolism of biological tissues. However, because of low sensitivity, many experiments that reveal higher-order structure and function have been limited due to inherent trade-offs between data acquisition time, signal-to-noise ratio, and resolution. This paper describes the further development of a statistical framework for MR image reconstruction which helps to mitigate these effects. Specifically, we advocate the collection of high-resolution multi-modal MR imaging data, and utilize the correlation between features in these different images to reduce noise while maintaining resolution. The proposed approach is illustrated with joint reconstruction of relaxometry and spectroscopic imaging data in a mouse model of stroke.