WO2009037711A2 - Procédés de reconstruction d'images pour des sondes irm - Google Patents

Procédés de reconstruction d'images pour des sondes irm Download PDF

Info

Publication number
WO2009037711A2
WO2009037711A2 PCT/IL2008/001261 IL2008001261W WO2009037711A2 WO 2009037711 A2 WO2009037711 A2 WO 2009037711A2 IL 2008001261 W IL2008001261 W IL 2008001261W WO 2009037711 A2 WO2009037711 A2 WO 2009037711A2
Authority
WO
WIPO (PCT)
Prior art keywords
functions
voxel
phase
phase encoding
data points
Prior art date
Application number
PCT/IL2008/001261
Other languages
English (en)
Other versions
WO2009037711A3 (fr
Inventor
Yuval Zur
Original Assignee
Topspin Medical (Israel) Ltd.
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 Topspin Medical (Israel) Ltd. filed Critical Topspin Medical (Israel) Ltd.
Publication of WO2009037711A2 publication Critical patent/WO2009037711A2/fr
Publication of WO2009037711A3 publication Critical patent/WO2009037711A3/fr

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/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
    • G01R33/3808—Magnet assemblies for single-sided MR wherein the magnet assembly is located on one side of a subject only; Magnet assemblies for inside-out MR, e.g. for MR in a borehole or in a blood vessel, or magnet assemblies for fringe-field MR
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/42—Details of probe positioning or probe attachment to the patient
    • A61B8/4209—Details of probe positioning or probe attachment to the patient by using holders, e.g. positioning frames
    • 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/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/32—Excitation or detection systems, e.g. using radio frequency signals
    • G01R33/34—Constructional details, e.g. resonators, specially adapted to MR
    • G01R33/341—Constructional details, e.g. resonators, specially adapted to MR comprising surface coils
    • 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/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
    • G01R33/383—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using permanent magnets
    • 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/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
    • G01R33/385—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using gradient magnetic field coils
    • 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
    • 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
    • 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/5617—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 RF refocusing, e.g. RARE
    • 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/56518—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to eddy currents, e.g. caused by switching of the gradient magnetic field
    • 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/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
    • 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/56572—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of a gradient magnetic field, e.g. non-linearity of a gradient magnetic field

Definitions

  • the present invention in some embodiments thereof, relates to self-contained MRI probes and methods of image reconstruction for them, and, more particularly, but not exclusively, to methods of image reconstruction for probes which have substantially nonuniform field gradients, and/or have limited signal to noise ratio (SNR), which is often the case for MRI probes small enough to use inside the body.
  • SNR signal to noise ratio
  • Conventional medical MRI systems have a field of view with very uniform static magnetic field in the bore of a large magnet, and use gradient coils that produce highly linear gradient fields in three dimensions, for gradient encoding of images.
  • Conventional MRI systems sometimes use RF receiver probes that can be used inside blood vessels, the rectum (for prostate imaging) and other body cavities, in order to improve signal to noise ratio (SNR) when imaging regions adjacent to the probe, as described, for example, by US patent 5,699,801 to Atalar, and in US patent 5,476,095 to Schnall et al.
  • SNR signal to noise ratio
  • US patent 5,304,930 to Crowley et al describes a larger self-contained "inside-out" MRI device, designed to be used outside the body, to image a part of the body.
  • the static magnetic field is not nearly uniform, but has a substantial gradient in one direction, for example the x-direction, which is used for gradient encoding of images.
  • This field gradient is required to be very linear, and the surfaces of constant field are required to be very flat, in the y-z plane.
  • K. P. Pruessmann et al "SENSE: Sensitivity Encoding for Fast MRI," Magnetic Resonance in Medicine 42:952-962 (1999), describes a method of using multiple RF receiving antennas in conventional MRI in a very uniform magnetic field and very linear gradient fields.
  • the different receiving antennas have different spatial patterns of sensitivity, which can be used to decrease acquisition time when used in conjunction with gradient phase encoding.
  • Pruessmann et al describe how to find a reconstruction matrix which is optimized for a given set of receiving antennas, in that it tranforms MRI data in phase space into voxel densities, for voxels which are as close as possible to a set of ideal voxel functions, but can be expressed as a linear combination of phase encoding functions, including the spatial sensitivity of the receiving coils.
  • An aspect of some embodiments of the invention concerns methods of MRI image reconstruction which provide improved accuracy when the gradient fields are mildly or very nonlinear functions of position, and/or provide improved SNR when SNR is marginal.
  • a method of reconstructing an image from a set of MRI data points in phase space acquired from an imaging slice, using phase encoding with one or more gradient fields that are each nonlinear as functions of position in a phase encoding direction comprising: a) multiplying a vector of the data points by a reconstruction matrix which differs from a fourier transform matrix at least due to taking into account the nonlinearity of the one or more gradient fields; and b) obtaining a vector of voxel densities for the image.
  • the reconstruction matrix is the product of an encoding matrix and the inverse of a correlation matrix of phase encoding functions for the phase encoding with the one or more nonlinear gradient fields, each element of the phase encoding matrix being an integral, over at least one phase encoding direction, of a product of one of the phase encoding functions and one of a set of ideal voxel functions, each of which is compact around a different grid point in at least one of the one or more phase encoding directions.
  • the ideal voxel functions are delta functions, boxcar functions, half-sin functions, or cubic spline functions of position in at least one of the one or more phase encoding directions.
  • the MRI data points consist of K points with phase encoding differing only in characteristics of a single gradient field pulse which varies in one phase encoding direction
  • the voxel densities consist of densities of N voxel functions of the one phase encoding direction
  • the reconstruction matrix is a K by N matrix
  • the vector of the voxel densities is the product of the vector of the K data points, and the reconstruction matrix for the one phase encoding direction.
  • the imaging slice has two phase encoding directions
  • the MRI data points are phase encoded by two gradient fields, one for each of the directions.
  • the MRI data points comprise an array of Kx by Ky points in two dimensional phase space
  • the voxel densities comprise an array of Nx by Ny densities for an array of two-dimensional voxel functions each of which is a product of one of a set of Nx voxel functions of x, a first of the two phase encoding directions, and one of a set of Ny voxel functions of y, a second of the two phase encoding directions
  • transforming the data points comprises: a) multiplying a matrix of the Kx by Ky data points by a Kx by Nx reconstruction matrix for the first of the two phase encoding directions; and b) multiplying the product by a Ky by Ny reconstruction matrix for the second of the two phase encoding directions, to obtain an Nx by Ny matrix of the voxel den
  • the one-dimensional voxel functions of x have a minimum of a measure of difference from a set of Nx ideal compact voxel functions localized around an array of grid points, for any set of Nx one-dimensional functions that can be expressed as linear combinations of Kx phase encoding functions of x for Kx evenly spaced amplitudes of the gradient field for x, the reconstruction matrix for x being a linear transformation from the phase encoding functions to the voxel functions of x, and wherein the one-dimensional voxel functions of y have a minimum of a measure of difference from a set of Ny ideal compact voxel functions localized around an array of grid points, for any set of Ny one- dimensional functions that can expressed as linear combinations of Ky phase encoding functions of y for Ky evenly spaced amplitudes of the gradient field for y, the reconstruction matrix for y being a linear transformation from the phase encoding functions to the voxel functions of y.
  • the MRI data points comprise K data points, each characterized by a different pair (kx,ky), the voxel functions comprise N voxel functions, each peaked at a different pair of values (x,y) of the two dimensions x and y, and transforming the data points comprises multiplying a vector of the K data points by a K by N reconstruction matrix for the two dimensions, to obtain a vector of the N voxel densities in the two dimensions.
  • a method of calculating an image reconstruction matrix for MRI data acquired from an imaging slice using phase encoding with one or more gradient fields that are each a nonlinear function of position in a phase encoding direction comprising: a) choosing a set of ideal voxel functions; b) calculating a set of phase encoding functions, one for each value of each of the one or more gradient fields at which the MRI data are acquired; and c) finding a set of voxel functions, each of which has a minimum of a measure of difference from one of the ideal voxel functions, subject to the constrain that the voxel function can be expressed as a linear combination of the phase encoding functions; wherein the reconstruction matrix is the linear transformation from the phase encoding functions to the voxel functions.
  • finding the set of voxel functions comprises: a) calculating an encoding matrix, each element of which is an integral over the one or more phase encoding directions of one of the phase encoding functions and one of the ideal voxel functions; b) calculating a correlation matrix of the phase encoding functions; c) multiplying the encoding matrix by the inverse of the correlation matrix, thereby obtaining the reconstruction matrix which is a linear transformation from the phase encoding functions to the voxel functions.
  • a method of reconstructing an image in an imaging slice from a set of MRI data points in phase space acquired using phase encoding by one or more gradient fields, at least one of which has a field gradient that is non-uniform in magnitude or direction or both, in the imaging slice comprising: a) transforming the set of data points to a set of voxel densities, each voxel being identified by a different set of values of the one of more gradient fields; b) for each voxel, determining a position in the imaging slice where the one or more gradient fields have the values in the set identifying that voxel; and c) reconstructing the image by calculating a density as a function of position in the imaging slice, using the voxel densities and the positions in the imaging slice of the voxels.
  • calculating a density as a function of position in the imaging slice comprises adjusting one or more of the voxel densities responsive to an estimated volume of the voxel calculated from a gradient of each of the one or more gradient fields at the position of the voxel.
  • a method of generating an MRI image of an imaging region with reduced noise comprising: a) applying to one or more slices of the imaging region NMR pulse sequences suitable for separating even parity and odd parity echoes, for a plurality of points in phase space; b) accumulating NMR signals for each point in phase space for each slice, and separating the accumulated signals into a set of even parity data points and a set of odd parity data points in phase space; c) identifying an expected phase of a signal part of each data point; d) subtracting from each data point a part that is expected to be substantially noise because it has phase orthogonal to the expected phase of the signal part; and e) using the sets of even parity and odd parity data points after subtracting the parts expected to be noise, to reconstruct an image with reduced noise.
  • the NMR pulse sequences comprise, for a point in phase space for a slice, a first NMR pulse sequence consisting of an excitation RF pulse, followed by one or more gradient pulses, followed by a train of refocusing pulses with substantially zero phase difference from the excitation pulse, and a second NMR pulse sequence consisting of an excitation pulse, followed by one or more gradient pulses, followed by a train of refocusing pulses with substantially 90 degree phase difference from the excitation pulse, and separating the signals into sets of even and odd parity data points comprises finding a sum of the accumulated signals from the first and second pulse sequences for each data point, and finding a difference of the accumulated signals from the first and second pulse sequences for each data point.
  • using the sets of even parity and odd parity data points to reconstruct an image comprises adding data of one parity to data of the other parity that has been flipped in each phase encoding direction, either in phase space or in position space.
  • identifying an expected phase of the signal part of a data point comprises finding phases of one or more data points for a plurality of different resonance frequencies, fitting the phases to a function of resonance frequency, and identifying the value of the fitted function of resonance frequency as the expected phase for all data points in phase space for that resonance frequency.
  • a method of MRI imaging comprising: a) applying an NMR pulse sequence to one or more imaging slices of an imaging region; b) acquiring a set of MRI data points from the one or more imaging slices; and c) reconstructing an image of the imaging region using any of the methods according to an exemplary embodiment of the invention.
  • an MRI imaging system comprising: a) a self-contained MRI probe comprising a magnet assembly that produces a static magnetic field in an imaging region, one or more sets of gradient coils that each produce a magnetic gradient field in the imaging region, and one or more RF antennas that together produce RF fields that excite nuclei in the imaging region and receive NMR signals from the excited nuclei; and b) a controller that directs the probe to produce pulse sequences of gradient fields and RF fields in the one or more slices of the imaging region, and that accumulates NMR signals to produce MRI data points in phase space, and that reconstructs an image from the MRI data points according to any of the methods according to an exemplary embodiment of the invention.
  • FIG. 1 is a flowchart showing a phase-cycling method of reconstructing an MRI image with reduced noise, according to an exemplary embodiment of the invention
  • FIG. 2 schematically shows an MRI pulse sequence that could be used for the phase- cycling method of FIG. 1 ;
  • FIG. 3 is flowchart for a method of finding a reconstruction matrix that can be used for MRI data acquired using phase encoding in a nonlinear gradient field, according to an exemplary embodiment of the invention
  • FIG. 4 schematically shows a one-dimensional density profile and a nonlinear gradient field profile, to illustrate the method of FIG. 3;
  • FIG's. 5A and 5B schematically show voxel functions found using the method of FIG. 3 with the gradient field profile of FIG. 4;
  • FIG. 6 schematically shows a plot of simulated MRI data points for the density and gradient field profiles of FIG. 4
  • FIG. 7 schematically shows a reconstructed density profile, found for the density profile and gradient field profile of FIG. 4, using a reconstruction matrix found using the method of FIG. 3, according to an exemplary embodiment of the invention
  • FIG. 8 schematically shows a two-dimensional object, to illustrate reconstructing two-dimensional MRI images using a reconstruction matrix found by the method of FIG. 3, according to an exemplary embodiment of the invention
  • FIG. 9 schematically shows density profiles reconstructed in two different ways from simulated MRI data of the object of FIG. 8, according to an exemplary embodiment of the invention.
  • FIG. 10 schematically shows contours of constant gradient field, for each of two gradient fields that are nonlinear functions of two gradient encoding directions, to illustrate an image reconstruction method according to an exemplary embodiment of the invention.
  • FIG. 11 schematically shows an MRI system which uses methods for image reconstruction according to an exemplary embodiment of the invention, including one or more of the methods of FIG's. 1, 3, and 10.
  • the present invention in some embodiments thereof, relates to self-contained MRI probes and methods of image reconstruction for them, and, more particularly, but not exclusively, to methods of image reconstruction for probes which have substantially nonuniform field gradients, and/or have limited signal to noise ratio (SNR), which is often the case for MRI probes small enough to use inside the body.
  • SNR signal to noise ratio
  • An aspect of some embodiments of the invention concerns a method of using phase cycling to improve SNR in MRI, by separating MRI data points of odd and even parity, correcting their phase, removing a part which is likely to be substantially all noise because of its phase, and combining the even and odd parity data to form an image.
  • the even and odd parity data is separated by running two pulse sequences for each point in phase space and each slice, for example: a first sequence consisting of an excitation pulse followed by one or more gradient pulses and a train of refocusing pulses that are 90 degrees out of phase with the excitation pulse, and a similar second sequence in which the refocusing pulses are in phase with the excitation pulse.
  • a best fit for example a linear fit with an offset, is found for the phase as a function of resonance frequency.
  • the phases for all data points, including those further out in phase space, are then reduced by this fitted function. After this phase correction, the dominant points in phase space will generally all be nearly real.
  • phase-corrected odd and even parity data is combined by first flipping one of them with respect to each of the gradient encoding axes, either in phase space or in actual space, and then adding the even and odd parity data points or images together.
  • An aspect of some embodiments of the invention concerns a method of image reconstruction for MRI data obtained using phase encoding in one or more phase encoding directions, with a gradient field that is a nonlinear function of the position in at least one gradient encoding direction.
  • the method uses a reconstruction matrix which takes into account the nonlinearities of the gradient fields.
  • the method is optionally more computationally intensive than the usual FFT methods that are used with highly linear gradients in conventional MRI, but provides more accurate images when one or more gradient field is a substantially nonlinear function of position in a phase encoding direction.
  • the reconstruction matrix is chosen so that the associated voxel functions will match as closely as possible ideal voxel functions, for example compact functions around a set of grid points, such as delta functions or boxcar functions, subject to the constraint that the voxel functions are linear combinations of phase encoding functions associated with the nonlinear field gradients.
  • the method is used in one phase encoding direction and one gradient field at a time, either to reconstruct a one-dimensional image, or to reconstruct a two-dimensional image with gradient fields and voxel functions that can be treated to good approximation as separable in the two gradient encoding directions.
  • the method is used in a fully two-dimensional way, to reconstruct a two- dimensional image.
  • An aspect of some embodiments of the invention concerns a method of calculating an image reconstruction matrix, for use with MRI data acquired using phase encoding with one or more gradient fields that are nonlinear functions of position in the phase encoding direction.
  • An aspect of some embodiments of the invention concerns a method of image reconstruction, for MRI data acquired using phase encoding with one or more gradient field that is a nonlinear function of position in the phase encoding direction, taking into account changes in the positions that different voxels would have if the gradient fields were linear functions of position.
  • the method can be used if a reconstruction matrix is used to reconstruct the image, as described above, or it can be used to correct distortions in images reconstructed by conventional methods such as FFT which do not take into account the nonlinearity of the gradient fields.
  • FIG. 1 is a flowchart showing a method of using phase-cycling to reduce noise in MRI images, according to an exemplary embodiment of the invention.
  • the method relies on the use of pulse sequences that can be used to separate MRI data into echoes of even parity, and odd parity.
  • pulse sequences that can be used to separate MRI data into echoes of even parity, and odd parity.
  • a first pulse sequence is applied to the imaging region, and data is accumulated.
  • a second pulse sequence is applied, and data is accumulated, in 104.
  • the first pulse sequence is a CPMG type of pulse sequence.
  • An excitation pulse 902 at frequency fi selecting a particular imaging slice, nominally with 90 degree flip angle, is applied first. This is followed by a gradient pulse 904, followed by a refocusing pulse 906 at frequency fi at a time ⁇ following the excitation pulse, at a phase that differs by 90 degrees from the phase of the excitation pulse.
  • an echo 908 occurs, and is accumulated as part of the MRI data for this value of the integrated strength of gradient pulse 906.
  • a train of further refocusing pulses 910 is applied at frequency fi, at time intervals of 2 ⁇ , and echoes 912 are received and accumulated half-way between successive refocusing pulses.
  • Train 910 continues for a time T, which may be adjusted to produce a desired degree of T 2 or diffusion weighting. Alternatively, T 2 or diffusion time Tj is found directly by measuring the decay rate of echoes in train 910, so the weighting does not depend on T.
  • Train 910 is followed by another series of pulses, at a different frequency f 2 to excite a different imaging slice, including an excitation pulse 914, a gradient pulse 916, and refocusing pulses 918, 90 degrees out of phase with the excitation pulse, with echoes 920 halfway between the refocusing pulses.
  • a second pulse sequence 922 is applied at frequency fi, with the same integrated strength of the gradient pulse, and the echoes are accumulated into an MRI data point, but this time the refocusing pulses are in phase with the excitation pulse.
  • a series of pulses, each extending for a time T, at different frequencies, is followed after a time T R by another first type of pulse sequence at frequency fi, this time with a different integrated strength of the gradient pulse, representing a different point in phase space, generally called k-space in the case of linear field gradients.
  • Data from the different imaging slices, from the different points in phase space, and from pulse sequences of the first or second type, is accumulated separately.
  • the time interval between the excitation pulse and the first refocusing pulse which is also the time interval between the first refocusing pulse and the first echo, need not be the same interval ⁇ as the interval between the other refocusing pulses in the train, but can be several times longer, as explained in the co-filed application "MRI Magnet and Coil Configurations.” This may allow more time for the gradient pulse, which can result in less heating of the gradient coil for the same integrated gradient field strength.
  • a second gradient pulse may be applied during the same longer time interval between the first refocusing pulse and the beginning of the train of refocusing pulses.
  • Pulse sequences of the first and second type need not be applied consecutively for the same point in phase space, but, for example, pulses of the first type maybe applied at all points in phase space, before moving on to pulses sequences of the second type., or any other pattern of interleaving may be used.
  • Odd parity data is data from echoes emitted by components of the excited nuclei that have undergone an odd number of reversals in phase due to refocusing pulses
  • even parity data is data from echoes emitted by components of the excited nuclei that have undergone an even number of reversals in phase sue to refocusing pulses.
  • the odd parity and even parity components of the data are optionally obtained by first adding together the data from the first and second types of pulse sequences, for a given point in phase space and a given slice, and then subtracting them.
  • the odd and even data is optionally separately corrected for phase errors.
  • phase errors are likely to be primarily due to a shift in the timing of the echo, due to errors in the timing of the pulses, and these phase errors will be linear with resonance frequency.
  • the real part of the data is taken, in order to reduce noise in the data, since the imaginary parts of phase corrected data points are likely to be only due only to noise.
  • the even and odd parity data is used to reconstruct a combined image incorporating all the data.
  • one parity of data is reflected, in all the phase encoding directions, before adding it to the other parity of data, because odd echoes are reversed in phase relative to even echoes. This is done, for example, by reconstructing separate images from the even parity and odd parity data, reflecting one of the images along each of the phase encoding axes, and then adding the two images together.
  • data of one of the parities is reflected in phase space along each of the phase encoding directions, and the data points of the two parities are then added together. The resulting data points are then used to reconstruct an image.
  • the image reconstruction can be done using any known method of reconstructing an image from phase encoded MRI data.
  • a method such as FFT may be used.
  • FFT Fast Fourier transform
  • other methods such as those described below, may yield more accurate images.
  • FIG. 3 is a flowchart 300, describing a method for finding a reconstruction matrix which can be used to convert a vector of phase encoded MRI data points into a vector of voxel densities, taking into account nonlinearities in the gradient fields that are used to phase encode the data.
  • the image reconstruction algorithm may be more computationally intensive then FFT, which is used when the gradient fields are highly linear, because it is much faster than multiplying an array of data points in phase space by an arbitrary reconstruction matrix, to find an array of voxel densities.
  • the reconstruction matrix may yield more accurate results than FFT, and other conventional methods of image reconstruction, especially with high nonlinear gradient fields, when FFT may fail completely. Calculating the reconstruction matrix may be computationally intensive, but only has to be done once for a given set of gradient field profiles.
  • the method will be described first for a one-dimensional case, with only a single phase encoding direction y, and voxel function depending only y.
  • the method may be used for more than one phase encoding direction, and two dimensional voxel functions, simply by replacing the single indexes for points in phase space, and grid points in configuration space, by pairs of indexes, representing two phase encoding directions, and replacing integrals over the single phase encoding direction y by integrals over both phase encoding directions.
  • Another way of using the reconstruction matrix for two-dimensional imaging which is much less computational intensive and which may be used when the two directions are separable to good approximation, will be described afterwards.
  • phase encoding function e ⁇ y) exp[i ⁇ -kGTh S -y].
  • the phase encoding functions are used to weight the spin density p(y) to produce the data points m ⁇ for each k. It should be understood that the spin density p(y) refers not the total density of hydrogen nuclei (or whatever species is being imaged), but the density weighted by Ti, T 2 , or diffusion, or any other weighting, whatever the case may be.
  • the weighting depends on parameters of the pulse sequence used to obtain the data points, or on analysis of the echoes that are received at different times during the pulse sequence, as is known in the art MRI.
  • the reconstruction methods to be described here start with the data points mk, regardless of how they are obtained, so can be used with any type of weighting.
  • ideal voxel functions I n (y) are chosen for each of N grid points indexed by n.
  • the ideal voxel functions may be compact functions centered around grid points y n , for example delta functions ⁇ (y - y n ), or boxcar functions of y - y n , or half-sin functions, or cubic splines.
  • a goal is to find actual voxel functions f n (y) which are as close as possible to the ideal voxel functions, given the constraint that they are linear combinations of the phase encoding functions Q ⁇ y).
  • the actual voxel functions depend on an N by K reconstruction matrix F n ⁇ which relates the data points m k to voxel densities V n of the reconstructed image:
  • a K by K correlation matrix C of the encoding functions is optionally found, defined by
  • the reconstruction matrix F is found at 310, optionally from
  • the voxel functions f n (y) satisfy the condition that for each n, is minimized, so the voxel functions are as close as possible to the ideal voxel functions (with this measure of the difference between them) subject to the constraint that they are expressed as linear combinations of the K phase encoding functions.
  • another criterion is used to choose F, for example the measure of closeness of the voxel functions to the ideal voxel functions uses absolute value instead of absolute value squared, or the measure is weighted in some way, or instead of minimizing the measure, it is only made adequately small by some criterion.
  • Eq. (9) as a measure of closeness which is to be minimized has the potential advantage that F can be expressed in the simple form given by Eq. (8).
  • solving Eq. (3) may give the same result as taking a fourier transform of the data points ⁇ ik to find the voxel densities V n .
  • This form of F is referred to herein, including in the claims, as a fourier transform matrix.
  • finding the voxel densities may be done with a considerably faster computation by using FFT software to calculate the fourier transform, rather than by using general matrix multiplication software which does not take advantage of the special relations between the elements of F given by Eq. (10).
  • the voxel functions are found using Eq. (5).
  • the voxel functions are sine functions of ⁇ -k max GT-(y - y n ), where k max is the maximum value of k for the data points.
  • FIG's. 4, 5, 6, and 7 illustrate one-dimensional image reconstruction using a reconstruction matrix F found by this procedure for the nonlinear gradient field profile 402 shown in FIG. 4, in an imaging region 404 extending from y m j n to y max .
  • the density profile 406 is flat in a central portion 408 of imaging region 404, and zero elsewhere.
  • the time interval T hS that the gradient field is turned on is chosen so that ⁇ i(y) varies by 2 ⁇ over the width of imaging region 404,
  • T Using a value of T of about this magnitude has the potential advantage of avoiding aliasing within the imaging region, which could occur if larger T were used, and not spending a large part of the time sampling data from parts of k space that do not contribute very much to the reconstructed image, which could occur if much smaller T were used.
  • the number N of voxels is 60, with grid points y n uniformly distributed in the imaging region, and delta functions ⁇ (y - y n ) used as the ideal voxel functions.
  • FIG. 7 shows a plot 700 of the reconstructed image density, as a superposition of the voxel functions f n (y), each one weighted by the voxel density V n .
  • the reconstructed image density resembles the actual density 406 shown in FIG. 4, but the edges are not as sharp because of the limited number of points in phase space used to represent the image.
  • the method outlined in FIG. 3 for finding the reconstruction matrix F can be used for two-dimensional imaging as well, but with each one-dimensional voxel function and ideal voxel function, replaced by two-dimensional functions, of y and z, for example, and the one-dimensional array of data points and encoding functions, indexed by k, replaced by two-dimensional encoding functions, and a two-dimensional array of data points, indexed by k y and k z .
  • the matrixes F, E, and C can be enormous, and take a long time to calculate.
  • the two-dimensional voxel functions f nyj nz(y,z) can each be taken to be the product of a y voxel function which depends only on y and the y grid index ny, and a z voxel function which depends only on z and the z grid index nz.
  • the two-dimensional encoding functions e ky,kz (y,z) may be approximated by a product of a y encoding function e (y) that depends on y and the index ky, and a z encoding function ⁇ ( Z ) that depends only on z and the index kz. Then the two-dimensional array of data points in phase space, m, is given by
  • Vfe P *" d y [ m dz p(y> z ) e ⁇ y )ky M e (z)fe M 02) where p(y,z) is the two-dimensional density.
  • the two-dimensional array of voxel densities V ny ,nz is given in matrix form by
  • Fy and Fz are separate one-dimensional reconstruction matrixes in y and z. They are each found according to the procedure described in FIG. 3, using respectively the y gradient field profile B (y) (y) which is taken to be independent of z, and the z gradient field profile B( Z )(z) which is taken to be independent of y.
  • FIG. 8 shows a simulated oval-shaped test object 802, of uniform density, in a square imaging region extending between z min and z ma ⁇ , and between y min and y max , used to illustrate two-dimensional image reconstruction using Eq. (12).
  • the y gradient field is assumed to have the same profile as in FIG. 4, and to be independent of z, while the z gradient field is assumed to have the same profile and to be independent of y.
  • the same number of grid points, and the same number of k values is used in each of y and z, and the voxel functions and reconstruction matrixes in y and z have the same forms as the one- dimensional voxels shown in FIG's.
  • reconstructed image 900 has a boundary that closely follows the boundary of the original object being imaged, but its edge is not as sharp, and it has some small oscillations near the edge, because of the limited range of ky and kz of the encoding functions used to approximate it.
  • Table 1 An examination of Table 1 shows that even with the limited number of data samples, most data samples mky ⁇ z are relatively small in amplitude and contribute little to the reconstructed image, with most of the contribution to the image coming from a few larger components of m k y ⁇ z with ky and kz smallest in absolute value.
  • data samples taken at larger absolute values of ky and kz use lower RF power, with the nominally 90 degree excitation RF pulses, and/or the nominally 180 degree refocusing RF pulses, having actual flip angles substantially less than 90 degrees and 180 degrees respectively. Although using this option reduces the SNR for those data samples, it has very little effect on the appearance and the SNR of the reconstructed image.
  • Table 2 lists the amplitudes of simulated data samples m ⁇ y , ⁇ z taken with reduced RF power for larger absolute values of ky and kz, rounded to the nearest 0.05.
  • FIG. 9 shows a plot of density for a reconstructed image 902, using these values of m Ky, ⁇ z - There is very little difference between image 902, with reduced RF power, and image 900.
  • the total RF power in the simulated data acquisition for image 802 was 12 times lower than for image 900.
  • the reduced flip angles used are chosen so that the RF power is approximately as low as it can be, while having little expected effect on the quality of the images.
  • the expected effect on the quality of the images, for a given set of flip angles is estimated by calculating simulated images having characteristics, such as contrast and scale of structure, similar to that in actual medical images, for the application that the probe is intended to be used for, for example looking for prostate cancer, or looking for plaque in arteries.
  • the voxel functions will not generally have peak values at a set of points on an equally spaced a rectangular grid, but, due to the nonlinearity and non-separability of the y and z gradient fields, the voxels will typically be localized on an array of curved lines of constant y gradient field and constant z gradient field, on an imaging slice of constant resonance frequency. This may be seen in plot 1006 in FIG. 10, which shows an imaging slice with isophase curves 1002 and 1004 for the y gradient and z gradient fields. To the extent that the gradient fields are small compared to a static magnetic field, the isophase curves on the imaging slice closely follow curves where the component of the gradient field parallel to the static magnetic field is constant.
  • the positions of each rectangular grid point labeled by (ny, nz) associated with a voxel density V ny , n z is adjusted to reflect the actual location of the grid point where the y and z gradient fields have a given value.
  • the element V ny nz is moved from its old assumed position (y ny , z nz ), in which the y coordinate depends only on ny and the z coordinate depends only on nz, and is assigned to a new position (y' ny ,nz, z' ny ,nz)-
  • the new grid points do not, in general, form a rectangular grid, with the y coordinates depending only on ny and the z coordinates depending only on nz, but the y and z coordinates of the points each depend on both ny and nz.
  • the new grid points are related to the old grid points by solving two simultaneous equations for each grid point (ny,nz), to find the two unknown coordinates y' ny,nz and z' ny nz of the new grid point:
  • Bz Vny,nz > Z n'y,nz ) ⁇ z ⁇ Z nz ) (14b)
  • B y (y) and B z (z) are the approximate values for the y-gradient and z-gradient fields used to find the separate y-direction and z-direction reconstruction matrixes used in Eq. (12).
  • the new grid point is the point where the actual gradient fields B y (y,z) and B z (y,z) have the same values as the approximate gradient fields B y (y) and B z (z) have at the old grid point.
  • Another method of two-dimensional image reconstruction, used in some embodiments of the invention, is similar to the method described above using Eq.
  • non-Cartesian coordinates Y and Z are used from the beginning, instead of y and z. Coordinates Y and Z are chosen so that the nominally y-gradient field B y (y,z) is a function only of Y, not of Z, and the nominally z-gradient field B z (y,z) is a function only of Z, not of Y. Then the method using Eq. (12), and described above, is carried through, but using Y and Z everywhere instead of y and z. The end result is a two-dimensional density profile depending on Y and Z. Optionally, the density profile is then converted to Cartesian coordinates y and z, and displayed.
  • the density is corrected to take into account the varying voxel volume, due to the nonlinearity and non-separability of the gradient fields in y and z.
  • this may be done automatically if a fully two-dimensional reconstruction matrix were used to find reconstruct the image, such a correction can also be done if separate one-dimensional reconstruction matrixes are used for y and z, as in Eq. (12), to correct for the non-separability of the y and z gradient fields, as well as if a method such as FFT is used, ignoring the nonlinearity of the field gradients, which may be feasible for field gradients that are only mildly nonlinear. As may be seen from FIG.
  • the volume of a voxel may depend both on the spacing between adjacent curves of constant y gradient field and constant z gradient field at the location of the voxel, and it will also, in general, depend on the angle that the curves of constant y gradient field make to the curves of constant z gradient field. If the resulting differences in voxel size are not taken into account, then the image will show artifacts, with apparently higher density in regions where the voxel size is bigger due to these effects.
  • FIG. 10 schematically shows an example of contours 1002 of constant B y , and contours 1004 of constant B z , as a function of y and z in an imaging slice 1006.
  • contours of constant Y and constant Z are also respectively contours of constant Y and constant Z, in the second two-dimensional image reconstruction method described above. Near the corners of the imaging slice, some contours 1002 intersect some contours 1004 twice, for example at point 1008 and at point 1010. This may cause imaging artifacts, unless at least one of the intersection points, for example point 1010, is effectively outside the imaging region, for example because the RF field is too low there, or because the RF antenna is too insensitive to signals from that point, or because there are not enough nuclei to be excited there.
  • FIG. 12 shows an MRI system 1100 which uses any of the methods described above, singly or in combination, for reconstructing images.
  • System 1 100 comprises an MRI probe 1102, optionally a self-contained MRI used inside the body for medical imaging, such as the prostate probe described in the co-filed applications "MRI Probe” and "MRI Magnet and Coil Configurations.”
  • MRI Probe optionally a self-contained MRI used inside the body for medical imaging, such as the prostate probe described in the co-filed applications "MRI Probe” and "MRI Magnet and Coil Configurations.”
  • Such probes are likely to use substantially nonlinear gradient fields, and/or have relatively low SNR, so would be especially likely to benefit from these image reconstruction methods, but the methods could be used in any MRI system using phase encoding, including a conventional MRI system using a bore magnet.
  • a controller 1104 controls the pulse sequences of probe 1 102 and receives and analyzes data from probe 1102. Alternatively, these functions could be performed by separate elements of system 1100.
  • Controller 1104 is, for example, a computer with a control interface and data interface to probe 1102. Controller 1104 is programmed, in its role as data analyzer, to reconstruct images from probe 1102, using any of the methods described above, and has a sufficient CPU, memory, and storage space to perform this task in a reasonably short time, optionally in real time.
  • an output device 1 106 such as a computer monitor or a printer, linked to controller 1104, provides output images from probe 1102, optionally in real time.
  • an image file for example an image file in some standard format, is stored, and/or transferred to another location, and used to provide an image later.
  • an input device 1 108 such as a keyboard or mouse linked to controller 1104, is used by an operator of the system, such as a physician performing a prostate examination and/or biopsy with the system, to instruct the controller when to initiate the imaging procedure, and optionally to set one or more free parameters of the imaging procedure.
  • an operator of the system such as a physician performing a prostate examination and/or biopsy with the system, to instruct the controller when to initiate the imaging procedure, and optionally to set one or more free parameters of the imaging procedure.
  • system 1 100 does not include probe 1102, but controller 1104 obtains data from an MRI probe and analyzes the data, according to any of the methods described above.
  • controller 1104 obtains data from an MRI probe and analyzes the data, according to any of the methods described above.
  • the term "about” refers to ⁇ 10 % .
  • compositions comprising, “comprising”, “includes”, “including”, “having” and their conjugates mean “including but not limited to”. This term encompasses the terms “consisting of and “consisting essentially of.
  • Consisting essentially of means that the composition or method may ⁇ include additional ingredients and/or steps, but only if the additional ingredients and/or steps do not materially alter the basic and novel characteristics of the claimed composition or method.
  • a compound or “at least one compound” may include a plurality of compounds, including mixtures thereof.
  • range format is merely for convenience and brevity and should not be construed as an inflexible limitation on the scope of the invention. Accordingly, the description of a range should be considered to have specifically disclosed all the possible subranges as well as individual numerical values within that range. For example, description of a range such as from 1 to 6 should be considered to have specifically disclosed subranges such as from 1 to 3, from 1 to 4, from 1 to 5, from 2 to 4, from 2 to 6, from 3 to 6 etc., as well as individual numbers within that range, for example, 1, 2, 3, 4, 5, and 6. This applies regardless of the breadth of the range.

Landscapes

  • Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Signal Processing (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Pathology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Vascular Medicine (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

L'invention concerne un procédé pour reconstruire une image à partir d'un ensemble de points de données IRM dans un espace de phase acquis à partir d'une tranche d'imagerie, au moyen d'un codage de phase par un ou plusieurs champs de gradient, chacun étant non linéaire en fonction des positions dans un sens de codage de phase. Selon l'invention, le procédé comprend les étapes suivantes : a) multiplication d'un vecteur de points de données par une matrice de reconstruction qui est différente d'une matrice de transformée de Fourier, au moins à cause de la non linéarité d'un ou plusieurs champs de gradient; et b) obtention d'un vecteur de densités d'éléments de volume pour l'image.
PCT/IL2008/001261 2007-09-20 2008-09-21 Procédés de reconstruction d'images pour des sondes irm WO2009037711A2 (fr)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
US96021307P 2007-09-20 2007-09-20
US96021207P 2007-09-20 2007-09-20
US60/960,212 2007-09-20
US60/960,213 2007-09-20

Publications (2)

Publication Number Publication Date
WO2009037711A2 true WO2009037711A2 (fr) 2009-03-26
WO2009037711A3 WO2009037711A3 (fr) 2009-08-13

Family

ID=40149595

Family Applications (3)

Application Number Title Priority Date Filing Date
PCT/IL2008/001261 WO2009037711A2 (fr) 2007-09-20 2008-09-21 Procédés de reconstruction d'images pour des sondes irm
PCT/IL2008/001260 WO2009037710A2 (fr) 2007-09-20 2008-09-21 Sonde d'irm
PCT/IL2008/001259 WO2009037709A2 (fr) 2007-09-20 2008-09-21 Aimant de sonde irm et configurations de bobine

Family Applications After (2)

Application Number Title Priority Date Filing Date
PCT/IL2008/001260 WO2009037710A2 (fr) 2007-09-20 2008-09-21 Sonde d'irm
PCT/IL2008/001259 WO2009037709A2 (fr) 2007-09-20 2008-09-21 Aimant de sonde irm et configurations de bobine

Country Status (1)

Country Link
WO (3) WO2009037711A2 (fr)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015164109A1 (fr) * 2014-04-24 2015-10-29 Dignity Health Système et procédé d'imagerie par résonance magnétique à blocs multiples en spirale
CN115604465A (zh) * 2022-12-14 2023-01-13 浙江荷湖科技有限公司（Cn） 基于相空间连续性的光场显微图像无损压缩方法及装置

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009130002A1 (fr) 2008-04-23 2009-10-29 Gsi Helmholtzzentrum Für Schwerionenforschung Gmbh Dispositif de production de champ magnétique
US9089314B2 (en) 2010-01-27 2015-07-28 Medtronic Cryocath Lp Partially compliant balloon device
US9119550B2 (en) * 2012-03-30 2015-09-01 Siemens Medical Solutions Usa, Inc. Magnetic resonance and ultrasound parametric image fusion
SG11202112987SA (en) 2019-06-25 2021-12-30 Promaxo Inc Systems and methods for image reconstruction in magnetic resonance imaging

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5886524A (en) * 1995-10-30 1999-03-23 Siemens Aktiengesellschaft Method of distortion correction for gradient non-linearities in nuclear magnetic resonance tomography apparatus
WO2004065977A1 (fr) * 2003-01-21 2004-08-05 Koninklijke Philips Electronics N.V. Procede de resonance magnetique avec gradients de champ magnetique non lineaires

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4862128A (en) * 1989-04-27 1989-08-29 The United States Of America As Represented By The Secretary Of The Army Field adjustable transverse flux sources
US5304930A (en) * 1993-02-01 1994-04-19 Panacea Medical Laboratories Remotely positioned MRI system
US5572132A (en) * 1995-08-15 1996-11-05 Pulyer; Yuly M. MRI probe for external imaging
US6704594B1 (en) * 2000-11-06 2004-03-09 Topspin Medical (Israel) Limited Magnetic resonance imaging device
US6400149B1 (en) * 2001-05-24 2002-06-04 Schlumberger Technology Corporation Nuclear magnetic resonance apparatus and method for generating an axisymmetric magnetic field having straight contour lines in the resonance region
EP1718194A4 (fr) * 2004-01-20 2010-05-12 Topspin Medical Israel Ltd Sonde irm destinee a l'imagerie de la prostate
US20060084861A1 (en) * 2004-10-18 2006-04-20 Topspin Medical (Isreal) Ltd. Magnet and coil configurations for MRI probes
US7372263B2 (en) * 2005-11-23 2008-05-13 Baker Hughes Incorporated Apparatus and method for measuring cased hole fluid flow with NMR
US8105233B2 (en) * 2007-10-24 2012-01-31 Tarek Ahmed Nabil Abou El Kheir Endoscopic system and method for therapeutic applications and obtaining 3-dimensional human vision simulated imaging with real dynamic convergence

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5886524A (en) * 1995-10-30 1999-03-23 Siemens Aktiengesellschaft Method of distortion correction for gradient non-linearities in nuclear magnetic resonance tomography apparatus
WO2004065977A1 (fr) * 2003-01-21 2004-08-05 Koninklijke Philips Electronics N.V. Procede de resonance magnetique avec gradients de champ magnetique non lineaires

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
CHING-MING LAI ET AL: "Reconstructing NMR images under magnetic fields with large inhomogeneities" JOURNAL OF PHYSICS E. SCIENTIFIC INSTRUMENTS, IOP PUBLISHING, BRISTOL, GB, vol. 15, no. 10, 1 October 1982 (1982-10-01), pages 1093-1100, XP020016649 ISSN: 0022-3735 *
CHING-MING LAI ET AL: "Reconstructing NMR images under nonlinear field gradients" JOURNAL OF PHYSICS E. SCIENTIFIC INSTRUMENTS, IOP PUBLISHING, BRISTOL, GB, vol. 16, no. 1, 1 January 1983 (1983-01-01), pages 34-38, XP020017001 ISSN: 0022-3735 -& CHING-MING LAI: "CORRIGENDUM; Reconstructing NMR images under nonlinear field gradients" JOURNAL OF PHYSICS E. SCIENTIFIC INSTRUMENTS, IOP PUBLISHING, BRISTOL, GB, vol. 16, no. 4, 1 April 1983 (1983-04-01), pages 355-355, XP020017070 ISSN: 0022-3735 *
LANGLOIS S ET AL: "MRI geometric distortion: A simple approach to correcting the effects of non-linear gradient fields" JOURNAL OF MAGNETIC RESONANCE IMAGING, SOCIETY FOR MAGNETIC RESONANCE IMAGING, OAK BROOK, IL, US, vol. 9, no. 6, 1 January 1999 (1999-01-01), pages 821-831, XP002399605 ISSN: 1053-1807 *
PRUESSMANN K P ET AL: "SENSE: SENSITIVITY ENCODING FOR FAST MRI" MAGNETIC RESONANCE IN MEDICINE, ACADEMIC PRESS, DULUTH, MN, US, vol. 42, no. 5, 1 November 1999 (1999-11-01), pages 952-962, XP000866655 ISSN: 0740-3194 cited in the application *
SEKIHARA K ET AL: "Image restoration from non-uniform magnetic field influence for direct Fourier NMR imaging" PHYSICS IN MEDICINE AND BIOLOGY, TAYLOR AND FRANCIS LTD. LONDON, GB, vol. 29, no. 1, 1 January 1984 (1984-01-01), pages 15-24, XP020022739 ISSN: 0031-9155 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015164109A1 (fr) * 2014-04-24 2015-10-29 Dignity Health Système et procédé d'imagerie par résonance magnétique à blocs multiples en spirale
US10222439B2 (en) 2014-04-24 2019-03-05 Dignity Health System and method for spiral multislab magnetic resonance imaging
CN115604465A (zh) * 2022-12-14 2023-01-13 浙江荷湖科技有限公司（Cn） 基于相空间连续性的光场显微图像无损压缩方法及装置

Also Published As

Publication number Publication date
WO2009037711A3 (fr) 2009-08-13
WO2009037709A2 (fr) 2009-03-26
WO2009037709A3 (fr) 2009-06-11
WO2009037710A3 (fr) 2009-05-22
WO2009037710A2 (fr) 2009-03-26

Similar Documents

Publication Publication Date Title
US7408345B2 (en) Generalized MRI reconstruction with correction for multiple image distortion
Wang et al. In vivo method for correcting transmit/receive nonuniformities with phased array coils
EP3004908B1 (fr) Irm parallèle avec correction de distorsion b0 et separation eau-graisse de type dixon au moyen d'une reconstruction sense regularisée
US6476606B2 (en) Method for parallel spatial encoded MRI and apparatus, systems and other methods related thereto
US7417427B2 (en) Magnetic resonance data acquisition method and apparatus
US7999545B2 (en) Method for obtaining amplitude and phase profiles of RF pulses for spatially selective excitation
US10180474B2 (en) Magnetic resonance imaging apparatus and quantitative magnetic susceptibility mapping method
US9638777B2 (en) Electric properties tomography imaging method and system
Sled et al. Understanding intensity non-uniformity in MRI
Sung et al. Measurement and characterization of RF nonuniformity over the heart at 3T using body coil transmission
US20040227514A1 (en) Synthetic images for a magnetic resonance imaging scanner using linear combination of source images to generate contrast and spatial navigation
US20110260726A1 (en) Phase labeling using sensitivity encoding: data acquisition and image reconstruction for geometric distortion correction in epi
Tam et al. Null space imaging: nonlinear magnetic encoding fields designed complementary to receiver coil sensitivities for improved acceleration in parallel imaging
US7235971B2 (en) Shimming of MRI scanner involving fat suppression and/or black blood preparation
JP6333293B2 (ja) 金属耐性ｍｒ画像化
US20160291106A1 (en) Mr image reconstruction using compressed sensing
US6323646B1 (en) Method and apparatus for producing diffusion weighted MR images
US20150077107A1 (en) Method for B0 Field Correction in Magnetic Resonance
US6239599B1 (en) Method and apparatus for identifying errors in magnetic resonance imaging examinations
WO2013181703A1 (fr) Méthode d'estimation de la vitesse d'absorption spécifique
CN106796274B (zh) 具有伪迹抑制的propeller-mr成像
WO2009037711A2 (fr) Procédés de reconstruction d'images pour des sondes irm
JP2006130285A (ja) 磁気共鳴イメージング装置
US20010054898A1 (en) Magnetic resonance imaging compensated for very rapid variations in static magnetic field
US10107883B2 (en) Ultrafast MRI system and method

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 08808062

Country of ref document: EP

Kind code of ref document: A2

NENP Non-entry into the national phase in:

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 08808062

Country of ref document: EP

Kind code of ref document: A2