WO2003003037A2 - Method and apparatus for enhanced multiple coil mr imaging - Google Patents

Method and apparatus for enhanced multiple coil mr imaging Download PDF

Info

Publication number
WO2003003037A2
WO2003003037A2 PCT/US2002/019275 US0219275W WO03003037A2 WO 2003003037 A2 WO2003003037 A2 WO 2003003037A2 US 0219275 W US0219275 W US 0219275W WO 03003037 A2 WO03003037 A2 WO 03003037A2
Authority
WO
WIPO (PCT)
Prior art keywords
signals
values
coils
signal
received signals
Prior art date
Application number
PCT/US2002/019275
Other languages
French (fr)
Other versions
WO2003003037A3 (en
Inventor
G. Randy Duensing
Steve Varosi
Scott B. King
Original Assignee
Mri Devices Corporation
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 Mri Devices Corporation filed Critical Mri Devices Corporation
Priority to EP02756225A priority Critical patent/EP1412769A2/en
Priority to AU2002322248A priority patent/AU2002322248A1/en
Publication of WO2003003037A2 publication Critical patent/WO2003003037A2/en
Publication of WO2003003037A3 publication Critical patent/WO2003003037A3/en

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/32Excitation or detection systems, e.g. using radio frequency signals
    • G01R33/34Constructional details, e.g. resonators, specially adapted to MR
    • G01R33/341Constructional details, e.g. resonators, specially adapted to MR comprising surface coils
    • G01R33/3415Constructional details, e.g. resonators, specially adapted to MR comprising surface coils comprising arrays of sub-coils, i.e. phased-array coils with flexible receiver channels
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image 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/5611Parallel 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5608Data 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

Definitions

  • the subject invention was made with government support under a research project supported by the National Institutes of Health (NIH), Department of Health and Human Services, Grant Number 5R44 RR11034-03. The government has certain rights in this invention.
  • the subject invention relates to the field of medical technology and can provide for an improved method of creating composite images from a plurality of individual signals.
  • the subject invention is particularly advantageous in the field of magnetic resonance imaging (MRI) where many individual images can be used to create a single composite image.
  • MRI magnetic resonance imaging
  • each receiver can be used to produce an individual image of the subject such that if there is n receivers there will be n images. The n images can then be processed to produce a single composite image.
  • RE radio frequency
  • Some systems also incorporate measurement and use of the noise variances of each coil.
  • Each of the n individual channel gains can then be adjusted after acquisition to produce equal noise variance individual images.
  • the SOS algorithm can then be applied. This additional procedure tends to improve the signal-to-noise ratio (SNR) of the process but may still fail to optimize the SNR of the resultant composite image. This results in an equation:
  • the subject invention pertains to a method and apparatus for improved processing of electrical signals.
  • a specific embodiment of the subject invention can be used with MRI devices.
  • the subject method and apparatus can be used to combine a plurality of individual images into a single composite image.
  • the composite image can have reduced distortion and/or increased signal to noise ratio.
  • the subject method and apparatus can be installed as an aftermarket addition to existing MRI devices in order to take advantages of the method herein described.
  • the subject invention can be incorporated into new MRI devices and/or systems.
  • Figure 1 shows a schematic representation of a specific embodiment of the subject invention.
  • Figure 2 shows a schematic representation of a specific embodiment of the subject invention.
  • Figure 3 shows a clinical image superimposed with percent gains of the optimal SNR with respect to the embodiments shown in Figure 2 pver standard SOS SNR.
  • Figure 4A shows a flowchart for standard SOS processing.
  • Figure 4B shows a flowchart for a specific embodiment of the subject invention.
  • Figure 4C shows a flowchart for a specific embodiment of the subject invention.
  • Figure 4D shows a flowchart for a specific embodiment of the subject invention.
  • the subject invention pertains to a method and apparatus for enhanced multiple coil imaging.
  • the subject invention is advantageous for use in imaging devices, such as MRIs where multiple images can be combined to form a single composite image.
  • the subject method and apparatus utilize a novel process of converting from the original signal vector in the time domain to allow the subject invention to be installed in-line with current MRI devices.
  • An MRI device which could be used with the subject invention, typically has multiple RF coils and receivers where each coil can produce an image of the subject.
  • a preferred first step in the use of a multi -receiver MRI system is to produce n data sets which could be used to produce n images, one from each of the n receivers employed by the device.
  • a next step can be to produce a single composite image using the n individual images.
  • This composite imaging is commonly formed through the manipulation of the n individual images in a sum-of-squares (SOS) algorithm.
  • SOS sum-of-squares
  • the value of each pixel in the composite image is the sum of the squares of the corresponding values of the pixels from each of the n individual images.
  • the sum of the value of each pixel in the composite image is the sum of the magnitude squared of the corresponding pixels from each of the individual images.
  • Some systems also incorporate measurement and use of the noise variances of each coil .
  • Each of the n individual channel gains can then be adjusted after acquisition to produce equal noise variance individual images.
  • the SOS algorithm can then be applied. This additional procedure tends to improve the signal-to- noise ratio (SNR) of the process but may still fail to optimize the SNR of the resultant composite image. This results in an equation:
  • the SOS algorithm is optimal if noise covariance matrix is the identity matrix.
  • the noise covariance matrix can be referred to as N, a nxn Hermitian symmetric matrix.
  • N a priori knowledge (i.e. no information about the coil receive profiles).
  • Optimal then considers only the knowledge of the noise covariance matrix and assumes that the relative pixel value at each location is approximately the relative reception ability of each coil. It can be shown that the SOS algorithm is optimal if the noise covariance matrix is the identity matrix. This method can be summarized by the following simple equation:
  • the process of converting from the original signal vector to a new signal vector can be a linear process corresponding to, for example, multiplication by constant values and addition of vector elements.
  • the process can occur before or after Fourier Transformation. That is, the process can occur in the time domain or in the image domain.
  • the subject method and apparatus can operate in the time domain to convert an original signal vector, which is typically in a correlated noise basis, to a new signal vector, which can be in an uncorrelated noise basis.
  • standard reconstruction algorithms for producing composite images can be used. These algorithms are already in place on many multi-channel MRI systems.
  • the subject method and apparatus can incorporate the processing of each signal with the coil so that the processing does not show to a user, such as to create the effect of "pseudo coils" having 0 noise covariance.
  • the subject invention can be implemented via software and/or hardware.
  • Figure 1 shows a schematic representation of a specific embodiment of the subject invention.
  • This embodiment can be placed after preamplif ⁇ cation and can also be placed within a specific coil package.
  • the device parameters can be gathered from average noise covariance data obtained, for example, from tests on typical subjects.
  • the noise parameters can be obtained internal to the device using noise correlators whose outputs could be used to set gains or attenuations to produce the desired result.
  • Alternative embodiments can be placed before preamplifications.
  • a coil array can have some degree of symmetry. With respect to coil arrays having some degree of symmetry, there can exist eigenvalues with degenerate eigenvectors. Generally, two eigenvectors with the same (or similar) eigenvalue can represent similar effective imaging profiles, often having some shift or rotation in the patterns. Accordingly, channels with the same eigenvalues can be added with some particular phase, resulting in little loss in SNR and partially parallel acquisition (PPA) capability. This can allow for maximizing the info ⁇ nation content per channel.
  • PPA partially parallel acquisition
  • the number of channels may be reduced from n to m, where m ⁇ n, channels with little loss in performance.
  • such a reduction in the number of channels can be implemented with respect to quadrature volume coil channels. For example, if two volume coils or modes have fields which are uniform and perpendicular, then the circularly polarized addition of these fields, for example by adding the fields 90 degrees out of phase, can allow all of the signal into one channel and no signal into the other channel.
  • FIGS. 4A-4B a comparison of the standard method for producing a single composite image from the output signals of N coils to an embodiment of the subject invention which incorporates a decorrelator into which N signals are inputted and M signals are outputted, where M ⁇ N, such that a single composite image is produced from the M outputted signals after processing.
  • Figures 4C and 4D show flowcharts for other embodiments of the subject invention.
  • Figure 4 A shows a flowchart for the standard method of producing a single composite image from the output signals of N coils in a MRI system.
  • FIG. 4B shows a flowchart for a specific embodiment of the subject invention similar to the standard method of Figure 4A with the N signals outputted from the preamplifier inputted to a decorrelator which outputs M signals, where M ⁇ N. Accordingly, the remaining processing shown in Figure 4B is for M signals, rather than N.
  • Figure 4C shows a flowchart for another specific embodiment of the subject invention which can implement the "optimal" reconstruction utilizing a software implementation.
  • the software implementing the optimal reconstruction is based on the previously measured noise covariance matrix U, such that processing of the N signals outputted after processing by the N image domain matrices produces the output image.
  • This approach would be a more efficient algorithm than prior art would provide.
  • This example describes a 4-channel hardware RF combiner network which realizes a basis change to reduced noise correlation.
  • the RF combiner network can realize a basis change to minimal noise correlation.
  • Such a basis change to minimal noise correlation can optimize SNR achievable using the Sum-of-Squares reconstruction method.
  • the network can be realized using passive RF components.
  • a specific embodiment was tested specifically using a 4-channel head coil on a 1.5T scanner. SNR gains of over 30% can be achieved in the periphery, with corresponding loss in uniformity. Noise correlation can be greatly reduced when using the combiner. Results from the measurements agree closely with software techniques.
  • a head coil design has demonstrated greater than 20% peripheral SNR gain by using "optimal" reconstruction offline from raw data.
  • "Optimal" can be defined as reconstructing with the highest SNR using the signal plus noise as an estimate of the true signal, but with the full noise covariance matrix taken into account.
  • the standard Sum- of-Squares reconstruction can be optimized for SNR by employing the noise correlation matrix. However, there are practicality issues and can be a lack of any substantial gain in image quality (i.e. 10% or less).
  • the Sum of Squares (with normalized variances) can be optimal of the noise covariance matrix is diagonal (no correlation between channels).
  • a hardware combiner circuit has been developed to achieve this SNR gain for a specific coil on a typical MRI scanner without software modification.
  • N K ' - ' . 5 . 5 1 . 5 . 7 0 1 . 7
  • the combiner circuit was easily inserted and removed without disturbing the coil to enable accurate comparisons.
  • SENSE and SMASH reconstruction algorithms use the noise correlation routinely, so the use of a hardware combiner solely to achieve better SNR is limited.

Abstract

The subject invention pertains to a method and apparatus for enhanced multiple coil imaging. The subject invention is advantageous for use in imaging devices, such as MRIs where multiple images can be combined to form a single composite image. In one specific embodiment, the subject method and apparatus utilize a novel process of converting from the original signal vector in the time domain to allow the subject invention to be installed in-line with current MRI devices.

Description

DESCRIPTION
METHOD AND APPARATUS FOR ENHANCED MULTIPLE COIL IMAGING
The subject invention was made with government support under a research project supported by the National Institutes of Health (NIH), Department of Health and Human Services, Grant Number 5R44 RR11034-03. The government has certain rights in this invention.
Cross-Reference to Related Application
The present application claims the benefit of U.S. Provisional Application No. 60/299,012, filed June 18, 2001, which is hereby incorporated by reference herein in its entirety, including any, figures, tables, or drawings.
Background of the Invention
The subject invention relates to the field of medical technology and can provide for an improved method of creating composite images from a plurality of individual signals. The subject invention is particularly advantageous in the field of magnetic resonance imaging (MRI) where many individual images can be used to create a single composite image.
In the early stages of MRI development, typical MRI systems utilized a single receiver channel and radio frequency (RE) coil. In order to improve performance, multi- coil systems employing multiple RF coils and receivers can now be utilized. During operation of these multi-receiver systems, each receiver can be used to produce an individual image of the subject such that if there is n receivers there will be n images. The n images can then be processed to produce a single composite image.
Many current systems incorporate a sum-of-squares (SOS) algorithm, where the value of each pixel in the composite image is the square-root of the sum of the squares of the corresponding values of the pixels from each of the n individual images. Where the pixel values are complex, the value of each pixel in the composite image is the square- root of the sum of the magnitude squared of the corresponding pixels from each of the individual images. In mathematical terms, if n coils produce n signals s=(sι, s2, s3, ..., sn) corresponding to the pixel values from a given location, the composite signal pixel is given by the following equation:
Figure imgf000003_0001
Some systems also incorporate measurement and use of the noise variances of each coil. Each of the n individual channel gains can then be adjusted after acquisition to produce equal noise variance individual images. Following this procedure, the SOS algorithm can then be applied. This additional procedure tends to improve the signal-to-noise ratio (SNR) of the process but may still fail to optimize the SNR of the resultant composite image. This results in an equation:
S' -[Diag(N)Yl - S
It can be shown that the SOS algorithm is optimal if the noise covariance matrix is the identity matrix. In order to further optimize the SNR of the resultant composite image, it would be helpful to have knowledge of the noise covariance matrix. Optimal SNR reconstruction in the presence of noise covariance can be summarized by the following simple equation:
U.S. Patent Nos.4,885,541 and 4,946,121 discuss algorithms relating to equations which are similar in form. Typically this method is applied in the image domain, after acquisition and Fourier transformation into separate images.
Brief Summary of the Invention
The subject invention pertains to a method and apparatus for improved processing of electrical signals. A specific embodiment of the subject invention can be used with MRI devices. The subject method and apparatus can be used to combine a plurality of individual images into a single composite image. The composite image can have reduced distortion and/or increased signal to noise ratio. In one embodiment, the subject method and apparatus can be installed as an aftermarket addition to existing MRI devices in order to take advantages of the method herein described. In another embodiment, the subject invention can be incorporated into new MRI devices and/or systems.
Brief Description of the Drawings Figure 1 shows a schematic representation of a specific embodiment of the subject invention.
Figure 2 shows a schematic representation of a specific embodiment of the subject invention.
Figure 3 shows a clinical image superimposed with percent gains of the optimal SNR with respect to the embodiments shown in Figure 2 pver standard SOS SNR. Figure 4A shows a flowchart for standard SOS processing. Figure 4B shows a flowchart for a specific embodiment of the subject invention. Figure 4C shows a flowchart for a specific embodiment of the subject invention.
Figure 4D shows a flowchart for a specific embodiment of the subject invention.
Detailed Description of the Invention The subject invention pertains to a method and apparatus for enhanced multiple coil imaging. The subject invention is advantageous for use in imaging devices, such as MRIs where multiple images can be combined to form a single composite image. In one specific embodiment, the subject method and apparatus utilize a novel process of converting from the original signal vector in the time domain to allow the subject invention to be installed in-line with current MRI devices. An MRI device, which could be used with the subject invention, typically has multiple RF coils and receivers where each coil can produce an image of the subject. Thus, a preferred first step in the use of a multi -receiver MRI system is to produce n data sets which could be used to produce n images, one from each of the n receivers employed by the device. In typical MRI devices, a next step can be to produce a single composite image using the n individual images. This composite imaging is commonly formed through the manipulation of the n individual images in a sum-of-squares (SOS) algorithm. In this process, the value of each pixel in the composite image is the sum of the squares of the corresponding values of the pixels from each of the n individual images. Where the pixel values are complex, the sum of the value of each pixel in the composite image is the sum of the magnitude squared of the corresponding pixels from each of the individual images. In mathematical terms, if n coils produce n signals s=(sι, s2, S3, ..., sn) corresponding to the pixel values from a given location, the composite signal pixel is given by the following equation:
Figure imgf000005_0001
Some systems also incorporate measurement and use of the noise variances of each coil . Each of the n individual channel gains can then be adjusted after acquisition to produce equal noise variance individual images. Following this procedure, the SOS algorithm can then be applied. This additional procedure tends to improve the signal-to- noise ratio (SNR) of the process but may still fail to optimize the SNR of the resultant composite image. This results in an equation:
Figure imgf000005_0002
It can be shown that the SOS algorithm is optimal if noise covariance matrix is the identity matrix. In order to further optimize the SNR of the resultant composite image, it would be helpful to have knowledge of the relative receiver profiles of all the coils and knowledge of the noise covariance matrix. For illustrative purposes, the noise covariance matrix can be referred to as N, a nxn Hermitian symmetric matrix. One can also consider optimal given no a priori knowledge (i.e. no information about the coil receive profiles). Optimal then considers only the knowledge of the noise covariance matrix and assumes that the relative pixel value at each location is approximately the relative reception ability of each coil. It can be shown that the SOS algorithm is optimal if the noise covariance matrix is the identity matrix. This method can be summarized by the following simple equation:
Figure imgf000005_0003
U.S. PatentNos.4,885,541 and 4,946,121, which are hereby incorporated herein by reference, discuss algorithms relating to equations which are similar in form. Typically this method is applied in the image domain, after acquisition and Fourier transformation into separate images. In order to further optimize the SNR of the resultant composite image, it would be helpful to have knowledge of the relative receiver profiles of all the coils and knowledge of the noise covariance matrix. As the matrix N"1 is Hermitian symmetric, N"1 can therefore be expressed in an alternate form as N = KύK.
This allows the equation above to be rewritten as:
-Js' - [K'KY1 - s = >/V - K^K'" s = ^(K''ts)' - K'
This is now in the fonn:
lsl-s
Therefore, it is as if the conventional SOS algorithm is perfonned on a new vector, S . Viewed another way, this optimal equation is equivalent to a Sum of Squares operation after a basis transformation to an uncorrelated basis. The new basis can be considered channels that are noise eigenmodes of the original coil. In a specific embodiment of the subject invention, the process of converting from the original signal vector to a new signal vector can be a linear process corresponding to, for example, multiplication by constant values and addition of vector elements. When only multiplication by constant values and addition of vector elements are used, the process can occur before or after Fourier Transformation. That is, the process can occur in the time domain or in the image domain.
Accordingly, the subject method and apparatus can operate in the time domain to convert an original signal vector, which is typically in a correlated noise basis, to a new signal vector, which can be in an uncorrelated noise basis. In a specific embodiment of the subject invention, standard reconstruction algorithms for producing composite images can be used. These algorithms are already in place on many multi-channel MRI systems. In addition, it may be possible to apply other algorithms for spatial encoding such as SMASH and SENSE with better efficiency.
In a specific embodiment, the subject method and apparatus can incorporate the processing of each signal with the coil so that the processing does not show to a user, such as to create the effect of "pseudo coils" having 0 noise covariance. The subject invention can be implemented via software and/or hardware.
Figure 1 shows a schematic representation of a specific embodiment of the subject invention. This embodiment can be placed after preamplifϊcation and can also be placed within a specific coil package. The device parameters can be gathered from average noise covariance data obtained, for example, from tests on typical subjects. Alternatively, the noise parameters can be obtained internal to the device using noise correlators whose outputs could be used to set gains or attenuations to produce the desired result. Alternative embodiments can be placed before preamplifications.
The subject method can allow for a reduction in the number of channels and, in a specific embodiment, be generalized to an optimal reduction in the number of channels. In some embodiments, a coil array can have some degree of symmetry. With respect to coil arrays having some degree of symmetry, there can exist eigenvalues with degenerate eigenvectors. Generally, two eigenvectors with the same (or similar) eigenvalue can represent similar effective imaging profiles, often having some shift or rotation in the patterns. Accordingly, channels with the same eigenvalues can be added with some particular phase, resulting in little loss in SNR and partially parallel acquisition (PPA) capability. This can allow for maximizing the infoπnation content per channel. In a specific embodiment, adding the signals from channels with identical (or similar) eigenvalues, the number of channels may be reduced from n to m, where m < n, channels with little loss in performance. In a specific embodiment, such a reduction in the number of channels can be implemented with respect to quadrature volume coil channels. For example, if two volume coils or modes have fields which are uniform and perpendicular, then the circularly polarized addition of these fields, for example by adding the fields 90 degrees out of phase, can allow all of the signal into one channel and no signal into the other channel. In this example, the use of a phased addition associated with the uniform modes both having the same eigenvalves and being related through a rotation of 90 degrees allows the conversion of two coils to one channel to be accomplished. Referring to Figures 4A-4B, a comparison of the standard method for producing a single composite image from the output signals of N coils to an embodiment of the subject invention which incorporates a decorrelator into which N signals are inputted and M signals are outputted, where M < N, such that a single composite image is produced from the M outputted signals after processing. Figures 4C and 4D show flowcharts for other embodiments of the subject invention. Figure 4 A shows a flowchart for the standard method of producing a single composite image from the output signals of N coils in a MRI system. The N signals are preamplifϊed and then mixed to lower frequencies. The N lower frequency signals are then sampled by A/D converters to produce N digital signals. 2D Fourier Transforms are applied to the N digital signals followed by multiplication by N image domain matrices to produce a single composite image. Figure 4B shows a flowchart for a specific embodiment of the subject invention similar to the standard method of Figure 4A with the N signals outputted from the preamplifier inputted to a decorrelator which outputs M signals, where M < N. Accordingly, the remaining processing shown in Figure 4B is for M signals, rather than N. Figure 4C shows a flowchart for another specific embodiment of the subject invention which can implement the "optimal" reconstruction utilizing a software implementation. In this embodiment, the software implementing the optimal reconstruction is based on the previously measured noise covariance matrix U, such that processing of the N signals outputted after processing by the N image domain matrices produces the output image. Figure 4D shows the flowchart for another specific embodiment in accordance with the subject invention where, after preamplification and optionl additional amplification, the N signals from the N coils undergo direct digitization and are inputted to a DSP for 2D Fourier Transform and matrix operations using the K matrix, where N is the noise covariance matrix and N=KlK to produce an output image. This approach would be a more efficient algorithm than prior art would provide.
Example 1
This example describes a 4-channel hardware RF combiner network which realizes a basis change to reduced noise correlation. The RF combiner network can realize a basis change to minimal noise correlation. Such a basis change to minimal noise correlation can optimize SNR achievable using the Sum-of-Squares reconstruction method. The network can be realized using passive RF components. A specific embodiment was tested specifically using a 4-channel head coil on a 1.5T scanner. SNR gains of over 30% can be achieved in the periphery, with corresponding loss in uniformity. Noise correlation can be greatly reduced when using the combiner. Results from the measurements agree closely with software techniques.
A head coil design has demonstrated greater than 20% peripheral SNR gain by using "optimal" reconstruction offline from raw data. "Optimal" can be defined as reconstructing with the highest SNR using the signal plus noise as an estimate of the true signal, but with the full noise covariance matrix taken into account. The standard Sum- of-Squares reconstruction can be optimized for SNR by employing the noise correlation matrix. However, there are practicality issues and can be a lack of any substantial gain in image quality (i.e. 10% or less). The Sum of Squares (with normalized variances) can be optimal of the noise covariance matrix is diagonal (no correlation between channels). To achieve this gain in clinical practice, a hardware combiner circuit has been developed to achieve this SNR gain for a specific coil on a typical MRI scanner without software modification.
The standard Sum-of-Squares (SoS) reconstruction from an n-channel array for each pixel is ύ sos = S S where S = tfj • • ■ Sll ) . The optimal reconstruction using only signal data is Sop, = s'N" s where N is the noise correlation matrix. To yield this optimal result from a standard SoS operation, a signal basis change can be employed. Since in general N' = N, there exists a K constructed via eigenvalue/vector decomposition such that N = K'K and so Sopt — s'K K S . Thus we can define
$ = K S as our new signal vector and ^opt = S S which is just a SoS operation.
In hardware, multiplication of the 4 element signal vector s by a matrix is equivalent to using a 4-in-to-4-out network which in general is complicated by phase shifts and gain scaling, and must almost certainly be placed after preamplification. The 4-channel head coil used in this example has 4 sets of opposing parallel-combined loops and was described in King, S.B. et al. 9th ISMRM proceedings, p. 1090, 2001, which is incorporated herein by reference. Its symmetry yields an all-real noise correlation matrix and corresponding signal basis change matrix of the form shown below: 1 . 5 .5 . 1 - . 7 1 0 . 7
. 5 1 . 1 - .5 . 7 1 0 - . 7
N = K ' - ' . 5 . 5 1 . 5 . 7 0 1 . 7
. 5 - . 5 . 5 1 .7 0 - 1 .7
The basis change matrix diagonalizes N, so in the new basis there is no noise correlation, which maximizes SNR. The ± 0.7 entries in K are just 3dB attenuations, resulting in the very simple circuit realization of Figure 2. Each 4-way splitter output is 3dB lower than each 2-way output so no additional attenuations are needed. The -1 polarity inversions are realized using baluns. All components are standard RF parts which utilize ferrite cores, so the circuit was kept far out of the bore. Images were made using a 1.5T scanner. Figure 3 shows a clinical image superimposed with percent gains of the hardware-optimal SNR over regular SoS SNR.
The combiner circuit was easily inserted and removed without disturbing the coil to enable accurate comparisons.
Clearly the SNR gain is mainly in the periphery and the resulting image suffers in uniformity between the center and the peripheiy. Hardware results are also in close agreement with software-optimal reconstruction using normal coil data. Noise correlation is significantly decreased when using the hardware-optimal combiner, and the matrix is nearly diagonal as expected:
Noise Correlation Matrices
Normal Coil w Hardware Combiner
1.00 0.50 0.40 0.11 1.00 0.05 0.12 0.05
0.50 1.00 0.10 0.57 0.05 1.00 0.07 0.15
0.40 0.10 1.00 0.48 0.12 0.07 1.00 0.10
0.11 0.57 0.48 1.00 0.05 0.15 0.10 1.00
The results of this example illustrate that a hardware De-correlator can be used to approximate "optimal" reconstruction using software sums of squares algorithm. As the symmetry of noise covariance entries is more important than particular values, this device can show robustness across different loading conditions.
SENSE and SMASH reconstruction algorithms use the noise correlation routinely, so the use of a hardware combiner solely to achieve better SNR is limited.
However, this and similar combiners may find use in reducing computational time and improving numerical stability in the various new reconstruction algorithms being developed.
It should be understood that the examples and embodiments described herein are for illustrative purposes only and that various modifications or changes in light thereof will be suggested to persons skilled in the art and are to be included within the spirit and purview of this application and the scope of the appended claims.

Claims

ClaimsWe claim:
1. A method of processing magnetic resonance imaging signals from a plurality of magnetic resonance imaging coils, comprising: deteπnining a noise covariance matrix, N, of a plurality of magnetic resonance imaging coils; receiving a corresponding plurality of signals,
Figure imgf000012_0001
s2, s3, ..., sn) from the plurality of coils, wherein the plurality of signals represent a corresponding plurality of pixel values for a location;
calculating a composite pixel value for the location,
Figure imgf000012_0002
where S' =
(K'"s)' , S = (K'"s) andN = K'K .
2. The method according to claim 1, wherein the noise covariance matrix, N, of the plurality of magnetic resonance imaging coils is a Hermitian symmetric matrix.
3. The method according to claim 1, wherein (s = §ι, s2, . . . ,sn) is produced by inputting s=(sι, s2, ..., sn) into a circuit, wherein the output of the circuit is
§ι = a\, si, + bi s2 + . . . + W] sn §2 = a2, si, + b2 s2 + . . . + w2 sn sn = an, sι, + bn s2 + . . . + wn sn wherein ai, a2, . . . , an, bj, b2, . . .,bn , Wi, w2, . . ., wn are constants.
4. The method according to claim 3, wherein ai, a2, . . . , an, bla b2, . . .bn , wi, w2, . . ., wn are values of K' such that ai, a2, . . . , an, bi, b2, . . .bn , wj, w2, . . ., wn are equal to
Figure imgf000012_0003
, respectively.
5. The method according to claim 4, wherein K is constructed via eigenvalue/vector decomposition of N.
6. The method according to claim 5, wherein when the eigenvalue/vector decomposition of N yields one or more eigenvalues with substantially degenerate eigenvectors, the values of K' are adjusted such that the output of the circuit is s = (§ι, s2, . . . ,sm) where m < n.
7. The method according to claim 6, wherein substantially all of the image information in received signals s=(sι, s2, ..., sn) is in output signals s = (§ι, s2, . . .
8. The method according to claim 6, wherein the adjustment of the values of
K' results in the received signals with substantially the same eigenvalues being added together with a phase.
9. The method according to claim 8, wherein the received signals are received from quadrature volume coils having fields which are substantially uniform and substantially peipendicular, wherein the values of K' are adjusted to accomplish the circularly polarized addition of the two received signals, such that substantially all of the image information in the two received signals is in one output signal.
10. The method according to claim 4, further comprising: preamplifying the received signals prior to inputting s=(sls s2, ..., sn) into the circuit.
11. The method according to claim 10, further comprising: mixing output signal S to lower frequencies; sampling S with a lower frequency by A/D converters to produce a digital S signal; applying a 2D Fourier Transform to the digital S signal; and processing the S signal after 2D Fourier Transform applied with image domain matrices to produce a plurality of pixel values, S, for a location, wherein a composite pixel value for the location, S' • S , utilizes S after processing with image domain matrices.
PCT/US2002/019275 2001-06-18 2002-06-18 Method and apparatus for enhanced multiple coil mr imaging WO2003003037A2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
EP02756225A EP1412769A2 (en) 2001-06-18 2002-06-18 Method and apparatus for enhanced multiple coil mr imaging
AU2002322248A AU2002322248A1 (en) 2001-06-18 2002-06-18 Method and apparatus for enhanced multiple coil mr imaging

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US29901201P 2001-06-18 2001-06-18
US60/299,012 2001-06-18

Publications (2)

Publication Number Publication Date
WO2003003037A2 true WO2003003037A2 (en) 2003-01-09
WO2003003037A3 WO2003003037A3 (en) 2003-02-27

Family

ID=23152950

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2002/019275 WO2003003037A2 (en) 2001-06-18 2002-06-18 Method and apparatus for enhanced multiple coil mr imaging

Country Status (4)

Country Link
US (2) US7057387B2 (en)
EP (1) EP1412769A2 (en)
AU (1) AU2002322248A1 (en)
WO (1) WO2003003037A2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006122013A2 (en) * 2005-05-06 2006-11-16 Invivo Corporation Method and apparatus for adaptive channel reduction for parallel imaging
WO2008010221A1 (en) * 2006-07-17 2008-01-24 Mri Microsystems Ltd. Method and system for signal processing

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7057387B2 (en) * 2001-06-18 2006-06-06 Invivo Corporation Method and apparatus for enhanced multiple coil imaging
WO2005012931A1 (en) * 2003-07-25 2005-02-10 National Research Council Of Canada Stacked coil array for magnetic resonance experiments
WO2005047914A1 (en) * 2003-11-12 2005-05-26 Philips Intellectual Property & Standards Gmbh Parallel mr imaging method
US7279893B1 (en) * 2006-04-20 2007-10-09 General Electric Company Receiver channel data combining in parallel mr imaging
US20100013472A1 (en) * 2006-09-13 2010-01-21 Martin Buehrer Method and apparatus for coil array compression
WO2008134891A1 (en) * 2007-05-03 2008-11-13 National Research Counsil Of Canada Rf based spatially selective excitation in mri
US8076939B2 (en) * 2008-04-10 2011-12-13 The General Hospital Corporation Method for fast magnetic resonance radiofrequency coil transmission profile mapping
CN102257405B (en) * 2008-11-12 2015-11-25 拜耳医疗保健公司 Orthogonal Endo-rectal coil and the interfacing equipment for this orthogonal Endo-rectal coil
US8581589B2 (en) * 2009-02-02 2013-11-12 The General Hospital Corporation System and method for variable mode-mixing in magnetic resonance imaging
US8970217B1 (en) 2010-04-14 2015-03-03 Hypres, Inc. System and method for noise reduction in magnetic resonance imaging
CA2803252C (en) 2010-07-01 2018-10-30 Medrad, Inc. Multi-channel endorectal coils and interface devices therefor
DE102010041212B4 (en) * 2010-09-22 2012-07-26 Siemens Aktiengesellschaft Compensation of echo-time-independent phase or magnitude components in recorded MR image data
MX2014002574A (en) * 2011-09-07 2014-06-05 Koninkl Philips Nv Mode splitter/combiner for noise figure minimization and control of power in rf coil arrays.
EP3698157A1 (en) * 2017-10-18 2020-08-26 Institut du Cerveau et de la Moelle Epiniere-ICM Optimised generation of mri images by a multi-antenna mri system
DE102020208611A1 (en) * 2020-07-09 2022-01-13 Siemens Healthcare Gmbh Device and method for data compression in local coils

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001073463A1 (en) * 2000-03-24 2001-10-04 Koninklijke Philips Electronics N.V. Magnetic resonance imaging method with sub-sampling

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4885541A (en) 1988-08-19 1989-12-05 General Electric Company Apparatus and method for enhanced multiple coil nuclear magnetic resonance (NMR) imaging
US4947121A (en) 1988-08-19 1990-08-07 General Electric Company Apparatus and method for enhanced multiple coil nuclear magnetic resonance (NMR) imaging
DE69527122T2 (en) 1994-04-08 2003-01-09 Koninkl Philips Electronics Nv METHOD AND DEVICE FOR GENERATING IMAGES BY MEANS OF MAGNETIC RESONANCE
EP0803738B1 (en) * 1996-04-24 2003-04-02 Philips Corporate Intellectual Property GmbH Method of synthesizing an image for producing a combination image from input images
US6898582B2 (en) * 1998-12-30 2005-05-24 Algodyne, Ltd. Method and apparatus for extracting low SNR transient signals from noise
US6160398A (en) 1999-07-02 2000-12-12 Vista Clara, Inc. Adaptive reconstruction of phased array NMR imagery
EP1307759A1 (en) * 2000-07-31 2003-05-07 Koninklijke Philips Electronics N.V. Magnetic resonance method for forming a fast dynamic imaging
US6771067B2 (en) * 2001-04-03 2004-08-03 The United States Of America As Represented By The Department Of Health And Human Services Ghost artifact cancellation using phased array processing
US6486671B1 (en) * 2001-05-14 2002-11-26 Ge Medical Systems Global Technologies Company Llc MRI image quality improvement using matrix regularization
US6675106B1 (en) * 2001-06-01 2004-01-06 Sandia Corporation Method of multivariate spectral analysis
US6584413B1 (en) * 2001-06-01 2003-06-24 Sandia Corporation Apparatus and system for multivariate spectral analysis
US7057387B2 (en) * 2001-06-18 2006-06-06 Invivo Corporation Method and apparatus for enhanced multiple coil imaging

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001073463A1 (en) * 2000-03-24 2001-10-04 Koninklijke Philips Electronics N.V. Magnetic resonance imaging method with sub-sampling

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
HAYES C E ET AL: "NOISE CORRELATIONS IN DATA SIMULTANEOUSLY ACQUIRED FROM MULTIPLE SURFACE COIL ARRAYS" MAGNETIC RESONANCE IN MEDICINE, ACADEMIC PRESS, DULUTH, MN, US, vol. 16, no. 2, 1 November 1990 (1990-11-01), pages 181-191, XP000175902 ISSN: 0740-3194 *
K.P.PRUESSMANN ET AL.: "Advances in Sensitivity Encoding with Arbitrary k-Space Trajectories" PROC. OF THE INTERNAT. SOCIETY FOR MAGN. RESON. IN MEDICINE, 9TH SCIENTIFIC MEETING AND EXHIBITION, EUROPEAN SOCIETY FOR MAGNETIC RESONANCE IN MEDICINE AND BIOLOGY, 18TH ANNUAL MEETING AND EXHIBITION, HELD JOINTLY, 21-27 APRIL 2001, GLASGOW, UK, vol. 1, page 4 XP002215281 *
ROEMER P B ET AL: "THE NMR PHASED ARRAY" MAGNETIC RESONANCE IN MEDICINE, ACADEMIC PRESS, DULUTH, MN, US, vol. 16, no. 2, 1 November 1990 (1990-11-01), pages 192-225, XP000175903 ISSN: 0740-3194 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006122013A2 (en) * 2005-05-06 2006-11-16 Invivo Corporation Method and apparatus for adaptive channel reduction for parallel imaging
WO2006122013A3 (en) * 2005-05-06 2006-12-28 Invivo Corp Method and apparatus for adaptive channel reduction for parallel imaging
US7576536B2 (en) 2005-05-06 2009-08-18 Invivo Corporation MRI method and apparatus for adaptive channel reduction in parallel imaging
WO2008010221A1 (en) * 2006-07-17 2008-01-24 Mri Microsystems Ltd. Method and system for signal processing

Also Published As

Publication number Publication date
AU2002322248A1 (en) 2003-03-03
US20060055405A1 (en) 2006-03-16
US7057387B2 (en) 2006-06-06
WO2003003037A3 (en) 2003-02-27
US20030038632A1 (en) 2003-02-27
EP1412769A2 (en) 2004-04-28
US7141971B2 (en) 2006-11-28

Similar Documents

Publication Publication Date Title
US7141971B2 (en) Method and apparatus for enhanced multiple coil imaging
US6377044B1 (en) Multi-mode receiver coils for MRI
US20080265890A1 (en) Simultaneous Mr-Excitation of Multiple Nuclei with a Single Rf Amplifier
US5389880A (en) Magnetic resonance imaging apparatus
EP0344293B1 (en) Nuclear magnetic resonance (nmr) imaging with multiple surface coils
Wright et al. Theory and application of array coils in MR spectroscopy
JP3825685B2 (en) Magnetic resonance imaging equipment using high frequency coils
US5578925A (en) Vertical field quadrature phased array coil system
EP0845683A1 (en) Multi-mode RF coil
US20050062472A1 (en) Mri tunable antenna and system
US4885541A (en) Apparatus and method for enhanced multiple coil nuclear magnetic resonance (NMR) imaging
King et al. Optimum SNR data compression in hardware using an Eigencoil array
US6914432B2 (en) Phased array coil assembly and method and system for employing the same
JP5641331B2 (en) RF coil array having coil elements with two preamplifiers
US20170234948A1 (en) Method and apparatus for hardware rf receiver channel reduction
US7235973B2 (en) Phased array coils utilizing selectable quadrature combination
US20120019250A1 (en) Noise matching in couplet antenna arrays
He et al. Four-channel magnetic resonance imaging receiver using frequency domain multiplexing
US10101418B2 (en) Flexible signal multiplexer for MRI receiving systems
JPH0542125A (en) Nuclear magnetic resonance device
Varosi et al. A hardware combiner to achieve ‘optimal’SNR using sum-of-squares reconstruction
EP0266398A1 (en) Reduced noise nmr localization system
Minakov et al. The Mode Switching Matrix of the Receiving Channel with a Different Type of Polarization for Advanced Magnetic Resonance Imaging for Early Cancer Diagnosis
KR20220115219A (en) High frequency coil for magnetic resonance imaging system
Bankson et al. A frequency independent time-domain-multiplexed receiver for MR imaging and spectroscopy

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A2

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NO NZ OM PH PL PT RO RU SD SE SG SI SK SL TJ TM TN TR TT TZ UA UG US UZ VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A2

Designated state(s): GH GM KE LS MW MZ SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

AK Designated states

Kind code of ref document: A3

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NO NZ OM PH PL PT RO RU SD SE SG SI SK SL TJ TM TN TR TT TZ UA UG US UZ VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A3

Designated state(s): GH GM KE LS MW MZ SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
WWE Wipo information: entry into national phase

Ref document number: 2002756225

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 2002756225

Country of ref document: EP

REG Reference to national code

Ref country code: DE

Ref legal event code: 8642

NENP Non-entry into the national phase

Ref country code: JP

WWW Wipo information: withdrawn in national office

Country of ref document: JP