To generate a reconstructed MR image 36, iterative reconstruction of multi-gradient echo MR imaging data is performed. Iterative reconstruction uses a model 30 that links the multi-gradient echo MR imaging data and the reconstructed MR image. The model includes a parameterized magnetic field variation component 32. While performing the iterative reconstruction, the parameter of the parameterized field variation component of the model optimizes a cost function 40 that depends on the partial derivative of the reconstructed MR image with respect to the parameterization field variation component parameter of the model. Updated for. The image is further processed to generate an R2 * map 50, a SWI image 52, or a QMS map 54.再構成されたＭＲ画像３６を生成するために、複数勾配エコーＭＲ撮像データの反復再構成が実施される。反復再構成は、複数勾配エコーＭＲ撮像データと再構成されたＭＲ画像とをリンクするモデル３０を使用する。モデルは、パラメータ化磁場変動成分３２を含む。反復再構成を実施する間に、モデルのパラメータ化磁場変動成分のパラメータは、モデルのパラメータ化磁場変動成分のパラメータに関して再構成されたＭＲ画像の偏導関数に依存するコスト関数４０を最適化するために更新される。画像は、更に、Ｒ２＊マップ５０、ＳＷＩ画像５２、又はＱＭＳマップ５４を生成するために処理される。
