WO2001058352A1 - Dispositif d'imagerie par resonance magnetique - Google Patents

Dispositif d'imagerie par resonance magnetique Download PDF

Info

Publication number
WO2001058352A1
WO2001058352A1 PCT/JP2001/000832 JP0100832W WO0158352A1 WO 2001058352 A1 WO2001058352 A1 WO 2001058352A1 JP 0100832 W JP0100832 W JP 0100832W WO 0158352 A1 WO0158352 A1 WO 0158352A1
Authority
WO
WIPO (PCT)
Prior art keywords
measurement
space
magnetic field
origin
measured
Prior art date
Application number
PCT/JP2001/000832
Other languages
English (en)
Japanese (ja)
Inventor
Takayuki Abe
Shigeru Watanabe
Original Assignee
Hitachi Medical 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
Priority claimed from JP2000029817A external-priority patent/JP3847512B2/ja
Priority claimed from JP2000097632A external-priority patent/JP3847519B2/ja
Application filed by Hitachi Medical Corporation filed Critical Hitachi Medical Corporation
Priority to US10/203,260 priority Critical patent/US6611144B2/en
Publication of WO2001058352A1 publication Critical patent/WO2001058352A1/fr

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5601—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution involving use of a contrast agent for contrast manipulation, e.g. a paramagnetic, super-paramagnetic, ferromagnetic or hyperpolarised contrast agent

Definitions

  • the present invention relates to a magnetic resonance imaging apparatus (hereinafter, referred to as an MRI apparatus) for obtaining a tomographic image of a desired part of a subject using a nuclear magnetic resonance (hereinafter, abbreviated as “MRJ”) phenomenon, and particularly to a vascular system.
  • MRI apparatus magnetic resonance imaging apparatus
  • MRJ nuclear magnetic resonance
  • the present invention relates to an MRI apparatus capable of securing a desired rendering range and image quality in a minimum time required for rendering a travel.
  • An MRI device measures the density distribution, relaxation time distribution, etc. of nuclear spins (hereinafter simply referred to as spins) at a desired examination site in a subject using NMR phenomena, and uses the measured data to determine the Are displayed as images.
  • spins nuclear spins
  • MRA MR angiography
  • 'A common method of using a contrast agent is to combine a T1 shortened contrast agent such as M-DTPA with a short TR sequence in the gradient echo system.
  • This method includes a T1 shortened contrast agent Since the blood flow spin has a T1 shorter than that of the surrounding tissue, saturation is unlikely to occur when a high-frequency magnetic field excitation is applied for a short repetition time TR, and a relatively high signal is emitted from other tissues.
  • This method is used to visualize blood vessels filled with blood containing a contrast agent with high contrast to other tissues Measurement of Volume data including blood vessels while the contrast agent stays in the target blood vessel (Specifically, three-dimensional measurement) is performed, and the obtained three-dimensional images are superimposed to perform a projection process to draw a blood flow. Three-dimensional to get information Based on the three-dimensional Daladier cement echo method to get over data Is used.
  • (1) injection method of a contrast agent, and (2) imaging time and timing are important.
  • the contrast agent must be injected into the blood vessel to be imaged so as to maintain a stable and high concentration. For this reason, rapid injection using an automatic injector is generally used.
  • the imaging timing so that the concentration of the contrast agent in the artery becomes high during data collection.
  • the contrast agent concentration peaks in the center of k-space (low frequency region), which controls the contrast of the image, and the timing is set according to the data collection method used and the sequence of the sequence. I do.
  • the data collection method mainly consists of a sequential order in which measurement is performed from one high-frequency side of k-space to the other high-frequency side through a low-frequency region, and measurement is alternately performed from the low-frequency region of k-space toward both high-frequency ends.
  • There is a centric order and centric trick order is generally used.
  • Three-dimensional measurement ( ⁇ In the centric order, one of the phase-encoding-donorap and the slice-en code loop is the outer loop, the other is the inner loop, and either or both are controlled by centric, order.
  • the centric order in this case is not a true centric order because the distance from the origin in k-space to the measurement point (sampling point) fluctuates as shown in Fig. 1 (b), In some cases, and arteriovenous separation was insufficient. '.
  • the optimal imaging timing force is shifted by S ⁇
  • the image quality deteriorates due to the acquisition of low-frequency information when the agent is thin.
  • the measurement time is too early, the data in the low-frequency region will be sampled during the time when the signal from the blood vessel is extremely low, while the data in the high-frequency region will have a high signal from the blood vessel. Since the signal is sampled during the time zone, a ringing artifact without a DC component is generated.
  • the imaging time becomes longer because the measurement is performed in the high-frequency region in the vertical and horizontal directions of the k space with respect to the origin.
  • an object of the present invention is to provide an MRI apparatus capable of delineating the entire target blood vessel with high contrast in a short time while reducing the influence on the image quality due to the shift of the optimal imaging timing. It is another object of the present invention to provide an MRI apparatus which is hardly affected by body motion and can separate and depict arteries and veins with MRA. It is another object of the present invention to provide a data collection method suitable for MRA. Disclosure of the invention
  • the measurement points in the k space are divided into two groups, and in the first group to be measured first, the distance from the origin in the k space gradually approaches U In the other group of measurements, the measurement order of the measurement points is changed from the low-frequency component so that the distance gradually increases.
  • Adopt a data collection method that controls toward high frequency components.
  • the MRI apparatus of the present invention comprises: a static magnetic field generating means for generating a static magnetic field in a space where a subject is placed; and a gradient magnetic field generating means for applying a gradient magnetic field in the Rice, phase encoding, and reading directions to the space.
  • a transmission system for irradiating a high-frequency magnetic field to cause nuclear magnetic resonance to an atomic nucleus of the biological tissue of the subject; a receiving system for detecting an echo signal emitted by the nuclear magnetic resonance;
  • a gradient magnetic field generation means a control system for controlling the transmission system and the reception system, a signal processing system for performing an image rain synthesis operation using an echo signal detected by the reception system, and a means for displaying the obtained image.
  • the control system executes a three-dimensional sequence for giving a slice code and a phase encode, and at this time, a total of k space defined by the number of slice encodes and the number of vertical encodes is provided.
  • the points are divided into two groups.In the first group, the distance from the origin in k-space to the measurement point gradually decreases in the order of measurement.In the second group, the distance from the origin in k-space to the measurement point is measured.
  • the gradient magnetic field generating means is controlled in the slice direction and the phase encoder direction so that the values gradually increase in the order of measurement.
  • the method of dividing the measurement points in k-space into two groups is that at least one group contains the measurement points from the low-frequency region to the high-frequency region, and the measurement point of the other group contains at least the measurement points in the low-frequency region. I just need. Also, the number of measurement points actually measured among the measurement points belonging to the two groups may be the same or different. That is, one of the forces of the two groups may include a non-measurement point (a point that is not measured).
  • the measurement points are divided into two groups according to the k-space region.
  • the control system attaches the slice encode and the vertical encode.
  • the k-space defined by the number of slice encodes and the number of phase encodes is divided into two, and the distance from the origin in k-space to the measurement point in k-space is
  • the gradient magnetic field generating means in the slice direction and the phase encode direction is controlled so that the distance gradually decreases in the order of measurement, and in the other region, the distance from the origin on the k-space gradually increases in the order of measurement.
  • the measurement points in the k space are divided into two groups having a complex conjugate relationship with each other.
  • the control system executes the three-dimensional sequence for applying the slice encoding and the phase encoding, and at this time, sets the measurement point of the lc space defined by the slice encoding number and the phase encoding number to the origin.
  • the gradient magnetic field generating means in the slice direction and the cross-encode direction are controlled so as to gradually increase in the order of measurement.
  • the adjacent measurement points it is preferable to divide the adjacent measurement points so that they belong to different groups.
  • the measurement points that are partially adjacent to each other near the origin must belong to the same group. Therefore, in this specification, “division of adjacent measurement points so that they belong to different groups” means that the condition that the complex conjugate relationship is satisfied and the condition that adjacent measurement points belong to different groups is maximal. A state that is satisfied.
  • control system does not measure all of the measurement points in one of the two divided areas but measures a small number of measurement points compared to the other area.
  • the control for measuring is performed.
  • control is performed to measure a smaller number of measurement points than in the other group.
  • a predetermined region of a subject is selected and excited, and a gradient magnetic field that encodes at least two directions is applied.
  • the step of measuring one signal is repeated a plurality of times while changing the strength of the gradient magnetic field.
  • the measurement space defined by the two-direction encoder gradient magnetic field strength is divided into two, Measurement is performed sequentially on the two divided areas, and at that time, in the first measurement area, the measurement is performed so that the distance from the origin of the measurement space to the measurement point gradually decreases in the measurement order, and the area measured later In, measurement is performed such that the distance from the origin of the measurement space to the measurement point gradually increases in the order of measurement.
  • the point closest to the k-space from the previous measurement point is measured as the next measurement point.
  • the three-dimensional image data acquisition method of the present invention includes a step of selecting and exciting a predetermined region of the subject, applying a gradient magnetic field that encodes in at least two directions, and measuring an echo signal generated from the region.
  • the measurement points in the measurement space defined by the two-direction encode gradient magnetic field intensities are shared with the origin and are mutually complex. Adjacent measurement points in a conjugate relationship are divided into first and second groups so that they belong to different groups, and measurements are sequentially performed on the first and second groups.
  • the distance from the origin of the measurement space to the measurement ⁇ is measured so as to gradually decrease in the order of measurement
  • the distance from the origin of the measurement space to the measurement point is measured. So that it gradually increases in the order of measurement To measure.
  • Fig. 1 (a) measurement without distance fluctuation from the origin can be performed, and the time when the lowest frequency component is measured and the target blood By matching the point at which the signal intensity of the tube peaks with the contrast agent, the target blood vessel can be drawn with high contrast. In addition, even if there is a slight deviation between the point at which the low-frequency component is measured and the peak of the signal strength, the low-frequency component can be reliably measured, and there is no image deterioration.
  • Figures 1 (b) to 1 (d) show the conventional Centric orderer, elliptic force centric order, and sequential channel order. Indicates the change in distance from the k-space origin. -According to a preferred aspect of the data collection method of the present invention, in the first group, of all the measurement points, only the measurement points of the minus part are measured, and in the second group, all the measurement points are measured.
  • one group can collect the data that was not measured without measuring some of the measurement points.
  • the measurement of unnecessary data having a low signal intensity can be eliminated, and a good image can be obtained.
  • FIG. 1 is a schematic diagram illustrating a data collection method adopted by an MRI apparatus according to the present invention and a conventional data collection method.
  • FIG. 2 is a block diagram illustrating an overall configuration of the MRI apparatus to which the present invention is applied.
  • FIG. 4 is a diagram showing an embodiment of a pulse sequence of a contrast MRA measurement performed by the MRI apparatus of the present invention.
  • FIG. 4 is a diagram schematically showing an embodiment of a k-space data collection order according to the present invention.
  • FIG. 5 is a view for explaining MRA imaging by the MRI apparatus of the present invention
  • FIG. 6 is a view showing a result of simulation of MRA imaging by the MRI apparatus of the present invention and MRA imaging by the conventional method
  • FIG. 8 is a diagram schematically illustrating another embodiment of the k-space data collection order according to the present invention.
  • FIG. 8 is a diagram schematically illustrating another embodiment of the k-space data collection order according to the present invention.
  • FIG. 10 is a diagram for explaining MRA imaging by the MRI apparatus of the present invention
  • FIG. 11 is a diagram schematically showing another embodiment of the k-space data collection order according to the present invention
  • FIG. FIG. 10 illustrates an image reconstruction method to which the data acquisition method of FIG. 10 is applied
  • FIG. 13 is a diagram illustrating a simulation for evaluating MRA imaging by the MRI apparatus of the present invention
  • FIG. 15 is a figure which shows typically the other example of the data collection order of k space by this invention.
  • FIG. 2 is a block diagram showing the overall configuration of the MRI apparatus according to the present invention.
  • This MKI device obtains a tomographic image of the subject using the ⁇ R phenomenon, and includes a static magnetic field generating magnet 2, a magnetic field gradient generating system 3, a sequencer 4, a transmitting system 5 ', and a receiving system 6. , A signal processing system 7, and a central processing unit (CPU) 8.
  • CPU central processing unit
  • the static magnetic field generating magnet 2 generates a uniform static magnetic field around the subject 1 in the direction of its body axis or in a direction orthogonal to the body axis, and has a certain space around the subject 1.
  • a permanent magnet type, normal conduction type or superconducting type magnetic field generating means is provided.
  • the magnetic field gradient generating system 3 includes a gradient magnetic field coil 9 wound in three directions of X, ⁇ , and Z, a gradient magnetic field IT source 10 for driving each gradient magnetic field coil, and a force.
  • a gradient magnetic field power supplies 10 of the respective coils By driving the gradient magnetic field power supplies 10 of the respective coils in accordance with the instructions from 4, the gradient magnetic fields Gx, Gy, Gz in the three axial directions of X, ⁇ , and Z are applied to the subject 1.
  • this gradient magnetic field By applying this gradient magnetic field, a specific slice or slab of the subject 1 can be selectively excited, and the position of the measurement point (sampling point) in the measurement space (k space) and the measurement order must be specified. Can be.
  • the sequencer 4 operates under the control of the CPU 8 and sends various commands necessary for data collection of tomographic images of the subject 1 to the magnetic field gradient generating system 3, the transmitting system 5, and the receiving system 6.
  • the operation timing of the magnetic field gradient generation system 3, the transmission system 5, and the reception system 6 controlled by the sequencer 4 is called a pulse sequence.
  • a sequence for three-dimensional blood flow imaging is adopted as one of the pulse sequences. The control of the sequencer 4 will be described later in detail.
  • the transmission system 5 irradiates a high-frequency magnetic field to cause nuclear magnetic resonance in nuclei of atoms constituting the biological tissue of the subject 1 by high-frequency pulses sent from the sequencer 4, and modulates with the high-frequency oscillator 11. 12 and high-frequency amplifier 13 and high frequency on the transmitting side And a wave coil 14a.
  • the high-frequency pulse output from the high-frequency oscillator 11 is amplitude-modulated by the modulator 12 according to the instruction of the sequencer 4, and the high-frequency pulse subjected to the amplitude modulation is amplified by the high-frequency amplifier 13 and then placed close to the subject 1.
  • the electromagnetic wave is applied to the subject 1 by supplying it to the high-frequency coil 14a.
  • the receiving system 6 detects the echo signal (fan R signal) emitted by nuclear magnetic resonance of the nucleus of the living body a ⁇ of the subject 1, and includes a high-frequency coil 14b and an amplifier 15 on the receiving side and a quadrature phase detector. It consists of a container 16, A / D conversion ⁇ 3 ⁇ 4 17, and power.
  • the electromagnetic (response signal) of the response of the subject 1 due to the electromagnetic wave radiated from the high-frequency coil 14a on the transmitting side is detected by the high-frequency coil 14b arranged close to the subject 1.
  • the detected echo signal is input to the A / D converter 17 via the amplifier 15 and the quadrature detector 16 and is converted into a digital value. Further, the quadrature detector is output at a timing according to an instruction from the sequencer 4.
  • the collected data is sampled into two series and is sent to the signal processing system 7.
  • the signal processing system 7 includes a CPU 8, a recording device such as a magnetic disk 18 and a magnetic tape 19, and a display 20 such as a CRT.
  • the CPU 8 performs processing such as Fourier transform and correction coefficient calculation image reconstruction on the echo signal (digitized data) sent from the receiving system 6 and performs an appropriate operation on the signal intensity distribution of an arbitrary cross section or a plurality of signals. The obtained distribution is imaged and displayed on the display 20 as a tomographic image.
  • the high-frequency coils 14a and 14b on the transmitting side and the receiving side and the gradient magnetic field coil 9 are installed in the magnetic field space of the static magnetic field generating magnet 2 arranged in the space around the subject 1. .
  • the sequencer 4 operates according to a predetermined pulse sequence, here a three-dimensional MM sequence, according to the operation timing of the magnetic field gradient generation system 3, the transmission system 5, and the reception system 6. Control the ringing.
  • This pulse sequence is pre-installed as a program in the memo provided in the CPU 8), and can be executed by the user selecting as appropriate according to the purpose of photographing, similarly to other pulse sequences. That is, the sequence 4 in which the MRA using the contrast agent is selected via the input device of the CPU 8 is controlled by the CPU 8 to execute a three-dimensional MRA sequence.
  • This panless sequence is a sequence based on the gradient echo method, as shown in FIG. 3, for example, and is common to three-dimensional MRA sequences. That is, the high-frequency magnetic field panelless RF is simultaneously applied with the region selection gradient magnetic field Gs to excite the region (slab) including the target blood vessel, and then the gradient magnetic field pulse Gel in the slice direction and the gradient in the directional code direction. The magnetic field pulse Ge2 is applied, then the readout gradient magnetic field Gr is applied, and the polarity is inverted to measure the echo signal. From the high-frequency magnetic field pulse RF to the echo signal measurement, the magnetic field strength of the gradient magnetic field Gel in the slice direction and the gradient magnetic field Ge2 in the phase encode direction are changed at a predetermined repetition time TR to obtain three-dimensional data.
  • the number of encoders in the slice direction and the phase encoder direction determines the image resolution in the direction A, and is set in advance in consideration of the measurement time and the like. For example, the number of encodes in the phase encoding direction is set to 128, 256, and the slice direction is set to 10 to 30.
  • the k-space (ky-kz space, where the slice direction is the z direction and the phase encoding direction is the y direction) is defined by the number of encoders in the slice direction and the phase encode direction. That is, in the sequence of FIG. 3, the signals measured when the gradient magnetic field intensity in the slice direction is a certain value Gel (Gz) and the gradient magnetic field intensity in the phase encoder direction is a certain value Ge2 (Gy) are Gz and Gy. Is located at the grid point (ky, kz) in k-space corresponding to.
  • the three-dimensional MRA sequence itself shown in FIG. 3 is general in MRA, but in this sequence adopted in the present embodiment, the data collection method uses the conventional centric orderer or elliptic force sensor. This is different from trick ordering.
  • the ky-kz space is divided into two along the ky axis or the kz axis, and in one area where measurement is first started, the area is determined from the point at a large distance from the origin 0 in the k space. Start the measurement, and then perform sampling control from the high-frequency component to the low-frequency component so that the sampling point gradually approaches the origin 0.
  • the sampling point is controlled from the low frequency component to the high frequency component so that the distance from the origin 0 gradually increases from or near the origin 0.
  • this data collection method is as follows: 1) One of the two divided areas starts from the point farthest from the origin, and thereafter, determines the subsequent sampling points so that the distance from the origin gradually decreases. In the other area, start from the origin or the point closest to the origin, and then determine the subsequent sampling points so that the distance from the origin gradually increases. 2) Make the distance between sampling points adjacent in time the shortest. Is defined by the AND of the two conditions. .
  • the condition of 2) is not essential, but by making sampling points adjacent to each other as close as possible, it is possible to reduce the artifat.
  • the condition in 2) may be to select not the distance between the two sampling points but the force with the same ky value, for example.
  • the optimum sampling point may be determined in consideration of not only the relationship between two sampling points but also the relationship between a plurality of sampling points in the next measurement and further subsequent measurements.
  • Fig. 4 shows a simplified example of the data collection method described above, showing the data collection order in a 5 * 9 Matritus k-space where the number of slice codes is 5 and the number of phase codes is 9.
  • circled numbers indicate the data collection order. This divided into two k-space in kz axis, toward the origin (No. 25) a point in a high frequency region indicated with the lower region (E ⁇ C region) "S tart" from (No. 1) in numerical order Measure In the upper ⁇ area (C ⁇ E area), measurement is performed from the origin to the high frequency area in the numerical order up to number 45. Note that the distance ⁇ z between two adjacent points on the coordinates of the k space is 1 ZFOVz.
  • Timing imaging is performed by the test injection method.
  • a small amount of a contrast medium (about 1 to 2 ml) is first injected for test injection to obtain a time signal curve at the target site as shown in FIG.
  • the arrival time tl of the contrast agent is measured from this curve, and the timing for performing the main imaging is determined based on the result.
  • the timing imaging method sets R0I at a specific site in the monitor area for the arrival of the contrast agent, captures the signal change of the same site, and exceeds the set threshold.
  • a method of automatically starting imaging at a point in time or a method called fluoroscopy, in which a target blood vessel is observed in real time by repeating short-time imaging and display, and imaging is started when an appropriate signal rise is obtained. It is also possible to employ these methods. However, the test injection method is preferable because the timing can be accurately measured by using a contrast agent prior to the main imaging.
  • the main imaging is performed as shown in FIG. 5 (b).
  • the main imaging may be performed only after the injection of the contrast agent, but preferably, images before and after the injection of the contrast agent are taken.
  • the imaging before and after the imaging is performed continuously for the same slice or slab position under the same conditions.
  • the imaging sequence is a sequence based on the short TR 3D gradient echo method as shown in Fig. 3.
  • a gradient magnetic field for rephasing the dephasing due to the flow that is, Gradient Moment Nulling may be added, but this is not essential. Instead, it is preferable to use a rather simple gradient echo for shortening TR / TE.
  • the imaging time T is determined when the pulse sequence repetition time TR, matrix size (slice chain code number and phase end number), and addition number are determined, so that the target blood vessel obtained by the above timing imaging reaches the contrast agent of interest.
  • imaging start time t2 imaging is started after injection of contrast agent so that data measurement in the low-frequency region of the ky-kz space is performed when the contrast agent reaches the target blood vessel. Set the time up to).
  • the sequencer 4 uses the gradient magnetic field in the slice direction in the region to be measured first (for example, the E ⁇ C region). And the gradient magnetic field noise in the direction of the phase encoder are controlled so that high-frequency components and low-frequency components are measured in order, and in the subsequent measurement area (for example, C ⁇ E area), low-frequency components and high-frequency components are measured. Control to measure the thigh. In this case, as described above, control is performed so that the distance from the origin to the sampling point gradually decreases in the first region, and the distance from the origin gradually increases in the subsequent region.
  • Figure 6 shows the results of simulating the difference in arteriovenous separation due to the different data collection methods.
  • FOV 320
  • TR 10 ms
  • number of phase encodes 160
  • number of slice encodes 16
  • image matrix 256 *
  • slice 5
  • the arteriovenous separation is expressed as the ratio of the signal strength of the artery to the signal strength of the vein.
  • the data collection method of the present invention has an increased arteriovenous signal ratio as compared to sequential ordering and elliptical centric ordering. Therefore, even if there is a vein near the artery that is confusing, it is possible to draw only the artery with high contrast.
  • the three-dimensional data of only the blood vessel can be obtained.
  • the difference processing is performed, for example, by performing a complex difference between slices at the same slice position in three dimensions.
  • the difference may be a difference between absolute values.
  • the method of removing tissue other than blood vessels by performing differential processing between images before and after contrast in this manner is called 3D MR-DSA (Digital Subtraction Angiography) and is a known method, and is not essential in the present invention. In particular, it is suitable for delineating thin blood vessels in which it is difficult to obtain sufficient contrast with tissues other than blood vessels.
  • the three-dimensional data after the difference processing can be viewed three-dimensionally by projecting it in any direction, such as coronal section, sagittal section, and transverse axis.
  • a projection method a known maximum intensity projection method or the like can be employed.
  • a sequence based on the gradient echo method is exemplified as a three-dimensional MRA sequence, but an EPI (Echo Planer Imaging) method for measuring a plurality of echo signals with one excitation, a split-type EPI, and the like are also employed. can do.
  • EPI Echo Planer Imaging
  • the k space may be divided by the force ky axis which indicates the case where the k space is divided by the kz axis.
  • the areas may be asymmetric.
  • the imaging time as a whole can be reduced.
  • the number of measurement points in the area to be measured first is reduced, or the number of measurement points in the area to be measured later is reduced, so that the imaging time is shorter than in Fig. 5 (b).
  • the measurement method shown in Fig. 5 (d) is effective when the distance between the target blood vessel and the nearby vein is short, and the time when the contrast agent reaches the target blood vessel after reaching the target blood vessel is short.
  • the data in the area with a small number of measurement points may be estimated from the data in the area where the unmeasured data (the data in the shaded area in Fig. 7) is measured, or May be. Also in this embodiment, in addition to the above-described effects, it is possible to selectively obtain a high-contrast artery image while reducing the influence of the shift in the three-dimensional measurement timing.
  • data collection in k-space is not limited to rectangular matrices as shown in Figs. 4 and 7, and data within a circle (ellipse) centered on the origin (No. 15) as shown in Fig. 8 It is also possible.
  • the number of slice encodes is 5, and the number of phase encodes is 9.
  • both ky and kz are high-frequency components. The area (outside the circle) is collected, but data concentrically arranged around the origin (No. 15) is collected.
  • the order of data collection is as follows: in the lower half area, which is measured first, from the center to the center in order from the side farther from the center, and in the upper half area, the measurement is in the direction away from the center. I do.
  • the embodiment in which the region of the k space is divided and the measurement points are divided into two groups according to the region has been described. It is also possible to divide them into groups.
  • FIG. 9 shows a k-space of an 8 * 8 matrix in which the number of slice encodes is 8 and the number of phase encodes is 8, as an example of the data collection method according to the second embodiment, which is zero-purified.
  • Lattice points belonging to these two groups have a complex conjugate relationship with each other, and adjacent lattice points belong to different groups. However, in order to satisfy the complex conjugate relation, adjacent lattice points belong to one group near the origin. '
  • the first of these two groups which measures first, starts the measurement from the point in k-space where the distance from the origin 0 force is large, and then gradually increases the sampling point from the high-frequency component so as to approach the origin 0 gradually. Sampling control is performed toward low frequency components. Also the second In the group, on the contrary, the sampling point is controlled from the low-frequency component to the high-frequency component such that the distance from the origin 0 gradually increases from or near the origin 0.
  • circled numbers indicate the data collection order. No hierarchy is the measurement point of the same number, indicating that it may be measured from either of them ⁇
  • the measurement starts from the grid point (number 1) farthest from the origin (the grid point with number 33), that is, the highest frequency component, and then the grid point of number 2 Measurements are sequentially performed up to the origin, such as the grid point of number 3.
  • measurement of the second group is performed.
  • measurement is started from the lattice point (No. 34) closest to the origin, that is, the low-frequency component, and measurement is performed in order from the origin in order.
  • the subject is placed in the measurement space in the static magnetic field magnet, an imaging region including a target blood vessel is determined, and timing imaging is performed.
  • the timing imaging is performed by, for example, a test injection method. That is, a small amount of a contrast medium (about 1 to 2 ml) was test-injected, and as shown in Fig. 10, the time at the target site " ⁇ sign curve was obtained, from which the arrival time of the syrup (signal intensity) Measure the tl and determine the timing to perform the main imaging based on the result.
  • 'Timing imaging perform the main imaging as shown in Fig. 10 (b). May perform only imaging after the injection of the contrast medium, but preferably images before and after the injection of the contrast medium.Imaging before and after the imaging is performed under the same conditions under the same slice or slab position. Is performed continuously.
  • the imaging sequence is a sequence based on the short TR 3D gradient echo method as shown in Fig. 3.
  • a gradient magnetic field that is, Gradient Moment Nulling may be added for rephasing the phase due to the flow, but this is not essential. Instead, it is preferable to use a simple gradient echo for shortening TR / TE.
  • the imaging time T is determined. Based on tl, imaging start time t2 (from the injection of the contrast agent to the start of imaging, so that data measurement in the low-frequency region of the ky-kz space is performed when the contrast agent reaches the target blood vessel. Time).
  • imaging measurement of the first group is first started, and then measurement of the second group is performed.
  • the sequencer 4 controls both the gradient magnetic field pulse in the slice direction and the gradient magnetic field in the phase encode code in the first group to be measured first so that the high frequency component and the low frequency component are sequentially measured, and then the measurement is performed thereafter.
  • control is performed so that measurement is performed in order from low-frequency components to high-frequency components. In this way, three-dimensional image data after contrast is obtained.
  • the data collection method shows a case where all the measurement points belonging to both the first group and the second group are measured, as shown in Fig. 10 (c).
  • the first group may employ a data collection method in which measurement of a predetermined high-frequency component is omitted and low-frequency components are measured in a short time.
  • An example of such a data collection method is shown in FIG. ⁇ 11 also illustrates an 8 * 8 matrix k-space in which the number of slice switches is 8, and the number of phase switches is 8.
  • the k-space is divided into two groups under the same conditions as in the embodiment shown in Fig. 9, but here the measurement is performed first. Measure only the components.
  • grid points present in 4 * 4 madridas in the low frequency region are measured among grid points in the k space.
  • the grid point with the greatest distance of the origin is used as the starting point, and the grid point of number 2 and the grid point of number 3 are measured up to the origin.
  • the grid point (number 10) adjacent to the origin From the origin, and measure all grid points belonging to the second group up to the highest frequency component in the order away from the origin.
  • the data of the high frequency region that is not measured in the first group can be estimated based on the complex characteristics of the first group and the second group.
  • a method of estimating unmeasured data a method based on a known half Fourier reconstruction method can be employed.
  • FIG. 12 schematically shows these processes.
  • the measured data is subjected to one-dimensional Fourier transform in the frequency encoder direction (k x direction) to obtain the actual measured data in the three-dimensional hybrid space.
  • the three-dimensional estimated data is obtained from the actual measurement data, and the hybrid spatial data is obtained by combining the actual measurement data and the estimation data.
  • Three-dimensional image data is obtained by subjecting the hybrid space data to two-dimensional Fourier transform.
  • the spatial resolution is not degraded.
  • the low-frequency component in k-space is measured when the contrast agent reaches the target blood vessel and the signal intensity of the blood flowing through the target blood vessel becomes the highest. This means that artery images can be drawn with high contrast.
  • the concentration of the contrast agent rises sharply.
  • FIG. 13 and 14 show the results of simulating differences in arteriovenous separation due to differences in imaging methods (data collection methods).
  • a simulated artery and vein was used, and a contrast medium was flowed into it at a flow rate of 40 cm, an arterial venous return time of 7 seconds, and an injection speed of 2 cc / s.
  • FIG. 13 shows the signal strength under the above conditions.
  • the peak of the signal from the artery is seen first, and the peak of the signal from the vein appears later.
  • FIG. 14 (a) is an image obtained by the imaging method of the present invention
  • FIG. 14 (b) is an image obtained by elibutanore centric ordering.
  • veins as well as arteries are imaged, and arteriovenous separation is not complete, whereas the imaging method of the present invention allows only arteries to be drawn with bulk contrast.
  • data collection in k-space is not limited to the rectangular matrix shown in Fig. 9 and Fig. 11, and data in a circle (ellipse) centered on the origin as shown in Fig. 15 can also be collected. It is possible.
  • measurement starts at a point in k-space that is farther from the origin 0, and then gradually moves the sampling point from the high-frequency component to the low-frequency component so as to approach the origin 0 gradually. Sampling control.
  • the sampling point is controlled from the low-frequency component to the high-frequency component such that the distance from the origin 0 gradually increases from or near the origin 0.
  • the grid points (measurement points) in the ky-kz space are divided into two groups, and one of the groups measured first has a distance from the origin in the k space.
  • the sampling point is controlled from the high-frequency component to the low-frequency component so that the sampling point gradually approaches, and in the other group, the sampling point is shifted from the low-frequency component to the high-frequency component so that the distance gradually increases.
  • Sampling control is performed for each component, so that the effect of high contrast can be achieved while reducing the effect of the imaging timing shift.
  • a vein-separated surface image can be obtained.
  • the imaging time can be reduced by reducing the number of measurement points in one of the two divided areas.

Landscapes

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

Abstract

Alors qu'une mesure ARM de contraste est exécutée, on fait en sorte que la concentration de l'agent de contraste atteigne son sommet lorsque la composante basse fréquence est mesurée par régulation de l'ordre de mesure de l'espace k, considérant la distance depuis l'origine. Les points de mesure dans l'espace k sont divisés en deux groupes: lorsque la concentration de l'agent de contraste dans un vaisseau cible devient élevée, la mesure du premier groupe est entamée. Pour le premier groupe, les points d'échantillon sont commandés par échantillonnage à partir de la composante haute fréquence vers la composante basse fréquence de manière que la distance depuis l'origine diminue graduellement. Pour l'autre groupe à mesurer consécutivement, les points d'échantillon de l'autre côté de l'origine sont commandés par échantillonnage à partir de la composante basse fréquence vers la composante haute fréquence, de manière que la distance depuis l'origine augmente graduellement. L'influence de l'erreur de la synchronisation d'imagerie dans la mesure ARM de contraste est réduite, et tout le vaisseau est représenté dans un contraste élevé. L'image d'une artère peut être représentée sélectivement.
PCT/JP2001/000832 2000-02-07 2001-02-07 Dispositif d'imagerie par resonance magnetique WO2001058352A1 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US10/203,260 US6611144B2 (en) 2000-02-07 2001-02-07 Magnetic resonance imaging device

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2000029817A JP3847512B2 (ja) 2000-02-07 2000-02-07 磁気共鳴イメージング装置
JP2000-29817 2000-02-07
JP2000-97632 2000-03-31
JP2000097632A JP3847519B2 (ja) 2000-03-31 2000-03-31 磁気共鳴イメージング装置

Publications (1)

Publication Number Publication Date
WO2001058352A1 true WO2001058352A1 (fr) 2001-08-16

Family

ID=26585003

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2001/000832 WO2001058352A1 (fr) 2000-02-07 2001-02-07 Dispositif d'imagerie par resonance magnetique

Country Status (1)

Country Link
WO (1) WO2001058352A1 (fr)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0523314A (ja) * 1991-07-19 1993-02-02 Hitachi Medical Corp Ｍｒｉ装置
JPH10277002A (ja) * 1997-04-07 1998-10-20 Ge Yokogawa Medical Syst Ltd Ｍｒ連続撮像方法およびｍｒｉ装置
JP2000350715A (ja) * 1999-05-28 2000-12-19 Siemens Ag 磁気共鳴により時間および位置分解された三次元のデータレコードを取得する方法およびこの方法を実施するための装置
JP2001070283A (ja) * 1999-09-08 2001-03-21 Toshiba Corp Ｍｒイメージングのデータ収集方法及びｍｒｉ装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0523314A (ja) * 1991-07-19 1993-02-02 Hitachi Medical Corp Ｍｒｉ装置
JPH10277002A (ja) * 1997-04-07 1998-10-20 Ge Yokogawa Medical Syst Ltd Ｍｒ連続撮像方法およびｍｒｉ装置
JP2000350715A (ja) * 1999-05-28 2000-12-19 Siemens Ag 磁気共鳴により時間および位置分解された三次元のデータレコードを取得する方法およびこの方法を実施するための装置
JP2001070283A (ja) * 1999-09-08 2001-03-21 Toshiba Corp Ｍｒイメージングのデータ収集方法及びｍｒｉ装置

Similar Documents

Publication Publication Date Title
JP3847512B2 (ja) 磁気共鳴イメージング装置
US6760611B1 (en) Magnetic resonance imaging method and device therefor
US8121668B2 (en) System and method for providing flow-enhanced signal intensity during a functional MRI process
US9700220B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
US9341694B2 (en) Method and magnetic resonance system for distortion correction in magnetic resonance imaging
US8587306B2 (en) Magnetic resonance imaging apparatus and multi-contrast acquiring method
US8466679B2 (en) Magnetic resonance imaging apparatus and method configured for susceptibility-emphasized imaging with improved signal-to-noise ratio
JP5719968B2 (ja) Ｍｒデータを収集する方法及び装置
US9588205B2 (en) Simultaneous and dynamic determination of longitudinal and transversal relaxation times of a nuclear spin system
US10302729B2 (en) Method and magnetic resonance apparatus for speed-compensated diffusion-based diffusion imaging
JP2004526491A (ja) 拡散強調された磁気共鳴画像化データの取得方法および装置
US20040113613A1 (en) Flow imaging using balanced phase contrast steady state free precession magnetic resonance imaging
JPH0263437A (ja) ３次元核磁気共鳴造影図を作る方法と装置
US7307420B2 (en) MRI method for simultaneous phase contrast angiography and invasive device tracking
WO2011007691A1 (fr) Dispositif d'imagerie par résonance magnétique et procédé d'imagerie par résonance magnétique
US6198960B1 (en) Flip angle modulated magnetic resonance angiography
WO2004080302A1 (fr) Systeme d'imagerie par resonance magnétique
US7167740B2 (en) Measuring method in magnetic resonance imaging device and magnetic resonance imaging device
US7660618B2 (en) Magnetic resonance imaging device
JPH07327960A (ja) 磁気共鳴イメージング装置
JP2012505708A (ja) 磁気共鳴イメージング法における到達時刻マッピング法
US11327138B2 (en) Method for compensating eddy currents when creating measurement data by means of magnetic resonance
JP3847519B2 (ja) 磁気共鳴イメージング装置
WO2016170863A1 (fr) Dispositif d'imagerie par résonance magnétique et procédé d'imagerie par résonance magnétique
JP3903276B2 (ja) 磁気共鳴撮像装置

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): US

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE TR

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: 10203260

Country of ref document: US

122 Ep: pct application non-entry in european phase