A method for diffusion magnetic resonance imaging may be provided. The method may comprise steps of performing sampling on an object at N diffusion weighted directions to acquire undersampled but complementary k-space data, combining the complementary data from different directions to obtain a full sampled k-space data, performing initial reconstruction based on common information among images at the N diffusion weighted directions, and performing joint iterative regularized reconstruction to k-space data in all diffusion weighted directions based on the initial images to obtain the desired final images. Due to the utilization of common information at the N diffusion weighted directions, the acquisition efficiency may be enhanced and image resolution and SNR of acquired images may be improved accordingly.