[go: up one dir, main page]

EP2288941A1 - Étiquetage de phase utilisant le codage de sensibilité : acquisition de données et reconstruction d'image pour correction de distorsion géométrique en imagerie écho-planaire - Google Patents

Étiquetage de phase utilisant le codage de sensibilité : acquisition de données et reconstruction d'image pour correction de distorsion géométrique en imagerie écho-planaire

Info

Publication number
EP2288941A1
EP2288941A1 EP09739969A EP09739969A EP2288941A1 EP 2288941 A1 EP2288941 A1 EP 2288941A1 EP 09739969 A EP09739969 A EP 09739969A EP 09739969 A EP09739969 A EP 09739969A EP 2288941 A1 EP2288941 A1 EP 2288941A1
Authority
EP
European Patent Office
Prior art keywords
images
space
image
map
phase
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
EP09739969A
Other languages
German (de)
English (en)
Other versions
EP2288941A4 (fr
Inventor
Udomchai Techavipoo
Thomas P. Leist
Song Lai
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Thomas Jefferson University
Original Assignee
Thomas Jefferson University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Thomas Jefferson University filed Critical Thomas Jefferson University
Publication of EP2288941A1 publication Critical patent/EP2288941A1/fr
Publication of EP2288941A4 publication Critical patent/EP2288941A4/fr
Withdrawn legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56563—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5611—Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5616—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using gradient refocusing, e.g. EPI
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56341—Diffusion imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56509—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/4806—Functional imaging of brain activation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56554—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by acquiring plural, differently encoded echo signals after one RF excitation, e.g. correction for readout gradients of alternating polarity in EPI

Definitions

  • the technical field generally relates to magnetic resonance imaging. More specifically, the invention relates to systems and methods for correcting geometric distortion in echo planar imaging (EPI) with phase labeling using sensitivity encoding.
  • EPI echo planar imaging
  • Magnetic resonance imaging uses a magnetic field and radio frequency (RF) energy pulses as a non-invasive method for analyzing objects.
  • MRI is used extensively in medical imaging.
  • an object or patient is placed in an external magnetic field.
  • the nuclear magnetic moments of the nuclei in the patient are excited at specific spin precession frequencies that are proportional to the external magnetic field.
  • Radio frequency (RF) signals resulting from the precession of the spins are collected using receiver coils.
  • the magnetic fields are altered using magnetic gradients, and the RF signals are collected that represent different anatomical regions of the patient.
  • the collected signals are combined to produce a volumetric image of the nuclear spin density of the patient.
  • RF pulses radio frequency (RF) pulses
  • the RF pulses have a unique frequency chosen to excite a particular nuclei, such as hydrogen, in the patient.
  • the RF pulses excite the hydrogen nuclei, thereby increasing the energy state of the nuclei.
  • the hydrogen nuclei relax and release RF energy as a free induction decay signal.
  • the free induction decay signal can be transformed into an echo.
  • the echoes are observed, measured, and converted into anatomical images.
  • the RF pulses may have a wide frequency spectrum to excite nuclei over a large range of resonant frequencies. Similarly, the RF pulses may have a narrow frequency spectrum to excite nuclei over a more narrow range of resonant frequencies. Combined, or composite RF pulses appear as a series of RF pulses and may be used to excite nuclei over different ranges of resonant frequencies. Composite RF pulses may be transmitted to excite multiple ranges of resonant frequencies allowing for simultaneous collection of received signals from multiple areas of interest, such as multiple areas of the patient anatomy.
  • a narrow RF pulse When a narrow RF pulse is applied, only those nuclei that precess at that particular frequency will be excited, and subsequently relax and release RF energy as a free induction decay signal.
  • the nuclei "resonate" to that particular frequency. For example, if a magnetic gradient caused hydrogen nuclei to precess at rates from 42.56 MHz at the low end of the gradient to 42.6 MHz at the high end, and the gradient were set up such that the high end was located at the patient's forehead and the low end at the patient's chin, then a 42,56 MHz RF pulse would excite the hydrogen nuclei in a slice near the chin, and a 42.6 MHz pulse would excite the hydrogen nuclei in a slice near the forehead.
  • the second dimension of the planar patient image is derived using a phase-encoding gradient.
  • a phase-encoding gradient is briefly applied in the y- direction in order to cause the magnetic vectors of nuclei along different portions of the gradient to have a different phase advance.
  • the sequence of RF and gradient pulses is then repeated to collect all the data necessary to produce an image.
  • the magnitude of the phase-encoding gradient is stepped as the number of repetitions progresses. For example, the phase-encoding gradient may be evenly incremented after each repetition.
  • the number of repetitions of the pulse sequence is dependent upon the type of image desired and can be any integer, such as 512. 1024, and the like.
  • the polarity of the phase encoding gradient may also be reversed to collect additional RF signal data. For example, when the number of repetitions is 512, 256 of the repetitions may be positive, and the other 256 repetitions may have a negative polarity phase encoding gradient of the same magnitude.
  • the x-axis gradient is the "frequency encoding gradient" or "read gradient.”
  • This gradient causes the relaxing nuclei to precess at different rates, so that the nuclei near one end of the gradient begin to precess at one higher rate, and those at the other end of the gradient precess at an even higher rate. As these nuclei relax, the nuclei at the high end of the gradient give off the highest frequency RF signals.
  • the read gradient is applied when the RF signals are to be measured.
  • the second and third dimensions of the image are acquired using a fast Fourier transformation (FFT).
  • the FFT decomposes the received RF signal into a sum of sine waves. Each of the sine waves have a different frequency, phases, and amplitudes. For example:
  • either the x-axis or y-axis may be set up as the slice-select axis depending upon the desired image orientation and the anatomical structure of the patient to be examined.
  • the x-axis may be utilized as the slice-select axis to acquire sagittal images
  • the y-axis may be utilized as the slice-select axis to acquire coronal images.
  • the slice select gradient, phase-encoding gradient, and read gradient are orthogonal.
  • the result of the MRI scan in the frequency domain representation (k-space) is converted to image display in the time domain data after a 2D or 3D fast Fourier transformation (FFT).
  • FFT fast Fourier transformation
  • the read gradient is along the x-axis and the phase-encoding gradient is along the y-axis.
  • an additional phase-encoding gradient is along z-axis to acquire data in a third dimension.
  • the k-space is an array of numbers/data whose Fourier transform is the MRI image.
  • the k-space data is arranged in an inhomogeneous distribution such that the data at the center of the k-space map contains low spatial frequency data. That is, the low spatial frequency data is the general spatial shape of the patient being scanned.
  • the data at the edges of the k-space map contains high spatial frequency data including the spatial edges and details of the patient.
  • the main magnetic field, the gradient magnetic fields, and the frequency composition of the RF excitation pulses may all cause quicker de-phasing if any of these elements are nonuniform.
  • Echo-planar imaging is a fast MRI data acquisition technique that has been widely used in many imaging applications for both structural and functional studies. Echo- planar imaging (EPl) techniques are capable of acquiring an entire MRI image in only a fraction of a second.
  • EPI is sensitive to magnetic field inhomogeneity caused by imperfect magnetic field shimming and tissue susceptibility differences. Tissue susceptibility differences result from nuclei resonating at different frequencies depending upon the material in which the nuclei are located. Magnetic field inhomogeneity creates local magnetic field gradients superimposed on the applied spatial encoding gradients. Therefore, spatial locations are decoded incorrectly, resulting in geometric distortion.
  • a number of distortion correction algorithms have been attempted previously.
  • conventional distortion correction techniques involve two steps, namely acquiring field inhomogeneity information and applying the field inhomogeneity information to correct the distortion.
  • To obtain the field inhomogeneity information single or multiple scans are performed with varying scan parameters, such as the echo time, the direction or polarity of the phase and frequency encoding gradients, the location of the k-space data, and the number of echoes or the number of encoding directions.
  • the field inhomogeneity information may be represented in terms of field maps, phase shifts, amplitude and phase errors, point spread functions, or relationships between corrected and distorted coordinates.
  • field inhomogeneity information is generally derived from a given patient position in the magnet and is valid only for that particular patient position. Therefore, conventional scans for field inhomogeneity information have to be acquired immediately before or after the EPI measurements to reduce patient position inconsistencies. For lengthy or repeated EPI measurements, such as DTI, fMRI, or dynamic contrast agent studies, patient motion during or between the scans would invalidate the patient position consistency requirement for applying the field inhomogeneity information for geometric distortion correction. Further, it is also impractical to sacrifice the temporal resolution by inserting reference scans for obtaining field inhomogeneity information in between the dynamic points. For these reasons, EPI images in many applications are often used without geometric distortion correction.
  • Previous attempts to address geometric distortion correction have involved obtaining field maps from the EPI measurements for every dynamic time point.
  • field maps are directly estimated from k-space echo displacements of each individual region (or pixel) in the spatial domain. This method can be applied to gradient-echo and asymmetric spin-echo images, but not to symmetric spin-echo images. Further, it is computationally intensive because the echo displacements are estimated pixelwise.
  • the EPI trajectories are modified such that multiple low-resolution scans with different echo times are embedded within each dynamic point. For example, the center of the k-space is collected twice. Field maps obtained using this approach have a relatively lower spatial resolution.
  • the present invention provides a system, method, and computer program product for providing geometric distortion correction in echo planar imaging (EPI) that generates corrected images without the shortcomings of previous techniques.
  • a system, method, and computer program product in accordance with the present invention employs p_hase labeling using sensitivity encoding (PLUS), that utilizes the EPI measurement data themselves to correct geometric distortion, without acquiring separate scans for field maps and coil sensitivity maps.
  • the system and method of the present invention integrates the technique of rjhase labeling for additional coordinate encoding (PLACE) for mapping field inhomogeneity with the methods of simulated phase evolution rewinding (SPHERE) that applies the PLACE-derived field maps to correct the distortion.
  • PLACE sensitivity encoding
  • SPHERE simulated phase evolution rewinding
  • the PLACE technique requires at least two images to generate a field map. Instead of acquiring phase images with different echo times as in other techniques, in PLACE, the phase images are acquired with different pre- phase -encoding gradients. Without changing other parameters of the pulse sequence, PLACE varies the area of the pre-phase-encoding gradient by adding/subtracting an area equal to a multiple (N) of the area of a phase-encoding pulse (blip). This manipulation shifts the k- space data down/up N lines from the original k-space data. Field maps from PLACE possess the same distorted spatial domain as the distorted images and are easy to pass to the SPHERE calculations.
  • SPHERE simulates the k-space data by re-phase-encoding each spatial location of a distorted image with additional reverse local phase error generated from a given field map. Finally, the simulated k-space data is inverse Fourier transformed to generate a corrected image with reduced artifacts and improved signal-to-noise ratios (SNRs).
  • SNRs signal-to-noise ratios
  • the present invention employs phase labeling using sensitivity encoding (PLUS), that utilizes the EPI measurement data themselves to correct geometric distortion, without acquiring separate scans for field maps and coil sensitivity maps.
  • the PLUS technique of the present invention incorporates a parallel imaging technique using phased array coils and their sensitivity maps to produce multiple images from a single measurement scan.
  • the trajectories of k-space data are modified such that the phase differences of the images after parallel imaging reconstruction contain local phase shifts. These phase shifts are applied afterwards to correct the distortion in the average of the images reconstructed previously using parallel imaging.
  • ghosting artifacts are also removed by alternately assigning k-space lines with the same readout directions into groups before applying parallel imaging reconstruction.
  • the measurement data themselves are also used to derive sensitivity maps for parallel imaging.
  • FIGURE IA illustrates a modified EPI trajectory for PLUS.
  • FIGURE 1 B illustrates the pulse sequence diagram for the EPI trajectory of
  • FIGURE IA is a diagrammatic representation of FIGURE IA.
  • FIGURES 1C- IF show k-space data rearranged into four groups according to EPI trajectories.
  • FIGURE 2 shows a block diagram illustration of a phase labeling using sensitivity encoding system in accordance with the present invention.
  • FIGURES 3A-3B show a comparison of artifacts in gradient-echo images acquired with a modified EPI trajectory and with a conventional EPI trajectory.
  • FIGURES 4A-4F illustrate the generation of a coil sensitivity map in accordance with the present invention.
  • FIGURES 5A-5D shows averaged unfolded images from a PLUS EPI trajectory and a conventional EPI trajectory.
  • FIGURES 6A-6F show a comparison of phase-shift maps estimated using PLUS and phase-shift maps estimated using PLACE, and the difference between the two phase-shift maps.
  • FIGURES 7A-7D show a comparison of average unfolded images and the geometric distortion correction using the phase-shift maps using PLUS and geometric distortion correction using the phase-shift maps using PLACE.
  • FIGURES 8A-8D show the geometric distortion correction results of conventional
  • FIGURES 9A-9D show 128-shot conventional gradient-echo images and Tl- weighted MP-RAGE images.
  • FIGURES 1 OA-IOB illustrate two EPI trajectories for PLUS.
  • FIGURE 11 illustrate an MRI system in accordance with the present invention. DETAILED DESCRIPTION OF THE INVENTION
  • the present invention includes a system, method, and computer program product for providing geometric distortion correction in echo planar imaging (EPI) that generates corrected images without the shortcomings of previous techniques.
  • a system, method, and computer program product in accordance with the present invention employs phase labeling using sensitivity encoding (PLUS), that utilizes the EPI measurement data themselves to correct geometric distortion, without acquiring separate scans for field maps and coil sensitivity maps.
  • the system and method of the present invention integrates the technique of phase labeling for additional coordinate encoding (PLACE) for mapping field inhomogeneity with the methods of simulated phase evolution rewinding (SPHERE) that applies the PLACE-derived field maps to correct the distortion.
  • PLACE sensitivity encoding
  • SPHERE simulated phase evolution rewinding
  • the system, method, and computer program product of the present invention corrects geometric distortion caused by magnetic field inhomogeneity in echo planar imaging (EPI) by using local phase shifts derived directly from the EPI measurement itself, without the need for extra field map scans or coil sensitivity maps.
  • EPI echo planar imaging
  • phase labeling for additional coordinate encoding requires extra scans that have different pre-phase-encoding areas compared to the EPI measurement scan.
  • phase labeling using sensitivity encoding (PLUS) techniques in accordance with the present invention skips lines in each of the k-space data sets and acquires the remaining lines in a time multiplexed fashion, resulting in a new k-space trajectory.
  • the k-space trajectory for PLUS shown in FIGURE IA may be modified in a number of ways. One example is, where k-space lines may be rearranged into two and three groups, respectively. Additional example trajectories are shown in FIGURE 10 and in the accompanying discussion to FIGURE 10.
  • a PLUS trajectory (shown in FIGURE IA as thin dotted lines, for example reference numeral 1 10) and its pulse sequence diagram 190 are presented in FIGURE IA and FIGURE IB, respectively.
  • the acquisition order of the second and the third lines of the k-space data are switched.
  • the acquisition order of the second and third lines are switched. This may be illustrated by starting at PLUS trajectory acquisition starting point 112 and moving along the PLUS trajectory 1 10 through acquisition points 114, 1 16, 1 18, 120, 122.
  • Some of the pulses (blips) shown in FIGURE IB are doubled as represented by double triangles 115, 117, 119, 121.
  • the k-space data shown in FIGURE IA are then rearranged into four group Gl, G2, G3, G4 by alternately assigning a k-space line to a group. That is, the thick lines (solid Ll, L5, dash-dash L2, L6, dash-dot L3, Ll, and dotted lines L4, L8) represent lines of the k-space data in FIGURE IA and their assignments to groups Gl, G2, G3, G4.
  • the four EPI trajectories shown as thin dotted lines 140, 150, 160, 170 presented in FIGURES 1C-1F are equivalent to the trajectories used in PLACE that have -I 5 0, 0, and +1 additional pulses (blips), respectively, in the pre-phase gradients.
  • the missing k-space lines in each group Gl, G2, G3, G4 will be filled up using parallel imaging reconstruction as described below. This reconstruction is equivalent to unfolding the images in the spatial domain.
  • the phase differences between the unfolded images from FIGURES 1C and ID and from FIGURES IE and IF are calculated and then averaged to obtain a phase-shift map.
  • SPHERE simulated phase evolution rewinding
  • a block diagram of the steps in carrying out the PLUS method in accordance with the present invention is shown in FIGURE 2.
  • FIGURE 1 1 illustrates a magnetic resonance imaging (MRI) system in accordance with the present invention configured to carry out the PLUS method.
  • the MRI system includes magnetic coils 1105 that provide a constant main magnetic field.
  • the main magnetic field strength may be 0.5 Tesla, 1 Tesla, 1.5 Tesla, 2 Tesla, and the like.
  • Magnetic gradients are generated by gradient amplifiers, including the x-gradient 1 1 10, y-gradient 1 115, and z-gradient 1 120.
  • the x-gradient 1 1 10, y-gradient 1115, and z-gradient 1120 produce waveforms that modify he main magnetic field by generating gradient magnetic fields through gradient coils 1125.
  • a radio frequency (RF) pulse generator 1130 creates RF pulses and transmits the pulses using RF coil 1 135.
  • the transmitted pulses excite the nuclei of the patient or object being examined.
  • Phased array coil 1135 may also receive the relaxation RF signals from the nuclei in the object as the excited nuclei precess.
  • the nuclei may be hydrogen nuclei, or other atoms.
  • a spectrometer 1140 processes the received RF emission signals from the nuclei and further processes the received signals using the PLUS techniques of the present invention.
  • the images may be further processed using computer system 1 145, and the images or imaging information may be displayed on a display device 1150.
  • the sensitivity encoding (SENSE) process is used for parallel imaging reconstruction in accordance with the present invention.
  • Other parallel imaging algorithms may be applied as well.
  • the SENSE reconstruction performs essentially two steps, namely calculating coil sensitivity maps and applying the maps to unfold images.
  • the SENSE reconstruction process is modified so that phases of unfolded images contain additional phases of a coil sensitivity map. These additional phases will be canceled out when calculating a phase-shift map.
  • the modified SENSE reconstruction processes described are independent from the MRl scanner internal image reconstruction software and the coil sensitivity maps are calculated directly from the measurement data.
  • a coil sensitivity map is calculated, and then the map is applied to unfold the images.
  • the complex-valued reference image si c is given by:
  • p and c* are the complex-valued images of the object and the k th coil sensitivity, respectively; a and ⁇ are their phases, respectively; n* is the zero-mean, complex- valued white noise from the A* coil.
  • SOS square root of the sum of the squares
  • Equation [1] the number of coils in the phased array coil.
  • can be estimated by masking the object region and fitting the image with a two-dimensional polynomial.
  • the order of the polynomial may be the 5 lh order.
  • the object phases a have to be removed from Equation [3] as well.
  • the phase differences between the images from the N" 1 coil and from the other coils may be utilized.
  • the object phases a may be removed and images containing phase differences ⁇ k - ⁇ N may be created as follows:
  • Equation [0044] where the double primes denote the modifications of the measurement and the noise.
  • the phases of the image in Equation [4] are smooth.
  • the phase differences ⁇ k - ⁇ N can be estimated by low-pass filtering the image in Equation [4] and calculating the image phases.
  • An alternative approach is to find the phase differences from Equation [1], instead of Equation [3]. This approach can avoid the noise amplification from dividing the measurements with the sum of the squares (SOS) image. 2) Image unfolding.
  • the coil sensitivity map may be used to unfold the images.
  • the equation for a folded image is as follows.
  • Equation [5] Equation [5]
  • the number of folds has to be less than the number of coils so that v can be estimated.
  • the method of the present invention estimates v using least squares.
  • the magnitude of the unfolded image ⁇ u ⁇ can be found by multiplying the magnitude of the result
  • phase-shift map A ⁇ in accordance with the present invention may be calculated from two images that have different additional blip areas in the pre-phase-encoding gradients. That is:
  • a ⁇ J— ATg[RPR[I 1 )RPR ⁇ I j )' ) , Equation [7] [0049]
  • Arg(.) represents the phase angles of a complex-valued image
  • RPRQ is an operation to remove a phase ramp in the phase-encoding direction (or the y-direction) of the image that is created by adding blips to the k-space data
  • the indices i and/ denote the numbers of blips added to the k-space data.
  • the relationship between this phase-shift map and the field map ABo is given by:
  • phase ramp is equal to 2 ⁇ nylN y , where n is the number of additional blips and N y is the number of pixels in the y-direction.
  • the phase ramp may be removed using circular shifts in the k y direction.
  • a phase-shift map in accordance with the present invention can be calculated by:
  • each image pair in the image products of this equation is derived from the k-space data in which the lines are scanned in the same direction. For example, the k-space lines of the images u ⁇ and « 2 are scanned from left to right, while those of the images uj and «4 are scanned from right to left. Therefore, the phases accumulated along the readout direction are canceled out.
  • Simulated phase evolution rewinding is used to apply a phase-shift map to correct a distorted image.
  • the distorted image has to be upsampled before applying SPHERE in order to avoid k-space aliasing artifacts caused by scaling the spatial domain.
  • upsamp ⁇ ng to four times is sufficient to avoid k-spac ⁇ aliasing artifacts.
  • higher upsarnpling may also be employed.
  • a x-k y space data h after rewinding the phase error using a phase-shift map A ⁇ is given by:
  • n is an index for the k y direction
  • p and q are indices of the x and y directions, respectively
  • N y is the total number of the original pixels in the y direction.
  • the components of h(p, ⁇ ) where q ⁇ ⁇ (R- ⁇ )N y l2 and q > (R- ⁇ )N y !2 are removed afterwards to prevent k-space aliasing.
  • inverse Fourier transformation is applied in the k y direction and downsampling is applied to reduce the number of pixels to match the original matrix size.
  • a method of acquiring data and images in accordance with the present invention may be performed using a 3 Tesla MRl scanner with an eight-channel sensitivity encoding (SENSE) head coil to receive imaging signals from the object under study.
  • the head coil may employ a phased array configuration of coils.
  • a pulse sequence was developed and used to acquire k-space data for phase labeling using sensitivity encoding (PLUS) as illustrated in FIGURE IB.
  • PLUS sensitivity encoding
  • the phase encoding gradient pulse (blip) area is doubled at some of the phase encoding points 115, 1 17, 1 19, 121 as indicated by double triangles.
  • the k-space data is then rearranged into four groups Gl, G2, G3, G4.
  • the four groups Gl, G2, G3, G4 represent data in common EPI trajectories 140, 150, 160, 170 with different numbers of additional blips added to the pre-phase encoding gradients.
  • the numbers of additional blips are -1, 0, 0, and 1, respectively.
  • k-space lines are traversing in the same readout direction.
  • the modified EPI trajectory was based on a general gradient echo EPI sequence for fMRI studies.
  • the pulse sequence 190 used to acquire k-space data for phase labeling using sensitivity encoding may include the following scan parameters.
  • TR may be set to 2500 ms, and TE set to 35 ms.
  • the FOV was 384 mm with 3x3 mm2 in-plane resolution for an imaging phantom and was 410 mm with 3.2x3.2 mm2 in-plane resolution for a human brain.
  • Slice thickness was 4 mm.
  • the magnitude and phase images were recorded from the individual coils of the eight-channel sensitivity encoding (SENSE) head coil.
  • SENSE eight-channel sensitivity encoding
  • Additional imaging sequences were also performed for comparison purposes, including common gradient-echo EPI data (without k-space trajectory modification), gradient-echo data using PLACE (5 dynamics with ⁇ 2 additional blips in pre-phase-encoding gradients), 128-shot conventional gradient-echo data (without using EPI readout), and high resolution Tl- weighted images using an MP-RAGE sequence. Post acquisition processing was performed on the acquired signal data.
  • FIGURE 3 is a comparison of artifacts in the gradient-echo images acquired with the modified EPI trajectory for the PLUS image 301 in FIGURE 3A) and with a conventional EPI trajectory 351 in FIGURE 3B.
  • These images 301, 351 are the square roots of the sums of the squares of the images from the multi-channel phased array coils used to receive the imaging signals from the object (imaging phantom) under study. The image from each coil is reconstructed using a 2D Fourier transformation.
  • the PLUS trajectory creates additional ghosting artifacts 31 1, 321 at ⁇ N/4 pixels (where N is the number of pixels in phase-encoding direction and the origin is at the center of the image 301).
  • FIGURES 4A-4F demonstrate the generation of a coil sensitivity map.
  • An individual image 401 from one coil of the phased array is shown in FIGURE 4A.
  • the upper and lower right parts 41 1, 421 of the object in FIGURE 4A are affected by the interference of the N/4 ghosts.
  • the individual image 401 from the phased array coil is divided by the square roots of the sums of the squares image 301 in FIGURE 3 A, and the resulting raw sensitivity magnitude image 441 is shown in FIGURE 4B. After the division operation is performed, the magnitudes of the quotient show some remaining interference patterns 412, 422 as seen in FIGURE 4B. However, these patterns 412, 422 disappear after polynomial fitting as shown in the image 451 in FIGURE 4C.
  • phase differences between the combined images from two coils of the phased array evolve slower as shown in image 471 in FIGURE 4E, and they become very smooth after low-pass filtering as shown in image 481 in FIGURE 4F.
  • the discontinuity of the shades of colors at the upper left parts 482 of the phantom is caused by phase wraparound. This discontinuity does not affect the parallel imaging reconstruction performed in accordance with the present invention.
  • FIGURE 5 A Averaged unfolded images and regular EPI images of a phantom and a brain are shown in FIGURE 5.
  • the phantom image 501 in FIGURE 5 A is generated from the same Ik- space data used in FIGURE 3A but includes the PLUS method described above. The ghosting artifacts and their interferences exhibited in FIGURE 3 completely disappear.
  • This averaged unfolded image 501 in FIGURE 5 A is comparable with the conventional EPI image shown in FIGURE 5B. The same comparison is shown for the brain images 521, 571 in FIGURES 5C and 5D.
  • FIGURE 6 Phase-shift maps for the phantom and brain images 501, 521 , 551, 571 in FIGURE 5 are shown in FIGURE 6,
  • the maps 601, 621 estimated using PLUS are shown in FIGURES 6A and 6D and the maps 651 , 671 estimated using PLACE are shown in FIGURES 6B and 6E.
  • the differences 680, 690 between the maps are also shown in FIGURES 6C and 6F.
  • the field inhomogeneity in the phantom is more pronounced than the field inhomogeneity in the brain for this particular slice.
  • the range of the colorbar shades is ⁇ /4 for the phantom images 601 , 651 and ⁇ /8 for the brain images 621 , 671.
  • phase-shift maps estimated using PLUS are slightly different from those estimated using PLACE. The differences between the two estimations are small.
  • the mean and the standard deviation are -2.94xlO "3 ⁇ 5.69x10 "3 ⁇ radian (-0.19 ⁇ 0.36 pixels) for the phantom and are 7.19 ⁇ lO "3 ⁇ 7.33 ⁇ l ⁇ "3 7t radian (0.46 + 0.47 pixels) for the brain.
  • phase-shift maps may be evaluated by applying them to correct the geometric distortion in the images in FIGURE 5 using a simulated phase evolution rewinding (SPHERE) algorithm. This evaluation of the phase-shift maps may be compared using the following image correction techniques.
  • SPHERE simulated phase evolution rewinding
  • phase-shift maps 601, 621 from PLUS in FIGURES 6A and 6D are applied to the averaged unfolded images 501, 521in FIGURES 5A and 5C.
  • the corrected images 701, 721 are shown in FIGURES 7A and 7C.
  • These images 701, 721 are typical images resulting by performing the PLUS method on the averaged unfolded images in accordance with the present invention. Results in the other categories below are shown for comparison purposes.
  • phase-shift maps 651 , 671 from PLACE in FIGURES 6B and 6E are applied to the averaged unfolded images 501 , 521in FIGURES 5A and 5C.
  • the corrected images 751 , 771 are shown in FIGURES 7B and 7D.
  • phase-shift maps 601, 621 from PLUS in FIGURES 6A and 6D are applied to the regular EPI images 551, 571 in FIGURES 5B and 5D.
  • the corrected images 801 , 821 are shown in FIGURES 8A and 8C.
  • phase-shift maps 651 , 671 from PLACE in FIGURES 6B and 6E are applied to the regular EPI images 551, 571 in FIGURES 5B and 5D.
  • the corrected images 851, 871 are shown in FIGURES 8B and 8D.
  • These images 851 , 871 are typical images resulting by performing the PLACE method in concert with the SPHERE correction.
  • FIGURE 9A a 128-shot conventional gradient echo image 901 of a phantom is shown
  • FIGURE 9C a 128-shot conventional gradient echo image 921 of a brain is shown
  • FIGURE 9B a Tl -weighted MP-RAGE image 951 of a phantom is shown
  • FIGURE 9D a Tl -weighted MP-RAGE image 971 of a brain is shown.
  • the corrected images 701 , 751 in FIGURES 7A and 7B are comparable. That is, the correction afforded by the PLUS techniques of the present invention compared to the PLACE techniques.
  • phantom images 701 , 751 in FIGURES 7A and 7B the distorted geometry of the checkerboard inside the phantom is corrected, and the correction is comparable to the reference images 901 , 951 in FIGURES 9A and 9B.
  • the upper half of the checkerboard appears noisier than the lower half.
  • the corrected brain images 721, 771 in FIGURES 7C and 7D are very comparable, with the exception of a faint skull structure 772 in FIGURE 7D that appears more contiguous than that of the corresponding skull structure 722 shown in FIGURE 1C.
  • the corrected images 801 , 851 in FIGURE 8 are also very similar and resemble the reference images 901, 951 in FIGURE 9.
  • the upper right corner 802 of the checkerboard in FIGURE 8A is slightly distorted compared to the upper right corner 852 in FIGURE 8B.
  • the brain images 821, 871 in FIGURES 8C and 8D are almost indistinguishable, except for a small part of the skull structure at the top of the image 871 in FIGURE 8D that has slightly improved image quality than the image 821 in FIGURE 8C.
  • a small area in the right hemisphere, anterior to the corpus callosum, in FIGURE 8C is also slightly brighter than the corresponding structure in FIGURE 8D,
  • the edges around the squares are sharper.
  • the difference in image quality may be observed, if at all, only at the faint skull structure around the brain.
  • the phase-shift maps used in FIGURE 8 were derived from four additional scans.
  • the phase-shift maps used in FIGURE 7 were from the PLUS technique and needed no additional scans.
  • the system and method of the present Invention employs phase labeling using sensitivity encoding (PLUS) to correct geometric distortion in EPI without using extra scans to map field inhomogeneity and coil sensitivity.
  • PLUS sensitivity encoding
  • the system and method of the present invention modifies a regular EPI trajectory and utilizes parallel imaging to create multiple images that are later used to reconstruct phase-shift maps for distortion correction. Since the phase-shift maps and the coil sensitivity are derived from the EPI measurements themselves, the method of the present invention eliminates problem arising from patient position inconsistencies between field mapping and coil sensitivity scans and the EPI measurements. In addition, with the method of the present invention, no interpolation is required since voxel positions and geometric distortions in the phase-shift maps, the coil sensitivity maps, and the measurements are similar. Therefore, the method of the present invention enables self- sufficient geometric correction of each dynamic point. Furthermore, the phase labeling using sensitivity encoding (PLUS) techniques of the present invention may be applied to both spin- echo and gradient-echo images.
  • PLUS sensitivity encoding
  • the method of the present invention eliminates ghosting artifacts by alternately assigning k-space lines that are scanned in the same direction into groups.
  • the unfolded image from each group is free from ghosting artifacts and their interferences on the object under study.
  • the average of the unfolded image is comparable to the image that is acquired using a regular EPI trajectory.
  • this averaged unfolded image has a signal-to- noise ratio (SNR) equal to ⁇ !g p of the SNR of the image with a full FOV, where g p is a geometry factor (g factor) at a pixel position p and it is greater than or equal to 1.
  • SNR signal-to- noise ratio
  • EPI trajectories for PLUS techniques may be created in many ways.
  • FIGURES 1OA and 1OB two additional EPI trajectories 1010, 1020 are shown.
  • the k-space lines in FIGURES 1OA and 1OB may be rearranged into 2 and 3 groups, respectively, according to the line order and direction shown as line styles (thick solid, thick dashed, and thick dash-dot lines).
  • lines Ll 1, L31, L51, L71 would be grouped together, and lines L21, L41, L61 , L81 would be grouped together.
  • lines Ll I l, L41 1, L71 1 would be grouped together, lines L21 1, L51 1 would be grouped together, and lines L31 1, L61 1 would be grouped together.
  • These groups create images that are folded 2 and 3 times, respectively, according to the number of groups.
  • the images after unfolding will be equivalent to the images with 0 blips added into the pre-phase-encoding gradient for the group L 11 , L31 , L51 , L71 and -2 additional blips added into the pre-phase-encoding gradient for the group L21, L41, L ⁇ l, L81.
  • the images after unfolding will be equivalent to the images with 0 blips added into the pre-phase-encoding gradient for the group Ll 1 1, L41 1 , L711, and -1 blips added into the pre-phase-encoding gradient for the group L311, L611, and +1 additional blips for the group L211 , L51 1, respectively.
  • the phase labeling using sensitivity encoding (PLUS) techniques of the present invention eliminate phase errors by reading trajectory lines in the same direction.
  • the alternative EPI trajectories illustrated in FIGURE 1OA includes lines Ll 1, L31, L51, L71 in one group that were read from left to right and lines L21, L41, L61, L81 in the other group were read from right to left. So, although the numbers of folds for the trajectories in FIGURES 1OA and 10 B are less than that of the trajectory in FIGURE 1, these trajectories have less potential for improvements using PLUS techniques.
  • the phases accumulated along the readout direction will not be canceled out when calculating the phase differences. As a consequence, the phase-shift map will likely contain positive and negative phase errors on the left and the right sides of the map.
  • each group contains k-space lines that were read in both directions. Therefore, ghosting artifacts may not be eliminated, resulting in a noisier phase- shift map. Therefore, additional ghosting artifact suppression techniques may improve the measurement accuracy when alternative trajectory lines are employed.
  • the alternative trajectory lines shown in FIGURES 1OA and 1OB scan k-space without backtracking as in previous work. Other alternative trajectories may be employed with the system and method of the present invention as well, depending upon the object under study, the type of acquisition performed, and the scan parameters chosen. If backtracking the k-space is employed, however, this may reduce the efficiency of the acquisition.
  • imaging sequences may use parallel imaging with a reduction factor less than or equal to one. Since the PLUS techniques utilize parallel imaging with a high reduction factor, the product of the reduction factors will ensure that the parallel imaging limit is not exceeded. That is, the number of coils or the noise level in the images will be sufficient for accurate phase-shift calculations.
  • the signal-to-noise ratio (SNR) of the images after distortion correction is limited by geometry factors described above. Even though the image resolution in the final PLUS images is comparable to the typical resolution used in the field, to efficiently employ PLUS techniques, a larger field-of-view (FOV) may be recommended.
  • FOV field-of-view
  • strong blip gradients may be used. These stronger blip gradients may also increase the eddy current distortion in the image, so eddy- current compensation should be properly calibrated.
  • the present invention presents a significant advancement over previous geometric distortion correction techniques by removing the need for extra field map scans.
  • the EPI scans and other fast MRI scans in accordance with the present invention maximize the benefit of shortened scan times while reducing artifacts.

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Vascular Medicine (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

La présente invention concerne un étiquetage de phase utilisant un système et procédé de codage de sensibilité en vue de corriger une distorsion géométrique causée par une inhomogénéité de champ magnétique en imagerie écho-planaire (IEP). Selon l'invention, des déphasages locaux obtenus directement de la mesure IEP elle-même sont utilisés sans avoir besoin de balayages de carte de champ ou de cartes de sensibilité de bobine supplémentaires. Le système et le procédé font appel à une imagerie parallèle et à une modification de trajectoire de l'espace k pour produire de multiples images à partir d'une acquisition unique. La mesure IEP est également utilisée pour obtenir des cartes de sensibilité en vue d'une reconstruction en imagerie parallèle. Les déphasages obtenus sont rétrospectivement appliqués à la mesure IEP en vue de corriger une distorsion géométrique dans la mesure elle-même.
EP09739969A 2008-05-02 2009-05-01 Étiquetage de phase utilisant le codage de sensibilité : acquisition de données et reconstruction d'image pour correction de distorsion géométrique en imagerie écho-planaire Withdrawn EP2288941A4 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US5005208P 2008-05-02 2008-05-02
PCT/US2009/042593 WO2009135167A1 (fr) 2008-05-02 2009-05-01 Étiquetage de phase utilisant le codage de sensibilité : acquisition de données et reconstruction d'image pour correction de distorsion géométrique en imagerie écho-planaire

Publications (2)

Publication Number Publication Date
EP2288941A1 true EP2288941A1 (fr) 2011-03-02
EP2288941A4 EP2288941A4 (fr) 2011-12-21

Family

ID=41255455

Family Applications (1)

Application Number Title Priority Date Filing Date
EP09739969A Withdrawn EP2288941A4 (fr) 2008-05-02 2009-05-01 Étiquetage de phase utilisant le codage de sensibilité : acquisition de données et reconstruction d'image pour correction de distorsion géométrique en imagerie écho-planaire

Country Status (2)

Country Link
EP (1) EP2288941A4 (fr)
WO (1) WO2009135167A1 (fr)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130076356A1 (en) * 2011-09-26 2013-03-28 Vladimir Jellus Magnetic resonance imaging method and apparatus to correct distortions due to inhomogeneities of the basic magnetic field

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102011077197B4 (de) * 2011-06-08 2013-05-16 Siemens Aktiengesellschaft Verzeichnungskorrektur bei einer Magnetresonanz-Bildgebung
EP3004908B1 (fr) * 2013-06-06 2021-04-28 Koninklijke Philips N.V. Irm parallèle avec correction de distorsion b0 et separation eau-graisse de type dixon au moyen d'une reconstruction sense regularisée
US10302726B2 (en) 2014-07-03 2019-05-28 Koninklijke Philips N.V. Image reconstruction for MRI using multiplexed sensitivity encoding

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5498963A (en) * 1994-08-05 1996-03-12 General Electric Company Correction of geometric distortion in MRI phase and phase difference images
US6144772A (en) * 1998-01-29 2000-11-07 Canon Kabushiki Kaisha Variable compression encoding of digitized images
EP1242832A1 (fr) * 1999-12-17 2002-09-25 Maxygen, Inc. Procedes de detection parallele de compositions ayant des caracteristiques desirees par spectroscopie irm
US6841998B1 (en) * 2001-04-06 2005-01-11 Mark Griswold Magnetic resonance imaging method and apparatus employing partial parallel acquisition, wherein each coil produces a complete k-space datasheet
US7945305B2 (en) * 2005-04-14 2011-05-17 The Board Of Trustees Of The University Of Illinois Adaptive acquisition and reconstruction of dynamic MR images
EP2041595A2 (fr) * 2006-07-06 2009-04-01 Koninklijke Philips Electronics N.V. Dispositif et procede d'imagerie par resonance magnetique

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130076356A1 (en) * 2011-09-26 2013-03-28 Vladimir Jellus Magnetic resonance imaging method and apparatus to correct distortions due to inhomogeneities of the basic magnetic field

Also Published As

Publication number Publication date
EP2288941A4 (fr) 2011-12-21
WO2009135167A1 (fr) 2009-11-05

Similar Documents

Publication Publication Date Title
US20110260726A1 (en) Phase labeling using sensitivity encoding: data acquisition and image reconstruction for geometric distortion correction in epi
Zaitsev et al. Point spread function mapping with parallel imaging techniques and high acceleration factors: Fast, robust, and flexible method for echo‐planar imaging distortion correction
In et al. Highly accelerated PSF-mapping for EPI distortion correction with improved fidelity
US7372269B2 (en) Magnetic resonance imaging method and apparatus
US7619411B2 (en) Generalized method for MRI chemical species separation using arbitrary k-space trajectories
CN104685368B (zh) 用于耐金属mr成像的方法和设备
US8085041B2 (en) Three-point method and system for fast and robust field mapping for EPI geometric distortion correction
CN104067137B (zh) 具有对流动伪影的抑制的mr成像
CN105308469A (zh) 具有B0失真校正的并行MR成像以及使用规则化的SENSE重建的多点Dixon水-脂肪分离
JP2004275198A (ja) 核磁気共鳴イメージング装置、および核磁気共鳴イメージング方法
CN104204838A (zh) 具有Dixon类型水/脂肪分离和关于主磁场的不均匀性的先验知识的MRI
JP7128972B2 (ja) 歪み補正を伴うｅｐｉ ｍｒ画像
US11009577B2 (en) System and method for Nyquist ghost correction in medical imaging
US10746832B2 (en) DTI with correction of motion-induced diffusion gradient inconsistency
Chung et al. An improved PSF mapping method for EPI distortion correction in human brain at ultra high field (7T)
EP2288941A1 (fr) Étiquetage de phase utilisant le codage de sensibilité : acquisition de données et reconstruction d'image pour correction de distorsion géométrique en imagerie écho-planaire
US11474178B2 (en) Method for generating a magnetic resonance image
Rodriguez et al. 3D FRONSAC with PSF reconstruction
Wild et al. Normalisation of metabolite images in 1H NMR spectroscopic imaging
Hsu et al. Correction for susceptibility-induced distortion in echo-planar imaging using field maps and model-based point spread function
JP7559763B2 (ja) ウェーブ符号化を用いたパラレルｍｒイメージング
Börnert 2D-RF-pulse-encoded curved-slice imaging
Techavipoo et al. Geometric distortion correction in EPI by phase labeling using sensitivity encoding (PLUS)
WO2019016245A1 (fr) Imagerie rm multi-tranche parallèle utilisant un calcul de moyenne de signal
Techavipoo et al. A robust approach for correction of susceptibility artifacts in EPI

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20101122

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO SE SI SK TR

AX Request for extension of the european patent

Extension state: AL BA RS

DAX Request for extension of the european patent (deleted)
RIC1 Information provided on ipc code assigned before grant

Ipc: G01R 33/561 20060101ALI20111013BHEP

Ipc: G01R 33/565 20060101ALI20111013BHEP

Ipc: G01V 3/00 20060101AFI20111013BHEP

A4 Supplementary search report drawn up and despatched

Effective date: 20111117

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN

18D Application deemed to be withdrawn

Effective date: 20111201