WO2018088096A1 - 磁気共鳴イメージング装置及び計算画像生成方法 - Google Patents

磁気共鳴イメージング装置及び計算画像生成方法 Download PDF

Info

Publication number
WO2018088096A1
WO2018088096A1 PCT/JP2017/036823 JP2017036823W WO2018088096A1 WO 2018088096 A1 WO2018088096 A1 WO 2018088096A1 JP 2017036823 W JP2017036823 W JP 2017036823W WO 2018088096 A1 WO2018088096 A1 WO 2018088096A1
Authority
WO
WIPO (PCT)
Prior art keywords
parameter
magnetic resonance
subject
imaging apparatus
resonance imaging
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
PCT/JP2017/036823
Other languages
English (en)
French (fr)
Japanese (ja)
Inventor
陽 谷口
久晃 越智
亨 白猪
悦久 五月女
俊 横沢
知樹 雨宮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi Ltd filed Critical Hitachi Ltd
Priority to US16/342,111 priority Critical patent/US10866296B2/en
Publication of WO2018088096A1 publication Critical patent/WO2018088096A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • 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/50—NMR imaging systems based on the determination of relaxation times, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo sequences
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5618—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using both RF and gradient refocusing, e.g. GRASE

Definitions

  • the present invention relates to a magnetic resonance imaging technique, and more particularly to a technique for estimating an object parameter by calculation.
  • a magnetic resonance imaging (MRI) apparatus uses a signal intensity or phase information of a nuclear magnetic resonance signal obtained from an atom constituting a tissue of a subject, mainly a hydrogen nucleus, to obtain a nuclear density (proton density) image of the tissue, An image of a moving part such as a blood flow is acquired.
  • the signal intensity and phase of the nuclear magnetic resonance signal are determined by the conditions at the time of imaging and the characteristics of the apparatus and the tissue of the subject.
  • a technique has been widely used in MRI apparatuses for obtaining and imaging, as parameters, characteristics that can be analyzed with respect to the relationship with the nuclear magnetic resonance signal among the characteristics of the apparatus and the tissue of the subject. Yes.
  • One such technique is to take a plurality of images with different imaging parameters and obtain subject parameters and device parameters by calculation for each pixel.
  • the imaging parameters are the repetition time, the set intensity of the high-frequency magnetic field, the phase of the high-frequency magnetic field, and the subject parameters are the longitudinal relaxation time, the transverse relaxation time, the spin density, the resonance frequency, the diffusion coefficient, and the irradiation of the high-frequency magnetic field.
  • intensity distribution Such as intensity distribution.
  • Device parameters include magnetic field strength, reception sensitivity distribution, etc., but also depend on the subject.
  • An image having the obtained object parameter value as a pixel value is called a calculation image or map.
  • Patent Document 1 a procedure for acquiring maps such as relaxation time, frequency, irradiation intensity of a high-frequency magnetic field, spin density, resonance frequency, etc. as a subject parameter or an apparatus parameter is a gradient echo (GE) -based high-speed method.
  • GE gradient echo
  • An RF-soiled GE that is an imaging sequence is disclosed as an example.
  • T2 longitudinal relaxation time
  • T2 * is a transverse relaxation time including the effect of static magnetic field inhomogeneity
  • T2 is a spin echo (SE: Spin Echo) imaging sequence.
  • SE Spin Echo
  • T2 is a more important parameter than T2 * in grasping the state of the tissue, and a T2 weighted image is more widely used than a T2 * weighted image in general clinical examinations.
  • a T2-weighted image can be obtained by using an SE-based imaging sequence.
  • an SE-based imaging sequence has a longer imaging time than a GE-based imaging sequence. Further, if it is attempted to acquire not only T2 but also T2 * and other parameters, a longer shooting time is required.
  • the present invention has been made in view of the above circumstances, and an object thereof is to obtain a plurality of subject parameter maps at high speed.
  • an imaging sequence that generates both gradient echoes and spin echoes is used, and the values of one or more parameters are calculated using the gradient echoes.
  • the value of another parameter is calculated using echo.
  • the parameter value obtained by calculating the other parameter may be used.
  • the MRI apparatus of the present invention includes a measurement unit that applies a high-frequency magnetic field and a gradient magnetic field to a subject and measures a nuclear magnetic resonance signal emitted from the subject, a control unit that controls the measurement unit according to a pulse sequence, and the measurement Using a nuclear magnetic resonance signal acquired by the unit and a signal function of the pulse sequence, a parameter calculation unit that calculates a parameter value of an object parameter related to the characteristics of the object, and the control unit includes the As the pulse sequence, the measurement unit is controlled using a pulse sequence that measures at least two types of nuclear magnetic resonance signals after one excitation high-frequency magnetic field is applied, and the parameter calculation unit includes the two types of nuclear magnetic resonance.
  • One of the signals is used to calculate parameter values of one or more subject parameters including the first subject parameter, and the two types of nuclear magnetic resonance signals are calculated. Using the other, to calculate a parameter value for one or more subject parameters including different second analyte parameters of said first object parameters.
  • the calculation image generation method of the present invention is characterized in that the characteristics of the subject are obtained using echo signals obtained by executing a pulse sequence including gradient echo measurement and subsequent spin echo measurement a plurality of times while changing the value of imaging parameters.
  • a method for generating a calculated image of an object parameter related to an object, wherein two or more objects including a first object parameter using the gradient echo obtained by a plurality of imaging and the signal function of the pulse sequence The parameter value of the parameter is calculated, and the parameter value of the second parameter is calculated using the spin echo obtained by a plurality of imaging operations and the signal function after the gradient echo measurement.
  • the present invention it is possible to obtain a plurality of parameters at high speed by calculating the parameters in stages using a plurality of nuclear magnetic resonance signals obtained in a single imaging sequence.
  • the spin echo for T2 calculation can be measured by using the waiting time after measuring the gradient echo. Therefore, it becomes possible to acquire the T2 map while minimizing the extension of the photographing time.
  • Functional block diagram of a computer. 6 is a flow showing the operation of the MRI apparatus of the embodiment.
  • (A), (b) is a figure which shows the image obtained by the parameter set of FIG. 6, respectively, (c), (d) is a calculation image obtained from the image of (a), respectively, The figure which shows the calculation image obtained from an image.
  • FIG. 1 is a block diagram showing a schematic configuration of the MRI apparatus 100 of the present embodiment.
  • the MRI apparatus 100 irradiates a high-frequency magnetic field and a nuclear magnetism with a magnet 101 that generates a static magnetic field, a gradient coil 102 that generates a gradient magnetic field, a sequencer 104, a gradient magnetic field power source 105, a high-frequency magnetic field generator 106, and the like.
  • a transmission / reception coil 107 that detects a resonance signal, a receiver 108, a calculator 109, a display 111, and a storage device 112 are provided.
  • a transmission coil and a reception coil may be provided separately.
  • the magnet 101, the gradient magnetic field coil 102, the gradient magnetic field power source 105, the sequencer 104, the high frequency magnetic field generator 106, the transmission / reception coil 107, and the receiver 108 are collectively referred to as a measurement unit 110.
  • a subject (for example, a living body) 103 is placed on a bed (table) in a static magnetic field generated by a magnet 101.
  • the sequencer 104 sends commands to the gradient magnetic field power source 105 and the high frequency magnetic field generator 106 to generate a gradient magnetic field and a high frequency magnetic field, respectively.
  • the high frequency magnetic field is applied to the subject 103 through the transmission / reception coil 107.
  • a nuclear magnetic resonance signal generated from the subject 103 is received by the transmission / reception coil 107 and detected by the receiver 108.
  • the sequencer 104 sets a nuclear magnetic resonance frequency (detection reference frequency f0) as a reference for detection.
  • the detected signal is sent to the computer 109, where signal processing such as image reconstruction is performed.
  • the result is displayed on the display 111. If necessary, the detected signal and measurement conditions can be stored in the storage device 112.
  • the sequencer 104 normally performs control so that each device operates at a timing and intensity programmed in advance.
  • a program that particularly describes a high-frequency magnetic field, a gradient magnetic field, and the timing and intensity of signal reception is called a pulse sequence (imaging sequence).
  • imaging sequence a program that particularly describes a high-frequency magnetic field, a gradient magnetic field, and the timing and intensity of signal reception.
  • the MRI apparatus of this embodiment stores an imaging sequence that generates both gradient echoes and spin echoes.
  • the computer 109 includes a CPU and a memory, functions as a control unit that controls the operation of each unit described above, and functions as a calculation unit that performs various signal processing and calculations. Specifically, the echo signal is measured by operating the measurement unit according to the pulse sequence. The obtained echo signal is subjected to various signal processing to obtain a desired image. The image includes a calculated image in which the subject parameter is a pixel value. Programs and algorithms for control and computation are stored in the storage device 112, and the functions of the computer 109 are realized by the CPU of the computer 109 loading and executing the program stored in the storage device 112. Is done. Note that some of the functions of the computer 109 may be realized by hardware such as PLD (programmable logic device).
  • PLD programmable logic device
  • FIG. 2 shows a configuration example of the computer 109 for realizing the above-described processing.
  • the computer 109 includes a control unit 210 that controls the entire apparatus including the measurement unit and the calculation unit, and a calculation unit 230.
  • the calculation unit 230 includes an image reconstruction unit 231, a signal function generation unit 233, a parameter estimation unit 235, and an image generation unit 237.
  • the signal function generation unit 233, the parameter estimation unit 235, and the image generation unit 237 which are functional units as parameter calculation units among the functional units of the calculation unit 230, are computers provided independently from the MRI apparatus 100. Therefore, it may be realized in a computer capable of transmitting and receiving data to and from the computer 109 of the MRI apparatus 100.
  • FIG. 3 shows an outline of a procedure for creating a calculation image.
  • a combination of a plurality of shooting conditions is determined in advance and stored in the storage device 112 (S301).
  • the control unit 210 sets one of the plurality of imaging conditions, controls the measurement unit 110, executes a predetermined pulse sequence, and measures an echo signal, here, a gradient echo and a spin echo.
  • Shooting is performed (S302).
  • S303 For each echo signal, when the number of measurements necessary for image reconstruction is completed, shooting is performed under different shooting conditions (S303). Shooting is repeated until shooting of all combinations of planned shooting conditions is completed.
  • the computer 109 (signal function generating unit 233) generates the signal function (S304).
  • the computer 109 (image reconstruction unit 231) reconstructs an image (GE image and SE image) for each of the two types of echo signals obtained by a plurality of photographing, and the parameter estimation unit 235 uses the pixel values of these images.
  • the parameter is estimated using the signal function generated by the signal function generation unit 233 (S305). Note that some of the calculations performed by the computer 109 may be performed in parallel with the shooting before the shooting is completed. If the shooting sequence is the same, the signal function is the same even if the shooting conditions are different. Therefore, by storing the generated signal function, it is not necessary to generate the signal function every time shooting is performed, and the same signal function can be used repeatedly.
  • the parameter estimation includes estimation using a gradient echo (estimation of a first parameter) and estimation using a spin echo (estimation of a second parameter), and different types of parameters are calculated in each process.
  • the image generation unit 237 creates an image, that is, a calculated image (parameter image) having pixel values as values for all or some of the calculated parameters (S306).
  • the image generation unit 237 creates a display image including the calculation image as a display image or further including the calculation image, and displays the display image on the display 111 (S307).
  • the imaging conditions are parameters (imaging parameters) that can be arbitrarily set by the user during execution of the imaging sequence. Specifically, the repetition time (TR), the echo time (TE), the set intensity of the high-frequency magnetic field (flip angle ( Flip Angle (FA)), high-frequency magnetic field phase increment ( ⁇ ), and the like. In the present embodiment, a plurality of combinations of these values are prepared.
  • an imaging sequence that generates both gradient echoes and spin echoes after the application of one excitation pulse is used as the imaging sequence.
  • an imaging sequence for example, a sequence combining an RF-soiled GE sequence and an SE sequence (hereinafter referred to as a GE-SE sequence) can be used.
  • the RF-soiled GE sequence can perform 3D imaging at high speed, and the pixel values of the images obtained by this imaging sequence are mainly relaxation parameters T1, T2 *, spin density ⁇ , and apparatus parameters. It depends on B1 and Sc.
  • FIG. 4A An example of the GE-SE sequence is shown in FIG. 4A.
  • RF, A / D, Gs, Gp, and Gr represent a high-frequency magnetic field, signal reception, a slice gradient magnetic field, a phase encoding gradient magnetic field, and a readout gradient magnetic field, respectively.
  • This figure shows a case where the axis of the phase encoding gradient magnetic field Gp is uniaxial, but in the case of a 3D-sequence, a biaxial phase encoding gradient magnetic field is used.
  • This GE-SE sequence uses the waiting time after the echo measurement of the RF-soiled GE sequence to generate a spin echo by an inverted RF pulse.
  • the sequence up to the gradient echo measurement is shown in FIG. It is the same as the spoiled GE sequence. That is, first, the application of the slice gradient magnetic field pulse 401 and the irradiation of the radio frequency magnetic field (RF) pulse 402 excite magnetization of a certain slice in the target object. Next, after applying a slice rephase gradient magnetic field pulse 403, a phase encoding gradient magnetic field pulse 404, and a dephase readout gradient magnetic field 405, a nuclear magnetic resonance signal (gradient echo, first echo) while applying a readout gradient magnetic field pulse 406 is applied. ) 407 is measured. Finally, a phase encoding gradient magnetic field pulse 408 for dephase is applied.
  • RF radio frequency magnetic field
  • the application of the slice gradient magnetic field pulse 409-1 and irradiation of the inversion pulse 410-1 reverse the magnetization in the slice.
  • the nuclear magnetic resonance signal spin echo, second echo
  • a phase encoding gradient magnetic field pulse 414-1 for dephase is applied.
  • the same sequence as the sequence from the application of the slice gradient magnetic field pulse 409-1 to the application of the phase encoding gradient magnetic field pulse 414-1 for dephasing is repeated as many times as necessary.
  • a total of four inversion pulses (410-1 to 410-4) are applied, and a total of four spien echoes from the second echo to the fifth echo are measured (413- 1 to 413-4).
  • the procedure from the RF pulse 402 irradiation to the last spin echo measurement is repeated at a repetition time TR, and each echo from the first echo to the fifth echo is measured a plurality of times.
  • the intensity (phase encoding amount kp) of the phase encoding gradient magnetic field pulse (404, 408, 411-1 to 411-4, 414-1 to 414-4) is changed at each repetition, and the phase increment value of the RF pulse 402 is changed.
  • the number of inversion pulses in the GE-SE sequence is arbitrary. However, when the repetition time is as short as several tens of milliseconds, an even number is desirable. This is because, like RF-soiled GE, the waiting time required until the next excitation is shortened to shorten the imaging time. In other words, RF-soiled GE is designed so that the magnetization excited by the excitation pulse has some longitudinal magnetization component facing the direction of the static magnetic field, and the waiting time until the next excitation pulse can be shortened. ing.
  • the inversion pulse is irradiated, the longitudinal magnetization is reversed, and when it is irradiated again, the direction of the original static magnetic field is restored. As described above, by making the inversion pulse an even number of times, the longitudinal magnetization can be set in the direction of the static magnetic field before the next excitation pulse as in the RF-soiled GE in this sequence.
  • RF-soiled GE that does not use an inversion pulse is also used.
  • the imaging parameter set is a combination of imaging parameters such as FA (flip angle), TR (repetition time), TE (echo time), inversion pulse interval, and ⁇ (RF phase increment value) as predetermined parameter values.
  • FA flip angle
  • TR repetition time
  • TE echo time
  • ⁇ RF phase increment value
  • the parameter value is determined in consideration of the type of subject parameter to be calculated.
  • the FA is set to about 5 to 60 degrees in a normal RF-soiled GE, but is made as small as possible in the present embodiment. This is due to the following reason.
  • the estimated values of T1 and T2 are also influenced by the protein concentration, and T1 and T2 different from the case where T1 and T2 are individually measured by a normal method are obtained.
  • the maximum value of TR is determined in consideration of the imaging sequence used and the number of inversion pulses. For example, in the GE-SE sequence, since an inversion pulse is added to the RF-soiled GE, the TR becomes longer than that of the RF-soiled GE. For this reason, it is better that the shooting parameter set using the GE-SE sequence has a long TR. Note that if the interval between the inversion pulses and TR of the GE-SE sequence is increased, the estimation accuracy of the long T2 is improved, but the imaging time is increased.
  • the phase increment ⁇ of the RF pulse changes the phase of the RF pulse in order to make the influence of the transverse relaxation different.
  • signals having different influences of the transverse relaxation can be obtained.
  • the effect of lateral relaxation is eliminated by setting the phase increment to about 20 degrees, so the increment is changed within a range below that.
  • the TE and the inversion pulse interval are determined in consideration of device restrictions and SAR. These imaging parameters may be fixed values.
  • the number of shooting parameter sets is equal to or greater than the number of parameters estimated in parameter estimation described later (the number of unknowns). In this embodiment, since the number of unknowns is 4 (T1, T2, B1, Sc), the shooting parameter set is 4 or more. As the number of parameter sets, that is, the number of images obtained by shooting the parameter set increases, the estimation accuracy improves, but the shooting time increases accordingly.
  • FIG. 1 An example of a parameter set selected to minimize noise based on the error propagation law is shown in FIG.
  • six shooting parameter sets P1 to P each consisting of combinations of 10 degrees, 40 degrees FA, 2 degrees, 5 degrees, 7 degrees, 8 degrees, 22 degrees and TR of 10 ms, 30 ms, 40 ms, are set.
  • P6 is determined.
  • the gradient echo TEs are all set to 3 ms.
  • the imaging sequence is P3 with the longest TR, the GE-SE sequence shown in FIG. 4A, and the rest of the imaging sequence is the RF-soiled GE shown in FIG.
  • the interval between inversion pulses of the GE-SE sequence is 8 ms.
  • the measurement unit 110 captures a plurality of images using the plurality of imaging parameter sets described above under the control of the control unit 210. That is, a plurality of times of imaging are performed while changing the imaging parameter set, and a plurality of gradient images (GE images) and a plurality of spin echo images (SE images) are obtained.
  • GE images gradient images
  • SE images spin echo images
  • echoes obtained by the GE-SE sequence in FIG. 4A are arranged in k space as shown in FIG. 4B for each echo number, and an image is reconstructed by performing two-dimensional inverse Fourier transform. . For example, a first echo image is reconstructed from the first echo, and a second echo image is reconstructed from the second echo.
  • the first echo image is a GE image
  • the second and subsequent echo images are SE images.
  • one GE image is obtained by one imaging.
  • 6 GE images and 4 SE images are obtained.
  • the calculation unit 230 calculates the subject parameter and the apparatus parameter using the plurality of images acquired as described above.
  • processing of the calculation unit 230 related to parameter calculation will be described with reference to FIG.
  • FIG. 8 is a diagram showing the relationship between each process and input and output data.
  • the signal function generation unit 233 When the shooting sequence is determined, the signal function generation unit 233 generates a function (signal function) representing the signal intensity of each pixel obtained by the shooting sequence by numerical simulation.
  • a signal function of a GE-SE sequence is generated.
  • the signal function fs is a function of imaging parameters (FA, TR, TE, ⁇ ), apparatus parameters, and subject parameters, and is expressed as follows.
  • T1, T2, and ⁇ are the longitudinal relaxation time, lateral relaxation time, and spin density of the subject parameter, respectively, and B1 and Sc are parameters that depend on the characteristics of the device and the properties of the subject (herein referred to as device parameters).
  • device parameters RF irradiation intensity and receiving coil sensitivity.
  • B1 is a coefficient of FA at the time of photographing, it is converted into a product form with FA.
  • ⁇ and Sc act on the signal intensity as proportional coefficients, they are put out of the function, and TE is also put out of the function because it affects the signal intensity in the form of an exponential function.
  • equation (1) can be rewritten as the equation in the second stage.
  • a function f that is the basis of fs is created by numerical simulation. That is, the spin density ⁇ and B1 and Sc of the imaging target are set to 1 and TE is set to 0, and arbitrary values are set for the subject parameters T1 and T2, and the imaging parameters FA, TR, and ⁇ are set for these.
  • the signal is calculated by numerical simulation with comprehensive changes.
  • the range of parameter values to be changed includes the range of imaging parameters used for actual imaging (FIG. 3: S302) and the range of T1 and T2 of the subject.
  • An example of each parameter value of the imaging parameter and the subject parameter is shown below. The number after the parameter is the number to be changed, and the numerical value after “:” is the parameter value.
  • An imaging parameter set (173400 sets in the above example) 810 composed of all combinations of these imaging parameters and subject parameters is constructed, and each signal value is calculated by computer simulation.
  • an object model in which spins are arranged on lattice points, an imaging sequence, imaging parameters, and apparatus parameters are input, and a Bloch equation, which is a basic equation of the magnetic resonance phenomenon, is solved and a magnetic resonance signal is output.
  • the object model is given as a spatial distribution of spins ( ⁇ , M 0 , T1, T2).
  • ⁇ is a magnetic rotation ratio
  • M 0 is thermal equilibrium magnetization (spin density)
  • T1 and T2 are longitudinal relaxation time and transverse relaxation time, respectively.
  • the Bloch equation is a first-order linear ordinary differential equation and is expressed by the following equation.
  • (x, y, z) represents a three-dimensional orthogonal coordinate system, and z is equal to the direction of the static magnetic field (intensity is B 0 ).
  • (Mx, My, Mz) is spin
  • Gx, Gy, Gz are gradient magnetic field strengths in the subscript direction
  • H 1 is high-frequency magnetic field strength
  • ⁇ f 0 is the frequency of the rotating coordinate system.
  • a signal function f is created by interpolation, and fs 820 is created according to equation (1).
  • interpolation linear interpolation or spline interpolation of about 1st to 3rd order can be used.
  • the signal function fs can also be expressed by the following equation (3) obtained by modifying the equation (1).
  • FIG. 9 shows a part of the signal function f created as described above.
  • the signal function f obtained by the above simulation is a five-dimensional function with TR, FA, ⁇ , T1, and T2 as variables.
  • the parameter estimation is performed by using a plurality of imaging parameter sets 530 (S802), and using the echo-specific images 840 and 850 and the imaging sequence signal function 820, the subject of the actual imaging.
  • T2 * is estimated using a GE image
  • T2 is estimated using an SE image (first parameter estimation)
  • T2 is estimated using an SE image (second parameter estimation).
  • T1, T2, B1, and a are estimated using a plurality of GE images and the signal function 820 generated in the signal function generation S801. More specifically, the unknown values T1 and T2, B1, and a are estimated by fitting the signal value I for each pixel to the signal function fs (Equation (3)) generated in S801.
  • the function fitting can be performed by the least square method represented by Expression (4).
  • I is a pixel value of an image captured with a predetermined imaging parameter set (FA, ⁇ , TR, TE), and ⁇ is T1, T2, and B1 estimated by the signal function of Expression (3).
  • A is the sum of the residuals between the value calculated by substituting a and the pixel value of the image.
  • T1 and T2, B1, and a are estimated so that the total sum ⁇ of the residuals is minimized.
  • T2 estimated here is T2 * because a gradient echo image is used as the original image.
  • the estimated T1 and T2 *, B1, and a are output as parameters 860. Since the parameter is calculated as a numerical value for each pixel, it becomes a map of each parameter, that is, a calculated image.
  • T2 is estimated using a plurality of spin echo images 850.
  • Expression (5) representing T2 attenuation of the MR signal is used as a signal function.
  • T2 and a ′ are estimated by fitting the value I to equation (5).
  • the TEs of the four images SE1 to SE4 are, for example, 8 ms, 16 ms, 24 ms, and 32 ms, respectively, when the inversion pulse interval is 8 ms.
  • T1 calculated in the first parameter estimation S803 can be used as T1.
  • the function fitting can be performed by a least square method that minimizes a residual between a pixel value of an image obtained by photographing and a value calculated from the signal function (5).
  • T2 estimated in this way is different from T2 * obtained in the first parameter estimation, and is a “true lateral relaxation time” that is not affected by the static magnetic field inhomogeneity.
  • the parameter estimation unit 235 outputs T2 as a parameter.
  • “a ′” is the parameter a estimated in the first parameter estimation S803, T1, T2 (T2 *), and B1 as shown in the equation (6). It is expressed using. Therefore, in the second parameter estimation S804, only T2 may be estimated after giving in advance the values of a, T1, T2 (T2 *), and B1 calculated in S803 as the correct value of “a ′”. Thereby, the number of unknowns is reduced from two to one, so that the estimation accuracy of T2 is improved.
  • the above is parameter estimation processing S305 in FIG.
  • the calculated image may be generated for all the estimated parameters or for some parameters.
  • a T1 weighted image or a T2 / T2 * weighted image may be generated using a T1 image or a T2 / T2 * image.
  • the image generation unit 237 displays the generated calculation image or emphasized image on the display 111 in various display forms.
  • the display form may be, for example, a calculation image in which pixel values are expressed in black and white shading, or may be displayed in color.
  • the calculation image may be displayed alone or in parallel with the proton density image obtained in the photographing (S302). It is also possible to display parameter values and value ranges of specific parts as numerical values.
  • the GE echo and the SE echo can be obtained without extending the entire imaging time by performing the imaging sequence for acquiring the spin echo using the waiting time of the GE system sequence.
  • Two types of echoes can be acquired.
  • the number of spin echoes acquired in the GE-SE sequence is four, but this number is arbitrary as long as it is one or more.
  • the number of spin echoes may be 1 or more.
  • the number of spin echoes needs to be two or more.
  • the number of spin echoes is larger. Practically, it is desirable that the number is 3 or more.
  • the imaging sequence for measuring one gradient echo is adopted in the parameter set with a short TR (for example, P1, P2, and P4 to P6 in FIG. 6), but the first echo is used in the RF-soiled GE in FIG.
  • a multi-echo sequence for measuring second, third,..., Gradient echoes can be used.
  • the same phase encoding as that of the first echo to the echoes after the second echo, it is possible to acquire images of echoes having different TEs. Thereby, the estimation accuracy of the subject parameter can be increased as in the second modification.
  • ⁇ Modification 4> B1, T1, T2 *, and a are calculated in the first parameter estimation S803, and a (a ′) and T2 are calculated in the second parameter estimation S804, but T2 * and T2 Other parameters may be estimated by any process. However, by using the parameter estimated in one process for the other estimation process, the number of unknowns can be reduced, the number of necessary images (for example, the number of spin echoes) can be reduced, and the imaging time can be shortened.
  • DESCRIPTION OF SYMBOLS 100 MRI apparatus, 101: Magnet which generate

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Animal Behavior & Ethology (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Signal Processing (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
PCT/JP2017/036823 2016-11-14 2017-10-11 磁気共鳴イメージング装置及び計算画像生成方法 Ceased WO2018088096A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US16/342,111 US10866296B2 (en) 2016-11-14 2017-10-11 Magnetic resonance imaging apparatus and calculation image generation method using pulse sequence

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2016221863A JP6721489B2 (ja) 2016-11-14 2016-11-14 磁気共鳴イメージング装置及び計算画像生成方法
JP2016-221863 2016-11-14

Publications (1)

Publication Number Publication Date
WO2018088096A1 true WO2018088096A1 (ja) 2018-05-17

Family

ID=62110406

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2017/036823 Ceased WO2018088096A1 (ja) 2016-11-14 2017-10-11 磁気共鳴イメージング装置及び計算画像生成方法

Country Status (3)

Country Link
US (1) US10866296B2 (enExample)
JP (1) JP6721489B2 (enExample)
WO (1) WO2018088096A1 (enExample)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114690101A (zh) * 2022-03-30 2022-07-01 西门子数字医疗科技(上海)有限公司 磁共振成像中的参数定量检测方法、装置及磁共振扫描仪

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7161377B2 (ja) * 2018-11-06 2022-10-26 富士フイルムヘルスケア株式会社 画像処理装置、及び、それを含む磁気共鳴イメージング装置及び磁気共鳴イメージングシステム
US11796618B2 (en) * 2019-07-12 2023-10-24 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for magnetic resonance imaging
CN115184852B (zh) * 2022-07-18 2025-07-25 沈阳工业大学 一种磁共振成像的方法及系统
JP2024064249A (ja) * 2022-10-27 2024-05-14 富士フイルムヘルスケア株式会社 磁気共鳴イメージング装置、画像処理装置、及び画像処理方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01221150A (ja) * 1988-02-29 1989-09-04 Shimadzu Corp Ｍｒｉのデータ収集法
JPH04170938A (ja) * 1990-11-01 1992-06-18 Toshiba Corp 磁気共鳴イメージング方法
JPH0654816A (ja) * 1992-08-06 1994-03-01 Hitachi Ltd 核磁気共鳴を用いた検査方法
JPH0779949A (ja) * 1993-09-14 1995-03-28 Toshiba Corp 磁気共鳴映像装置
JPH0966043A (ja) * 1995-08-31 1997-03-11 Shimadzu Corp Ｍｒイメージング装置
JPH11322A (ja) * 1997-04-11 1999-01-06 Univ Pennsylvania 磁気共鳴映像法（ｍｒｉ）における横方向緩和速度への可逆性寄与分測定方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3445689A1 (de) * 1984-12-14 1986-06-19 Max-Planck-Gesellschaft zur Förderung der Wissenschaften e.V., 3400 Göttingen Verfahren und einrichtung zur ortsaufgeloesten untersuchung einer probe mittels magnetischer resonanz von spinmomenten
US4896113A (en) * 1988-11-25 1990-01-23 General Electric Company Use of repeated gradient echoes for noise reduction and improved NMR imaging
US5270654A (en) * 1991-07-05 1993-12-14 Feinberg David A Ultra-fast multi-section MRI using gradient and spin echo (grase) imaging
DE4432575C2 (de) 1993-09-14 2003-04-10 Toshiba Kawasaki Kk Verfahren zur Bildgebung der Gehirnfunktion mittels einer Kernspinresonanz-Vorrichtung und hieran angepasste Kernspinresonanz-Vorrichtung
US5594336A (en) * 1995-06-02 1997-01-14 Picker International, Inc. Three point technique using spin and gradient echoes for water/fat separation
US20080272779A1 (en) * 2005-06-01 2008-11-06 Koninklijke Philips Electronics, N.V. Determination of Distribution Information of a Contrast Agent by Mr Molecular Imaging
JP5449903B2 (ja) 2009-07-29 2014-03-19 株式会社日立メディコ 磁気共鳴イメージング装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01221150A (ja) * 1988-02-29 1989-09-04 Shimadzu Corp Ｍｒｉのデータ収集法
JPH04170938A (ja) * 1990-11-01 1992-06-18 Toshiba Corp 磁気共鳴イメージング方法
JPH0654816A (ja) * 1992-08-06 1994-03-01 Hitachi Ltd 核磁気共鳴を用いた検査方法
JPH0779949A (ja) * 1993-09-14 1995-03-28 Toshiba Corp 磁気共鳴映像装置
JPH0966043A (ja) * 1995-08-31 1997-03-11 Shimadzu Corp Ｍｒイメージング装置
JPH11322A (ja) * 1997-04-11 1999-01-06 Univ Pennsylvania 磁気共鳴映像法（ｍｒｉ）における横方向緩和速度への可逆性寄与分測定方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114690101A (zh) * 2022-03-30 2022-07-01 西门子数字医疗科技(上海)有限公司 磁共振成像中的参数定量检测方法、装置及磁共振扫描仪
US12306277B2 (en) 2022-03-30 2025-05-20 Siemens Healthineers Ag Method and apparatus for quantitative detection of parameters in MRI, and magnetic resonance scanner

Also Published As

Publication number Publication date
US20190250231A1 (en) 2019-08-15
US10866296B2 (en) 2020-12-15
JP6721489B2 (ja) 2020-07-15
JP2018078959A (ja) 2018-05-24

Similar Documents

Publication Publication Date Title
JP5449903B2 (ja) 磁気共鳴イメージング装置
JP6371554B2 (ja) 磁気共鳴イメージング装置
US11122988B2 (en) Magnetic resonance imaging apparatus
US8754645B2 (en) Method for spatially resolved determination of an MR parameter
JP6721489B2 (ja) 磁気共鳴イメージング装置及び計算画像生成方法
WO2012077543A1 (ja) 磁気共鳴イメージング装置及びコントラスト強調画像取得方法
US20150177350A1 (en) Methods and Systems for Improved Magnetic Resonance Acquisition
US11435423B2 (en) Image processing apparatus, magnetic resonance imaging apparatus including the same, and magnetic resonance imaging system
Malik et al. Phase relaxed localized excitation pulses for inner volume fast spin echo imaging
JP2018078959A5 (enExample)
JP7020930B2 (ja) 磁気共鳴イメージング装置、磁気共鳴イメージングシステム及びパラメータ推定方法
JP5684888B2 (ja) 磁気共鳴イメージング装置
JP2013192578A (ja) 磁気共鳴イメージング装置及びプログラム
JP6487554B2 (ja) 磁気共鳴イメージング装置
JP6114846B2 (ja) 磁気共鳴イメージング装置
JP7377231B2 (ja) 磁気共鳴イメージング装置、及び、画像処理方法
JP7448418B2 (ja) 磁気共鳴イメージング装置、及び、磁気共鳴イメージング方法
JP5881869B2 (ja) 磁気共鳴イメージング装置
JP2024064249A (ja) 磁気共鳴イメージング装置、画像処理装置、及び画像処理方法
JP2016140417A (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: 17870199

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 17870199

Country of ref document: EP

Kind code of ref document: A1