An iterative reconstruction is performed of multiple gradient echo MR imaging data to generate a reconstructed MR image (36). The iterative reconstruction uses a model (30) that links the MR imaging data and the reconstructed MR image. The model includes a parameterized magnetic field fluctuation component (32). During the performing of the iterative reconstruction, parameters of the parameterized magnetic field fluctuation component of the model are updated to optimize a cost function (40) dependent on partial derivatives of the reconstructed MR image with respect to the parameters of the parameterized magnetic field fluctuation component of the model. The image may be further processed to generate an R2* map (50), an SWI image (52), or a QSM map (54).L'invention concerne une reconstruction itérative effectuée à partir de multiples données d'imagerie RM à écho de gradient pour générer une image MR reconstruite (36). La reconstruction itérative utilise un modèle (30) qui relie les données d'imagerie RM et l'image RM reconstruite. Le modèle comprend une composante de fluctuation de champ magnétique (32) paramétré. Pendant la réalisation de la reconstruction itérative, des paramètres de la composante de fluctuation de champ magnétique paramétrée du modèle sont mis à jour pour optimiser une fonction de coût (40) en fonction des dérivées partielles de l'image RM reconstruite par rapport aux paramètres de la composante de fluctuation de champ magnétique paramétré du modèle. L'image peut en outre être traitée pour générer une carte R2* (50), une image SWI (52) ou une carte QSM (54).
