US4651098A - Method for imaging nuclear magnetic resonance signals by using non-linear magnetic field gradient - Google Patents

Method for imaging nuclear magnetic resonance signals by using non-linear magnetic field gradient Download PDF

Info

Publication number
US4651098A
US4651098A US06/656,932 US65693284A US4651098A US 4651098 A US4651098 A US 4651098A US 65693284 A US65693284 A US 65693284A US 4651098 A US4651098 A US 4651098A
Authority
US
United States
Prior art keywords
magnetic field
target
specifying
signals
nmr
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
US06/656,932
Other languages
English (en)
Inventor
Yoshifumi Yamada
Kunio Tanaka
Zenwemon Abe
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Utsunomiya University
Original Assignee
Utsunomiya University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Utsunomiya University filed Critical Utsunomiya University
Assigned to UTSUNOMIYA UNIVERSITY reassignment UTSUNOMIYA UNIVERSITY ASSIGNMENT OF ASSIGNORS INTEREST. Assignors: ABE, ZENWEMON, TANAKA, KUNIO, YAMADA, YOSHIFUMI
Application granted granted Critical
Publication of US4651098A publication Critical patent/US4651098A/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console

Definitions

  • This invention relates to a method for imaging nuclear magnetic resonance (to be referred to as "NMR", hereinafter) signals by using a non-linear magnetic field gradient which signals carry information on distribution of nuclear magnetic substance, such as hydrogen (H), fluorine (F), sodium (Na), carbon (C), phosphorus (P), and the like elements, in a body being measured (to be referred to as a "target”, hereinafter).
  • NMR nuclear magnetic resonance
  • F fluorine
  • Na sodium
  • C carbon
  • P phosphorus
  • an NMR computed tomography (to be referred to as the "NMR-CT", hereinafter), using a linear magnetic field gradient is known. More specifically, a conventional method for NMR imaging is based on the fact that when the NMR signals are detected under the presence of a linear magnetic field gradient, the spectrum of the NMR signals becomes proportionate to a projected image of the spin distribution in the target, the projected image being taken in a direction perpendicular to the direction of the magnetic field gradient.
  • the conventional method forms an image of NMR-CT by applying the same algorithm as that of the x-ray computed tomography (to be referred to as the "x-ray CT", hereinafter) to a plurality of such projected images detected by varying the direction of the linear magnetic field gradient to a number of orientations.
  • Another conventional method for NMR imaging is based on the fact that when free precession of nuclear magnetization is effected under the presence of a linear magnetic field gradient, the phase of the nuclear magnetization varies in proportion to the spatial position in the target, and that quantities proportional to one-dimensional Fourier transformation (FT) of the spin density distribution can be obtained by systematically varying the timing and the value of the magnetic field gradient for effecting the free precession of the nuclear magnetization.
  • FT one-dimensional Fourier transformation
  • the conventional method for imaging NMR signals using a non-linear magnetic field gradient has a shortcoming in that a specifying magnetic field to be used therein is required to have a strong non-linear gradient to specify a small region for measurement and to clearly define the boundary of such small region, and coils with comparatively large ampere-turns are necessary for the generation and scanning manipulation of the specifying magnetic field with such strong non-linear gradient.
  • the requirement for coils with large ampere-turns is a kind of limitation to the formation of a large imaging device, especially a large magnetic field generator.
  • the above-mentioned conventional method for imaging NMR signals by using the non-linear magnetic field gradient has another shortcoming in its slowness in producing two- or three-dimensional images.
  • the specifying magnetic field with the non-linear gradient specifies a small region for each direct measurement, so that each application of the high-frequency (to be referred to as "HF", hereinafter) pulse magnetic field for detection of nuclear magnetization yields measured data from only a small region, and it takes a long time to collect sufficient data for imaging the entire target on a two- or three-dimensional basis.
  • the methods using the non-linear magnetic field may be left behind as far as the speed of collecting the measured data is concerned.
  • an object of the present invention is to obviate the above-mentioned shortcomings of the prior art by providing an improved method for imaging NMR signals by using a non-linear magnetic field gradient without necessitating any excessive increase of the intensity of the specifying magnetic field with a non-linear gradient.
  • the method for imaging NMR signals by using a non-linear magnetic field gradient according to the invention eliminates the need of an excessive increase of ampere-turns for exciting the nuclear magnetization, produces measured data from an entire measuring region for each application of the HF pulse magnetic field for detection of the nuclear magnetization, and ensures efficient imaging in a short period of time with an excellent S/N ratio.
  • Another object of the invention is to provide a method for imaging NMR signals by using a non-linear magnetic field gradient, in which method NMR signals carrying information covering an entire imaging region can be obtained at one time by using a specifying magnetic field with a non-linear gradient including a change of magnetic field intensity along a quadratic curve, and a large number of NMR signals carrying measured data are obtained from each scanning point in response to scanning of the specifying magnetic field by either varying the current in a scanner coil or mechanically moving a generator coil thereof, and the NMR signals thus obtained are processed by a computer so as to produce an image of information of spin density distribution in the imaging region.
  • a further object of the invention is to provide a method for imaging NMR signals by using a non-linear magnetic field gradient, in which method sufficient phase change of the nuclear spin necessary for identifying individual points in an imaging region cn be obtained by scanning over a wide region with a comparatively weak specifying magnetic field with a non-linear gradient, so as to facilitate the use of a weak static magnetic field and to eliminate the need of a strong linear magnetic field gradient in the conventional NMR-CT for improving the resolution without reducing the S/N ratio.
  • a specifying magnetic field with a non-linear gradient of magnetic field intensity is applied to the target for a certain period of time so as to spatially encode the phase of the nuclear magnetization in a certain locality of the target in a non-linear fashion by causing free precession of the nuclear magnetization.
  • the specifying magnetic field has a magnetic field intensity varying in a non-linear manner.
  • the specifying magnetic field is spatially scanned by shifting the center of the non-linear gradient thereof over a center scan distance by a specifying field scanner.
  • the above encoding is successively repeated while effecting the spatial scanning of the specifying magnetic field at each encoding, and the encoded NMR signals from the target are detected by a receiver coil.
  • a plurality of the encoded NMR signals thus detected have different degrees of contribution from the nuclear magnetizations at different localities of the target.
  • Such encoded NMR signals are processed by a computer so as to produce an image of the information carried thereby, i.e., the information relating to the distribution of nuclear magnetic substance in the target.
  • FIG. 1 is a block diagram which schematically shows fundamental formation of a device for imaging NMR signals carrying internal information of a target by using a method according to the present invention
  • FIG. 2 is a schematic diagram showing formation and arrangement of various elements for generating a variety of magnetic fields for producing the NMR signals;
  • the curves (a) and (b) of FIG. 3 show the manner in which pulse-like magnetic fields are applied to a target in a method according to the present invention.
  • the curves (a) and (b) of FIG. 4 show the result of a numerical simulation of the reproduction of nuclear spin density distribution by a method for imaging NMR signals according to the present invention.
  • U1 through U4 are constituent units of a device suitable for carrying out the method according to the present invention
  • A is a target
  • B is an electromagnet for producing a homogeneous static magnetic field
  • C1 is a scanner coil for scanning a specifying magnetic field
  • C2 is a generator coil of the specifying magnetic field
  • D is a transmitter coil of HF magnetic field
  • T is a transmitter
  • R is a receiver
  • 1 is a radio frequency (RF) oscillator
  • 2 is a program pulse generator
  • 3 is a wave-shaping gate
  • 4 is an RF power amplifier
  • 5 is a preamplifier
  • 6 is a main amplifier
  • 7 is a phase-sensitive detector
  • 8 is an A/D (analog-to-digital) converter
  • 9 and 10 are D/A (digital-to-analog) converters
  • 11 and 12 are D.C. power amplifiers
  • 13 is a highly stabilized D.C. power source
  • 14 is a
  • FIG. 1 Overall arrangement of a device suitable for imaging the information on the inside of a target by a method for imaging NMR signals according to the invention is shown in FIG. 1, and the disposition of major constituent elements of the device is shown in FIG. 2.
  • a pair of electromagnets B facing each other with a spacing therebetween are energized by a highly stabilized D.C. power source 13, so as to produce a homogeneous static magnetic field H o .
  • the specifying magnetic field generator coil C2 and the scanner coil C1 of the specifying magnetic field are disposed in the proximity of the electromagnets B and energized by a specific magnetic field scanning circuit U3a, so as to superpose a specifying magnetic field ⁇ H s of suitable configuration onto the static magnetic field H o .
  • various other means may be used for producing the homogeneous static magnetic field H o , such as air-core coils, superconductive magnets or permanent magnets.
  • the coils C1 and C2 for generating and scanning the specifying magnetic field ⁇ H s are disposed along the z-axis, but it is also possible to dispose them along the x-axis or y-axis depending on the manner in which the specifying magnetic field is generated and scanned in the specific application.
  • a target A is placed in the composite magnetic field (H o + ⁇ H s ) thus formed by the superposition of the homogeneous static magnetic field and the specifying magnetic field.
  • the transmitter coil D is disposed around the target A and energized by the transmitter T of the second unit U2, so as to superpose a high-frequency magnetic field onto the composite magnetic field.
  • the frequency of the high-frequency magnetic field thus superposed coincides with the nuclear magnetic resonance frequency of the specific nucleus being measured within the target A, resonant energy absorption from the thus applied high-frequency magnetic field occurs, and the NMR alternating magnetic field in response to such energy absorption is picked up by the receiver coil E.
  • the electromotive force induced in the receiver coil E is applied to the receiver R of the second unit U2 and detected there as NMR signals.
  • the axes of the transmitter coil D and the receiver coil E are disposed perpendicular to the direction of the homogeneous static magnetic field H o produced by the static magnetic field generator electromagnets B.
  • the nuclear magnetic substance in the target A has a Larmor frequency which is determined by the intensity of the static magnetic field at the target A with a one-to-one relationship.
  • the high-frequency magnetic field of the Larmor frequency is applied in the form of a pulse having a suitable duration and a suitable magnitude, a detectable signal with diminishing amplitude is produced by nuclear spin excited in the target A even after removing such pulse-like high-frequency magnetic field, which diminising signal has a Larmor frequency depending on the static magnetic field surrounding the target A when the high-frequency magnetic field is removed therefrom.
  • the nuclear spins in the target A after the removal of the high-frequency magnetic field do free precession at a Larmor frequency which is determined by the surrounding static magnetic field with one-to-one relationship. Accordingly, if a specifying magnetic field ⁇ H s with a non-linear variation of magnetic field intensity is applied to the target A by the specifying magnetic field generator coil C2 for a certain period of time immediately after the removal of the high-frequency magnetic field, the phase of the nuclear spins at a region of the target A is shifted depending on the intensity of the specifying magnetic field ⁇ H s at that region.
  • information which varies depending on the position in the target A can be converted into phase information, and signals in the form of an integral of contributions from such phase information can be obtained.
  • the specifying magnetic field for producing the characterizing nuclear spin phase information from spatial information concerning the inside of a target can assume any of various configurations; namely, a bar-shaped specifying magnetic field having cylindrical contour surfaces of the magnetic field intensity, a star-like specifying magnetic field having saddle-shaped contour surfaces of the magnetic field intensity, or a spherical or ellipsoidal specifying magnetic field having spherical or ellipsoidal contour surfaces of the magnetic field intensity.
  • the change of magnetic field intensity at a point is proportional to the square of the distanace from the central position to that point.
  • the generator coil C2 of the specifying magnetic field in the first unit U1 of the illustrated embodiment can be made of parallel wires, a rectangular coil, or a circular coil, provided that the desired specifying magnetic field is generated thereby.
  • the specifying magnetic field is electrically scanned by the scanner coil C1 of the first unit U1, without affecting the configuration of the specific magnetic field to any significant extent. In practice, such scanning of the specifying magnetic field may be easily effected by shifting the generator coil C2 or moving the target A itself, instead of the scanning by the above scanner coil C1.
  • the generator coil C2 of the specifying magnetic field may be provided with accessories, such as a compensator coil means for providing compensation for systematic magnetic field intensity offset which may accompany the generation of the specifying magnetic field, a linear gradient generator coil means for generating a linear magnetic field in a different spatial direction, or a non-linear gradient generator coil means for generating a planar homogeneous magnetic field region.
  • a compensator coil means for providing compensation for systematic magnetic field intensity offset which may accompany the generation of the specifying magnetic field
  • a linear gradient generator coil means for generating a linear magnetic field in a different spatial direction
  • a non-linear gradient generator coil means for generating a planar homogeneous magnetic field region.
  • the first unit U1 comprising the homogeneous static magnetic field generator magnets B, the transmitter coil D and the receiver coil E is an essential NMR apparatus having those elements which are selectively driven by the combined action of the second unit U2 and the third unit U3 under the control of the forth unit U4.
  • the fourth unit U4 not only effects the overall control of the device of FIG. 1, but also processes the measured data and provides display of the outcome thereof.
  • a high-frequency magnetic field with the certain frequency spectrum is applied to the target A by the transmitter coil D, so that nuclear magnetization is excited either universally throughout the entire regions of the target A to be imaged or selectively on that selected or planar region of the target A, which selected or planar region is formed by superposing a specifying magnetic field onto the target A by the generator coil C2 of the first unit U1, the specifying magnetic field having a linear magnetic field gradient or planar non-linear magnetic field gradient with a planar region of homogeneous magnetic field intensity.
  • the nuclear magnetization thus excited in the target A has a magnetization component in a direction perpendicular to the homogeneous static magnetic field H o , and such magnetization component is maximized when the high-frequency magnetic field pulse for the excitation satisfies the conditions of the so-called 90° pulse.
  • the specifying magnetic field is applied to the target A by the generator coil C2 of the first unit U1.
  • the nuclear magnetization component perpendicular to the homogeneous static magnetic field H o is caused to do free precession at that Larmor frequency which corresponds to the magnetic field intensity of the specifying magnetic field ⁇ H s at the position of such nuclear magnetization component.
  • NMR signals are generated at various points in the target A, and the NMR signals are synthesized at the receiver coil E for producing a comparatively weak composite high-frequency current in the receiver coil E.
  • the high-frequency current is amplified by a preamplifier 5 and a main amplifier 6 of the second unit U2.
  • the amplifier current is applied to a phase-sensitive detector 7 for detection by using the oscillating frequency of the RF oscillator 1 as a reference signal, and the so-called free induction decay (to be referred to as the "FID", hereinafter) signal is obtained.
  • the specifying mangetic field ⁇ H s is applied to the target A for a time period ⁇ .
  • the nuclear magnetizations at different regions of the target A do free precessions at different Larmor frequencies depending on the magnetic field intensity of the specifying magnetic field ⁇ H, so that the amplitude of the FID signal at the end of the time period ⁇ is proportional to the composite or synthesized sum of the non-linear phase movement of individual nuclear magnetization at each region relative to the spatial coordinates, which non-linear phase movement can be in the form of a quadratic function.
  • the application of the specifying magnetic field ⁇ H s is ceased, and immediately thereafter a static magnetic field having a linear magnetic field gradient in a certain direction is either applied to the target A or not applied.
  • the FID signals due to the homogeneous static magnetic field H o or spin echo signals due to the above-mentioned applications of magnetic field in pulse-form are converted into digital signals by the A/D converter 8, and the thus converted digital signals are delivered to the control computer 14 for processing.
  • the above operation is repeated while shifting the specifying magnetic field ⁇ H s by a certain amount at a time by the scanner coil C1 of the specifying magnetic field, and the FID signals or the spin echo signals are successively converted into digital signals and applied to the control computer 14 as measured data signals.
  • the control computer 14 carries out frequency analysis and imaging operation on such measured data signals, and the image thus obtained is shown on the peripheral display 15.
  • control computer 14 provides signals for generating and scanning the specifying magnetic field, which signals are applied to the generator coil C2 and the scanner coil C1 after being processed by the A/D converters 9, 10 and the D.C. power amplifiers 11, 12. Thereby, the computer 14 controls the generation and scanning of the specifying magnetic field ⁇ H s and the generation and stipulation of the linear magnetic field gradient or the planar non-linear magnetic field gradient.
  • the highly stabilized D.C. power source 13 actuates the static electromagnets B under the control of the control computer 14 so as to produce the homogeneous static magnetic field H o .
  • a first embodiment to be described is an NMR imaging method which uses bar-shaped specifying magnetic field which has contour lines of magnetic field intensity extending in parallel to the direction of a homogeneous static magnetic field H o , the contour lines of the magnetic field intensity passing a closed circle on a plane such as x-y plane perpendicular to the homogeneous static magnetic field H o , as disclosed by the inventors in their Japanese Patent Laid-open Publication No. 133,192/1979 and their Japanese Patent Application No. 154,491/1982. If the coordinates of the center of such bar-shaped specifying magnetic field are given by (x', y'), the intensity deviation of the specifying magnetic field at a measuring point (x, y) in a target A is given by
  • h is a constant depending on the number of turns of, the dimensions of, and the current in the generator coil C2 for the specifying magnetic field
  • h o is a homogeneous magnetic field offset superposed on the homogeneous static magnetic field H o .
  • an operation of applying pulse-like magnetic field is effected: Namely, as the bar-shaped specifying magnetic field is scanned in the x-axis and y-axis directions, each time the central position (x', y') of the bar-shaped specifying magnetic field is passed the above-mentioned 90° pulse magnetic field is applied by the transmitter coil D. Immediately after the 90° pulse magnetic field, the bar-shaped specifying magnetic field given by the equation (1) is applied for a time period of ⁇ , as shown in FIG. 3. Thereafter, a static magnetic field having a linear gradient in the z-axis direction is applied by using a coil means therefor which is built in the generator coil C2 for the specifying magnetic field.
  • the curve (a) of FIG. 3 schematically shows the time relation between the specifying magnetic field and the linear magnetic field pulse, and also FID signals generated in response to the application of the magnetic fields.
  • the curve (b) of FIG. 3 shows the timing relating to scanning current pulses to be applied to the scanner coil C1 of the specifying magnetic field.
  • the center of the bar-shaped specifying magnetic field is shifted in the x-axis direction by varying the amplitude of the current pulse of the curve (b).
  • two phase-sensitive detectors 7 are provided, and one of them uses the high-frequency output signals from the RF oscillator 1 as reference signals thereof without any modification, while the other one of them shifts the phase of the high-frequency output signals from the RF oscillator 1 by 90° so as to use the thus phase-shifted HF signals as the reference signals therof.
  • the two detectors 7 produce a pair of FID signals which have a phase difference of 90°.
  • the pair of FID signals can be collectively dealt with as a complex FID signal, by dealing one of the two FID signals as a real part and the other one of the two FID signals as an imaginary part.
  • ⁇ (x,y,z) is a spin density distribution function, which can be treated as the nuclear magnetization distribution function under the thermal equilibrium conditions
  • g is a magnitude of the linear magnetic field gradient in the z direction
  • h o ' is a magnetic field offset superposed on the linear magnetic field gradient in the z-axis direction.
  • T 2 is a spin-spin relaxation time and ⁇ is the magnetogyric ratio.
  • the spin distribution in the z-axis direction can be separated from the Fourier transformation of those FID signals which respond to the application of the static magnetic field with the linear magnetic field gradient in the z-axis direction. Accordingly, since the spin distribution in the z-axis direction can be separated by applying the equation (3) to the FID signals obtained in response to the linear magnetic field gradient, three-dimensional spin distribution can be reproduced by reproducing the x-axis direction and y-axis direction spin distributions for each angular frequency ⁇ .
  • the integral term of the equation (3) is in the form of an integral of the product of the spin density distribution ⁇ (x,y,z) with the phase term e -j ⁇ h ⁇ (x'-x).spsp.2 + (y'-y).spsp.2 ⁇ over the entire region of the variables x, y, and such form of integral represents an integral of the nuclear magnetization at a point (x, y) after rendering a phase shift proportional to the square of the distance from the central position (x', y') of the bar-shaped specifying magnetic field to that point.
  • the phase shift varies in a non-linear manner or in line with a quadratic function.
  • the application of the bar-shaped specifying magnetic field ⁇ H s for a time period of ⁇ results in the encoding of phase of the nuclear magnetization at the point (x, y) in a manner proportional to the square of the distance from the central position (x', y') of the bar-shaped specifying magnetic field to that point.
  • the FID signals for each of the central points (x', y') in response to the pulse series shown in the curve (a) of FIG. 3 are measured after the moment of switching to the linear magnetic field gradient.
  • the details of the means for producing the specifying magnetic field are the same as those disclosed by the inventors in the above-mentioned Japanese Patent Publications.
  • the method for imaging NMR signal according to the invention has been described by referring to three-dimensional imaging of the spin density distribution with the use of the bar-shaped specifying magnetic field, but there are cases where two-dimensional imaging such as two-dimensional tomography is sufficient. Another embodiment of the invention for two-dimensional imaging will be described now.
  • a generator coil C2 has a built-in coil means for producing a linear magnetic field gradient in the y-axis direction. Under the presence of the thus produced linear magnetic field gradient in the y-axis direction, a selective 90° pulse which is wave-shaped so as to have a narrow frequency spectrum is applied to the target A. As a result, nuclear magnetization is excited in the target A only on a slice-like region thereof parallel to the x-z plane, and the thus excited nuclear magnetization has the same Larmor frequency as the frequency spectrum of the selective 90° pulse.
  • the means for the linear magnetic field gradient is switched to the generator coil C2 proper so as to produce the above-mentioned bar-shaped specifying magnetic field extending in the z-axis direction.
  • a linear magnetic field gradient in the z-axis direction is applied to the target A by another coil means built in the generator coil C2 for the specifying magnetic field.
  • Complex FID signals are produced in response to the pulse-like magnetic field application to the target A while shifting the bar-shaped specifying magnetic field in the x-axis direction. If the y-axis direction thickness of the slice-like region of the target A in parallel to the x-z plane is selected to be sufficiently thin, the y-component in the integral of the equation (2) becomes a constant. Accordingly, a two-dimensional tomogram for each selected slice-like region of the target A can be produced by effecting calculation of the equation (5) for two directions, i.e., x-axis and z-axis directions. In this case, the scanning of the bar-shaped specifying magnetic field only in the x-axis direction is sufficient, and scanning in the y-axis direction is of course unnecessary.
  • slice-like regions in parallel to the y-z plane can be selected for the two-dimensional tomography by replacing the x-axis direction with the y-axis direction in the above description of the two-dimensional tomography in parallel to the x-z plane, including the coil means for producing the linear magnetic field gradient which coil means are built in the generator coil C2 for the specifying magnetic field.
  • the above embodiments use the complex FID signals under the constant magnetic field gradient as shown in FIG. 3. It is also possible to use those complex spin echo signals which carry the same information such as the frequency spectrum and the phase.
  • the curve (a) of FIG. 3 after a certain time period t w from the moment of switching the bar-shaped specifying magnetic field to the linear magnetic field gradient in the z-axis direction, either the direction of the magnetic field gradient is reversed or the strong non-selective 180° pulse is applied, and the complex spin echo signals can be observed after another time period t w from the magnetic field gradient reversal or the application of the 180° pulse.
  • the complex spin echo signal has two components, one treated as a real portion and the other one treated as an imaginary portion, which are outputs from two phase-sensitive detectors 7 using two reference frequencies with a phase difference of 90° therebetween.
  • the specifying magnetic field has a profile of magnetic field intensity which varies in accordance with a quadratic function.
  • the bar-shaped specifying magnetic field is one example.
  • the spherical or ellipsoidal specifying magnetic field having spherical or ellipsoidal closed contours of magnetic field intensity, or the star-like specifying magnetic field having saddle-like open contours of magnetic field intensity can be also used in the method of the invention.
  • the central position of any of the above-mentioned specifying magnetic field can be scanned by an electric means.
  • the star-like specifying magnetic field can be easily produced by using annular coil pairs or the like.
  • the spherical or ellipsoidal specifying magnetic field can be produced by using a combination of the coils for producing the bar-like specifying magnetic field and the coils for the above-mentioned bar-shaped specifying magnetic field.
  • x, y, z are spatial coordinates in the target A
  • x', y', z' are the coordinates of the central point of the specifying magnetic field
  • h o is offset superposed on the homogeneous static magnetic field
  • h 1 and h 2 are constants depending on the shape, number of turns, and dimensions of and the current in the coils for producing the specifying magnetic field.
  • h 1 h 2 (>0); in the case of the ellipsoidal specifying magnetic field, h 1 >0, h 2 >0 (h 1 ⁇ h 2 ); and in the case of the star-like specifying magnetic field, h 1 >0 and h 2 ⁇ 0, or h 1 ⁇ 0 and h 2 >0.
  • imaging of the three-dimensional spin density distribution can be effected by the three-dimensional scanning of central positions (x', y', z') of the above-mentioned specifying magnetic field followed by the calculation of the convolution of the sampled data of the complex FID signals for each position at time ⁇ through Fourier transformation.
  • the complex FID signals under the presence of the homogeneous static magnetic field H o . More particularly, in the latter case, after the application of the 90° pulse the specifying magnetic field is applied for a time period ⁇ , in a manner similar to the pulse sequence shown in the curve (a) of FIG. 3, and then the free precession takes place under the presence of the homogeneous static magnetic field H o , so that complex FID signals under the presence of the homogeneous static magnetic field H o after the time ⁇ are detected as functions of time and Fourier transformation are taken from the complex FID signals, for selection and use of the frequency components rendering the maximum values.
  • the value of the Fourier transformation is proportional to the integral of the equation (8).
  • the linear magnetic field gradient of FIG. 3 is replaced with the homogeneous static magnetic field and the scanning is effected for the coordinates x', y', and z'.
  • the processing of the signal in this case is by using an algorithm which is substantially the same as that for the processing of the data by the equations (3) through (6) in the x-axis, y-axis and z-axis directions.
  • a method for imaging NMR signals carrying information on the inside of a target images the information covering the entire inside regions of the target very efficiently in a short time period by using a non-linear magnetic field gradient of comparatively weak intensity.
  • FIG. 4 shows the result of a numerical simulation of imaging NMR signals according to the method of the invention for a simple case of one-dimensional imaging.
  • the numerical simulation relates to the numerical calculation of one-dimensional spin density distribution in the x-axis direction.
  • the curve (a) of FIG. 4 shows reproduced spin density distribution, which is obtained by applying the above-mentioned calculations and operations of FID signals generated for a spin distribution as shown in the curve (b) of the figure.
  • the method for imaging NMR signals for slice-like regions is not restricted to the use of the bar-shaped specifying magnetic field, but the method similarly applies to cases where the spherical or ellipsoidal specifying magnetic field or the star-like specifying magnetic field produced by circular coil pairs is used.

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
US06/656,932 1983-10-11 1984-10-02 Method for imaging nuclear magnetic resonance signals by using non-linear magnetic field gradient Expired - Lifetime US4651098A (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP58189484A JPS6080747A (ja) 1983-10-11 1983-10-11 非線形磁場勾配を用いた核磁気共鳴映像化方法
JP58-189484 1983-10-11

Publications (1)

Publication Number Publication Date
US4651098A true US4651098A (en) 1987-03-17

Family

ID=16242029

Family Applications (1)

Application Number Title Priority Date Filing Date
US06/656,932 Expired - Lifetime US4651098A (en) 1983-10-11 1984-10-02 Method for imaging nuclear magnetic resonance signals by using non-linear magnetic field gradient

Country Status (3)

Country Link
US (1) US4651098A (ja)
JP (1) JPS6080747A (ja)
GB (1) GB2149921B (ja)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4838274A (en) * 1987-09-18 1989-06-13 Air Products And Chemicals, Inc. Perfluoro-crown ethers in fluorine magnetic resonance imaging
US5122748A (en) * 1990-08-13 1992-06-16 The Trustees Of Columbia University In The City Of New York Method and apparatus for spatial localization of magnetic resonance signals
DE4309958C1 (de) * 1993-03-26 1994-09-29 Markus Von Dr Kienlin Verfahren und Vorrichtung zur ortsauflösenden Magnetresonanzuntersuchung eines Meßobjekts
US5530354A (en) * 1994-07-29 1996-06-25 Medical Advances, Inc. Non-monotonic gradient coil system for magnetic resonance imaging
WO2004065977A1 (en) * 2003-01-21 2004-08-05 Koninklijke Philips Electronics N.V. Magnetic resonance method with non-linear magnetic field gradients
EP1359430A3 (en) * 2002-05-01 2005-04-06 GE Medical Systems Global Technology Company LLC Spatial encoding mr data of a moving subject using a higher-order gradient field
US20120025822A1 (en) * 2010-07-23 2012-02-02 Walter Witschey Method of MR (=magnetic resonance) with spatial encoding to generate an image of spectroscopic data
US20120133362A1 (en) * 2010-11-25 2012-05-31 Patrick Gross Magnetic resonance method and system for phase correction of magnetic resonance signals originating in mixed tissue
US20150362571A1 (en) * 2012-12-28 2015-12-17 Centre National De La Recherche Scientifique Method for analysis by nuclear magnetic resonance of a sample including a species to be characterized and a reference species

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62176442A (ja) * 1986-01-29 1987-08-03 横河メディカルシステム株式会社 核磁気共鳴断層撮影装置用スキヤンコントロ−ラ
JPS63216551A (ja) * 1987-03-06 1988-09-08 株式会社日立製作所 核磁気共鳴イメ−ジング装置
DE102005051021A1 (de) 2005-10-25 2007-04-26 Universitätsklinikum Freiburg Apparaturen und Verfahren zur kernspintomographischen Aufnahme mit lokalen Magnetfeldgradienten in Verbindung mit lokalen Empfangsspulen

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4319190A (en) * 1980-03-06 1982-03-09 Bell Telephone Laboratories, Incorporated Nuclear magnetic resonance imaging in space and frequency coordinates
US4442404A (en) * 1978-12-19 1984-04-10 Bergmann Wilfried H Method and means for the noninvasive, local, in-vivo examination of endogeneous tissue, organs, bones, nerves and circulating blood on account of spin-echo techniques
US4558277A (en) * 1982-03-13 1985-12-10 Bruker Medizintechnik Gmbh Method for measuring the nuclear magnetic resonance
US4573015A (en) * 1982-08-31 1986-02-25 Asahikawa Medical College Method of measuring internal information from a target by using nuclear magnetic resonance

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4442404A (en) * 1978-12-19 1984-04-10 Bergmann Wilfried H Method and means for the noninvasive, local, in-vivo examination of endogeneous tissue, organs, bones, nerves and circulating blood on account of spin-echo techniques
US4319190A (en) * 1980-03-06 1982-03-09 Bell Telephone Laboratories, Incorporated Nuclear magnetic resonance imaging in space and frequency coordinates
US4558277A (en) * 1982-03-13 1985-12-10 Bruker Medizintechnik Gmbh Method for measuring the nuclear magnetic resonance
US4573015A (en) * 1982-08-31 1986-02-25 Asahikawa Medical College Method of measuring internal information from a target by using nuclear magnetic resonance

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4838274A (en) * 1987-09-18 1989-06-13 Air Products And Chemicals, Inc. Perfluoro-crown ethers in fluorine magnetic resonance imaging
US5122748A (en) * 1990-08-13 1992-06-16 The Trustees Of Columbia University In The City Of New York Method and apparatus for spatial localization of magnetic resonance signals
DE4309958C1 (de) * 1993-03-26 1994-09-29 Markus Von Dr Kienlin Verfahren und Vorrichtung zur ortsauflösenden Magnetresonanzuntersuchung eines Meßobjekts
US5506504A (en) * 1993-03-26 1996-04-09 Markus Von Kienlin Method and apparatus for conducting a spatially resolving magnetic resonance examination of a test subject
US5530354A (en) * 1994-07-29 1996-06-25 Medical Advances, Inc. Non-monotonic gradient coil system for magnetic resonance imaging
US7558613B2 (en) 2002-05-01 2009-07-07 General Electric Company Spatial encoding MR data of a moving subject using a higher-order gradient field
EP1359430A3 (en) * 2002-05-01 2005-04-06 GE Medical Systems Global Technology Company LLC Spatial encoding mr data of a moving subject using a higher-order gradient field
US20060036161A1 (en) * 2002-05-01 2006-02-16 Hinks Richard S Spatial encoding mr data of a moving subject using a higher-order gradient
US20060061359A1 (en) * 2003-01-21 2006-03-23 Koninklijke Philips Electronics N.V. Magnetic resonance method with non-linear magnetic field gradients
US7208949B2 (en) 2003-01-21 2007-04-24 Koninklijke Philips Electronics N. V. Magnetic resonance method with non-linear magnetic field gradients
WO2004065977A1 (en) * 2003-01-21 2004-08-05 Koninklijke Philips Electronics N.V. Magnetic resonance method with non-linear magnetic field gradients
US20120025822A1 (en) * 2010-07-23 2012-02-02 Walter Witschey Method of MR (=magnetic resonance) with spatial encoding to generate an image of spectroscopic data
US8791700B2 (en) * 2010-07-23 2014-07-29 Universitaetsklinikum Freiburg Method of MR (=magnetic resonance) with spatial encoding to generate an image or spectroscopic data
US20120133362A1 (en) * 2010-11-25 2012-05-31 Patrick Gross Magnetic resonance method and system for phase correction of magnetic resonance signals originating in mixed tissue
US8866477B2 (en) * 2010-11-25 2014-10-21 Siemens Aktiengesellschaft Magnetic resonance method and system for phase correction of magnetic resonance signals originating in mixed tissue
US20150362571A1 (en) * 2012-12-28 2015-12-17 Centre National De La Recherche Scientifique Method for analysis by nuclear magnetic resonance of a sample including a species to be characterized and a reference species
US10241169B2 (en) * 2012-12-28 2019-03-26 Centre National De La Recherche Scientifque Method for analysis by nuclear magnetic resonance of a sample including a species to be characterized and a reference species

Also Published As

Publication number Publication date
JPS6080747A (ja) 1985-05-08
GB8425216D0 (en) 1984-11-14
GB2149921B (en) 1987-02-11
JPH0263009B2 (ja) 1990-12-27
GB2149921A (en) 1985-06-19

Similar Documents

Publication Publication Date Title
US5365172A (en) Methods and apparatus for MRI
US4516075A (en) NMR scanner with motion zeugmatography
US7382129B2 (en) 4 dimensional magnetic resonance imaging
US4431968A (en) Method of three-dimensional NMR imaging using selective excitation
Kwiat et al. A decoupled coil detector array for fast image acquisition in magnetic resonance imaging
US4307343A (en) Moving gradient zeugmatography
EP0042254B2 (en) Nuclear resonance apparatus including means for rotating the gradient of a magnetic field
US4573015A (en) Method of measuring internal information from a target by using nuclear magnetic resonance
JPS58151545A (ja) 核磁気共鳴画像形成方法と装置
US4651098A (en) Method for imaging nuclear magnetic resonance signals by using non-linear magnetic field gradient
US4910460A (en) Method and apparatus for mapping eddy currents in magnetic resonance imaging
US4654591A (en) NMR flow imaging using bi-phasic excitation field gradients
US4594550A (en) Method of scanning specifying magnetic field for nuclear magnetic resonance imaging
US4697149A (en) NMR flow imaging using a composite excitation field and magnetic field gradient sequence
US4876508A (en) Method and apparatus for NMR imaging
JP2003500138A (ja) Ｍｒｉコイル配列からの平行したデータ収集のための方法と装置
EP0246327A1 (en) Method and apparatus for nmr imaging
Bottomley Nuclear magnetic resonance: Beyond physical imaging: A powerful new diagnostic tool that uses magnetic fields and radio waves creates pictures of the body's internal chemistry
US11885864B2 (en) Magnetic resonance imaging apparatus and method of controlling the same
JPH09154831A (ja) Ｍｒイメージング方法及びｍｒｉ装置
JPH0374100B2 (ja)
CA1203282A (en) Nmr scanner with motion zeugmatography
JPH0316133B2 (ja)
JPH0245037A (ja) 磁気共鳴イメージング装置
JPS60201242A (ja) Νｍｒ断層像撮影装置

Legal Events

Date Code Title Description
AS Assignment

Owner name: UTSUNOMIYA UNIVERSITY, 350, MINE-MACHI, UTSUNOMIYA

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST.;ASSIGNORS:YAMADA, YOSHIFUMI;TANAKA, KUNIO;ABE, ZENWEMON;REEL/FRAME:004369/0205;SIGNING DATES FROM 19840912 TO 19840914

STCF Information on status: patent grant

Free format text: PATENTED CASE

FEPP Fee payment procedure

Free format text: PAYOR NUMBER ASSIGNED (ORIGINAL EVENT CODE: ASPN); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY

FPAY Fee payment

Year of fee payment: 4

FEPP Fee payment procedure

Free format text: PAT HLDR NO LONGER CLAIMS SMALL ENT STAT AS NONPROFIT ORG (ORIGINAL EVENT CODE: LSM3); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY

FPAY Fee payment

Year of fee payment: 8

FEPP Fee payment procedure

Free format text: PAT HOLDER CLAIMS SMALL ENTITY STATUS - SMALL BUSINESS (ORIGINAL EVENT CODE: SM02); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY

REFU Refund

Free format text: REFUND - PAYMENT OF MAINTENANCE FEE, 12TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: R185); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY

FPAY Fee payment

Year of fee payment: 12