Correction of magnetic resonance images using simulated magnetic resonance images

A computer program product comprising machine executable instructions for execution by a processor controlling a medical imaging system, wherein execution of the machine executable instructions causes the processor to:

• receive measured magnetic resonance imaging data descriptive of a first region of interest of a subject;
• receive a B0 map, and at least one of a transverse relaxation map and a longitudinal relaxation map, preferably both a T1 map, a T2 map and/or a T2* map, and a magnetization map each descriptive of a second region of interest of the subject, wherein the first region of interest of the subject and the second region of interest of the subject both comprise an overlapping region of interest;
• receive pulse sequence commands configured for controlling a magnetic resonance imaging system to acquire the measured magnetic resonance imaging data;
• calculate a simulated magnetic resonance image of the overlapping region of interest using at least the B0 map, and at least one of a transverse relaxation map and a longitudinal relaxation map, preferably both the T1 map, the T2 map, and/or T2*map the magnetization map, and the pulse sequence commands as input to a Bloch equation model;
• and
• reconstruct a corrected magnetic resonance image from the measured magnetic resonance imaging data for the overlapping region of interest by solving an inverse problem, wherein the inverse problem comprises an optimization of a cost function and a regularization term formed from the simulated magnetic resonance image.

A method of operating a medical imaging system, wherein the method comprises:

• receiving measured magnetic resonance imaging data descriptive of a first region of interest of a subject;
• receiving B0 map, and at least one of a transverse relaxation map and a longitudinal relaxation map, preferably both a T1 map, a T2 map and/or a T2* map, and a magnetization map each descriptive of a second region of interest of the subject, wherein the first region of interest of the subject and the second region of interest of the subject both comprise an overlapping region of interest;
• receiving pulse sequence commands configured for controlling a magnetic resonance imaging system to acquire the measured magnetic resonance imaging data;
• calculating a simulated magnetic resonance image of the overlapping region of interest using at least the B0 map, and at least one of a transverse relaxation map and a longitudinal relaxation map, preferably both the T1 map, the T2 map, and/or a T2* map, the magnetization map, and the pulse sequence commands as input to a Bloch equation model;
• and
• reconstructing a corrected magnetic resonance image from the measured magnetic resonance imaging data for the overlapping region of interest by solving an inverse problem, wherein the inverse problem comprises an optimization of a cost function and a regularization term formed from the simulated magnetic resonance image.