本发明提供一种在使用MRI计算由生物体组织间的磁化率差所引起的局部磁场分布时，运算时间短且计算高SNR局部频率分布的方法。该方法中，将使用MRI在至少两个以上的不同回波时间测量到的多回波复数图像变换为低分辨率图像。从低分辨率多回波复数图像的相位分布中分离出由全局性磁场变化引起的全局性频率分布和包括接收和发送相位等的偏移相位分布。使计算出的全局性频率分布和偏移相位分布高分辨率化。根据测量到的多回波复数图像和高分辨率全局性频率分布和高分辨率偏移相位分布，来计算各回波的局部频率分布。对各回波的局部频率分布进行加权平均来计算最终的局部频率分布。