WO2012005137A1 - 磁気共鳴イメージング装置及びｒｆパルス制御方法 - Google Patents

磁気共鳴イメージング装置及びｒｆパルス制御方法 Download PDF

Info

Publication number
WO2012005137A1
WO2012005137A1 PCT/JP2011/064737 JP2011064737W WO2012005137A1 WO 2012005137 A1 WO2012005137 A1 WO 2012005137A1 JP 2011064737 W JP2011064737 W JP 2011064737W WO 2012005137 A1 WO2012005137 A1 WO 2012005137A1
Authority
WO
WIPO (PCT)
Prior art keywords
sequence
pulse
flip angle
magnetic resonance
resonance imaging
Prior art date
Application number
PCT/JP2011/064737
Other languages
English (en)
French (fr)
Inventor
中西　健二
板垣　博幸
Original Assignee
株式会社 日立メディコ
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 株式会社 日立メディコ filed Critical 株式会社 日立メディコ
Priority to US13/808,659 priority Critical patent/US9594136B2/en
Priority to JP2012523823A priority patent/JP5942265B2/ja
Publication of WO2012005137A1 publication Critical patent/WO2012005137A1/ja

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
    • G01R33/543—Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5613—Generating steady state signals, e.g. low flip angle sequences [FLASH]
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems

Definitions

  • the present invention relates to a magnetic resonance imaging (hereinafter referred to as “MRI”) apparatus, and more particularly to application control of a high frequency magnetic field pulse that leads to a steady state free precession motion state in a short time.
  • MRI magnetic resonance imaging
  • MRI equipment measures nuclear magnetic resonance (NMR) signals generated by the nuclear spins that make up the body of a subject, especially the human body, and the shape and function of the head, abdomen, limbs, etc. in two or three dimensions It is a device that automatically images.
  • NMR nuclear magnetic resonance
  • the NMR signal is given different phase encoding depending on the gradient magnetic field, frequency-encoded, and measured as time series data.
  • the measured NMR signal is reconstructed into an image by two-dimensional or three-dimensional Fourier transform.
  • a pulse sequence (hereinafter referred to as an SSFP sequence) that measures an echo signal by making the magnetization a steady state free precession motion (SSFP; Steady State Free Free Precession, hereinafter referred to as a steady state). )It has been known.
  • SSFP sequence high-frequency magnetic field pulses (hereinafter referred to as RF pulses) are irradiated with a very short repetition time TR of about several ms to excite nuclear spins (hereinafter abbreviated as spins).
  • spins nuclear spins
  • the spin converges to a steady state in which excitation and relaxation are balanced after a time of about the longitudinal relaxation time (T1) has elapsed. While the echo signal measured in the steady state has a high signal-to-noise ratio (SNR; Signal to Noize Ratio), it is necessary to irradiate the RF pulse continuously until the steady state is reached.
  • SNR signal-to-noise ratio
  • a typical application is cine imaging.
  • transient states the state in the middle of the spin reaching the steady state from the thermal equilibrium state is called a transient state.
  • excitation and relaxation are not balanced, and the spin is unstable, so the intensity of the echo signal often oscillates. This vibration causes artifacts such as ghosting and blurring in the reconstructed image. Therefore, when measuring the echo signal in the transient state, a start-up sequence is executed as a preparation sequence before the data acquisition sequence for measuring the echo signal for image reconstruction, and the echo signal at the time of echo data collection is Suppresses vibration.
  • ⁇ / 2 method for example, Patent Document 1 as one method of the start-up sequence.
  • an RF pulse having a flip angle that is half the flip angle ( ⁇ ) of the RF pulse for echo data collection is irradiated before TR / 2 before the irradiation time of the first RF pulse.
  • on-resonance spin for spins having the same resonance frequency as the RF pulse irradiation frequency
  • vibration can be suppressed, but the spin having a resonance frequency different from the RF pulse irradiation frequency.
  • off-resonance spin vibration cannot be sufficiently suppressed.
  • Linear-Flip-Angle method for example, Non-Patent Document 1
  • the Linear Flip Angle method is a method of irradiating an RF pulse at a certain time interval TR while linearly increasing the flip angle as shown in the equation (1).
  • ⁇ (n) is the absolute value of the flip angle of the nth RF pulse
  • N is the total number of RF pulses in the start-up sequence
  • a ⁇ / N. That is, the flip angle ⁇ (n) specified by the equation (1) is a value on a straight line with an offset (intercept) of 0 (zero) and an inclination of a.
  • the flip angle is ⁇ (n) ⁇ ( ⁇ 1) n .
  • vibrations can be suppressed for off-resonance spins compared to the ⁇ / 2 method, but vibrations occur for on-resonance spins. This vibration can be suppressed by setting the total number of RF pulses to about several tens of times.
  • the echo signal oscillates not only for off-resonance spins but also for on-resonance spins.
  • the echo signal is transmitted with the vibration remaining. Must be measured.
  • artifacts such as ghosting and blurring may occur in the image.
  • the present invention has been made in view of the above problems, and even if the total number of RF pulses of the start-up sequence leading to a steady state is reduced, the vibration of the echo signal is reduced, and a ghost or blur is added to the image.
  • the purpose is to prevent artifacts such as rings from occurring.
  • the RF pulse sequence of the start-up sequence has a flip angle that increases monotonously with an offset.
  • the sum between two adjacent terms in a monotonically increasing number sequence is the flip angle.
  • the MRI apparatus of the present invention includes a start-up sequence that suppresses spin oscillation in a transient state by an RF pulse sequence in which the flip angle monotonously increases, and an SSFP sequence that measures an echo signal in a steady state.
  • a calculation processing unit that generates an imaging sequence and a measurement control unit that controls measurement of an echo signal from the subject based on the imaging sequence, and the RF pulse sequence has a flip angle that monotonously increases with an offset It is characterized by being.
  • the RF pulse control method of the present invention includes the steps of setting the number of RF pulses in the RF pulse sequence and the flip angle of the RF pulse in the SSFP sequence, and the number of RF pulses in the set RF pulse sequence and the RF in the SSFP sequence. And a flip angle calculation step for obtaining an RF pulse sequence by obtaining a flip angle that monotonously increases with an offset based on the flip angle of the pulse.
  • the vibration of the echo signal is reduced, and ghost, blurring, etc. Artifacts do not occur, and high-quality images can be acquired.
  • FIG. 6 is a diagram illustrating a monotonically increasing number sequence (a) of the first embodiment and a flip angle (b) of an RF pulse sequence based on the monotonically increasing number sequence.
  • the figure which shows the flip angle of the RF pulse sequence of Start-up sequence when N is 5 times.
  • the figure which shows the behavior of the spin of on resonance. (a) shows the case where the Start-up sequence of FIG. 3 is used, and (b) shows the case where the conventional Linear Flip Angle method is used.
  • FIG. 3 is a diagram illustrating functional blocks of an arithmetic processing unit according to the first embodiment.
  • 3 is a flowchart illustrating an operation flow of functional blocks of an arithmetic processing unit according to the first embodiment. The figure which shows an example of the input screen of imaging conditions.
  • FIG. 10 is a diagram showing a monotonically increasing number sequence (a) of Example 2 and a flip angle (b) of an RF pulse sequence based on the monotonically increasing number sequence.
  • FIG. 11 is a diagram showing a result of numerical calculation of vibrations in a transient state for off-resonance spins when the start-up sequence of FIG. 10 is used.
  • FIG. 1 is a block diagram showing the overall configuration of an embodiment of an MRI apparatus according to the present invention.
  • This MRI apparatus uses a NMR phenomenon to obtain a tomographic image of a subject 101.
  • a static magnetic field generating magnet 102, a gradient magnetic field coil 103, a gradient magnetic field power supply 109, and an RF transmission coil 104, an RF transmitter 110, an RF receiver coil 105, a signal detector 106, a signal processor 107, a measurement controller 111, an overall controller 108, a display / operation unit 113, and a subject 101 are mounted.
  • a bed 112 for taking the top plate into and out of the static magnetic field generating magnet 102.
  • the static magnetic field generating magnet 102 generates a uniform static magnetic field in the direction perpendicular to the body axis of the subject 101 in the vertical magnetic field method and in the body axis direction in the horizontal magnetic field method.
  • a permanent magnet type, normal conducting type or superconducting type static magnetic field generating source is arranged around the.
  • the gradient magnetic field coil 103 is a coil wound in the three-axis directions of X, Y, and Z that are the real space coordinate system (stationary coordinate system) of the MRI apparatus, and each gradient magnetic field coil is a gradient magnetic field that drives it.
  • a current is supplied to the power source 109.
  • the gradient magnetic field power supply 109 of each gradient coil is driven according to a command from the measurement control unit 111 described later, and supplies a current to each gradient coil.
  • gradient magnetic fields Gx, Gy, and Gz are generated in the three-axis directions of X, Y, and Z.
  • a slice gradient magnetic field pulse (Gs) is applied in a direction orthogonal to the slice plane (imaging cross section) to set a slice plane for the subject 101, orthogonal to the slice plane and orthogonal to each other.
  • Phase encoding gradient magnetic field pulse (Gp) and frequency encoding (leadout) gradient magnetic field pulse (Gf) are applied in the remaining two directions, and position information in each direction is encoded in the NMR signal (echo signal). .
  • the RF transmission coil 104 is a coil that irradiates the subject 101 with an RF pulse, and is connected to the RF transmission unit 110 and supplied with a high-frequency pulse current. As a result, an NMR phenomenon is induced in the spins of atoms constituting the living tissue of the subject 101.
  • the RF transmission unit 110 is driven in accordance with a command from the measurement control unit 111 described later, and the RF transmission coil 104 is arranged in proximity to the subject 101 after the high frequency pulse is amplitude-modulated and amplified. , The subject 101 is irradiated with an RF pulse.
  • the RF receiving coil 105 is a coil that receives an echo signal emitted by the NMR phenomenon of spin that constitutes the living tissue of the subject 101.
  • the received echo signal is connected to the signal detecting unit 106 and is received by the signal detecting unit 106. Sent.
  • the signal detection unit 106 performs detection processing of the echo signal received by the RF receiving coil 105. Specifically, the echo signal of the response of the subject 101 induced by the RF pulse irradiated from the RF transmission coil 104 is received by the RF receiving coil 105 disposed in the vicinity of the subject 101, and measurement control described later is performed. In accordance with a command from the unit 111, the signal detection unit 106 amplifies the received echo signal, divides the signal into two orthogonal signals by quadrature detection, and samples each by a predetermined number (for example, 128, 256, 512, etc.) Each sampling signal is A / D converted into a digital quantity and sent to a signal processing unit 107 described later. Therefore, the echo signal is obtained as time-series digital data (hereinafter referred to as echo data) composed of a predetermined number of sampling data.
  • echo data time-series digital data
  • the signal processing unit 107 performs various processes on the echo data, and sends the processed echo data to the measurement control unit 111.
  • the measurement control unit 111 mainly transmits various commands for collecting echo data necessary for reconstruction of the tomographic image of the subject 101 to the gradient magnetic field power source 109, the RF transmission unit 110, and the signal detection unit 106. And a control unit for controlling them. Specifically, the measurement control unit 111 operates under the control of the overall control unit 108 described later, and controls the gradient magnetic field power source 109, the RF transmission unit 110, and the signal detection unit 106 based on a predetermined pulse sequence. The echo necessary for reconstructing the image of the imaging region of the subject 101 is repeatedly executed by irradiating the subject 101 with the RF pulse and applying the gradient magnetic field pulse and detecting the echo signal from the subject 101. Control data collection.
  • the application amount of the phase encoding gradient magnetic field is changed in the case of two-dimensional imaging, and the application amount of the slice encoding gradient magnetic field is further changed in the case of three-dimensional imaging.
  • Values such as 128, 256, and 512 are normally selected as the number of phase encodings, and values such as 16, 32, and 64 are normally selected as the number of slice encodings.
  • the overall control unit 108 controls the measurement control unit 111 and controls various data processing and processing result display and storage, and includes an arithmetic processing unit 114 having a CPU and a memory, an optical disc, And a storage unit 115 such as a magnetic disk.
  • the measurement control unit 111 is controlled to execute the collection of echo data, and when the echo data is input from the measurement control unit 111, the arithmetic processing unit 114 converts the encoded information applied to the echo data. Based on this, it is stored in an area corresponding to the K space in the memory.
  • a group of echo data stored in an area corresponding to the K space in the memory is also referred to as K space data.
  • the arithmetic processing unit 114 performs processing such as signal processing and image reconstruction by Fourier transform on the K space data, and displays the resulting image of the subject 101 on the display / operation unit 113 described later. And is recorded in the storage unit 115.
  • the display / operation unit 113 includes a display unit for displaying the reconstructed image of the subject 101, a trackball or a mouse and a keyboard for inputting various control information of the MRI apparatus and control information for processing performed by the overall control unit 108. Etc., and an operation unit.
  • the operation unit is disposed in the vicinity of the display unit, and an operator interactively controls various processes of the MRI apparatus through the operation unit while looking at the display unit.
  • the radionuclide to be imaged by the MRI apparatus is a hydrogen nucleus (proton) which is the main constituent material of the subject, as is widely used in clinical practice.
  • proton the main constituent material of the subject
  • the form or function of the human head, abdomen, limbs, etc. is imaged two-dimensionally or three-dimensionally.
  • each flip angle of the RF pulse sequence constituting the Start-up sequence is controlled so that the offset (intercept) becomes a value on a straight line other than 0 (zero).
  • this embodiment will be described in detail with reference to FIGS.
  • ⁇ represents the flip angle of the SSFP sequence at the time of echo data collection
  • N represents the number of RF pulses in the start-up sequence.
  • ⁇ is the limit that the longitudinal relaxation time T1 and the transverse relaxation time T2 are sufficiently longer than the repetition time TR, and the relaxation process can be ignored.
  • the on resonance This means the absolute value of the angle that the spin makes with the Mz axis.
  • the absolute value ⁇ (n) of the flip angle of the nth RF pulse is calculated from ⁇ (n) according to the following equation (3).
  • the flip angle of the RF pulse sequence of the start-up sequence is ⁇ / 10, ⁇ 3 ⁇ / 10, 5 ⁇ / 10, ⁇ 7 ⁇ / 10, 9 ⁇ / 10 as shown in FIG. It becomes.
  • ⁇ 90 deg., 9 deg., ⁇ 27 deg., 45 deg., ⁇ 63 deg., 81 deg.
  • the flip angle of the conventional Linear Flip Angle method expressed by Equation (1) is 18 deg., ⁇ 36 deg., 54 deg., ⁇ 72 geg., 90 deg.
  • the gradient magnetic field pulse waveform is omitted.
  • the pulse sequence described in Patent Document 2 can be used, and thus detailed description thereof is omitted here.
  • Fig. 4 shows the behavior of on-resonance spin when the Start-up sequence is used.
  • (a) shows the case where the start-up sequence of this embodiment is used, and for comparison,
  • (b) shows the case where the conventional Linear ⁇ Flip Angle method is used.
  • the spin behavior on the Mx-Mz plane when executing the Start-up sequence is shown, and the numerical value in the figure represents the number of times of RF pulse irradiation.
  • the flip angle of the RF pulse in this embodiment shown in the figure is 9 deg., -27 deg., 45 deg., -63 deg., 81 deg., and in the case of the conventional Linear ⁇ Flip Angle method shown in (b)
  • the flip angle of the RF pulse is 18 deg., ⁇ 36 deg., 54 deg., ⁇ 72 geg., 90 deg.
  • the on-resonance spin applied with the start-up sequence of this example moves alternately with the Mz axis every time an RF pulse is applied, but at the end of the start-up sequence, the Mz axis Is exactly ⁇ / 2, so that the vibration of the echo signal can be suppressed.
  • the on-resonance spin to which the conventional Linear-Flip-Angle start-up sequence is applied moves alternately across the Mz axis each time an RF pulse is applied, and at the end of the start-up sequence, 401 and Since it vibrates between 402, the echo signal will also vibrate.
  • Fig. 5 shows the numerical calculation results of the vibration in the transient state for the off-resonance spin.
  • the spin whose vertical axis is zero is defined as “on resonance”, and the other spins are defined as “off resonance” (the same applies to FIG. 11 described later).
  • (a) shows the case where the start-up sequence of this embodiment is used, and for comparison,
  • (b) shows the case where the conventional Linear ⁇ Flip Angle method is used.
  • the horizontal axis represents the number of times of RF pulse irradiation
  • the vertical axis represents the offset frequency between repetition times TR
  • the absolute value of the transverse magnetization is displayed in a grace case.
  • Fig. 6 shows the result of numerical calculation of the flip angle dependence of the vibration in the transient state for the on-resonance spin using the Start-up sequence.
  • (a) shows the case where the start-up sequence of this embodiment is used, and for comparison,
  • (c) shows the case where the conventional LinearLineFlip Angle method is used.
  • the horizontal axis indicates the number of RF pulses
  • the vertical axis indicates the flip angle
  • the absolute value of the transverse magnetization is displayed in gray scale as in FIG.
  • the parameters used for the numerical calculation are the same as in FIG.
  • Figures (b) and (d) show the absolute values of transverse magnetization at the number of RF pulses for spins with flip angles FA of 45, 90, and 135 degrees, respectively, from figures (a) and (c).
  • FIG. 5B shows that no vibration occurs in the start-up sequence of this embodiment. This tendency is the same for an arbitrary flip angle.
  • the vibration of the echo signal can be suppressed at an arbitrary flip angle, but the conventional Linear Flip Angle In the method, it can be understood that the echo signal vibrates at any flip angle. Therefore, by using the start-up sequence of the present embodiment, it is possible to suppress vibration even when a flip angle shifts due to uneven irradiation.
  • FIG. 6 shows the numerical calculation results for the on-resonance spin, but the off-resonance spin is similar to that in FIG. It is possible to sufficiently suppress the vibration with respect to the spin at.
  • the first term of ⁇ (n) is 0, but vibration can be suppressed even with an arbitrary value.
  • FIG. 7 shows each function of the arithmetic processing unit 114 that obtains the flip angle of each RF pulse constituting the start-up sequence of the above-described embodiment and starts the SSFP sequence including the start-up sequence.
  • the arithmetic processing unit 114 of the present embodiment includes an imaging condition setting unit 701, an RF pulse sequence setting unit 702, a pulse sequence setting unit 703, and an imaging control unit 704.
  • the imaging condition setting unit 701 is connected to the display / operation unit 113 and the storage unit 115, and causes the display unit to display an input screen for accepting input of imaging condition setting (ie, imaging parameter value).
  • imaging condition setting ie, imaging parameter value
  • FIG. 9 shows an example of the input screen.
  • This input screen includes at least the flip angle ⁇ (901) in the SSFP sequence that is the data acquisition sequence, and the RF pulse number N (902) in the RF pulse sequence in the start-up sequence that is the preparation sequence before collecting the echo data.
  • a region for receiving input of the value is displayed for each imaging condition including.
  • the operator inputs the value of each imaging condition via the operation unit.
  • the imaging condition setting unit 701 stores various imaging condition values input by the operator in the memory as the various imaging condition setting values.
  • predetermined values stored in the storage unit 115 may be used.
  • the RF pulse sequence setting unit 702 is connected to the imaging condition setting unit 701 and calculates and sets the flip angle of the RF pulse sequence of the start-up sequence based on the input imaging conditions.
  • the RF pulse sequence setting unit 702 uses the values of the flip angle ⁇ and the number of RF pulses N input by the operator and stored in the memory or read from the storage unit 115, based on the formula (2).
  • An equality number sequence ⁇ (n) is created, and the absolute value ⁇ (n) of the flip angle of the nth RF pulse is calculated based on the equation (3) using the generated difference number sequence ⁇ (n).
  • the absolute value ⁇ (n) of the flip angle of the nth RF pulse may be directly calculated based on the equation (4).
  • the pulse sequence generation unit 703 is connected to the RF pulse sequence setting unit 702 and the storage unit 115, and sets application timings of RF pulses and gradient magnetic field pulses based on imaging conditions input by the operator and stored in the memory.
  • a pal sequence is generated for specific data. Since the pulse sequence of the present embodiment includes a start-up sequence as a preparation sequence and an SSFP sequence as a data collection sequence, the pulse sequence setting unit 703 continuously executes these two pulse sequences. Specific data for each is generated. Specifically, the start-up sequence data specifically defines the RF pulse sequence of the flip angle represented by the equation (5) in which the irradiation phase is added to the flip angle set by the RF pulse sequence setting unit 702. This is data on RF pulses and gradient magnetic field pulses.
  • the SSFP sequence data is data on an RF pulse and a gradient magnetic field pulse that define the measurement of an echo signal for image reconstruction by irradiating an RF pulse with a flip angle ⁇ .
  • the imaging control unit 704 is connected to the pulse sequence generation unit 703, the storage unit 115, and the measurement control unit 111, notifies the measurement control unit 111 of specific data of the pulse sequence generated by the pulse sequence generation unit 703, and performs measurement control.
  • the unit 111 is caused to execute the pulse sequence to start imaging.
  • the imaging condition setting unit 701 displays an input screen for accepting an input of imaging condition setting on the display unit.
  • the operator inputs and sets imaging conditions via this input screen.
  • the imaging condition setting unit 701 reads the imaging condition value stored in advance in the storage unit 115.
  • the imaging condition setting unit 701 stores the values of various imaging conditions that are input or read in the memory of the arithmetic processing unit 114.
  • step 802 the RF pulse sequence setting unit 702 calculates the flip angle of the RF pulse sequence of the start-up sequence based on the value of the imaging condition stored in the memory in step 801, and stores it in the memory of the arithmetic processing unit 114.
  • Step 803 the pulse sequence generation unit 703, based on the imaging condition value stored in the memory in Step 801 and the flip angle of the RF pulse sequence of the Start-up sequence stored in the memory in Step 802, Specific data for defining the flip angle of the RF pulse, the gradient magnetic field pulse shape, the application timing thereof, and the like are obtained, and specific data for continuously executing the Start-up sequence and the SSFP sequence is generated.
  • step 804 the imaging control unit 704 notifies the measurement control unit 111 of specific data of the pulse sequence generated by the pulse sequence generation unit 703 in step 803, and causes the measurement control unit 111 to execute the generated pulse sequence. To start imaging.
  • the MRI apparatus and the RF pulse control method of the present embodiment linearly increase the flip angle of each RF pulse constituting the RF pulse sequence of the start-up sequence with an appropriate offset. Thereby, it is possible to suppress vibration of the on-resonance spin. Therefore, even if the total number of start-up sequence RF pulses that lead to a steady state is reduced, the vibration of the echo signal is reduced, and artifacts such as ghosting and blurring do not occur in the image, and a high-quality image is acquired. be able to. That is, echo data collection can be started with a small number of idle shots.
  • a start-up sequence is performed using a number sequence that smoothly changes at the rising edge and the contact point with the ⁇ / 2 line, that is, a number sequence in which the difference between adjacent two terms monotonously increases and then decreases monotonically. Determine each flip angle of the RF pulse sequence.
  • this embodiment will be described in detail with reference to FIGS.
  • the number sequence ⁇ is created based on the straight line connecting the two points 0 and ⁇ / 2.
  • it may be a monotonically increasing number sequence. Therefore, as a result of numerical calculation by changing the monotonically increasing number sequence variously, the number sequence that smoothly changes at the rise and the contact point with the ⁇ / 2 line, that is, the difference between adjacent binomials increases monotonously after the monotonic increase.
  • the inventor has discovered that by using the decreasing number sequence to determine each flip angle of the RF pulse sequence of the Start-up sequence, the oscillation of the echo signal is suppressed for a wider offset frequency.
  • ⁇ (n) ⁇ / 4 ⁇ (1-cos (n ⁇ / N)) (6)
  • This number sequence is shown in FIG.
  • ⁇ (n) ⁇ / 2 ⁇ (1- ⁇ ⁇ cos (n ⁇ / N ⁇ ))
  • n 1, 2,. , N (7)
  • ⁇ 2 (1 + cos ( ⁇ / N)) / 2
  • sin ( ⁇ ) sin ( ⁇ / N) / 2 ⁇
  • cos ( ⁇ ) (1 + cos ( ⁇ / N)) / 2 ⁇ .
  • FIG. 11 shows the result of numerical calculation of the vibration of the echo signal in the transient state for the off-resonance spin (in FIG. 11, the vertical axis is a spin other than zero).
  • the coordinate axes and the method of numerical calculation are the same as in FIG. Comparing FIG. 5 (a) with FIG. 11, it can be seen that FIG. 11 can suppress the vibration with respect to a wider frequency offset. Note that FIG. 11 shows that vibration is also suppressed for spins of on-resonance (spins whose vertical axis is zero in FIG. 11).
  • each function of the arithmetic processing unit 114 that obtains the flip angle of each RF pulse in the start-up sequence of the above-described embodiment and starts the SSFP sequence including the start-up sequence will be described.
  • Each function of the arithmetic processing unit 114 in the present embodiment is the same as that of the first embodiment shown in FIG. 7, but the processing contents of the RF pulse sequence setting unit 702 are different.
  • the processing contents of the RF pulse sequence setting unit 702 are different.
  • only different portions will be described, and description of the same portions will be omitted.
  • the RF pulse sequence setting unit 702 of the arithmetic processing unit 114 of the present embodiment creates an arithmetic sequence ⁇ (n) based on the equation (6), and uses the generated arithmetic sequence ⁇ (n) (3 ) To calculate the absolute value ⁇ (n) of the flip angle of the nth RF pulse.
  • the absolute value ⁇ (n) of the flip angle of the nth RF pulse may be directly calculated based on the equation (7).
  • processing flow of this embodiment which is performed in cooperation with the functional units of the arithmetic processing unit 114 of this embodiment, is the same as the processing flow of the first embodiment shown in FIG.
  • the MRI apparatus and the RF pulse control method of the present embodiment use the rising angle and the flip angle calculated from the equation (3) based on the numerical sequence that smoothly changes at the point of contact with the ⁇ / 2 line.
  • a start-up sequence comprising a plurality of RF pulse sequences is used.
  • 1 subject 2 static magnetic field generation system, 3 gradient magnetic field generation system, 4 sequencer, 5 transmission system, 6 reception system, 7 signal processing system, 8 central processing unit (CPU), 9 gradient magnetic field coil, 10 gradient magnetic field power supply, 11 High frequency transmitter, 12 modulator, 13 high frequency amplifier, 14a high frequency coil (transmitting coil), 14b high frequency coil (receiving coil), 15 signal amplifier, 16 quadrature phase detector, 17 A / D converter, 18 magnetic disk, 19 optical disc, 20 display, 21 ROM, 22 RAM, 23 trackball or mouse, 24 keyboard

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

　過渡状態での計測において、ゴーストやブラーリング等のアーチファクトを低減するために、Start-upシーケンスのRFパルス系列を、オフセットを有して単調増加するフリップ角とする。例えば、単調増加数列の隣り合う二項間の和をフリップ角とする。具体的には、RFパルス系列のRFパルス数と、SSFPシーケンスにおけるRFパルスのフリップ角とを設定し、設定されたRFパルス系列のRFパルス数とSSFPシーケンスにおけるRFパルスのフリップ角とに基づいて、オフセットを有して単調増加するフリップ角を求めてRFパルス系列とする。

Description

磁気共鳴イメージング装置及びＲＦパルス制御方法
　本発明は、磁気共鳴イメージング(以下、「MRI」という)装置に関し、特に、短時間で定常状態自由歳差運動状態に導く高周波磁場パルスの印加制御に関する。
　MRI装置は、被検体、特に人体の組織を構成する原子核スピンが発生する核磁気共鳴(NMR)信号を計測し、その頭部、腹部、四肢等の形態や機能を2次元的に或いは3次元的に画像化する装置である。撮影においては、NMR信号には、傾斜磁場によって異なる位相エンコードが付与されるとともに周波数エンコードされて、時系列データとして計測される。計測されたNMR信号は、2次元又は3次元フーリエ変換されることにより画像に再構成される。
　上記MRI装置において、代表的な高速撮像法として、磁化を定常状態自由歳差運動(SSFP；Steady State Free Precession、以下、定常状態という)にしてエコー信号を計測するパルスシーケンス(以下、SSFPシーケンスという)が知られている。SSFPシーケンスでは、数ms程度の非常に短い繰り返し時間TRで高周波磁場パルス(以下、RFパルスという)を照射して原子核スピン(以下、スピンと省略)を励起する。また、励起されたスピンは、縦緩和と横緩和により、もとの熱平衡状態に回復しようとする。この繰り返しにより、スピンは、縦緩和時間(T1)程度の時間が経過した後に、励起と緩和がバランスした定常状態に収束する。定常状態で計測されたエコー信号は高い信号雑音比(SNR ；Signal to Noize Ratio)となる一方で、定常状態に至るまでの期間は連続的にRFパルスを照射することが必要であるため、心臓シネ撮像などへの適用が代表的である。
　これに対して、スピンが熱平衡状態から定常状態に到達するまでの途中の状態は、過渡状態と呼ばれる。過渡状態では、励起と緩和がバランスしておらず、スピンが不安定であることから、エコー信号の強度はしばしば振動する。この振動は、再構成画像においてゴーストやブラーリング等のアーチファクトの原因となる。そこで、過渡状態でエコー信号の計測を行う際には、画像再構成用のエコー信号を計測するデータ収集シーケンスの前に準備シーケンスとしてのStart-upシーケンスを実行し、エコーデータ収集時のエコー信号の振動を抑制する。
　Start-upシーケンスの一つの方法として、α/2法(例えば特許文献1)がある。この方法では、エコーデータ収集のRFパルスのフリップ角(α)の半分のフリップ角のRFパルスを、第一回目のRFパルスの照射時間よりもTR/2前に照射する。α/2法では、RFパルスの照射周波数と同じ共鳴周波数を持つスピン(以下、オンレゾナンスのスピン)については、振動を抑制することができるが、RFパルスの照射周波数と異なる共鳴周波数を持つスピン(以下、オフレゾナンスのスピン)については、十分に振動を抑制することができない。
　Start-upシーケンスの別の方法として、オフレゾナンスのスピンの振動をより効果的に抑制できるLinear Flip Angle法(例えば非特許文献1)が挙げられる。Linear Flip Angle法は、一定の時間間隔TRで、(1)式に示す様に、フリップ角を線形に増加させながらRFパルスを照射する方法である。
　　φ(n)＝a×n(n＝1,2,…,N)　　　　　(1)
　ここで、φ(n)はn番目のRFパルスのフリップ角の絶対値、NはStart-upシーケンスのRFパルスの総数であり、a＝α/Nである。即ち、(1)式で指定されるフリップ角φ(n)は、オフセット(切片)が0(ゼロ)で傾きがaの直線上の値となる。実際には、RFパルスは照射位相を180deg.ずらしながら照射するため、フリップ角はφ(n)×(－1)nとなる。Linear Flip Angle法では、オフレゾナンスのスピンについては、α/2法よりも振動を抑制可能であるが、オンレゾナンスのスピンについては振動が生じる。この振動は、RFパルスの総数を数10回程度とすれば抑制可能である。
特開平8-56932号公報 特開2004-329268号公報
V. S. Deshpande, et al., ISMRM 10, p1586 (2002)
　Linear Flip Angle法では、Start-upシーケンスのRFパルスの総数が少ない場合、オフレゾナンスのスピンだけではなく、オンレゾナンスのスピンについてもエコー信号が振動してしまう。特に、セントリックオーダーのK空間走査のように、数回程度の空打ちの直後から画像再構成に用いるエコー信号の計測を開始しなければならない場合においては、振動が残存した状態でエコー信号を計測しなければならない。その結果、画像にゴーストやブラーリング等のアーチファクトが生じる可能性がある。
　そこで、本発明は、上記課題を鑑みてなされたものであり、定常状態に導くStart-upシーケンスのRFパルスの総数を少なくしても、エコー信号の振動を低減して、画像にゴーストやブラーリング等のアーチファクトが生じないようにすることを目的とする。
　上記目的を達成するために、本発明は、Start-upシーケンスのRFパルス系列を、オフセットを有して単調増加するフリップ角とする。例えば、単調増加数列の隣り合う二項間の和をフリップ角とする。
　具体的には、本発明のMRI装置は、フリップ角が単調増加するRFパルス系列により過渡状態でのスピンの振動を抑制するStart-upシーケンスと、定常状態でエコー信号を計測するSSFPシーケンスとからなる撮像シーケンスを生成する演算処理部と、撮像シーケンスに基づいて被検体からのエコー信号の計測を制御する計測制御部と、を備え、RFパルス系列は、オフセットを有して単調増加するフリップ角であることを特徴とする。
　また、本発明のRFパルス制御方法は、RFパルス系列のRFパルス数と、SSFPシーケンスにおけるRFパルスのフリップ角とを設定するステップと、設定されたRFパルス系列のRFパルス数とSSFPシーケンスにおけるRFパルスのフリップ角とに基づいて、オフセットを有して単調増加するフリップ角を求めてRFパルス系列とするフリップ角算出ステップと、を有することを特徴とする。
　本発明のMRI装置及びRFパルス制御方法によれば、定常状態に導くStart-upシーケンスのRFパルスの総数を少なくしても、エコー信号の振動を低減して、画像にゴーストやブラーリング等のアーチファクトが生じず、高画質の画像を取得することができる。
本発明が適用されるMRI装置の構成の概要を示すブロック図。 実施例1の単調増加数列(a)と、この単調増加数列に基づくRFパルス系列のフリップ角(b)を示す図。 Nが5回の場合、Start-upシーケンスのRFパルス系列のフリップ角を示す図。 オンレゾナンスのスピンの振舞を示す図。(a)は図3のStart-upシーケンスを用いた場合を示し、(b)は従来のLinear Flip Angle法を用いた場合を示す。 オフレゾナンスのスピンの振舞を示す図。(a)は図3のStart-upシーケンスを用いた場合を示し、(b)は従来のLinear Flip Angle法を用いた場合を示す。 Start-upシーケンスを用いた場合のオンレゾナンススピンについて、過渡状態での振動のフリップ角依存性を数値計算した結果を示す図。(a)(b)は図3のStart-upシーケンスを用いた場合を示し、(c)(d)は従来のLinear Flip Angle法を用いた場合を示す。 実施例1の演算処理部の機能ブロックを示す図。 実施例1の演算処理部の機能ブロックの動作フローを示すフローチャート。 撮像条件の入力画面の一例を示す図。 実施例2の単調増加数列(a)と、この単調増加数列に基づくRFパルス系列のフリップ角(b)を示す図。 図10のStart-upシーケンスを用いた場合のオフレゾナンスのスピンについて過渡状態での振動を数値計算した結果を示す図。
　以下、添付図面に従って本発明のMRI装置の好ましい実施例について詳説する。なお、発明の実施例を説明するための全図において、同一機能を有するものは同一符号を付け、その繰り返しの説明は省略する。
　最初に、本発明に係るMRI装置を図1に基づいて説明する。図1は、本発明に係るMRI装置の一実施例の全体構成を示すブロック図である。
　このMRI装置は、NMR現象を利用して被検体101の断層画像を得るもので、図1に示すように、静磁場発生磁石102と、傾斜磁場コイル103及び傾斜磁場電源109と、RF送信コイル104及びRF送信部110と、RF受信コイル105及び信号検出部106と、信号処理部107と、計測制御部111と、全体制御部108と、表示・操作部113と、被検体101を搭載する天板を静磁場発生磁石102の内部に出し入れするベッド112と、を備えて構成される。
　静磁場発生磁石102は、垂直磁場方式であれば被検体101の体軸と直交する方向に、水平磁場方式であれば体軸方向に、それぞれ均一な静磁場を発生させるもので、被検体101の周りに永久磁石方式、常電導方式あるいは超電導方式の静磁場発生源が配置されている。
　傾斜磁場コイル103は、MRI装置の実空間座標系(静止座標系)であるX、Y、Zの3軸方向に巻かれたコイルであり、それぞれの傾斜磁場コイルは、それを駆動する傾斜磁場電源109に接続され電流が供給される。具体的には、各傾斜磁場コイルの傾斜磁場電源109は、それぞれ後述の計測制御部111からの命令に従って駆動されて、それぞれの傾斜磁場コイルに電流を供給する。これにより、X、Y、Zの3軸方向に傾斜磁場Gx、Gy、Gzが発生する。
2次元スライス面の撮像時には、スライス面(撮像断面)に直交する方向にスライス傾斜磁場パルス(Gs)が印加されて被検体101に対するスライス面が設定され、そのスライス面に直交して且つ互いに直交する残りの2つの方向に位相エンコード傾斜磁場パルス(Gp)と周波数エンコード(リードアウト)傾斜磁場パルス(Gf)が印加されて、NMR信号(エコー信号)にそれぞれの方向の位置情報がエンコードされる。
　RF送信コイル104は、被検体101にRFパルスを照射するコイルであり、RF送信部110に接続され高周波パルス電流が供給される。これにより、被検体101の生体組織を構成する原子のスピンにNMR現象が誘起される。具体的には、RF送信部110が、後述の計測制御部111からの命令に従って駆動されて、高周波パルスが振幅変調され、増幅された後に被検体101に近接して配置されたRF送信コイル104に供給されることにより、RFパルスが被検体101に照射される。
　RF受信コイル105は、被検体101の生体組織を構成するスピンのNMR現象により放出されるエコー信号を受信するコイルであり、信号検出部106に接続されて受信したエコー信号が信号検出部106に送られる。
　信号検出部106は、RF受信コイル105で受信されたエコー信号の検出処理を行う。具体的には、RF送信コイル104から照射されたRFパルスによって誘起された被検体101の応答のエコー信号が被検体101に近接して配置されたRF受信コイル105で受信され、後述の計測制御部111からの命令に従って、信号検出部106が、受信されたエコー信号を増幅し、直交位相検波により直交する二系統の信号に分割し、それぞれを所定数(例えば128、256、512等)サンプリングし、各サンプリング信号をA/D変換してディジタル量に変換し、後述の信号処理部107に送る。 従って、エコー信号は所定数のサンプリングデータからなる時系列のデジタルデータ(以下、エコーデータという)として得られる。
　信号処理部107は、エコーデータに対して各種処理を行い、処理したエコーデータを計測制御部111に送る。
　計測制御部111は、被検体101の断層画像の再構成に必要なエコーデータ収集のための種々の命令を、主に、傾斜磁場電源109と、RF送信部110と、信号検出部106に送信してこれらを制御する制御部である。具体的には、計測制御部111は、後述する全体制御部108の制御で動作し、ある所定のパルスシーケンスに基づいて、傾斜磁場電源109、RF送信部110及び信号検出部106を制御して、被検体101へのRFパルスの照射及び傾斜磁場パルスの印加と、被検体101からのエコー信号の検出と、を繰り返し実行し、被検体101の撮像領域についての画像の再構成に必要なエコーデータの収集を制御する。繰り返しの際には、2次元撮像の場合には位相エンコード傾斜磁場の印加量を、3次元撮像の場合には更にスライスエンコード傾斜磁場の印加量も、変えて行なう。位相エンコードの数は通常1枚の画像あたり128、256、512等の値が選ばれ、スライスエンコードの数は、通常16,32,64等の値が選ばれる。これらの制御により信号処理部107からのエコーデータを全体制御部108に出力する。
　全体制御部108は、計測制御部111の制御、及び、各種データ処理と処理結果の表示及び保存等の制御を行うものであって、CPU及びメモリを内部に有する演算処理部114と、光ディスク、磁気ディスク等の記憶部115とを有して成る。具体的には、計測制御部111を制御してエコーデータの収集を実行させ、計測制御部111からのエコーデータが入力されると、演算処理部114がそのエコーデータに印加されたエンコード情報に基づいて、メモリ内のK空間に相当する領域に記憶させる。メモリ内のK空間に相当する領域に記憶されたエコーデータ群をK空間データともいう。そして演算処理部114は、このK空間データに対して信号処理やフーリエ変換による画像再構成等の処理を実行し、その結果である被検体101の画像を、後述の表示・操作部113に表示させると共に記憶部115に記録させる。
　表示・操作部113は、再構成された被検体101の画像を表示する表示部と、MRI装置の各種制御情報や上記全体制御部108で行う処理の制御情報を入力するトラックボール又はマウス及びキーボード等の操作部と、から成る。この操作部は表示部に近接して配置され、操作者が表示部を見ながら操作部を介してインタラクティブにMRI装置の各種処理を制御する。
　現在MRI装置の撮像対象核種は、臨床で普及しているものとしては、被検体の主たる構成物質である水素原子核(プロトン)である。プロトン密度の空間分布や、励起状態の緩和時間の空間分布に関する情報を画像化することで、人体頭部、腹部、四肢等の形態または、機能を2次元もしくは3次元的に撮像する。
　次に、本発明のMRI装置及びRFパルス制御方法の実施例1について説明する。本実施例は、Start-upシーケンスを構成するRFパルス系列の各フリップ角を、オフセット(切片)が0(ゼロ)以外の直線上の値となるように、制御する。以下、図2～図9に基づいて本実施例を詳細に説明する。
　最初に、Start-upシーケンスにおける各RFパルスのフリップ角についての本実施例の制御の概要を説明する。
　図2のように、0とα/2の間を線形に増加するような数列Θ(n)を作成する。つまり、初項0、公差α/(2N)の等差数列を作成する。
　　　Θ(n)＝(α/(2N))×n　(n＝1, 2, …N)　　　　(2)
　ここで、αはエコーデータ収集時のSSFPシーケンスのフリップ角を、NはStart-upシーケンスのRFパルス数を表す。例えば、Nが5回の場合は、α/10、2α/10、3α/10、4α/10、5α/10(＝α/2)となる。このΘ(n)は、縦緩和時間T1と横緩和時間T2が繰り返し時間TRに較べて十分長く、緩和過程が無視できるような極限で、n番目のRFパルスを照射した後で、オンレゾナンスのスピンがMz軸とのなす角の絶対値を意味する。
　そこで、以下の(3)式に従い、Θ(n)からn番目のRFパルスのフリップ角の絶対値φ(n)を算出する。
　　　φ(n)＝Θ(1)　　　　　　 n＝1のとき、
　　　φ(n)＝Θ(n)＋Θ(n－1)　　2≦n≦Nのとき　　　　　(3)
　(2)式を(3)式に代入して解くと、下記(4)式のようになる。
　　　φ(n)＝a×n＋b(n＝1, 2, …, N)　　　(4)
　ここで、傾き：a＝α/N、オフセット(切片)：b＝－α/(2N)である。(4)式で表される直線の例としてN＝5の場合を図2(b)に示す。つまり、本実施例のStart-upシーケンスを構成するRFパルス系列の各RFパルスのフリップ角は、オフセット(切片)が0(ゼロ)でない直線上の値となり、この直線に沿ってフリップ角がRFパルスの照射毎に単調に増加する。フリップ角が単調に増加するという点では従来のLinear Flip Angle法と同様であるが、従来のLinear Flip Angle法でのフリップ角の絶対値は、(1)式のようにオフセットがないのに対して、本実施例のフリップ角制御では(4)式のようにオフセットbが存在する点が異なる。
　実際には、照射位相を180deg.ずらしながら各RFパルスを照射するため、n番目のRFパルスのフリップ角は(5)式となる。
　　　φ(n)＝｛a×n＋b｝×(－1)n      (n＝1,2,…,N)　　　(5)
　例えば、Nが5回の場合、Start-upシーケンスのRFパルス系列のフリップ角は、図3に示す様に、α/10,－3α/10,5α/10,－7α/10,9α/10となる。更にα＝90deg.の場合、9deg.,－27deg.,45deg.,－63deg.,81deg.となる。これに対して、(1)式で示す従来のLinear Flip Angle法のフリップ角は、同一条件で18deg.,－36deg.,54deg.,－72geg.,90deg.である。なお、図3では、傾斜磁場パルス波形を省略してある。傾斜磁場パルス波形の詳細は、例えば特許文献2に記載のパルスシーケンスを用いることができるので、ここでの詳細説明は省略する。
　Start-upシーケンスを用いた場合の、オンレゾナンスのスピンの振舞を図4に示す。(a)図に、本実施例のStart-upシーケンスを用いた場合を示し、比較のために、(b)図に従来のLinear Flip Angle法を用いた場合を示す。共に、縦緩和時間T1と横緩和時間T2が繰り返し時間TRに較べて十分長く、緩和過程を無視できるような極限を想定し、N＝5,α＝90degとした場合のRFパルス系列を有して成るStart-upシーケンスを実行する場合のMx－Mz平面でのスピンの挙動を示したものであり、図中の数値はRFパルスの照射回数を表す。(a)図に示す本実施例の場合のRFパルスのフリップ角は9deg.,－27deg.,45deg.,－63deg.,81deg.であり、(b)図の従来のLinear Flip Angle法の場合のRFパルスのフリップ角は、18deg.,－36deg.,54deg.,－72geg.,90deg.となる。
　図4に示す様に、本実施例のStart-upシーケンスを印加されるオンレゾナンスのスピンは、RFパルスの印加毎にMz軸を挟んで交互に運動するが、Start-upシーケンス終了時にはMz軸とのなす角がちょうどα/2となるので、エコー信号の振動を抑制することが可能である。これに対して従来のLinear Flip Angle法のStart-upシーケンスが印加されるオンレゾナンスのスピンは、RFパルスの印加毎にMz軸を挟んで交互に運動し、Start-upシーケンス終了時には、401と402の間を振動することになるのでエコー信号も振動してしまう。
　図5に、オフレゾナンスのスピンについて過渡状態での振動を数値計算した結果を示す。なお、図5では、縦軸がゼロのスピンが「オンレゾナンス」で、それ以外のスピンが「オフレゾナンス」と定義する(後述の図11でも同様)。(a)図に、本実施例のStart-upシーケンスを用いた場合を示し、比較のために、(b)図に従来のLinear Flip Angle法を用いた場合を示す。共に、横軸はRFパルスの照射回数を、縦軸は繰り返し時間TR間でのオフセット周波数を示し、横磁化の絶対値をグレースケースで表示した。数値計算に用いたパラメータは、T1＝1200ms、T2＝250ms、FA＝90deg.、TR＝4.0ms、N＝5である。
　(b)図に示す様に、従来のLinear Flip Angle法を用いた場合では、502の領域のようにRFパルスの照射周波数とのずれが小さいスピン(オフレゾナンススピン)が振動しているのに対して、本実施例のStart-upシーケンスを用いた場合は、(a)図に示す様に、オンレゾナンス付近でのスピン501について振動を十分に抑制していることがわかる。
　図6に、Start-upシーケンスを用いた場合のオンレゾナンスのスピンについて、過渡状態での振動のフリップ角依存性を数値計算した結果を示す。(a)(b)図に本実施例のStart-upシーケンスを用いた場合を示し、比較のために、(c)(d)図に従来のLinear Flip Angle法を用いた場合を示す。(a)(c)図の横軸はRFパルスの回数を、縦軸はフリップ角を示し、図5と同様に横磁化の絶対値をグレースケールで表示した。数値計算に用いたパラメータは図5と同じである。
　(b)(d)図は、(a)(c)図のそれぞれからフリップ角FAが45deg.,90deg.,135deg.であるスピンに関して、各RFパルスの回数における横磁化の絶対値を示す。(d)図より、Linear Flip Angle法では、FAが45deg.,90deg.,135deg.であるスピンが振動する。一方、(b)図より、本実施例のStart-upシーケンスでは、振動が生じないことがわかる。この傾向は、任意のフリップ角でも同様であり、(a)(c)図により、本実施例のStart-upシーケンスでは任意のフリップ角でエコー信号の振動を抑制できるが、従来のLinear Flip Angle法では、いずれのフリップ角でもエコー信号が振動していることが理解できる。したがって、本実施例のStart-upシーケンスを用いることで、照射不均一によりフリップ角にずれが生じた場合でも振動を抑制することが可能である。
　図6では簡単のために、オンレゾナンスのスピンについての数値計算結果を示したが、オフレゾナンスのスピンに関しては図5と同様、本実施例のStart-upシーケンスを用いた場合は、オンレゾナンス付近でのスピンについて振動を十分に抑制可能である。
　なお、本実施例では、Θ(n)の初項を0としたが、任意の値でも振動を抑制することが可能である。
　次に、上述した本実施例のStart-upシーケンスを構成する各RFパルスのフリップ角を求めて、該Start-upシーケンスを含むSSFPシーケンスを起動する演算処理部114の各機能を図7に示す機能ブロック図を用いて説明する。本実施例の演算処理部114は、撮像条件設定部701と、RFパルス系列設定部702と、パルスシーケンス設定部703と、撮像制御部704と、を有してなる。
　撮像条件設定部701は、表示・操作部113と記憶部115に接続され、表示部に撮像条件の設定(つまり撮像パラメータの値)入力を受け付ける入力画面を表示させる。操作者はこの入力画面を介して撮像条件を設定する。図9に入力画面の一例を示す。この入力画面には、少なくとも、データ収集シーケンスであるSSFPシーケンスにおけるフリップ角α(901)と、エコーデータ収集前の準備シーケンスであるStart-upシーケンスにおけるRFパルス系列のRFパルス数N(902)と、を含む撮像条件毎にその値の入力を受け付ける領域が表示される。操作者は、操作部を介して各撮像条件の値を入力する。そして、撮像条件設定部701は、操作者により入力された各種撮像条件の値を、該各種撮像条件の設定値として、メモリに記憶する。なお、フリップ角αとRFパルス数Nについては、記憶部115に記憶された所定の値を用いてもよい。
　RFパルス系列設定部702は、撮像条件設定部701に接続され、入力された撮像条件に基づいて、Start-upシーケンスのRFパルス系列のフリップ角を計算して設定する。RFパルス系列設定部702は、操作者が入力してメモリに記憶された、或いは、記憶部115から読み出した、フリップ角αとRFパルス数Nの値を用いて、(2)式に基づいて等差数列Θ(n)を作成し、この生成した等差数列Θ(n)を用いて(3)式に基づいてn番目のRFパルスのフリップ角の絶対値φ(n)を算出する。或いは、(4)式に基づいてn番目のRFパルスのフリップ角の絶対値φ(n)を直接算出してもよい。
　パルスシーケンス生成部703は、RFパルス系列設定部702と記憶部115に接続され、操作者が入力しメモリに記憶された撮像条件に基づいて、RFパルスや傾斜磁場パルスの印加タイミング等を設定する具体的データを求めてパルシーケンスを生成する。本実施例のパルスシーケンスは、準備シーケンスとしてのStart-upシーケンスと、データ収集シーケンスとしてのSSFPシーケンスと、から成るので、パルスシーケンス設定部703は、これらの2つのパルスシーケンスを連続して実行するための具体的データをそれぞれ生成する。具体的には、Start-upシーケンスのデータは、RFパルス系列設定部702で設定されたフリップ角に照射位相を付加した(5)式で表されるフリップ角のRFパルス系列を具体的に規定するRFパルスや傾斜磁場パルスについてのデータである。また、SSFPシーケンスのデータは、フリップ角αのRFパルスを照射して画像再構成用のエコー信号の計測を規定するRFパルスや傾斜磁場パルスについてのデータである。
　撮像制御部704は、パルスシーケンス生成部703と記憶部115と計測制御部111に接続され、パルスシーケンス生成部703で生成されたパルスシーケンスの具体的データを計測制御部111に通知し、計測制御部111にそのパルスシーケンスを実行させて撮像を開始する。
　次に、上記演算処理部114の各機能部が連携して行なう、本実施例の処理フローを図8に示すフローチャートに基づいて説明する。
　ステップ801で、撮像条件設定部701は、表示部に撮像条件の設定入力を受け付ける入力画面を表示する。操作者はこの入力画面を介して撮像条件を入力設定する。或いは、撮像条件設定部701は、記憶部115に事前に記憶された撮像条件の値を読み込む。そして、撮像条件設定部701は、入力された或いは読み込んだ各種撮像条件の値を演算処理部114のメモリに記憶する。
　ステップ802で、RFパルス系列設定部702は、ステップ801でメモリに記憶された撮像条件の値に基づいて、Start-upシーケンスのRFパルス系列のフリップ角を計算して演算処理部114のメモリに記憶する。
　ステップ803で、パルスシーケンス生成部703は、ステップ801でメモリに記憶された撮像条件の値と、ステップ802でメモリに記憶されたStart-upシーケンスのRFパルス系列のフリップ角と、に基づいて、RFパルスのフリップ角や傾斜磁場パルス形状及びそれらの印加タイミング等を規定する具体的データを求めて、Start-upシーケンスとSSFPシーケンスを連続して実行するための具体的データを生成する。
　ステップ804で、撮像制御部704は、ステップ803でパルスシーケンス生成部703が生成したパルスシーケンスの具体的データを計測制御部111に通知して、計測制御部111に生成されたパルスシーケンスを実行させて撮像を開始する。
　以上までが、本実施例の処理フローの説明である。
　以上説明したように、本実施例のMRI装置及びRFパルス制御方法は、Start-upシーケンスのRFパルス系列を構成する各RFパルスのフリップ角を、適切なオフセットと共に線形に増加させる。これにより、オンレゾナンススピンの振動を抑制することが可能である。従って、定常状態に導くStart-upシーケンスのRFパルスの総数を少なくしても、エコー信号の振動を低減して、画像にゴーストやブラーリング等のアーチファクトが生じず、高画質の画像を取得することができる。つまり、少ない空打ち回数でエコーデータ収集を開始することが可能となる。
　次に、本発明のMRI装置及びRFパルス制御方法の実施例2について説明する。本実施例は、立ち上がり、及び、α/2の直線との接点で滑らかに変化する数列、つまり、隣り合う二項間の差が単調増加の後単調減少する数列を用いて、Start-upシーケンスのRFパルス系列の各フリップ角を決定する。以下、図10，11に基づいて本実施例を詳細に説明する。
　前述の実施例1では、0とα/2の2点を繋ぐ直線をもとに数列Θを作成したが、一般には単調増加数列であればよい。そこで、単調増加数列を様々に変化させて数値計算した結果、立ち上がり、及び、α/2の直線との接点で滑らかに変化する数列、つまり、隣り合う二項間の差が単調増加の後単調減少する数列を用いて、Start-upシーケンスのRFパルス系列の各フリップ角を決定することにより、より広いオフセット周波数に対してエコー信号の振動が抑制されることを発明者は発見した。そのような数列の例として、
　　　Θ(n)＝α/4×(1－cos(nπ/N))　　　(6)
　が考えられる。この数列を図10(a)に示す。(6)式の数列について、(3)式に従い漸化式を解くと
　　　φ(n)＝α/2×(1-β×cos(nπ/N－γ))、n＝1,2,…,N　　　(7)
　となる。ここで、β2＝(1＋cos(π/N))/2、sin(γ)＝sin(π/N)/2β、cos(γ)＝(1＋cos(π/N))/2βである。(7)式で表される曲線の例としてN＝5の場合を図10(b)に示す。実際には、照射位相を180deg.ずらしながら各RFパルスを照射するため、n番目のRFパルスのフリップ角は(5)式のように表される。
　このようなフリップ角を有するRFパルス系列のStart-upシーケンスにより、実施例1と同様に、オンレゾナンスのスピンについては、過渡状態でのエコー信号の振動を抑制可能である。図11に、オフレゾナンスのスピン(図11で縦軸がゼロ以外のスピン)について過渡状態でのエコー信号の振動を数値計算した結果を示す。座標軸、及び、数値計算の方法は図5(a)と同様とする。図5(a)と図11を比較すると、図11ではより広い周波数オフセットに対して振動が抑制可能であることが分かる。なお、図11は、オンレゾナンスのスピン(図11で縦軸がゼロのスピン)についても振動が抑制されていることが示されている。
　次に、上述した本実施例のStart-upシーケンスにおける各RFパルスのフリップ角を求めて、該Start-upシーケンスを含むSSFPシーケンスを起動する演算処理部114の各機能を説明する。本実施例における演算処理部114の各機能は、図7に示した前述の実施例1と同様であるが、RFパルス系列設定部702の処理内容が異なる。以下異なる箇所のみを説明し、同一の箇所については説明を省略する。
　本実施例の演算処理部114のRFパルス系列設定部702は、(6)式に基づいて等差数列Θ(n)を作成し、この生成した等差数列Θ(n)を用いて(3)式に基づいてn番目のRFパルスのフリップ角の絶対値φ(n)を算出する。或いは、(7)式に基づいてn番目のRFパルスのフリップ角の絶対値φ(n)を直接算出してもよい。
　また、本実施例の演算処理部114の各機能部が連携して行なう、本実施例の処理フローは図8に示した前述の実施例1の処理フローと同様なので説明を省略する。
　以上説明したように、本実施例のMRI装置及びRFパルス制御方法は、立上り、及び、α/2の直線との接点で滑らかに変化する数列をもとに(3)式から算出したフリップ角のRFパルス系列を有して成るStart-upシーケンスを用いる。これにより、オンレゾナンススピンの振動を抑制するだけではなく、より広い周波数オフセットに関してオフレゾナンススピンの振動を抑制することが可能である。従って、実施例1と同様に、定常状態に導くStart-upシーケンスのRFパルスの総数を少なくしても、エコー信号の振動を低減して、画像にゴーストやブラーリング等のアーチファクトが生じず、高画質の画像を取得することができる。
　1　被検体、2　静磁場発生系、3　傾斜磁場発生系、4　シーケンサ、5　送信系、6　受信系、7　信号処理系、8　中央処理装置(CPU)、9　傾斜磁場コイル、10　傾斜磁場電源、11　高周波発信器、12　変調器、13　高周波増幅器、14a　高周波コイル(送信コイル)、14b　高周波コイル(受信コイル)、15　信号増幅器、16　直交位相検波器、17　A/D変換器、18　磁気ディスク、19　光ディスク、20　ディスプレイ、21　ROM、22　RAM、23　トラックボール又はマウス、24　キーボード

Claims (15)

  1. 　フリップ角が単調増加するRFパルス系列により過渡状態でのスピンの振動を抑制するStart-upシーケンスと、定常状態でエコー信号を計測するSSFPシーケンスとからなる撮像シーケンスを生成する演算処理部と、
    　前記撮像シーケンスに基づいて被検体からのエコー信号の計測を制御する計測制御部と、
    を備えた磁気共鳴イメージングであって、
    　前記RFパルス系列は、オフセットを有して単調増加するフリップ角であることを特徴とする磁気共鳴イメージング装置。
  2. 　請求項1記載の磁気共鳴イメージング装置において、
    　前記オフセットは、前記SSFPシーケンスのRFパルスのフリップ角(α)と、前記RFパルス系列のRFパルス数(N)に応じて変わる値であることを特徴とする磁気共鳴イメージング装置。
  3. 　請求項2記載の磁気共鳴イメージング装置において、
    　前記SSFPシーケンスのRFパルスのフリップ角(α)と、前記RFパルス系列のRFパルス数(N)の入力を受け付ける入力手段を備えていることを特徴とする磁気共鳴イメージング装置。
  4. 　請求項2記載の磁気共鳴イメージング装置において、
    　前記RFパルス系列のフリップ角は、単調増加数列の隣り合う二項間の和であることを特徴とする磁気共鳴イメージング装置。
  5. 　請求項4記載の磁気共鳴イメージング装置において、
    　前記単調増加数列は、ゼロと前記SSFPシーケンスのRFパルスのフリップ角(α)の半分の値との間の値をとることを特徴とする磁気共鳴イメージング装置。
  6. 　請求項2記載の磁気共鳴イメージング装置において、
    　前記RFパルス系列のフリップ角は、線形に増加することを特徴とする磁気共鳴イメージング装置。
  7. 　請求項6記載の磁気共鳴イメージング装置において、
    　前記RFパルス系列のフリップ角は、傾きがα／N、オフセットが－α／(2N)である直線に沿って増加することを特徴とする磁気共鳴イメージング装置。
  8. 　請求項2記載の磁気共鳴イメージング装置において、
    　前記RFパルス系列のフリップ角は、隣り合う二項間の差が単調増加の後単調減少する数列の隣り合う二項間の和であることを特徴とする磁気共鳴イメージング装置。
  9. 　請求項2記載の磁気共鳴イメージング装置において、
    　前記RFパルス系列のフリップ角は、立ち上がり、及び、α／2の直線との接点で滑らかに変化する数列の隣り合う二項間の和であることを特徴とする磁気共鳴イメージング装置。
  10. 　請求項2記載の磁気共鳴イメージング装置において、
    　前記RFパルス系列のフリップ角は、曲線に沿って増加することを特徴とする磁気共鳴イメージング装置。
  11. 　請求項10記載の磁気共鳴イメージング装置において、
    　前記RFパルス系列のn(n＝1,2,・・・,N)番目のフリップ角φ(n)は、
    　　φ(n)=α／2×(1－β×cos(nπ／N－γ))
    　　β2＝(1＋cos(π／N))／2
    　　sin(γ)＝sin(π／N)／2β
    　　cos(γ)＝(1＋cos(π／N))／2β
    　で表されることを特徴とする磁気共鳴イメージング装置。
  12. 　フリップ角が単調増加するRFパルス系列により過渡状態でのスピンの振動を抑制するStart-upシーケンスと、定常状態でエコー信号を計測するSSFPシーケンスとからなる撮像シーケンスにおける前記RFパルス系列のフリップ角を制御するRFパルス制御方法であって、
    　前記RFパルス系列のRFパルス数と、前記SSFPシーケンスにおけるRFパルスのフリップ角とを設定するステップと、
    　前記設定されたRFパルス系列のRFパルス数とSSFPシーケンスにおけるRFパルスのフリップ角とに基づいて、オフセットを有して単調増加するフリップ角を求めて前記RFパルス系列とするフリップ角算出ステップと、
    　前記求めたフリップ角のRFパルス系列を有する撮像シーケンスを実行するステップと、
    　を有することを特徴とするRFパルス制御方法。
  13. 　請求項12記載のRFパルス制御方法において、
    　前記フリップ角算出ステップは、単調増加数列の隣り合う二項間の和として、前記RFパルス系列の各RFパルスのフリップ角を求めることを特徴とするRFパルス制御方法。
  14. 　請求項13記載のRFパルス制御方法において、
    　前記単調増加数列は、初項がゼロでない、等差数列であることを特徴とするRFパルス制御方法。
  15. 　請求項13記載のRFパルス制御方法において、
    　前記単調増加数列は、隣り合う二項間の差が単調増加の後単調減少する数列であることを特徴とするRFパルス制御方法。
PCT/JP2011/064737 2010-07-07 2011-06-28 磁気共鳴イメージング装置及びｒｆパルス制御方法 WO2012005137A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US13/808,659 US9594136B2 (en) 2010-07-07 2011-06-28 Magnetic resonance imaging apparatus and RF pulse control method
JP2012523823A JP5942265B2 (ja) 2010-07-07 2011-06-28 磁気共鳴イメージング装置及びｒｆパルス制御方法

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2010-154409 2010-07-07
JP2010154409 2010-07-07

Publications (1)

Publication Number Publication Date
WO2012005137A1 true WO2012005137A1 (ja) 2012-01-12

Family

ID=45441120

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2011/064737 WO2012005137A1 (ja) 2010-07-07 2011-06-28 磁気共鳴イメージング装置及びｒｆパルス制御方法

Country Status (3)

Country Link
US (1) US9594136B2 (ja)
JP (1) JP5942265B2 (ja)
WO (1) WO2012005137A1 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20200281500A1 (en) * 2017-09-19 2020-09-10 University Of Cincinnati Systems and Methods for Detecting Small Physiological or Pathological Changes Using High Resolution Magnetic Resonance Imaging

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004329268A (ja) * 2003-04-30 2004-11-25 Hitachi Medical Corp 磁気共鳴イメージング装置
JP2006506115A (ja) * 2002-11-14 2006-02-23 コーニンクレッカ　フィリップス　エレクトロニクス　エヌ　ヴィ 平衡定常状態自由歳差運動ｍｒ撮影のための磁化準備シーケンス
JP2007068798A (ja) * 2005-09-08 2007-03-22 Ge Medical Systems Global Technology Co Llc Ｒｆパルス印加方法およびｍｒｉ装置
JP2008295925A (ja) * 2007-06-04 2008-12-11 Hitachi Medical Corp 核磁気共鳴イメージング装置
JP2009279202A (ja) * 2008-05-22 2009-12-03 Hitachi Medical Corp 核磁気共鳴イメージング装置

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003250775A (ja) * 2002-02-25 2003-09-09 Ge Medical Systems Global Technology Co Llc Ｍｒｉ装置およびｍｒａ撮影方法
DE102004005174B4 (de) * 2003-02-06 2014-11-20 Ge Medical Systems Global Technology Company, Llc MRI-Verfahren mit adiabatischer Vorbereitungssequenz
US7576535B2 (en) * 2003-12-08 2009-08-18 The Board Of Trustees Of The Leland Stanford Junior University Multi-compartment separation in magnetic resonance using transient steady-state free precession imaging
US7253620B1 (en) * 2004-03-08 2007-08-07 United States Of America, As Represented By The Secretary, Department Of Health And Human Services Spectrally selective suppression with steady-state free precession
JP2007068796A (ja) * 2005-09-08 2007-03-22 Ge Medical Systems Global Technology Co Llc Ｒｆパルス印加方法およびｍｒｉ装置
JP4864969B2 (ja) * 2006-06-16 2012-02-01 株式会社日立メディコ 磁気共鳴イメージング装置
US7479783B2 (en) * 2006-11-15 2009-01-20 Beth Israel Deaconess Medical Center, Inc. Echo train preparation for fast spin-echo acquisition
US8269494B2 (en) * 2008-02-29 2012-09-18 University Of Southern California Preparing the magnetization state of a sample for alternating repetition time steady state free precession magnetic resonance imaging
US8169218B2 (en) * 2008-04-18 2012-05-01 Beth Israel Deaconess Medical Center, Inc. Methods and apparatus for non-contrast enhanced pulmonary magnetic resonance imaging

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006506115A (ja) * 2002-11-14 2006-02-23 コーニンクレッカ　フィリップス　エレクトロニクス　エヌ　ヴィ 平衡定常状態自由歳差運動ｍｒ撮影のための磁化準備シーケンス
JP2004329268A (ja) * 2003-04-30 2004-11-25 Hitachi Medical Corp 磁気共鳴イメージング装置
JP2007068798A (ja) * 2005-09-08 2007-03-22 Ge Medical Systems Global Technology Co Llc Ｒｆパルス印加方法およびｍｒｉ装置
JP2008295925A (ja) * 2007-06-04 2008-12-11 Hitachi Medical Corp 核磁気共鳴イメージング装置
JP2009279202A (ja) * 2008-05-22 2009-12-03 Hitachi Medical Corp 核磁気共鳴イメージング装置

Also Published As

Publication number Publication date
US20130106417A1 (en) 2013-05-02
JPWO2012005137A1 (ja) 2013-09-02
JP5942265B2 (ja) 2016-06-29
US9594136B2 (en) 2017-03-14

Similar Documents

Publication Publication Date Title
US9687172B2 (en) System for motion corrected MR diffusion imaging
JP5942271B2 (ja) 磁気共鳴イメージング装置及び流体強調画像取得法
US8148981B2 (en) MRI apparatus and MRI method for SSFP with center frequency and 1st order gradient moments zeroed
JP5740307B2 (ja) 磁気共鳴イメージング装置及び傾斜磁場印加方法
JP6371554B2 (ja) 磁気共鳴イメージング装置
JP5484341B2 (ja) 磁気共鳴イメージング装置及び息止め撮像方法
WO2012077543A1 (ja) 磁気共鳴イメージング装置及びコントラスト強調画像取得方法
JPWO2004004563A1 (ja) 磁気共鳴イメージング装置及び渦電流補償導出方法
US9594135B2 (en) Magnetic resonance imaging apparatus and method for adjusting excitation region
JP5808659B2 (ja) 磁気共鳴イメージング装置及びＴ１ρイメージング法
WO2013002232A1 (ja) 磁気共鳴イメージング装置及びその傾斜磁場出力波形の測定方法
JP4230875B2 (ja) 磁気共鳴イメージング装置
JP5564213B2 (ja) 磁気共鳴イメージング装置
JP5942265B2 (ja) 磁気共鳴イメージング装置及びｒｆパルス制御方法
JP6203610B2 (ja) 磁気共鳴イメージング装置
JP2011078574A (ja) 磁気共鳴イメージング装置及び残留磁場抑制方法
JP5421600B2 (ja) 核磁気共鳴イメージング装置および核磁気共鳴イメージング装置の作動方法
JP5758230B2 (ja) 磁気共鳴イメージング装置及び反転ｒｆパルス位相制御方法
JP6157976B2 (ja) 磁気共鳴イメージング装置、及び方法
JP6118127B2 (ja) 磁気共鳴イメージング装置
JP2016131847A (ja) 磁気共鳴イメージング装置および磁気共鳴イメージング方法
JP2017123888A (ja) 磁気共鳴イメージング装置及びリフェーズ傾斜磁場印加方法
WO2012137606A1 (ja) 磁気共鳴イメージング装置及びラディアルサンプリング方法
JP5283213B2 (ja) 磁気共鳴イメージング装置
JP2012065736A (ja) 磁気共鳴イメージング装置及び実数成分画像取得方法

Legal Events

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

Ref document number: 11803472

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2012523823

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 13808659

Country of ref document: US

122 Ep: pct application non-entry in european phase

Ref document number: 11803472

Country of ref document: EP

Kind code of ref document: A1