Cross-Reference to Related Applications

[0001] This application claims priority from U.S. Provisional Patent Application 63/116,551. filed November 20, 2020, which is incorporated herein by reference in its entirety.

Technical Field

[0002] The present invention relates to Nuclear Magnetic Resonance (NMR) and Magnetic Resonance Imaging (MRI), and more particularly to mapping the generated magnetic field (Bo); and correcting inhomogeneity in the generated magnetic field (Bo).

Background Art

[0003] Nuclear Magnetic Resonance (NMR) and Magnetic Resonance Imaging (MRI) utilize strong magnets to generate a homogeneous magnetic field (Bo) across objects ranging from small sample tubes to the human body (Fig. 1). Although the magnets used are very homogeneous, small imperfections in the magnet result in spatial variations in magnetic field strength on the order of parts per million. Additionally, the field becomes even more inhomogeneous when the sample (e.g. tubes or people) are inserted into them. For example, for the human head, the interfaces between air (e.g. sinuses), bone (skull) and the soft tissue (brain) distort the magnetic field locally and lead to variations (i.e. inhomogeneity) in the Bo field strength across the object. Since the measured frequency from a specific location varies linearly with Bo strength at that location, a range in frequencies will be detected across the object. This results in unwanted broadening of the signals from the sample, thereby lowering the signal to noise ratio. Additionally, in MRI, this inhomogeneity also results in geometric distortions (stretching/compression) in the shape of the object being imaged.

[0004] To overcome these effects, passive iron shims or active Bo correction coils, also known as Bo shims or simply “shims” may be incorporated into the hardware used within the bore of the magnet. Active shims are designed such that each coil generates a unique spatial variation in BO field, with its overall strength and sign determined by the amount of current applied to the individual coils. When applied to the object, these coils provide small Bo correction fields; either increasing or decreasing the local Bo field strength, so as to make the Bo field across the object more spatially homogeneous. Thus, when the correct set of currents are applied to the individual active shims, the shims generate a spatially vary ing Bo field distribution matching, but opposite in sign, to the distortions generated by the object and any residual imperfections in the magnet.

[0005] Placement of the shims and/ or adj ustment of the current applied to the active shims is optimally achieved by first accurately characterizing the spatial variation in Bo field across the object by mapping the Bo field. Once this has been accomplished, the amount of current used to drive each of the shim coils can be calculated using calibrated images of the specific Bo fields generated by the individual shim coils. However, this process is dependent upon an accurate mapping of the existing Bo field distribution. Further, for MRI Applications or high throughput NMR applications, it is important to keep the duration of these measurements as short as possible to maximize the available study time for the primary measurements. Thus methods which can enhance the accuracy of the measurement while minimizing the time required to make it are valuable tools for both NMR and MRI applications.

Summary of the Embodiments

[0006] The system and methods of the various embodiments of the invention described below achieve rapid and highly accurate Bo mapping that may be used, for example, to determine the placement of passive shims and/or shim currents for active shims for the purpose of reducing inhomogeneity when using Magnetic Resonance Imaging (MRI) systems. While the below described embodiments refer to MRI systems, it is to be understood that the invention is applicable to other imaging systems known in the art, such as, without limitation, Nuclear Magnetic Resonance (NMR) systems.

[0007] In accordance with an embodiment of the invention, a method is presented for mapping and correcting the inhomogeneity of a magnetic field within an object using an Magnetic Resonance Imaging (MRI) system (or NMR system) where there is a single dominant resonance. The method includes acquiring at least three MRI images, each at different echo times (TE). At least two ΔTE images (ΔTE i=i ... N) are generated based on the at least three MRI images, wherein the subscripts I = 1 . . . N refer to images with sequentially increasing ΔTE times. Aliasing in the ΔTE ₁ image is permitted. The ΔTE times of ΔTE ₁ and ΔTE ₂ are set such that the alias points at which wrapping occurs in ΔTE 1 does not overlap with the alias points of ΔTE ₂ . Each ΔTE images is unwrapped. A final Bo map is set to the unwrapped ΔTEN image.

[0008] In accordance with related embodiments of the invention, the method may further include determining shim currents and/or passive shims based on the final Bo map, the shim currents for reducing Bo inhomogeneity. The shim currents and/or passive shims may be applied to the sample to reduce Bo inhomogeneity in the sample. The shim currents and/or passive shims may generate a Bo field opposite in sign to the inhomogeneity of the magnetic field.

[0009] In accordance with another related embodiment of the invention, the at least two ΔTE images may differ in their echo time (TE) such that the difference in resulting phase (ΔΦ ) of any two ΔTE images varies according to ΔΦ = ΔTE ^. Δv, where Δv is the difference in frequency at a voxel in the image from a central frequency, and ΔTE is the difference in echo time between any pair of images.

[0010] In accordance with yet another embodiment of the invention, unwrapping each ΔTE image to form a final Bo map include: unwrapping aliasing in the ΔTE ₁ image based on information in the next longest ΔTE image ( the ΔTE ₂ image) , to form an unwrapped ΔTE ₁ image; and using the unwrapped ΔTE ₁ image to unwrap the next longest ΔTE image (the ΔTE ₂ image). Unwrapping aliasing in the ΔTE ₁ and ΔTE ₂ images may include calculating a logistic temporal unwrapping parameter LTUP, where LTUP images the fractional difference measured between two ΔTE images into a set of values describing the extent of aliasing. Unwrapping each ΔTE image to form a final Bo map may further include, for i=3 to N, using the previously unwrapped ΔTE i-i image to unwrap the ΔTE i image. Using the unwrapped ΔTE ₁ image to unwrap the next longest ΔTE image (the ΔTE ₂ image) may include using a numerical estimate of the difference between the unwrapped value in the ΔTE ₁ image with values representing different degrees of aliasing in the ΔTE ₂ .

[0011] In accordance with still another embodiment of the invention, each of the MRI images may be acquired using gradient echo, asymmetric spin echo, and/or stimulated echo acquisitions. Acquiring each ΔTE image may include a single pattern matching algorithm to simultaneously unwrap each ΔTE images. Acquiring the at least three MRI images may include using simultaneous multi-slice methods, K-space undersampling in the phase encoding direction, and/or use of 2 and 3 dimensional readout trajectories, to decrease acquisition times.

[0012] In accordance with another embodiment of the invention, a system is presented for mapping and correcting the inhomogeneity of the magnetic field within an object using Magnetic Resonance Imaging (MRI) where there is a single dominant resonance. The system includes an MRI imaging device for acquiring at least three MRI images, each at different echo times (TE). The MRI imaging device includes a controller. The controller is configured to generate at least two ΔTE images (ΔTE i=i . N ) based on the at least three MRI images, wherein the subscripts i=l ... N refer to images with sequentially increasing ΔTE times; wherein aliasing in the ΔTE ₁ image is permitted, and wherein the ΔTE times of ΔTE ₁ and ΔTE ₂ are set such that the alias points at which wrapping occurs in ΔTE ₁ does not overlap with the alias points of ΔTE ₂ . The controller is further configured to unwrap each ΔTE image, and set the final Bo map to the unwrapped ΔTEN image.

[0013] In accordance with related embodiments of the system, the controller may be further configured to determine shim currents and/or passive shims based on the final Bo map, the shim currents or passive shims for reducing Bo inhomogeneity. The controller may be further configured to apply the shim currents and/or passive shims to the sample to reduce Bo inhomogeneity in the sample. The shim currents and/or passive shims may generate a Bo field opposite in sign to the inhomogeneity of the magnetic field.

[0014] In accordance with another related embodiment of the invention, the at least two ΔTE images may differ in their echo time (TE) such that the difference in resulting phase ( ΔΦ ) of any two MRI ΔTE images varies according to ΔΦ= ΔTE Δv , wΔhere v isΔ the difference in frequency at a voxel in the image from a central frequency, and ΔTE is the difference in echo time between any pair of images.

[0015] In accordance with a further related embodiment of the invention, wherein to unwrap each ΔTE image to form a final Bo map, the controller may be configured to: unwrap aliasing in the ΔTE ₁ image is based on information in the next longest ΔTE image ( the ΔTE ₂ image) , to form an unwrapped ΔTE ₁ image, and use the unwrapped ΔTE ₁ image to unwrap the next longest ΔT E image (the ΔTE ₂ image). To unwrap aliasing in the ΔTE ₁ and ΔTE ₂ images, the controller may be configured to include calculating a logistic temporal unwrapping parameter LTUP, where LTUP images the fractional difference measured between two ΔTE images into a set of values describing the extent of aliasing. To unwrap each ΔTE image to form a final Bo map, the controller may be configured, for i=3 to N, use the previously unwrapped ΔTE _i-1 image to unwrap the ΔTE i image. To use the unwrapped ΔTE ₁ image to unwrap the next longest ΔTE image (the ΔTE ₂ image), the controller may be further configured to use a numerical estimate of the difference between the unwrapped value in the ΔTEi image with values representing different degrees of aliasing in the ΔTE ₂ .

[0016] In accordance with still further embodiments of the invention, the MRI device may acquire each of the MRI images using gradient echo, asymmetric spin echo, and/or stimulated echo acquisitions. To unwrap each ΔTE image, the controller may be configured to use a single pattern matching algorithm to simultaneously unwrap each ΔTE images. The MRI device may be configured to acquire the at least three MRI images using simultaneous multi-slice methods, K-space undersampling in the phase encoding direction, and/or use of 2 and 3 dimensional readout trajectories, so as to decrease acquisition times.

Brief Description of the Drawings

[0017] The foregoing features of embodiments will be more readily understood by reference to the following detailed description, taken with reference to the accompanying drawings, in which:

[0018] Fig. 1 shows a schematic of a Magnetic Resonance Imaging (MRI) system (prior art);

[0019] Fig. 2 shows a conventional gradient echo pulse sequence for Bo mapping;

[0020] Figs. 3A and 3B show Aliasing and Uncertainty in Bo Mapping;

[0021] Fig. 4 shows an example of logistical temporal unwrapping, in accordance with embodiment of the invention;

[0022] Fig. 5 shows the effects of uncertainty on LTUP in the ΔTE = 2ms image, in accordance with an embodiment of the invention;

[0023] Fig. 6 shows the effects of uncertainty in the ΔTE = 2ms image and the ΔTE = 9ms image, in accordance with an embodiment of the invention;

[0024] Fig. 7 compares images and Bo maps acquired from a healthy subject using conventional methodology; and using the Logistical Temporal Unwrapping methodology in accordance with embodiments of the invention;

[0025] Fig. 8 is an illustrative flow diagram of a Logistical Temporal Unwrapping methodology, in accordance with an embodiment of the invention; and

[0026] Fig. 9 shows three wrapped ΔTE maps and the final unwrapped Bo map obtained from a spherical phantom, in accordance with various embodiments of the invention.

Detailed Description of Specific Embodiments

[0027] Definitions. As used in this description and the accompanying claims, the following terms shall have the meanings indicated, unless the context otherwise requires:

[0028] The term “image” refers to a dataset in which a parameter(s) such as, without limitation, phase, Bo field variation etc... , is provided as a two or three dimensional map reflecting the spatial variation of that parameter. [0029] In illustrative embodiments of the invention, a system and method is provided for mapping and correcting the inhomogeneity of a magnetic field within an object using Magnetic Resonance Imaging (MRI) or Nuclear Magnetic Resonance (NMR) systems. The system and method includes logistical temporal unwrapping that enables the use of longer ΔTE images to unwrap aliasing in the shortest ΔTE image, and utilizes the patterns of aliasing in all images to unwrap the most accurate longer ΔTE image. Aliasing in the shortest ΔTE image is advantageously allowed, minimizing the need to accurately predict the maximum inhomogeneity present. Additionally, the ΔTE values used can lengthen by factors larger than two, enabling greater accuracy with fewer ΔTE images and thus shorter acquisition times. Details are described below.

[0030] Fig. 1 shows a conventional Magnetic Resonance Imaging (MRI) system. The MRI system includes a magnet, gradient coil (used for spatial localization), shim coil (used to correct Bo inhomogeneity) and RF coil (used for excitation and detection of signals). A controller, which may include a computer or other processing elements is used for control of the system. Shim and RF amplifiers may be used to drive their respective coils in the magnet. While the below described embodiments refer to MRI systems, it is to be understood that the invention is applicable to NMR or other imaging systems known in the art as well.

[0031] Performance Considerations for B ₀ Mapping

[0032] Fig. 2 shows a conventional gradient echo pulse sequence for Bo mapping. The sequence generally includes: a slice selective RF pulse (RF, SLICE) followed by a variable delay, Δ, an incremented phase encoding gradient (PE) and a readout gradient (RO). The acquisition is successively looped over N slices, M TE times, and finally P phase encode steps.

[0033] Bo mapping is commonly achieved using gradient echo images acquired with two or more different echo times. For samples or tissues with a single dominant resonance (e.g. water) the phase of the acquired images is given by the product of the difference in frequency of the dominant resonance from the carrier frequency (Δ v ) and the echo time (TE ) and a spatially varying constant (Φ ₀ (r)) dependent on a number of parameters including RF coil system, reconstruction method and other instrumental factors.

[0034]

[0035] Thus the phase difference, (ΔΦ), between any two images is given by [2], such that the only unknown is the frequency offset v. Δ

[0036]

[0037] Thus, assuming a single dominant resonance such that the distribution of frequencies in the image are solely related to the distribution of B ₀ field present, the spatial variation in ΔΦ provides a direct measure of the inhomogeneity (Δ B ₀ ) present in the sample or tissue.

[0038] For a given level of phase noise (∈ _pfiase ) and constant signal level, the measured phase (ΔΦ _msd ) and true phase (AΦ _true ) and measured (Δ v _msd ) and true frequencies (Δ v _true ) are given by Eqs. 3 and Eq. 4.

[0039]

[0040] Δ

[0041] Thus the accuracy of the measurement generally increases with increasing ΔTE.

[0042] However, since the phase ^<p _msd is periodic, with periodicity 2π , such that for long ΔTEs (i.e. more accurate measurements) if the phase will alias such that the true phase and frequency is given by Eq. 5 and Eq. 6 with n being an integer describing the amount of wrapping.

[0043] Δ [0044]

[0045] Thus Δv _true may not be uniquely determined in the presence of wrapping and noise/uncertainty. Thus the most accurate (longest ΔTE) is prone to large systematic errors due to aliasing. For example, Fig 3 A shows the measured values (Δ v _msd ) along the vertical axis as a function of the true value (Av _true ) along the horizontal axis for a ΔTE=2ms image over the range of +500Hz. The wrapping points are at +250Hz. Displayed in Fig 3B (lines added for visualization) is an expansion of the region from -270 to -240Hz (box in Fig 3a) displaying Δv _msd for the ΔTE=2ms image and a ΔTE=2ms with an “error” of +10Hz. For values between -270 and -260 Hz, both images alias and their difference is - 10Hz. However from -260 to -250Hz, the ΔTE=2ms image with a +10Hz error term does not alias, while the ΔTE=2ms images does resulting a difference of +490Hz.

[0046] Methods of unwrapping

[0047] To overcome this limitation two broad methods, spatial and temporal unwrapping have been described to unalias the data. See, for example: Jenkinson, M., Fast, automated, N-dimensional phase-unwrapping algorithm. Magn Reson Med, 2003. 49(1): p. 193-7; Hetherington, H.P., et al., Robust fully automated shimming of the human brain for high- field 1H spectroscopic imaging. Magnetic resonance in medicine, 2006. 56(1): p. 26-33;

Robinson, S., H. Schodl, and S. Trattnig, A method for unwrapping highly wrapped multiecho phase images at very high field: UMPIRE. Magn Reson Med, 2014. 72(1): p. 80-92; Geiger, Y. and A. Tai, Optimal echo times for multi-gradient echo-based B0 field-mapping. NMR Biomed, 2020. 33(7): p. e4316; and Dagher, J., T. Reese, and A. Bilgin, High- resolution, large dynamic range field map estimation. Magn Reson Med, 2014. 71(1): p. 105-17, each of which is incorporated herein by reference in its entirety.

[0048] Spatial unwrapping relies on the spatial continuity of the data to identify discontinuities in the Bo data as indicators of aliasing. See Jenksinson et al. These methods perform acceptably under conditions of moderate wrapping (i.e. moderate values of ΔTE) , such that the accuracy of the measurements is limited by the ΔTE used. For temporal unwrapping, multiple ΔTE values are used, with the shortest ΔTE chosen to eliminate aliasing (see Hetherington et al.), i.e. Under this condition, the values in the shortest ΔTE are used to “correct” aliased values in longer ΔTE images. Unfortunately, if the maximum uncertainty in the shortest ΔTE image ∈ ^max , where k is a scaling value such that all “outliers” are included, exceeds the bandwidth of the longer ΔTE image (1/ΔTE) the extent of aliasing will be ambiguous. Thus the “longer” ΔTE image cannot use an arbitrarily long ΔTE. For a normal distribution of noise/uncertainty, with being the standard deviation of the noise/uncertainty in the frequency domain, k would reflect number of standard deviations to account for the maximum noise value present in the sample

[0049] To overcome this limitation previously, a multi-ΔTE acquisition used a “boot strap” approach where: 1) the initial ΔTE value is chosen to be sufficiently short to eliminate aliasing; 2) each subsequent ΔTE image uses a factor of 2 in increasing duration and 3) each subsequent ΔTE image is used as the reference image to unwrap the next longest ΔTE image. Thus as the accuracy increases and bandwidth decreases in each subsequent image, the uncertainty in the “reference” ΔTE image decreases preserving the ability to unwrap the data unambiguously. However this method is limited by 1) the requirement to capture all inhomogeneity present (i. e. preclude aliasing) in the shortest ΔTE image and 2) the number of images required to get to arbitrarily long ΔTE values for a desired level of accuracy. These limitations manifest as longer acquisition times (more ΔTE images) for both increasing initial inhomogeneity and final accuracy.

[0050] Logistical Temporal Unwrapping

[0051] In illustrative embodiments of the invention, logistical temporal unwrapping provides the following unique advantages: 1) aliasing in the shortest ΔTE image is allowed, minimizing the need to accurately predict the maximum inhomogeneity present; and 2) the ΔTE values used can lengthen by factors larger than 2, enabling greater accuracy with fewer ΔTE images and thus shorter acquisition times. This methodology achieves these advantages by 1) using the longer ΔTE images to unwrap aliasing in the shortest ΔTE image and 2) utilizing the patterns of aliasing in all images to unwrap the most accurate image, i.e. the longest ΔTE image.

[0052] Unwrapping Aliasing in the Shortest TE image [0053] Aliasing results in discontinuities in the measured phase at well-defined frequency values in the bandwidth of the measurement. By comparing the measured values in the shortest ΔTE image with that measured in the longer ΔTE images, aliasing can be detected. Specifically a logistic temporal unwrapping parameter (LTUP) can be calculated, which enables the detection of wrapping in the shortest ΔTE image.

[0054]

[0055] Fig. 4 displays an example of the process graphically using ΔTE values of 2 and 9ms (based on three MRI images with TE times of 4, 6 and 13 ms) with up to a factor of 2 in aliasing (+500Hz, maximal expected bandwidth) present in the 2ms image, in accordance with an embodiment of the invention. When the value of the LTUP is (LTUP=0,+2, or +4, i.e., even numbers) the measured value in the shortest ΔTE image is unaliased, whereas if the value is (LTUP=+l,+3,+5, i.e., odd numbers) the value is aliased and can be corrected. For example, in the past conventional methodology, to achieve a true bandwidth of +500Hz, and achieve reasonable accuracy 4 ΔTEs would be acquired (i.e. 1, 2, 4 and 8ms). Whereas, under the current method, only 2 ΔTEs would be acquired yielding equivalent maximal bandwidth, but with superior accuracy (9 vs 8ms) with fewer ΔTEs (2 versus 4).

[0056] The effects of noise/uncertainty in the measurements is most pronounced at the aliasing points, where a few Hz of “noise” can dramatically change the measured values. For example, for a true value of 250Hz, if ∈ ^max =10Hz in the ΔTE=2ms image, the “true” value will map into a set of values between +240Hz and +260Hz. However, due to aliasing this set of true values {240, 241, ...260} maps into the values {240, 241, ...250, -249, -248,... - 240}. This is graphically represented in Fig. 5, which shows the effects of uncertainty in the ΔTE=2ms image on LTUP, in accordance with an embodiment of the invention, by two lines 2ms(-10Hz) and 2ms(+10Hz) representing the extremes of the numerical range. Assuming ∈ ^max ~0 in the ΔTE=9ms image, the true value of 250Hz maps into a value of 27,8Hz in the ΔTE=9ms image. However, due to the uncertainty in the ΔTE=2ms image, the LTUP can take on not only non-integral values, e.g. (240-27.7)/55.6 = +3.82, but also very different values when the frequency is aliased, e.g. (-240/-27.7)/55.6 ~ -4.81. [0057] The non-integral values of LTUP can be addressed by calculating LTUP*, where LTUP*=round(LTUP). After correction for non-integral values, however LTUP*= -5 or +4. Similarly for frequencies within + ∈ ^max of -250Hz, the LTUP* can take on values of - 4 or +5. However, the set of LTUP* values over these two regions are distinct, such that the LTUP* can still be used to determine if the value has aliased or not. For example if the true frequency is 245Hz, the +10Hz limits will place the measured value between 235 and 255Hz, yielding value of 235 to 250 and -250 to -245Hz in ΔTE=2ms image with two possible LTUP* values, +4 and -5. This range of values in the ΔTE=2ms image could also have been generated by a true frequency of -255Hz, with LTUP* values of -4 and +5. Outside of the potential aliasing regions, the calculated LTUP* values for the extremes of uncertainty in the ΔTE=2ms image are coincident (Fig. 5) and no ambiguity for aliasing exists. Thus the LTUP*, even in the presence of noise provides unambiguous unwrapping over these regions.

[0058] Unwrapping the ΔTE _long Image

[0059] Once the ΔTE=2ms image has been unwrapped and is outside of ±∈ ₂ (the uncertainty/noise), of the alias point of the image + l/(2ΔTE _long ). cΔan also be unwrapped. Fig. 6 assists in showing the unwrapping the ΔTE _long Image, in accordance with an embodiment of the invention (Fig. 6 includes experimental uncertainty in both the ΔTE ₂ ms and 9ms images). Specifically, for each in the ΔTE=9ms image there are nine theoretical potential values in the ΔTE=9ms image, which are given by: forming the pairs

Thus the unwrapped value ΔTE=9ms image, is achieved by selecting the nth pair of values , minimizing E _short (n) where

[0060] However, if the value of is within ±∈ _long (the uncertainty/noise) of the alias point of the image ±l/(2ΔTE _long ), “false” aliasing due to noise may have occurred. If “false” aliasing due to noise does occur, ∈ _short (n) will be larger than ∈ _short , the uncertainty in the ΔT ∈ _short image. In this case the value is used for at this location to unwrap the ΔTE _long image.

[0061] Alternatively, the process of unwrapping can be more generally viewed from the perspective of pattern matching, i. e. the finding the best numerical agreement between pairs of measured values (i.e. and the n sets of theoretical potential pairs i.e. of values, where the number of pairs is dependent upon the extent of aliasing and maximum bandwidth.

[0062] Thus, in general, in the presence of noise or uncertainty in the ΔTEiong image, unambiguous unwrapping can be achieved if uncertainties in the frequency region about the alias points in the image are spanned by continuous sets of differing values from the ΔTEshort image. Continuity of these values in the ΔTEshort image is achieved when the specific ΔTE values for both the ΔTEshort and ΔTEiong images are chosen such that the alias points+e ^max do not overlap. Once the ΔTEshort image is unwrapped, the ΔTEiong image can be unwrapped.

[0063] As described, we have allowed for noise and uncertainty in the ΔTEshort image but not the ΔTEiong image(s). The presence of noise or uncertainty in the ΔTEiong image pixels with true values near the alias points of the ΔTEiong image can result in multiple LTUP* values. To unalias the ΔTEshort image the regions of uncertainty about the alias points in the short ΔTE image are extended by the noise in the ΔTEiong image, i.e. ∈ ^max = e^ ^x _t + eiong- As long as these bands, alias points + ∈ ^max , are within the continuous region of the ΔTEiong image, the ΔTEshort image can be unaliased as described previously (see Fig. 5) using the LTUP* Since the ΔTE _short image is continuous over these regions, the ΔTEiong image can be unwrapped as described previously by creating the set of paired values and selecting the nth pair such that is a minimum or when near the wrap points of the ΔTEiong image and is large.

[0064] Creating Very High Accuracy Maps with Additional ΔT E images.

[0065] As described, we have utilized 2 ΔTE images ΔTE _(on9 ) to generate a Bo map. However, the ΔTE times chosen are limited by the maximum bandwidth of the sample, the aliasing patterns and the uncertainty in the individual ΔTE maps. Thus in general arbitrarily long ΔTE times for very high accuracy may not be possible under conditions of significant uncertainty and large sample bandwidth. To overcome this limitation additional ΔTE times (ΔTE ₃ , ΔTE ₄ , ... ) can be acquired with increasing length. Each additional ΔTE image then forms a new “pair” of images, with the long ΔTE image from the previous pair (e.g. ΔTE ₂ ) now serving as the short ΔTE image and the new image (e. g. ΔTE ₃ ) serving as ΔTE _long . Since the new ΔTE _sftort image has already been unwrapped the requirement for non-overlapping alias point regions, can be relaxed, and the new ΔTE chosen must only satisfy the criteria that the maximal uncertainty in the new pair of images (e g for be less than the bandwidth of the new ΔTE image, (e.g. 1/ΔTE ₃ ). Thus the process becomes recursive for each ΔTE image. Notably, as described previously, the multiplicative factor used in increasing the length of subsequent ΔTE times is not constrained by a factor of 2, but rather the actual noise/ uncertainty in the images. This allows for very efficient increases in accuracy per additional ΔTE image acquired. At some point, the noise/uncertainty in the longer ΔTE image may not decrease inversely with increasing duration.

[0066] Fig. 7 compares scout images and Bo maps acquired from a healthy subject using conventional methodology using ΔTE times of 1, 2, 4 and 8ms; and using the Logistical Temporal Unwrapping methodology in accordance with embodiments of the invention, with ΔTE times of 2 and 9ms. The maps display every 4th slice in the image set, from a set of 33 slices. The images are plotted on a scale of +100Hz. The target region of interest for shimming is shown in red on the scout images (top row) and black on the Bo maps. Despite fewer ΔTE times, the Logistical Temporal Unwrapping methodology shows a small but significant reduction (10%, 20.2 versus 22.5) in standard deviation of Bo over the entire brain in comparison to the conventional method. This is consistent with the lengthening of the maximum echo time from 8ms to 9ms (12.5%).

[0067] Fig. 8 is an illustrative flow diagram of a Logistical Temporal Unwrapping methodology, in accordance with an embodiment of the invention. At step 801, the varying echo times (TE) are selected based on BW _max and accuracy, which may be determined, for example, by ΔTE long. The Bo mapping image data is then acquired, step 803. Illustratively, at least two ΔT E images (ΔT E i=i .. N) may be generated based on at least three MRI images, wherein the subscripts i = 1 ... N refer to images with sequentially increasing ΔTE times. Advantageously, aliasing in the ΔTE ₁ image is permitted, and the ΔTE times of ΔTE ₁ and ΔTE ₂ are set such that the alias points at which wrapping occurs in ΔTE ₁ does not overlap with the alias points of ΔTE ₂ . The phase differences in the ΔTE images are calculated in the ΔTE maps are calculated and converted to frequency, step 805. Each ΔTE image is then unwrapped. More specifically, as an example, LTUP* is calculated and the ΔTE ₁ (short image) is unwrapped, step 807. Then, the ΔTE ₂ (long image) is unwrapped using the unwrapped ΔTE ₁ (short image), step 809. If there are more ΔTE images, step 811, set the unwrapped ΔTE (short image) to the unwrapped ΔTE (long image), and the ΔTE (long image) to the next longest ΔTE image, step 813. If there are no more ΔTE images at step 811, set the final Bo map to the last unwrapped ΔTEN (long image), step 815.

The shim currents and/or passive shims may then be calculated to generate the Bo correction field, step 817, which may then be applied to improve the Bo homogeneity.

[0068] Illustratively, the following is a more specific algorithm for unwrapping and setting of the final B0 map, in accordance with an embodiment of the invention.

1. Select a set of evolution times, a. A minimum of 3 evolution times to generate at least 2 ΔTE images b. Select the individual ΔTE times in the first two ΔTE image such that the frequency regions about the alias points (2/ΔTE), including the maximal uncertainty do not overlap over the range of maximum expected frequencies (BW _max ) in the sample or tissue. c. Select additional ΔTE times such that where the subscripts short and long refer to images with sequentially increasing ΔTE times. d. Aliasing in the image is permitted over BW _max .

2. Acquire the data.

3. Calculate the phase in the ΔTE images and convert to frequency.

4. Unwrap the ΔTE _short image: a. Divide the bandwidth of the ΔTE _sfwrt image into two conceptual regions: i. Frequencies within of the alias points of the image ii. Frequencies outside of this region b. For frequencies within region (ii) i. Calculate the sets of theoretical values for LTUP* over as a function of extent of aliasing (n=0,±l,... .) ii. Calculate LTUP* from the measured data iii. Use the measured LTUP* to look up the value of n from the theoretical values in 4.b.i iv. c. For frequencies within region (i) i. Perform the same calculation as in 4.b.i over the the range of frequencies spanned by the alias point ii. Calculate LTUP* from the measured data iii. Use the measured LTUP* to look up the value of n from the set of theoretical values in 4.c.i. Each aliased region will have multiple but distinct LTUP* values for a given value of n. iv.

5. Unwrap the long ΔTE image: a. Using the measured value in the image calculate the potential pairs of frequency values for the b. c. d.

Unwrap additional longer ΔTE images if present by . . . a. Substituting the parameters for the image with those from the b. Substituting the parameters for the image with those from the c. Repeating the calculation described in 5

7. After all images have been unwrapped, use the values from the longest ΔT E image to generate the final map.

8. Using the final Bo map described in 7 calculate a set of shim currents and/or determine passive shimming for generating a correction Bo field of approximately equal amplitude but opposite sign.

9. Apply the calculated shim currents and/or passive shims to reduce the Bo inhomogeneity in the sample.

[0069] Fig. 9 shows three wrapped ΔTE maps and the final unwrapped Bo map obtained from a spherical phantom, in accordance with various embodiments of the invention. More particularly, the figure shows nine slices from a set of 33 slices acquired from the spherical phantom in which the shims were purposely offset to generate a significant amount of inhomogeneity. The bottom three rows show the wrapped maps from the: ΔTE =2 ms image (bottom row); 9 ms image; (2 ^nd row from bottom); and 21 ms image, plotted on their scales of +/-BW/2 where BW=1/ΔTE. The top row shows the final Bo map obtained from the unwrapped last ΔTE = 21 ms map, using the above-described methodology and embodiments of the invention.

[0070] Embodiments of the invention may be located implemented in part in any conventional computer programming language. For example, preferred embodiments may be implemented in a procedural programming language (e.g., “C”) or an object oriented programming language (e.g., “C++”, Python). Alternative embodiments of the invention may be implemented as pre-programmed hardware elements, other related components, or as a combination of hardware and software components.

[0071] Embodiments also can be implemented in part as a computer program product for use with a computer system - for example, the controller of the MRI system described above. Such implementation may include a series of computer instructions fixed either on a tangible medium, such as a computer readable medium (e.g., a diskette, CD-ROM, ROM, or fixed disk) or transmittable to a computer system, via a modem or other interface device, such as a communications adapter connected to a network over a medium. The medium may be either a tangible medium (e.g., optical or analog communications lines) or a medium implemented with wireless techniques (e g., microwave, infrared or other transmission techniques). The series of computer instructions embodies all or part of the functionality previously described herein with respect to the system. Those skilled in the art should appreciate that such computer instructions can be written in a number of programming languages for use with many computer architectures or operating systems. Furthermore, such instructions may be stored in any memory device, such as semiconductor, magnetic, optical or other memory devices, and may be transmitted using any communications technology, such as optical, infrared, microwave, or other transmission technologies. It is expected that such a computer program product may be distributed as a removable medium with accompanying printed or electronic documentation (e.g., shrink wrapped software), preloaded with a computer system (e.g., on system ROM or fixed disk), or distributed from a server or electronic bulletin board over the network (e.g., the Internet or World Wide Web). Of course, some embodiments of the invention may be implemented as a combination of both software (e.g., a computer program product) and hardware. Still other embodiments of the invention are implemented as entirely hardware, or entirely software (e.g., a computer program product).

[0072] Although various exemplary embodiments of the invention have been disclosed, it should be apparent to those skilled in the art that various changes and modifications can be made which will achieve some of the advantages of the invention without departing from the true scope of the invention.