The present invention provides a technique for improving precision of a magnetic susceptibility by reducing and calculating an artifact and noise with a simple process in the calculation of magnetic susceptibility by using an MRI. In this technique, firstly, a magnetic field distribution is obtained through a phase image from a complex image obtained with the MRI. Then, a smoothing process is repeatedly performed on the magnetic susceptibility distribution calculated from the magnetic field distribution under restricted conditions based on a relational expression between the magnetic field and the magnetic susceptibility. In detail, the magnetic susceptibility distribution is calculated by a minimization process for repetitively calculating an error function defined by a preset initial magnetic susceptibility and magnetic distribution to minimize the error. At this time, an update magnetic susceptibility distribution, which is calculated on the basis of the error function during the minimization process, is smoothed every time the update magnetic susceptibility distribution is calculated.La présente invention concerne une technique permettant d'améliorer la précision d'une susceptibilité magnétique par réduction et calcul d'un artefact et du bruit par un procédé simple dans le calcul de la susceptibilité magnétique en utilisant une IRM. Dans cette technique, premièrement, une distribution de champ magnétique est obtenue par l'intermédiaire d'une image de phase dérivée d'une image complexe obtenue avec l'IRM. Ensuite, un processus de lissage est mis en œuvre à plusieurs reprises sur la distribution de susceptibilité magnétique calculée à partir de la distribution de champ magnétique dans des conditions restreintes sur la base d'une expression relationnelle entre le champ magnétique et la susceptibilité magnétique. Plus précisément, la distribution de susceptibilité magnétique est calculée par un processus de minimisation pour calculer de manière répétitive une foncti
