B0 field inhomogeneity estimation using internal phase maps from long single-echo MRI acquisitions

Abstract

To use once acquisition for estimating B0 field inhomogeneity with high resolution, shorten a time for the acquisition, and avoid a factor of an error.SOLUTION: A magnetic resonance (MR) image may be created from MR data by receiving the MR data, applying a transform to the MR data, where a result of the application is an image space representation of the MR data, determining a wrapped phase map of the image space representation of the MR data, obtaining an unwrapped phase map based on the wrapped phase map, scaling the unwrapped phase map into a B0 field map, reconstructing the MR image based on the MR data, correcting the MR image based on the B0 field map, and outputting the MR image. The scaling may be free of accounting for effects on the MR data by artifact sources secondary to B0 field inhomogeneities.SELECTED DRAWING: Figure 1

Description

[Technical Field]

[0001] [Priority Claim] This application claims priority to U.S. Provisional Application No. 63 / 124,911, filed December 14, 2020, the contents of which are incorporated herein by reference.

[0002] A method for estimating B0 field inhomogeneity in magnetic resonance imaging (MRI) is disclosed. [Background technology]

[0003] MRI (Magnetic Resonance Imaging) is used to acquire images of patients, but the quality of the acquired images can be impaired by various artifacts. For example, inhomogeneity in the magnetic field (e.g., the B0 field) causes distortion and blurring artifacts. As the echo time (TE) of the image acquisition increases, the artifacts become more severe.

[0004] Image quality can be improved by correcting artifacts. For example, understanding the impact of B0 field inhomogeneity allows for the reduction of artifacts caused by inhomogeneity. In some cases, B0 field maps are acquired in addition to patient images and used to correct for artifacts. These approaches are prone to error and require overlapping full scans, increasing already long examination times. Reducing the resolution of the scan, and therefore the B0 field map, shortens the long examination time, resulting in insufficient error correction capabilities. Summary of the Invention

[0005] By way of preamble, the preferred embodiments described below include a method and system for estimating B0 field inhomogeneity that uses a single acquisition to estimate B0 field inhomogeneity with high resolution, reducing acquisition time and avoiding sources of error.

[0006] A first aspect is a method for generating an MR image from MR (magnetic resonance) data, the method including the steps of: a processor receiving the MR data; the processor applying a transform to the MR data, the result of the application being an image space representation of the MR data; the processor determining a wrapped phase map of the image space representation of the MR data; the processor obtaining an unwrapped phase map based on the wrapped phase map; the processor scaling the unwrapped phase map to a B0 field map, the scaling being based solely on the effect of B0 field inhomogeneity on the phase of the MR data; the processor reconstructing an MR image based on the MR data; the processor correcting the MR image based on the B0 field map; and the processor outputting the MR image.

[0007] The method can include a processor applying a virtual coil combination to an image space representation of the MR data. The wrapped phase map can be determined based on the virtual coil combination applied to the image space representation of the MR data. The method can include a processor determining an intensity image based on the image space representation of the MR data. The method can include a processor applying a sum-of-squares operation to the image space representation of the MR data. The determination of the intensity image can be based on the sum-of-squares operation applied to the image space representation of the MR data. The acquisition of the unwrapped phase map can be based on the intensity image. The method can include a processor generating a mask image based on the intensity image and applying the mask image to the B field map. The contribution of heating, motion, radio frequency pulse inhomogeneity, coil sensitivity, and eddy currents, or a combination thereof, to the phase of the MR data is not accounted for in the scaling. The MR data is received from a clinical acquisition, and the correction of the MR image is based on the B field map of only that clinical acquisition. The scaling of the unwrapped phase map can be based on an approximation of the B field map where ΔB = φ / 2πTE. The MR data can be acquired from a non-Cartesian acquisition scheme. The MR data can be acquired with an echo time of at least 20 ms, a magnetic field strength of about 3 T, or a combination thereof. The method can further include the processor applying a low-pass filter to the B0 field map to remove signals above a frequency range of B0 field inhomogeneity. The transform can be an inhomogeneity transform.

[0008] A second aspect is a medical imaging system including an MR (magnetic resonance) imaging device, a processor, and a memory, the memory storing executable instructions that function to receive MR data from the MR imaging device, apply a transform to the MR data, the result of which is an image space representation of the MR data, determine a wrapped phase map of the image space representation of the MR data, obtain an unwrapped phase map based on the wrapped phase map, scale the unwrapped phase map to a B field map, the scaling ignoring contributions to the MR data by artifact sources secondary to the B field map, reconstruct an MR image based on the MR data, correct the MR image based on the B field map, and output the MR image.

[0009] The memory of the system can store instructions operable to apply a virtual coil combination scheme to an image space representation of the MR data. A wrapped phase map is determined based on the virtual coil combination scheme applied to the image space representation of the MR data. The memory of the system can store instructions operable to determine a magnitude image based on the image space representation of the MR data. The memory of the system can store instructions operable to apply a sum-of-squares operation to the image space representation of the MR data. The magnitude image is determined based on the sum-of-squares operation applied to the image space representation of the MR data. An unwrapped phase map is obtained based on the magnitude image. The memory of the system can store instructions operable to generate a mask image based on the magnitude image and apply the mask image to the B0 field map.

[0010] A third aspect is a non-transitory computer-readable medium having stored thereon processor-executable process steps, wherein the processor performs the process steps in a system that receives clinical MR data from an MR imaging device, applies a transform to the clinical MR data, where the result of the application is an image space representation of the clinical MR data, determines a wrapped phase map of the image space representation of the MR data, obtains an unwrapped phase map based on the wrapped phase map, scales the unwrapped phase map to a B0 field map according to a linear relationship between B0 and echo time, reconstructs an MR image based on the clinical MR data, corrects the MR image based on the B0 field map, and outputs the MR image.

[0011] The present invention is defined by the claims, and nothing in this section should be construed as limiting the scope of those claims. Additional aspects and advantages of the present invention are described below in conjunction with the preferred embodiments and may be claimed individually or in any combination. [Brief explanation of the drawings]

[0012] The components and drawings are not necessarily to scale, emphasis instead being placed upon illustrating the principles of the invention. Moreover, in the drawings, like reference characters refer to like parts throughout the various views. [Figure 1] 1 illustrates an embodiment of a method for generating a MR (Magnetic Resonance) image. [Figure 2] 1 illustrates an embodiment of a flowchart for image reconstruction with correction based on a B0 field map. [Figure 3] 10 illustrates an example of a phase map and a B0 field map. [Figure 4] 1 illustrates various phase maps. [Figure 5] 1 illustrates various susceptibility weighted images. [Figure 6] 1 illustrates an embodiment of a medical imaging system. DETAILED DESCRIPTION OF THE INVENTION

[0013] Patient-induced inhomogeneities in the B0 field can cause artifacts, such as distortion and blurring, in image acquisitions, which are more pronounced for acquisitions with long echo times, such as in susceptibility-weighted imaging (SWI) and / or T2* gradient echo (GRE) imaging.

[0014] SWI is used for high-resolution cerebral venography (e.g., diagnosing traumatic brain injury) due to its higher sensitivity to venous blood. Long echo times (e.g., 20–40 ms) are used in SWI to enhance magnetic susceptibility information. However, SWI data acquired with long echo times may contain amplified off-resonance artifacts, such as geometric distortion and image blur. Correcting the off-resonance artifacts could broaden the adoption of SWI.

[0015] Furthermore, artifacts accumulate across the readout direction by affecting the imaging gradient over time. In non-Cartesian sampling patterns (used in SWI and / or T2*GRE), this accumulation increases the severity of artifacts in the data. Non-Cartesian sampling patterns allow for better spatial resolution and faster acquisition, for example, when used in compressed sensing reconstruction. Correcting the artifacts allows for wider adoption of non-Cartesian sampling patterns.

[0016] In some cases, artifacts due to inhomogeneities are corrected based on information about the B field. Some methods of correcting artifacts involve establishing a B field map (ΔB) using an additional low-resolution MR acquisition of the patient in addition to the primary clinical acquisition, which requires a total of two volumetric acquisitions using separate echo times.

[0017] To calculate ΔB0 from the phase maps (φ1 and φ2) at each voxel of two volumes acquired at different echo times (e.g., TE1 and TE2, where TE1 > TE2), use the following equation: φ1=φ0+2π×ΔB0×TE1 φ2=φ0+2π×ΔB0×TE2 Transformed ΔB0=(φ1-φ2) / (2π×(TE1-TE2))

[0018] The φ term adds effects on the phase of other sources beyond the B field, such as heating, motion, radiofrequency pulse inhomogeneity, coil sensitivity, and eddy currents. The B field map shows the B field inhomogeneity and allows for correction of artifacts caused by the inhomogeneity. However, the additional acquisition increases the total examination time for the patient. Furthermore, this method is prone to errors in the correction if the patient moves between the B field map acquisition and the clinical acquisition.

[0019] By approximating the B0 field map based on only one echo acquisition (e.g., only the acquisition from the patient), additional acquisitions are avoided and / or examination time is reduced. The B0 field map is approximated based on only phase images from one acquisition (e.g., a SWI acquisition). When two acquisitions are used to determine the B0 field map, the second acquisition has lower resolution. Estimating the B0 field map using a single acquisition benefits from high spatial resolution, robustness to scan-to-scan motion, and minimal impact on image reconstruction time.

[0020] In some cases, to avoid additional acquisitions, a ΔB field map is estimated from clinical acquisitions. For example, a gradient-based method or a grid search method can be used. However, these methods have problems such as poor B field estimation, poor correction performance, and long estimation times.

[0021] Another method for estimating a B0 field map from a clinical acquisition is based on the principle that the principal components of the observed phase (from a single acquisition) are linearly related to ΔB0 and TE. In some cases, the relative contribution of non-ΔB0 terms (e.g., secondary error sources to B0 field inhomogeneity, such as heating, motion, RF pulse inhomogeneity, coil sensitivity, and eddy currents) is not significant. For long echo times, the contribution of secondary sources to the phase map is considered negligible compared to the contribution caused by B0 inhomogeneity. For example, conditions such as a 3 T field strength and an echo time of 20 ms or longer acquisition are suitable for SWI acquisitions. An estimated B0 field map generated with 0.6 mm isotropy can be estimated within 3 minutes from a single acquisition and provides comparable performance in terms of artifact correction compared to a B0 field map determined from a second acquisition. Therefore, acquisitions used to reconstruct and generate patient images for diagnostic or clinical purposes are also used to estimate the B0 field map.

[0022] As mentioned above, the φ term adds the effects of artifacts caused by secondary sources such as heating, motion, radio frequency pulse inhomogeneity, coil sensitivity, and eddy currents. For example, eddy currents are weak in a properly designed and / or calibrated MR scanner. These secondary sources of artifacts are ignored for certain imaging parameters (e.g., long TE at intermediate field strengths), so the B field map is approximated as follows: ΔB0=φ / 2πTE

[0023] Additional sources of image artifacts, such as large (lumped) motion, may not be negligible and are corrected by other techniques beyond B field map correction. Image artifacts caused by RF field emission and reception are each individually limited and correctable. Overall, their contribution to the φ term is approximately 15-20 times lower than 2πTEΔB (the effect of B field inhomogeneity), e.g., 3 T with a TE of at least 20 ms. Therefore, the φ term can be removed and the approximation of the B field map can be performed accordingly.

[0024] Approximating the B0 field map without considering the φ term offers significant benefits to medical professionals, patients, and hospitals. For example, the approximated B0 field map can be obtained without additional acquisitions with performance comparable to that of a B0 field map determined from two acquisitions. In one example, a B0 field map is determined from two acquisitions lasting 2 minutes and 43 seconds, while a B0 field map is approximated from a single clinical acquisition lasting 2 minutes and 55 seconds. Furthermore, the approximated B0 field map can produce image corrections comparable to those obtained with a B0 field map determined from two acquisitions (an example is shown in Figure 5). In addition, the resolution of the approximated B0 field map is the same as the corresponding acquisition, allowing for higher-fidelity corrections compared to a two-acquisition B0 field map. Determining a high-resolution B0 field map using a two-acquisition approach requires even longer acquisition times. In practice, the tradeoff between resolution and acquisition time for determining a two-acquisition B0 field map leads to a B0 field map with very low resolution, which in practice results in insufficient artifact correction during image reconstruction. As disclosed herein, approximating the B0 field map using a single acquisition improves each of these aspects.

[0025] In one application, a medical image scanner and reconstruction workflow uses a single acquisition approach to provide retrospective B field approximation to correct past image acquisitions where B field correction was not performed during clinical acquisition. This application provides significant benefits for sequences that are sensitive to ΔB artifacts, where ΔB is not / was not acquired during clinical imaging, even though the acquisition is at a long echo time.

[0026] FIG. 1 illustrates a method for generating a magnetic resonance image. More, fewer, or different steps may be performed. In some cases, one or more of steps 109, 111, 113, and / or 121 may be omitted. The steps may be performed in a different order than shown. For example, step 107 and / or step 111 may proceed directly from step 103. In another example, step 119 may proceed directly from step 115. In another example, step 123 may proceed directly from step 119. In another example, step 121 may be performed as part of step 103 or may proceed directly from step 103. A processor coupled to a memory may be configured to perform one or more of the steps. For example, processor 603 and memory 605 of FIG. 6 may be configured to perform one or more of the steps of FIG. 1. Additional, different, or alternative aspects of the processes are described below in conjunction with FIG. 2.

[0027] In step 101, MR data is received. The MR data can be received from a medical imaging device, such as an MR imaging device. The MR data is generated in a clinical MR acquisition, rather than an additional acquisition performed to determine a B field map. The clinical acquisition is an acquisition with parameters suitable for imaging patient tissue. For example, the clinical acquisition has a longer echo time than a subsequent acquisition performed to observe the B field. In this manner, the B field map is approximated based solely on the clinically acquired MR data, rather than on an additional acquisition. The MR data is k-space data, such as k-space data 201 of FIG. 2.

[0028] In some examples, the MR data is collected according to SWI or other long echo time (e.g., 20 ms or longer) processes. The MR data is collected in a Cartesian or non-Cartesian coordinate system. In one example, the MR data is acquired with an echo time of at least 20 ms and / or a magnetic field strength of about 3 T. In another example, the echo time is 30 ms. In yet another example, the magnetic field strength is less than 3 T, e.g., 1.5 T. In another example, the magnetic field strength is greater than 3 T, e.g., 7 T.

[0029] In step 103, an adjoint transform is applied to the MR data (e.g., adjoint transform 203). The transform receives k-space MR data as input. The result of applying this transform is an image space representation of the MR data. In one example, the transform is a nonuniform fast Fourier transform (NUFFT). Other transforms used in reconstructing object or image space from k-space can be used.

[0030] In step 105, a virtual coil combination method is applied to the image space representation (transformation) of the MR data (e.g., combination 205). The virtual coil combination method can be applied when an array of coils is used to receive MR signals for the MR data. In some examples, the wrapped phase map determined in step 107 (e.g., wrapped phase map 207) is determined based on the virtual coil combination method. In other examples, another combination method is used to determine the wrapped phase map from the image space representation of the MR data. The virtual coil combination method describes the behavior of the MR coils used to acquire the MR data. In this way, the virtual coil combination method allows information about the phase of the magnetic resonance signal to be obtained from the image space representation. The wrapped phase map determined from the image space representation (e.g., by applying the virtual coil combination method or another combination method) collects or represents the phase information. The term "wrapped" in "wrapped phase" is used because the phase of the signal is constrained to a range of, for example, -π to +π. In the wrapped phase map, if the measured phase value is outside the range, the value is corrected to fall within the range. The result of wrapping the phase information within a range is shown by the sharp local contrast bands in the phase maps of Figures 2, 3 and 4.

[0031] When MR data is acquired using a single coil, the wrapped phase map is determined directly from the MR data, i.e., the wrapped phase map is obtained without applying any combination method, such as virtual coil combination.

[0032] In step 109, a sum-of-squares technique (e.g., technique 209) is applied to the image space representation of the MR data. In some examples, the magnitude image (e.g., magnitude image 211) determined in step 111 is based on the sum-of-squares technique or other technique applied to the image space representation of the MR data. The magnitude image represents or indicates the strength of the MR signals received from the voxels of the MR data acquisition. In this way, the magnitude image indicates another aspect of the image space representation of the MR data in addition to the wrapped phase map.

[0033] In step 113, a mask image (e.g., mask image 213) is generated based on the intensity image. The mask image removes portions of the imaging volume that are inappropriate or undesirable for clinical evaluation and / or modification of the reconstructed image. For example, if brain tissue or vasculature are of particular interest, the mask image removes portions of the MR data that represent objects unrelated to the target, such as the skull.

[0034] In step 115, an unwrapped phase map is generated (e.g., by unwrapping 215). The unwrapped phase map is obtained based at least on the wrapped phase map, and in some cases based on the wrapped phase map and an intensity image. For example, the wrapped phase map is unwrapped by solving a Poisson equation weighted by the intensity of the previously estimated image. By solving the Poisson equation, phase values ​​that are not constrained to a range of, for example, -π to +π can be determined. The unwrapping process determines phase values ​​that are not constrained to a specified range, thereby more accurately and in greater detail indicating the local and global variance of the phase. The unconstrained phase values ​​allow for accurate correction based on the resulting B0 field map.

[0035] In step 117, a mask image is applied to the phase-unwrapped phase map (e.g., during or as part of unwrapping 215). Applying the mask image to the unwrapped phase map removes unwanted or irrelevant portions of the unwrapped phase map. In this manner, the magnitude image (the basis of the mask image) can be used to focus on the phase information used to create the B field map.

[0036] In step 119, the unwrapped phase map is scaled to obtain a B0 field map (e.g., B0 field map 221). The scaling proceeds, for example, as described with respect to scaling 217. As described above, artifacts in the MR data are caused by B0 field inhomogeneities and secondary sources. However, by focusing only on the contribution of B0 field inhomogeneities to MR data artifacts, the B0 field map is determined based on a linear relationship between the phase values ​​for voxels of the MR data and the echo time of the clinical acquisition. Based on the echo time, each voxel is "scaled" according to the phase value for the voxel, giving a local value for the B0 field at the voxel. The value of the B0 field at the voxel can be determined according to B0 = φ / 2πTE. The B0 field map is determined by combining the local B0 values ​​for each voxel (or at least the voxels of interest).

[0037] In step 121, a filter is applied to the B0 field map (e.g., proceeding as described with respect to filter 219). The filter may be, for example, a low-pass filter. The filter removes contributions to the MR data due to susceptibility in the B0 field map while preserving contributions to the MR data due to B0 field inhomogeneity. For example, the susceptibility contributions are distributed in one region of the spectral domain, and the filter removes the susceptibility contributions while simultaneously removing signals outside that region. By removing the susceptibility contributions in the B0 field map, only contributions due to B0 inhomogeneity are removed from the reconstructed image, and the susceptibility is present in the corrected image.

[0038] In step 123, an MR image is reconstructed from the MR data (e.g., according to the reconstruction in section C of FIG. 2 and / or the reconstruction in reconstruction 223). The image reconstruction process proceeds from the same MR data, the k-space MR data, used to estimate the B field map. Transformations, such as non-Fourier operations, are used to bring the k-space data into image space. The reconstructed image is affected by artifacts from various sources, such as B field inhomogeneities. The artifacts are corrected by applying the B field map to the reconstructed image, as in step 125. The output of the correction is, for example, a correction volume 225. The local B field value at a voxel can be used to weight the magnetic intensity of the reconstructed image at the same voxel. Through this weighting, the observed values ​​at the voxel are scaled. (Increase / Decrease) This eliminates the effects of a non-uniform B0 field throughout the imaging volume.

[0039] In step 127, the MR images are output. The MR images are output to one or more destinations. For example, the MR images are sent to a storage device (such as memory 605 or other repository), an MR system, and / or a display (such as display 611). The reconstructed MR images may be reformatted, such as by reconstructing a volume and volume rendering an image from the reconstructed volume.

[0040] FIG. 2 shows a flowchart for reconstructing an image with correction based on a B field map. Each process in FIG. 1 also describes additional, different, or alternative aspects of the process in FIG. 2. In section A, a phase map is obtained from k-space data 201 from a clinical MR acquisition. In section B, the phase map is scaled, and in some cases, filtered, to approximate a B field map that represents B field inhomogeneities at the time of acquisition. In section C, an image (e.g., a three-dimensional image) is reconstructed from the MR data and corrected based on the approximated B field map.

[0041] K-space data 201 (e.g., MR data from a clinical acquisition) is received, e.g., as described with respect to step 101. Beginning in section A, an intensity image 211 and / or a wrapped phase map 207 is obtained by computing an adjoint 203 of the density-compensated k-space data 201, e.g., as described with respect to one or more of steps 103, 105, 107, 109, and / or 111. While an adjoint NUFFT 203 is shown, other transforms (e.g., discrete Fourier transform, fast Fourier transform, partial (fast) Fourier transform, etc.) may also be used. Furthermore, while a virtual coil combination 205 and a sum-of-squares operation 209 are shown, other combinations (e.g., adaptive coil combination or susceptibility map) and other operations on the coil channels in the MR data may also be used, respectively.

[0042] In section B, an estimated (approximated) B field map 221 is generated by unwrapping 215, scaling 217, and / or low-pass filtering 219. Through unwrapping 215, an unwrapped phase map is determined from wrapped phase map 207, for example, as described with respect to process 115. The wrapped phase map is unwrapped 215 by solving a Poisson equation weighted by the intensity of a previously estimated image 211. The magnitude image 211 and associated mask 213 estimated in adjoint 203 are used to weight and stabilize the phase unwrapping algorithm 215. The unwrapping process 215 of section B receives the wrapped phase map 207, the magnitude image 211, and in some instances the mask image 213 as inputs, although other unwrapping techniques 215 may be used. In one example, an unwrapping technique 215 is used that unwraps the wrapped phase map without using the magnitude image 211 and / or the mask image 213. For example, Laplace's method or the Speedy rEgion-Growing Algorithm for Unwrapping Estimated Phase method (SEGUE) can be applied to unwrap the wrapped phase map 207.

[0043] The unwrapped phase map is linearly scaled 217 based on the relationship between the voxel-specific B value, the voxel-specific phase value, and the echo time of the k-space data 201 acquisition. For example, the voxel-specific B field value can be determined according to B = φ / 2πTE. A low-pass filter 219, such as a Hanning filter, can be applied to the scaled B field map to filter out high frequencies. The low-pass filter 219 prevents suppression of susceptibility contributions to the MR data in the reconstructed and corrected image 225, which is subject to SWI. The resulting approximation of the B field map 221 is used to correct artifacts in images reconstructed from the k-space data 201.

[0044] In section C, an image is reconstructed 223 based on the k-space MR data 201. Reconstruction 223 proceeds as described with respect to step 123. In some examples, non-Fourier operations are used to reconstruct an MR image from the k-space MR data 201. However, the reconstructed image is susceptible to artifacts. The artifacts are corrected in image space by applying an estimated B field map 221, such as a B field map estimated according to one or more steps in FIG. 1. After correction, a corrected imaging volume 225 is acquired.

[0045] For SWI, an additional step 231 is performed on the corrected volume 225 to obtain an SW image. For example, based on the corrected volume 225, an additional phase image can be determined, e.g., as described with respect to steps 103, 105, and 107 of FIG. 1 and elements 203, 205, and 207 of FIG. 2. A filter is applied to this additional phase image to remove low frequencies (e.g., as described with respect to step 121 and element 219) and retain high-frequency contributions to the phase image due to magnetic susceptibility. A mask is created from the filtered additional phase image based on multiplication of voxel values ​​(e.g., 4-5 times) and applied to the intensity image of the corrected volume 225. Through this process, an SWI image is obtained.

[0046] A minimum intensity projection 227 is applied to a filtered further phase image determined from the corrected image volume 225 via processing step 231 to generate a corrected SW image 229 .

[0047] In one example, the processes of Sections A and B are performed to estimate the B field map 221 without considering contributions to image artifacts from sources beyond B field inhomogeneity. In this example, the B field map is approximated in less than 3 minutes, with adjoint calculations completed in 127 seconds, coil combination using a virtual coil combination method completed in 31 seconds, mask generation from the intensity image completed in 5 seconds, and phase unwrapping, filtering, and scaling completed in a total of 12 seconds.

[0048] Figures 3, 4, and 5 show examples and comparisons of phase maps, acquired B field maps, estimated B field maps, and images corrected according to the B field maps. In the examples in Figures 3, 4, and 5, MR data were acquired as SWI image volumes for two volunteers at 3T (using Siemens Healthineers PrismaFIT, Erlangen, Germany) using a 64-channel head / neck coil array and a 3D spherical stack-of-SPARKLING sampling pattern. The MR data were acquired with the following parameters: acquisition time (TA) = 5 min, 0.6 mm isotropic resolution, field of view (FOV) = 240 mm, number of slices = 208, TE = 20 ms, relaxation time (TR) = 50 ms, and Bw = 33 Hz / px. An additional reference B field map was collected with the following parameters: TA = 2 min 43 s, FOV = 240 mm, 2 mm isotropic resolution, TE = 4.92 ms, and TE = 7.38 ms. The reference B field map served as a baseline against which the performance of the single-collection B field maps could be compared.

[0049] Image reconstruction (e.g., for the image shown in Figure 4) was performed iteratively in 3D using soft-thresholding regularization in the wavelet domain. Corrected non-Fourier operations were used with pre-computed density compensation. Channels were recombined using sum-of-squares (e.g., as described with respect to steps 109 and / or 111) and virtual coils (e.g., as described with respect to steps 105 and / or 107). Reconstructions were performed uncorrected and corrected using either acquired or estimated field maps. SWI minimum projections were applied at 16 mm or greater. All post-processing was performed using a 2560-core Quadro P5000 GPU with 16 GB of GDDR5 VRAM.

[0050] In FIG. 3, images associated with two patients are shown through rows 1 and 2. Column A shows a phase map, column B shows a B phase map based on a two-acquisition technique, and columns C, D, and E show B field maps estimated from a single clinical acquisition. The B maps in columns C, D, and E are estimated according to different unwrapping algorithms, such as the multiple unwrapping techniques described with respect to step 115 and element 215. The estimated B map in column C, using an unwrapping technique based on at least the unwrapped phase map 207 and magnitude image 211, demonstrates improved error correction performance over the B maps in columns D and E.

[0051] In Figure 4, the phase maps in rows 1, 2, 3, and 4 correspond to the images in rows 1, 2, 3, and 4 of Figure 5. Column A in Figure 4 shows SW images without correction and corresponds to the images in column A in Figure 5. Column B in Figure 4 shows SW images corrected with the reference B field map, and column C shows SW images corrected with the estimated B field map.

[0052] In Figure 5, rows 1, 2, 3, and 4 show different SW images reconstructed from MR data. Column A shows images without artifact correction, and column B shows images corrected with a reference B field map (e.g., a B field map established from two acquisitions). Column C shows the absolute difference between the images in columns A and B, allowing the extent of correction relative to the baseline B field map to be observed. Column D shows images corrected using a B field map estimated from clinical acquisitions alone, and column E shows the absolute difference between the images in columns A and D, allowing the extent of correction relative to the estimated B field map to be observed. Arrows in Figure 5 highlight subtle SW image features present in the images in columns A, B, and D.

[0053] As can be seen from FIG. 3 (e.g., columns B, C, D, and E), unwrapping (e.g., performed as described with respect to step 115) provides greater consistency in the interior regions of the B field map compared to the collected B field map. However, the phase map shown in column B of FIG. 4 (e.g., when compared to the phase map shown in column C of FIG. 4) suggests that the B values ​​in the collected B field map are too high because the phase is negatively wrapped (starting from the center and decreasing to −π until the phase value wraps to +π). As a result, the effects of inhomogeneity are overcompensated. The estimated B field map used in the correction image shown in column C of FIG. 4 yields better compensation. The resulting SW image in FIG. 5 shows equivalent correction between the reference B field map (shown in column C of FIG. 5) and the proposed B field map approximation method (shown in column E of FIG. 5).

[0054] In the case of intermediate magnetic fields (e.g., 3 T) and long echo times (e.g., 20 ms or longer), a single phase image (determined from a single clinical acquisition) can be used to quickly and efficiently estimate a B0 field map, allowing for correction of off-resonance artifacts with performance comparable to that of corrections performed on a B0 field map obtained using additional ΔB0 acquisitions. A single-acquisition B0 field estimate may have a duration comparable to that of a two-acquisition technique. However, a single-acquisition B0 field estimate is advantageous because the B0 field is estimated offline (e.g., on a computer and without additional image acquisition) with greater robustness to inter-scan motion and higher resolution. In one example, applying a low-pass filter (e.g., as described with respect to step 121) during the adjoint call (e.g., as described with respect to step 103) can further reduce the time required to estimate the B0 field map.

[0055] 6 shows an example embodiment of a medical imaging system 601. The system 601 includes a processor 603, a memory 605, a network adapter 607, an image sensor (medical imaging device) 609, and a display 611. The image sensor 609, the display 611, the processor 603, and the memory 605 may be part of a medical imaging device, a computer, a server, a workstation, or other system for image processing of patient image data.

[0056] Additionally, different or fewer components may be provided. For example, the image sensor 609 may be remote from the system 601. Image reconstruction and / or B0 field map estimation may be applied as a standalone application in the system 601 or a local device, or as a service deployed in a network (e.g., cloud) architecture. As another example, user input devices (e.g., keyboards, buttons, sliders, dials, trackballs, mice, or other devices) may be provided for user manipulation of the patient image data.

[0057] The processor 603 may be a controller, control processor, general-purpose processor, digital signal processor, three-dimensional data processor, graphics processing unit, application-specific integrated circuit, field programmable gate array, artificial intelligence processor, digital circuit, analog circuit, combination thereof, or other device now known or later developed for processing image data. The processor 603 may be one device, multiple devices, or a network of devices. In the case of multiple devices, parallel or serial division of processing may be used. The processor 603 operates according to or is configured by stored instructions, hardware, and / or firmware to perform the various processes described herein. For example, the processes of FIGS. 1 and 2 may be stored as instructions and configured to configure the processor 603 to perform the processes of FIGS. 1 and 2.

[0058] The memory 605 can be an external storage device, RAM, ROM, a database, and / or local memory (e.g., a solid-state drive or hard drive). The same or another non-transitory computer-readable medium can be used for instructions and other data. The memory 605 can be implemented using a database management system (DBMS) and can reside in memory such as a hard disk, RAM, or removable media. Additionally or alternatively, the memory 605 can be internal to the processor 603 (e.g., cache). The memory 605 can store MR data, phase maps, intensity images, B0 field maps, reconstructed images, and / or computer program instructions. Data stored in the memory 605 can be accessed and retrieved by the processor 603 or another processor.

[0059] The instructions, methods, and / or techniques described herein for performing the MR image generation, MR image correction, and / or B0 field estimation processes are provided in a non-transitory computer-readable storage medium or memory, such as a cache, buffer, RAM, removable media, a hard drive, or other computer-readable storage medium (e.g., memory 605). Computer-readable storage media include various types of volatile and non-volatile storage media. The functions, processes, or techniques illustrated in the figures or described herein are performed in response to one or more sets of instructions stored in a computer-readable storage medium. These functions, processes, or techniques are independent of the particular type of instruction set, storage medium, processor, or processing architecture, and are instead performed by software, hardware, integrated circuits, firmware, microcode, and the like, operating alone or in combination.

[0060] In one embodiment, the instructions are stored on a removable media device for reading by a local or remote system. In another embodiment, the instructions are stored in a remote location, for example, for transfer over a computer network. In yet another embodiment, the instructions are stored on a particular computer, CPU, GPU, or system. Some components and method steps making up the systems illustrated in the figures are implemented in software, so the actual connections between system components (or process steps) will vary depending on the programming scheme of the embodiment.

[0061] The network adapter 607 communicates with one or more wired or wireless networks. Data can be sent and received between the system 601 and other computers on the network via the adapter 607. For example, MR data can be collected from a remote computer, remote storage, or remote MR imaging system via the adapter 607. In another example, MR images can be sent to a remote computer via the adapter 607 (e.g., for storage or further processing).

[0062] The medical imaging device 609 can be an MR imaging device. For example, the medical imaging device records data based on a magnetic field applied to an imaging volume. Although one medical imaging device 609 is shown, more than one medical imaging device 609 may be provided. The medical imaging device 609 receives a patient, such that the MR data contains information about the patient's magnetic resonance within the imaging volume. The MR data collected by the medical imaging device 609 is saved. For example, the image data is stored in memory 605 or in storage remote from the system 601.

[0063] The display 611 may be a CRT, LCD, projector, plasma, printer, tablet, smartphone, or other display device now known or that may be developed in the future for displaying output such as patient history data, etc. In some examples, the display 611 may present visual or audiovisual output.

[0064] While the present invention has been described through various embodiments, it will be understood that many changes and modifications can be made without departing from the scope of the present invention. The foregoing detailed description is intended to be interpreted as illustrative rather than limiting, and it will be understood that it is the following claims, including all equivalents, that are intended to define the spirit and scope of the invention.

Claims

1. 1. A method for generating an MR (magnetic resonance) image from MR data, comprising: a processor receiving MR data obtained from only one echo acquisition; the processor applying a transform to the MR data, the result of which is an image space representation of the MR data; the processor applying a virtual coil combination to the image space representation of the MR data; the processor determining a wrapped phase map of the image space representation of the MR data, the wrapped phase map being determined based on the virtual coil combination applied to the image space representation of the MR data; the processor obtaining an unwrapped phase map based on the wrapped phase map; the processor scaling the unwrapped phase map to a B field map, the scaling being based solely on the effect of B field inhomogeneity on phases in the MR data; the processor reconstructing an MR image based on the MR data; the processor correcting the MR image based on the B0 field map; and wherein the processor outputs the MR image.

2. The method of claim 1 , further comprising the processor determining an intensity image based on the image space representation of the MR data.

3. the processor further applying a sum-of-squares operation to the image space representation of the MR data; The method of claim 2 , wherein the determination of the intensity image is based on the sum-of-squares operation applied to the image space representation of the MR data.

4. The method of claim 2 , wherein the obtaining of the unwrapped phase map is based on the intensity image.

5. generating a mask image based on the intensity image by the processor; The method of claim 4 further comprising the processor applying the mask image to the unwrapped phase map.

6. The method described in claim 1, wherein the contribution of heating, motion, radio frequency pulse non-uniformity, coil sensitivity, and eddy currents, or a combination thereof, to the phase in the MR data is not represented in the scaling.

7. the MR data is received from a clinical acquisition; The method of claim 1 , wherein the correction of the MR image is based on the B 0 field map of the clinical acquisition only.

8. 2. The method of claim 1, wherein the scaling of the unwrapped phase map is based on an approximation of the B field map with ΔB=φ / 2πTE.

9. The method of claim 1 , wherein the MR data is obtained from a non-orthogonal acquisition.

10. The method of claim 1 , wherein the MR data is acquired with an echo time of at least 20 ms, a magnetic field strength of about 3 T, or a combination thereof.

11. The method of claim 1 , further comprising the processor applying a low pass filter to the B0 field map to remove signals above a frequency range of the B0 field inhomogeneity.

12. The method of claim 1 , wherein the transformation applied to the MR data is a non-uniform transformation.

13. 1. A medical imaging system including an MR (magnetic resonance) imaging device, a processor, and a memory for storing instructions, The instructions, when executed by the processor, receiving MR data obtained by only one echo acquisition from the MR imaging device; applying a transform to the MR data, the result of which is an image space representation of the MR data; applying a virtual coil combination to the image space representation of the MR data; determining a wrapped phase map of the image space representation of the MR data, the wrapped phase map being determined based on the virtual coil combination applied to the image space representation of the MR data; obtaining an unwrapped phase map based on the wrapped phase map; scaling the unwrapped phase map to a B field map, the scaling being based solely on the effect of B field inhomogeneities on phases in the MR data; reconstructing an MR image based on the MR data; correcting the MR image based on the B0 field map; and outputting the MR image.

14. The instructions stored in the memory, when executed by the processor, The medical imaging system of claim 13 , further operable to determine an intensity image based on the image space representation of the MR data.

15. The instructions stored in the memory, when executed by the processor, 15. The medical imaging system of claim 14, further operative to apply a sum-of-squares operation to the image space representation of the MR data, and wherein the intensity image is determined based on the sum-of-squares operation applied to the image space representation of the MR data.

16. The medical imaging system of claim 14 , wherein the unwrapped phase map is obtained based on the intensity image.

17. The instructions stored in the memory, when executed by the processor, generating a mask image based on the intensity image; The medical imaging system of claim 16 , further operative to apply the mask image to the unwrapped phase map.

18. A non-transitory computer-readable medium having processor-executable process steps stored thereon, comprising: A system for executing the process steps by a processor, receiving clinical MR data from an MR imaging device, the clinical MR data being obtained from only one echo acquisition; applying a transform to the clinical MR data, the result of the application being an image space representation of the clinical MR data; applying a virtual coil combination method to the image space representation of the clinical MR data; determining a wrapped phase map of the image space representation of the clinical MR data, the wrapped phase map being determined based on the virtual coil combination applied to the image space representation of the clinical MR data; obtaining an unwrapped phase map based on the wrapped phase map; scaling the unwrapped phase map to a B field map, the scaling being based solely on the effect of B field inhomogeneity on phases in the clinical MR data; reconstructing an MR image based on the clinical MR data; correcting the MR image based on the B0 field map; and a non-transitory computer-readable medium operative to output the MR image.

Images