PROBLEM TO BE SOLVED: To provide a method of calculating a higher-SNR local frequency distribution in a shorter calculation time when a local magnetic field distribution resulting from magnetic susceptibility differences between living tissues is calculated using MRI.SOLUTION: Multiecho complex images measured at at least two or more different echo times using MRI are converted into low-resolution images. Global frequency distribution resulting from global magnetic field changes and offset phase distribution including reception and transmission phases are separated from a phase distribution of a low-resolution multiecho complex image. The calculated global frequency distribution and offset phase distribution are provided with higher resolution. A local frequency distribution of each echo is calculated from the measured multiecho complex image, high-resolution global frequency distribution, and high-resolution offset phase distribution. The local frequency distributions of the echoes are averaged on a weighted basis to calculate the final local frequency distribution.SELECTED DRAWING: Figure 6【課題】MRIを用いて、生体組織間の磁化率差に起因した局所磁場分布を算出するにあたり、演算時間が短く,SNRが高い局所周波数分布を算出する方法を提供する。【解決手段】MRIを用いて少なくとも2つ以上の異なるエコー時間で計測したマルチエコー複素画像を低解像度画像に変換する。低解像度マルチエコー複素画像の位相分布から、大域的な磁場変化に起因する大域的周波数分布と、受信や送信位相などを含むオフセット位相分布を分離する。算出した大域的周波数分布とオフセット位相分布を高解像度化する。計測したマルチエコー複素画像と高解像大域的周波数分布と高解像オフセット位相分布とから各エコーの局所周波数分布を算出する。各エコーの局所周波数分布を加重平均して最終的な局所周波数分布を算出する。【選択図】図6