METHOD AND APPARATUS FOR DETERMINATION OF PHASE DISTRIBUTIONS IN MAGNETIC RESONANCE IMAGING

Abstract

In a method and apparatus for determination of phase distributions in MR imaging, a measured phase distribution of the region of interest is combined with at least one second phase value to form a combination-phase distribution, wherein the phase values of the combination-phase distribution are restricted to a defined presentation interval. A correction-phase distribution is generated, based on a known magnetic field distribution in the region of interest. The phase values thereof are not restricted to the defined presentation interval. A corrected combination-phase distribution is generated using the correction-phase distribution and the combination-phase distribution, in which the phase values are restricted to the defined presentation interval. An absolute combination-phase distribution is generated from the corrected combination-phase distribution, in which the phase values are not restricted to the defined presentation interval.

Claims

1. A method for generating magnetic resonance (MR) image data representing a region of interest of an examination subject, said method comprising: operating an MR data acquisition scanner in order to obtain a first measured phase distribution of MR signals acquired from said region of interest; providing said first measured phase distribution of the region of interest to a computer and providing said computer with an input designating at least one second phase value of the region of interest; in said computer, combining said first measured phase distribution and said at least one second phase value in order to generate a combination-phase distribution, wherein phase values of the combination-phase distribution are restricted to a predetermined presentation interval; in said computer, generating a correction-phase distribution, based on a known magnetic field distribution in said scanner in said region of interest, wherein phase values of said correction-phase distribution are not restricted to said predetermined presentation interval; in said computer, generating a corrected combination-phase distribution using the correction-phase distribution and the combination-phase distribution, wherein phase values of the corrected combination-phase distribution are restricted to the predetermined presentation interval; in said computer, generating an absolute combination-phase distribution from the corrected combination-phase distribution, wherein phase values of the absolute combination-phase distribution are not restricted to said predetermined presentation interval; and in said computer, using said phase values of said absolute combination-phase distribution to reconstruct image data from said MR signals, and using said image data to present a visualization of said region of interest at a display screen.

2. A method as claimed in claim 1 comprising generating a corrected absolute combination-phase distribution using the correction-phase distribution and the absolute combination-phase distribution.

3. A method as claimed in claim 1, comprising determining the correction-phase distribution .sub.korr(r) according to .sub.korr(r)=B.sub.0(r)T.sub.eff, wherein is the gyromagnetic ratio, B.sub.0(r) is a magnetic field distribution in the region of interest (21), r is the position vector in the region of interest (21), and T.sub.eff is an effective evolution time.

4. A method as claimed in claim 1 comprising generating said absolute combination-phase distribution based on a phase-unwrapping technique.

5. A method as claimed in claim 1 comprising combining said first measured phase distribution and said at least one second measured phase value as a determination of a phase-difference distribution.

6. A method as claimed in claim 1 comprising acquiring said first measured phase distribution by measurement of a phase distribution of an MR image of the region of interest after a first evolution period of nuclear spins in the region of interest.

7. A method as claimed in claim 6 comprising providing said at least one second phase value as a second measured phase distribution of the region of interest.

8. A method as claimed in claim 7 comprising providing the second measured phase distribution as a measurement of a phase distribution of an MR image after a second evolution period of said nuclear spins that is different from said first evolution period.

9. A method as claimed in claim 1 comprising providing said at least one second phase value as a definition of a global phase value.

10. A magnetic resonance (MR) apparatus comprising: an MR data acquisition scanner; a computer configured to operate said MR data acquisition scanner in order to obtain a first measured phase distribution of MR signals acquired from said region of interest; said computer being provided with an input designating at least one second phase value of the region of interest; said computer being configured to combine said first measured phase distribution and said at least one second phase value in order to generate a combination-phase distribution, wherein phase values of the combination-phase distribution are restricted to a predetermined presentation interval; said computer being configured to generate a correction-phase distribution, based on a known magnetic field distribution in said scanner in said region of interest, wherein phase values of said correction-phase distribution are not restricted to said predetermined presentation interval; said computer being configured to generate a corrected combination-phase distribution using the correction-phase distribution and the combination-phase distribution, wherein phase values of the corrected combination-phase distribution are restricted to the predetermined presentation interval; said computer being configured to generate an absolute combination-phase distribution from the corrected combination-phase distribution, wherein phase values of the absolute combination-phase distribution are not restricted to said predetermined presentation interval; and said computer being configured to use said phase values of said absolute combination-phase distribution to reconstruct image data from said MR signals, and to use said image data to present a visualization of said region of interest at a display screen.

11. A non-transitory, computer-readable data storage medium encoded with programming instructions, said storage medium being loaded into a computer system of a magnetic resonance (MR) apparatus comprising an MR data acquisition scanner, said programming instructions causing said computer system to: operate said MR data acquisition scanner in order to obtain a first measured phase distribution of MR signals acquired from said region of interest; receive an input designating at least one second phase value of the region of interest; combine said first measured phase distribution and said at least one second phase value in order to generate a combination-phase distribution, wherein phase values of the combination-phase distribution are restricted to a predetermined presentation interval; generate a correction-phase distribution, based on a known magnetic field distribution in said scanner in said region of interest, wherein phase values of said correction-phase distribution are not restricted to said predetermined presentation interval; generate a corrected combination-phase distribution using the correction-phase distribution and the combination-phase distribution, wherein phase values of the corrected combination-phase distribution are restricted to the predetermined presentation interval; generate an absolute combination-phase distribution from the corrected combination-phase distribution, wherein phase values of the absolute combination-phase distribution are not restricted to said predetermined presentation interval; and use said phase values of said absolute combination-phase distribution to reconstruct image data from said MR signals, and use said image data to present a visualization of said region of interest at a display screen.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0023] FIG. 1 schematically illustrates an MR system with which the method for the determination of a phase distribution of a region of interest can be carried out according to the invention.

[0024] FIG. 2a shows a two-dimensional representation of a magnetic field distribution B.sub.M in a region of interest, which is based on system-specific circumstances of the MR system.

[0025] FIG. 2b shows a line graph of the magnetic field distribution shown in FIG. 2a along a central line of the two-dimensional representation in FIG. 2a.

[0026] FIG. 3a shows a two-dimensional representation of a magnetic field distribution B.sub.P induced by an object in the region of interest.

[0027] FIG. 3b shows a line graph of the magnetic field distribution shown in FIG. 3a along a central line of the two-dimensional representation in FIG. 3a.

[0028] FIG. 4a shows a two-dimensional representation of a superimposition of the magnetic field distributions B.sub.0=B.sub.M+B.sub.P in FIGS. 2a and 3a.

[0029] FIG. 4b shows a line graph of the magnetic field distribution shown in FIG. 4a along a central line of the two-dimensional representation in FIG. 4a.

[0030] FIG. 5a shows a two-dimensional representation of a phase-difference distribution (r) of the region of interest, which is based on the magnetic field distribution shown in FIGS. 4a and 4b and wherein the phase values thereof are restricted to the interval [; +].

[0031] FIG. 5b shows a line graph of the phase-difference distribution (r) shown in FIG. 5a along a central line of the two-dimensional representation in FIG. 5a.

[0032] FIG. 6a shows a two-dimensional representation of a corrected phase-difference distribution of the region of interest, which is based on the phase-difference distribution shown in FIG. 5a and wherein the phase values thereof are restricted to the interval [; +].

[0033] FIG. 6b shows a line graph of the corrected phase-difference distribution of the region of interest shown in FIG. 6a along a central line of the two-dimensional representation in FIG. 6a.

[0034] FIG. 7 shows a two-dimensional representation of an absolute phase-difference distribution *(r) of the region of interest, which was generated by a phase-unwrapping technique from the corrected phase-difference distribution in FIG. 6 and wherein the phase values thereof are not restricted to the interval [; +].

[0035] FIG. 8 shows a two-dimensional representation of a spatial magnetic field distribution of the B.sub.0 field, which was calculated from the absolute phase-difference distribution in FIG. 7.

[0036] FIG. 9 is a flowchart with steps for carrying out the method for the determination of a phase distribution of a region of interest according to an exemplary embodiment of the invention.

DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0037] FIG. 1 schematically shows an MR system with which the method for the determination of a phase distribution of a region of interest 21 according to the invention can be carried out.

[0038] A person 12 under examination, or more generally an object under examination, is in moved into the tunnel of the system, wherein a region of interest 21 is examined by the MR system. The magnetic resonance system has a scanner 10 with a basic field magnet that generates a basic magnetic field B0. The person 12 under examination, arranged on a bed 11, is moved into the center of the scanner 10 in order to record spatially encoded magnetic signals from the region of interest 21. Radiation of radio-frequency (RF) pulse sequences and switching of magnetic field gradients causes certain nuclear spins in the person 12 to be deflected from the magnetization generated by the basic field B0, out of the equilibrium position. Currents induced in RF reception coils, as the excited spins return to the equilibrium position, are detected as magnetic resonance signals, which are then transformed into image data. The general mode of operation for creating MR images and the detection of magnetic resonance signals are known to those skilled in the art, so a more detailed description is not necessary herein.

[0039] The magnetic resonance system further has a control computer 13 that controls the MR scanner 10. The central control computer 13 is designed to implement the method described below for the determination of a phase distribution of a region of interest 21, and includes a gradient controller 14 that controls switching the magnetic field gradients, and an RF controller 15 that controls radiating the RF pulses for deflecting the nuclear spins from the equilibrium position. A memory 16 stores the imaging sequences required for recording the MR images and all programs required for the operation of the MR system. A recording controller 17 controls the image recording and hence, dependent on the chosen imaging sequences, the sequence of the magnetic field gradients and RF pulses and the receiving intervals of MR signals. Thus, the recording controller 17 also controls the gradient controller 14 and the RF controller 15. In a computing processor 20, MR images are calculated (reconstructed) and visualized on a display 18. An operator can operate the MR system via an input unit 19 of the MR system. The memory 16 can have imaging sequences and program modules, which carry out the method according to the invention when one of the modules shown is executed in the computing processor 20. The memory 16 thus stores control information (code) that can be executed by the control computer 13. The recording controller 17 is also designed such that it is able to assist in the method described below for the determination of a phase distribution of the region of interest 21.

[0040] According to the invention, the MR system in FIG. 1 is designed such that, when the control information is executed in the MR control computer 13, it determines at least one phase distribution of the region of interest 21 in order to generate an MR image of the region of interest 21 and, for the generation of an absolute phase combination distribution, the phase values thereof are not restricted to a defined presentation interval and take account of a known magnetic field distribution in the region of interest 21.

[0041] The basis of the method according to the invention is the separation of magnetically induced phase differences, in particular phase differences induced by fixed components (e.g. RF whole body coil, covering, gradient coil, etc.) with known magnetic field distributions and phase differences induced by non-fixed components (e.g. patient bench, RF local coils, positioning aids) with known positions and magnetic field distributions. The magnetically induced phase difference can be derived from a known spatial magnetic field distribution B.sub.M(r), which is, for example, determined on the installation of the MR system, during regular service procedures or even before delivery as a design parameter. The spatial distribution of the desired signal-phase difference (r)=(r).sub.ref can be derived from the complex MR signal with amplitude A at the location r:

S(r)=A(r)exp(i(r))

[0042] Wherein the following relationship also applies:

(r)=.sub.0(r)T.sub.eff+.sub.acq(r)+.sub.coil(r)=B.sub.0(r)T.sub.eff+.sub.acq(r)+.sub.coil(r)

[0043] Here .sub.0 (r) designates the deviation of the local precession frequency from a reference value (generally the current center frequency of the RF receiving system), the gyromagnetic ratio and T.sub.eff an effective evolution time that is dependent upon the measurement. For a gradient echo measurement, T.sub.eff=TE applies, for example, wherein the echo time TE defines the time interval between an RF excitation pulse and the time of data recording. At this point, the local field deviation B.sub.0(r) contains the magnetic field distribution B.sub.M(r) and object-induced magnetic field variations B.sub.P(r), but also further influencing variables such as the effective change in the precession frequency due to chemical displacement B.sub.CS(r). For example, in the case of MR elastography, the signal can be impressed in a targeted manner with a phase term .sub.acq(r) independent of T.sub.eff and which provides the desired phase contrast. Local coils and RF transmission effects can cause an additional phase term .sub.coil(r) to occur, which can be determined and taken into account by reference measurements.

[0044] Particularly at the edge of the imaging region, the signal phase (r) has rapid spatial variations as a result of the contribution B.sub.M(r)T.sub.eff. The method according to the invention now provides for the precise separation of this known contribution before the use of phase-unwrapping methods or filter techniques. Hence, the latter work on a benign (less strongly spatially fluctuating)

.sub.eff(r)=(r)B.sub.M(r)T.sub.eff.

[0045] This enables the ascertainment of the sought-for spatial parameter distributions with improved robustness or precision. If necessaryfor example for the determination of B0 mapsthe separated magnetic component can be re-added to the phase difference at the end of the calculation. Since the magnetic field is known as an absolute value, this can take place without restriction of the value range of the phase.

[0046] The reference phase .sub.ref can be defined arbitrarily or determined by the phase at a position in the image (or a mean over a defined range). The reference phase may also relate to a further measurement with changed settings. For the ascertainment of B0 maps, it is, for example, possible to perform two measurements with different effective evolution times T.sub.eff,1 and T.sub.eff,2. The signal phase difference is then obtained as (r)=.sub.1(r).sub.2(r), wherein the following relationships apply:

[00001]$S_1\left(r\right)=A_1\left(r\right).Math.\exp \left(i.Math..Math.{}_1\left(r\right)\right),.Math.S_2\left(r\right)=A_2\left(r\right).Math.\exp \left(i.Math..Math.{}_2\left(r\right)\right),.Math.\begin{array}{c}{}_1\left(r\right)=.Math..Math..Math.{}_0\left(r\right).Math.T_{\mathrm{eff},1}+{}_{\mathrm{acq},1}\left(r\right)+{}_{\mathrm{coil},1}\left(r\right)\\ =.Math..Math..Math..Math..Math.B_0\left(r\right).Math.T_{\mathrm{eff},1}+{}_{\mathrm{acq},1}\left(r\right)+{}_{\mathrm{coil},1}\left(r\right),\end{array}$$\begin{array}{c}{}_2\left(r\right)=.Math..Math..Math.{}_0\left(r\right).Math.T_{\mathrm{eff},2}+{}_{\mathrm{acq},2}\left(r\right)+{}_{\mathrm{coil},2}\left(r\right)\\ =.Math..Math..Math..Math..Math.B_0\left(r\right).Math.T_{\mathrm{eff},2}+{}_{\mathrm{acq},2}\left(r\right)+{}_{\mathrm{coil},2}\left(r\right).\end{array}$

[0047] If the two measurements are performed with an identical coil configuration .sub.coil,1(r)a=.sub.coil,2(r) and, apart from the evolution time, an identical course .sub.acq,1(r)=.sub.acq,2(r), the sought-for signal-phase difference is obtained from the formula

(r)=B.sub.0(r)(T.sub.eff,1T.sub.eff,2).

[0048] Here once again, the known, rapidly spatially varying contribution due to the basic magnetic field can be separated by means of the formula

.sub.eff(r)=B.sub.0(r)(T.sub.eff,1T.sub.eff,2)B.sub.M(r)(T.sub.eff,1T.sub.eff,2),

[0049] and thus a more benign spatial phase distribution is obtained for the further data processing.

[0050] To summarize, the method according to the invention reduces the spatial variation of a signal phase thus permitting a more robust and more precise determination of relative phase relationships in the image.

[0051] As an example, the method is described for recording absolute B0 maps. The aim of recording B0 maps is to determine B.sub.0(r)=B.sub.M(r)+B.sub.P(r), wherein, for purposes of simplicity, no further influencing variables (e.g. B.sub.SC(r))) are taken into account. To this end, at least two measurements with different effective evolution times T.sub.eff,1 and T.sub.eff,2 are recorded and the spatial distribution of the relative signal phases is determined by the formula

(r)=B.sub.0(r)(T.sub.eff,1T.sub.eff,2)=(B.sub.M(r)+B.sub.P(r))(T.sub.eff,1T.sub.eff,2).

[0052] In the prior art, the desired B0 information would be determined from the following relationship:

B.sub.0(r)=(r)/((T.sub.eff,1T.sub.eff,2)).

[0053] Since the phase difference only has a value range of and +, the directly accessible value range for the field map is restricted to +/((T.sub.eff,1T.sub.eff,2)). For example, measuring with T.sub.eff,1=4.8 ms and T.sub.eff,2=9.6 ms produces a value range of about +/2.4 T. With a basic field strength of 1.5 T, this corresponds to a relative field change of +/1.6 ppm. At the edge of the imaging region, the inhomogeneities of the basic field magnet can very possibly be a few tens of ppm. As a result, this order of magnitude can no longer be depicted directly with the measurement, so unique determination of the B0 amplitude is not directly possible.

[0054] The following two-dimensional example assumes a magnetically induced field distribution according to a spherical harmonic function of the type B.sub.M(x, y)=B.sub.M,0*cos(5*)*(r/r.sub.0).sup.5, wherein is the angle in the plane, r is the distance to the center, r.sub.0 is a reference radius and B.sub.M,0 is a scaling factor. This is superimposed with an object-induced field distribution according to a quadratic function B.sub.P(x, y)=B.sub.P,0*(r/r.sub.0).sup.2.

[0055] The actual physical fields have the following forms, wherein the units used for representation were selected arbitrarily and are identical in all the figures.

[0056] FIG. 2a shows a two-dimensional representation of a magnetic field distribution B.sub.M in a region of interest, which is based on system-specific circumstances of the MR system, such as, for example, the design of the basic field magnet. The edge regions of the circular representation of the magnetic field distribution B.sub.M shown have strong spatial variation with a high spatial frequency of the magnetic field strength.

[0057] FIG. 2b shows a line graph of the magnetic field distribution shown in FIG. 2a along a central line of the two-dimensional representation in FIG. 2a.

[0058] FIG. 3a shows a two-dimensional representation of a magnetic field distribution B.sub.P induced by an object in a region of interest. The edge regions of the circular representation of the magnetic field distribution B.sub.P shown have uniform magnetic field strength.

[0059] FIG. 3b shows a line graph of the magnetic field distribution shown in FIG. 3a along a central line of the two-dimensional representation in FIG. 3a.

[0060] FIG. 4a shows a two-dimensional representation of a superimposition of the magnetic field distributions B.sub.0=B.sub.M+B.sub.P in FIGS. 2a and 3a. In the superimposed magnetic field distribution B.sub.0=B.sub.M+B.sub.P, the influence of the magnetic field distribution B.sub.M predominates so that the edge regions of the circular representation also have strong spatial variation with a high spatial frequency of the magnetic field strength.

[0061] FIG. 4b shows a line graph of the magnetic field distribution shown in FIG. 4a along a central line of the two-dimensional representation in FIG. 4a.

[0062] A measurement of phase distribution by the MR system can only detect the superimposed field distribution B.sub.0(x, y) via the determination of a phase difference (r) and supplies data with a restricted value range.

[0063] Therefore, evaluation routines see the following spatial phase distribution, which is depicted in FIG. 5a.

[0064] FIG. 5a shows a two-dimensional representation of a phase-difference distribution (r) of the region of interest, which is based on the magnetic field distribution shown in FIGS. 4a and 4b and wherein the phase values thereof are restricted to the interval [; +].

[0065] FIG. 5b shows a line graph of the phase-difference distribution (r) shown in FIG. 5a along a central line of the two-dimensional representation in FIG. 5a.

[0066] A number of phase jumps can be identified in the edge regions of the corrected phase-difference distribution in FIGS. 5a and 5b. These phase jumps are caused by the restriction of the phase values to the interval [; +].

[0067] The problem of the ambiguous values in the edge region is clearly identifiable. In the prior art, these rapid spatial variations had to be detected by means of suitable image processing methods (e.g. phase-unwrapping) and assigned correctly andin particular in the case of noisy datathis can result in errors.

[0068] According to the invention, first, the known contribution .sub.M(r)=B.sub.M(r)(T.sub.eff,1T.sub.eff,2) of the magnet is separated from the determined phase difference (r). Herein, the magnetically induced phase is restricted to the value range between and +. This achieves a more benign spatial phase distribution, which is substantially object-induced, as shown in FIG. 6.

[0069] FIG. 6a shows a two-dimensional representation of a corrected phase-difference distribution of the region of interest, which is based on the phase-difference distribution shown in FIG. 5a and wherein the phase values thereof are restricted to the interval [; +].

[0070] FIG. 6b shows a line graph of the corrected phase-difference distribution shown in FIG. 6a of the region of interest along a central line of the two-dimensional representation in FIG. 6a.

[0071] It can be clearly seen that the problem of ambiguous values in the edge region is much less pronounced. As a result, it is possible to assign an absolute phase differencefor example by means of phase-unwrapping methodswith a high degree of probability even in the case of noisy data. In order to generate a continuous phase-difference distribution from the corrected phase-difference distribution in FIG. 6 without restriction to the interval [; +], in a next step a phase-unwrapping technique is applied to the corrected phase-difference distribution. It is also possible to use filter methods to generate a continuous phase-difference distribution, as is known from the prior art. As a result, the object-induced phase difference *(r) is obtained without restriction of the value range, as shown in FIG. 7.

[0072] FIG. 7 shows a two-dimensional representation of an absolute phase-difference distribution *(r) of the region of interest, which was generated by a phase-unwrapping technique from the corrected phase-difference distribution in FIG. 6 and wherein the phase values thereof are not restricted to the interval [; +].

[0073] In the last step, the previously separated magnetic component is taken into account once again. Herein, it is essential that the value range of this contribution is not restricted since the absolute value of B.sub.M(r) is known a priori. In this way, the spatial distribution of the absolute phase differences is obtained directly and, derived therefrom, the sought-for spatial B0 distribution as in shown FIG. 8 is obtained from the formula:

*(r)=*.sub.P(r)+B.sub.M(r)(T.sub.eff,1T.sub.eff,2)

[0074] FIG. 8 shows a two-dimensional representation of a spatial magnetic field distribution of the B.sub.0 field, which was calculated from the absolute phase-difference distribution in FIG. 7.

[0075] FIG. 9 is a flowchart with steps for carrying out the method for the determination of a phase distribution of a region of interest according to an exemplary embodiment of the invention.

[0076] The method starts in step S40. In step S41, a first measured phase distribution of the region of interest is provided. In step S42, at least one second phase value is provided. In step S43, the first measured phase distribution and the at least one second phase value are combined with one another to generate a combination-phase distribution. Herein, the phase values of the combination-phase distribution are restricted to a defined presentation interval. In step S44, a correction-phase distribution is generated, which is based on the known magnetic field distribution in the region of interest. Herein, the phase values of the correction-phase distribution are not restricted to the defined presentation interval. In step S45, a corrected combination-phase distribution is generated using the correction-phase distribution and the combination-phase distribution, wherein the phase values of the corrected combination-phase distribution are restricted to the defined presentation interval. In step S46, an absolute combination-phase distribution is generated from the corrected combination-phase distribution, wherein the phase values of the absolute combination-phase distribution are not restricted to the defined presentation interval. The method ends in step S47, which includes using said phase values of said absolute combination-phase distribution to reconstruct image data from said MR signals, and using said image data to present a visualization of said region of interest at a display screen.

[0077] In an exemplary embodiment, when recording with a number of coil elements, the method according to the invention is first applied to the phase data of each individual element. This is following by a conventional combination taking account of the phase information.

[0078] In another exemplary embodiment, first the data from all coil elements is combined in the conventional way taking account of the phase information. This is followed by the application of the method according to the invention to the combination-phase data.

[0079] In a further exemplary embodiment, more than two measurements are performed with different effective evolution times, wherein the method according to the invention is applied to any combinations of signal phases.

[0080] In an additional exemplary embodiment, when performing more than one measurement with different effective evolution times, the method according to the invention is applied to both the individual measurements (i.e. before the combination of the signal phases) and to the combined measurements.

[0081] Due to the improved robustness to the strong magnetic field variations in the edge region, the method according to the invention enables high-precision results even with coarser voxels within the region of interest. This permits, for example, accelerated recording of a large spatial region with reduced resolution, for example for recording B0 maps during adjustment procedures.

[0082] In addition to the fixed basic field magnet, further elements can contribute to in local variations of the B0 field (and hence in phase variation). If the influence of these elements on the spatial magnetic field distribution and the position of the elements relative to the basic field magnet are known, these components can also be taken into account during the separation according to the invention. Examples of fixed elements include covering, an RF whole-body coil, gradient coils and shim coils. In the case of gradient coils and shim coils, the variable, known settings of the additional fields generated thereby are taken into account with linear geometry or also higher order geometry. Examples of non-fixed elements with a known position are a patient bench, positioning aids and local RF coils.

[0083] In summary, the method for the determination of a phase distribution of a region of interest is carried out. A combination of phase distributions of the region of interest, which were determined in order to generate an MR image of the region of interest by means of an MR system, is performed, wherein a known magnetic field distribution in the region of interest is taken into account. According to the invention a component of a priori known spatial variations of the signal phase of measured data are separated in order to a generate a more benign spatial phase distribution, which is substantially induced by an object.

[0084] According to the invention, the method utilizes previous knowledge of a spatial distribution of a basic magnetic field, in particular rapidly spatially varying contributions in an edge region of a region of interest due to system-induced known magnetic field variations, such as, for example, magnetic field variations induced by the design of the basic field magnet and thus reduces errors in the determination of phase differences in the edge region of the region of interest.

[0085] Although modifications and changes may be suggested by those skilled in the art, it is the intention of the Applicant to embody within the patent warranted hereon all changes and modifications as reasonably and properly come within the scope of the Applicant's contribution to the art.