WO2013114927A1 - 磁気共鳴撮影装置 - Google Patents
磁気共鳴撮影装置 Download PDFInfo
- Publication number
- WO2013114927A1 WO2013114927A1 PCT/JP2013/050290 JP2013050290W WO2013114927A1 WO 2013114927 A1 WO2013114927 A1 WO 2013114927A1 JP 2013050290 W JP2013050290 W JP 2013050290W WO 2013114927 A1 WO2013114927 A1 WO 2013114927A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- temperature
- information
- magnetic resonance
- voxel
- imaging apparatus
- Prior art date
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/46—NMR spectroscopy
- G01R33/4625—Processing of acquired signals, e.g. elimination of phase errors, baseline fitting, chemometric analysis
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/01—Measuring temperature of body parts ; Diagnostic temperature sensing, e.g. for malignant or inflamed tissue
- A61B5/015—By temperature mapping of body part
-
- 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/4804—Spatially selective measurement of temperature or pH
-
- 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
- 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/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
Definitions
- the present invention relates to magnetic resonance imaging technology.
- the present invention relates to a magnetic resonance spectroscopy (MRS) technique for calculating a temperature image from a resonance frequency difference between water and a metabolite, and a magnetic resonance spectroscopy imaging (MRSI) technique.
- MRS magnetic resonance spectroscopy
- MRSI magnetic resonance spectroscopy imaging
- the magnetic resonance imaging apparatus is an apparatus that acquires physical and chemical information of a measurement target by inducing a magnetic resonance phenomenon by irradiating a measurement target placed in a static magnetic field with a high-frequency magnetic field having a specific frequency.
- Magnetic Resonance Imaging which is currently widely used, mainly uses the nuclear magnetic resonance phenomenon of hydrogen nuclei in water molecules to visualize differences in hydrogen nuclei density and relaxation time that vary depending on biological tissues. It is a method to do. As a result, tissue differences can be imaged, which is highly effective in diagnosing diseases.
- MRS and MRSI separate the nuclear magnetic resonance signals for each molecule based on the difference in the resonance frequency (chemical shift) due to the difference in the chemical bond of the molecule (metabolite), and the concentration and relaxation time for each molecular species. This is a method of measuring the difference between the two.
- MRS is a method of observing molecular species in a selected spatial region
- MRSI is a method of imaging for each molecular species.
- target nuclei include 1 H (proton), 31 P, 13 C, and 19 F.
- MRS / MRSI proton MRSI
- MRS / MRSI proton MRSI
- NAA N-acetylaspartic acid
- MRS / MRSI can be applied not only to the measurement of metabolite concentration, but also to in-vivo temperature measurement using the resonance frequency difference between water and metabolite. It is known that the resonance frequency of water shifts with the temperature, and the shift amount has a temperature coefficient of ⁇ 0.01 ppm / ° C. (for example, Non-Patent Document 1). On the other hand, it is known that the resonance frequency of metabolites such as NAA does not change in the temperature range under the biological environment. There is a technique for measuring the temperature in a living body from the difference in resonance frequency between water and a metabolite using these characteristics (for example, see Non-Patent Document 2).
- the calculation of the temperature information is performed by fitting with a model function as follows, for example. First, water and metabolites (here, NAA is taken as an example) are measured individually or simultaneously. Then, the spectrum is obtained by performing Fourier transform in the time direction. The measured water and NAA spectrum peak regions (spectrum peaks) are fitted using, for example, the Lorentz function of the following equation (1).
- ⁇ is the frequency
- L i is the signal intensity
- ⁇ i is the resonance frequency of the target substance
- a i is the half width of the spectrum peak
- I i is the height of the spectrum peak
- ⁇ i is the phase
- c is a constant. Term.
- the fitting result depends on the spectrum quality such as the half width of the peak and the signal-to-noise ratio, and the calculated temperature varies.
- the operator cannot grasp the accuracy of the calculated temperature. Since the calculated temperature accuracy (reliability) is essential information for diagnosis, the temperature is calculated by measuring a plurality of times and the accuracy of the temperature calculated from the standard deviation thereof is confirmed. For this reason, it takes a lot of time to measure the temperature of one subject and the burden on the subject is also great.
- the present invention has been made in view of the above circumstances, and an object of the present invention is to provide a technique for acquiring in-vivo temperature information and its accuracy information, which can be completed in a short time and has less burden on the subject.
- a static magnetic field generation unit that generates a static magnetic field in a space in which the subject is placed, a high-frequency magnetic field irradiation unit that irradiates the subject with a high-frequency magnetic field, and a gradient magnetic field application that applies a gradient magnetic field to the subject
- a control unit a receiving unit that receives a nuclear magnetic resonance signal generated from the subject, a gradient magnetic field applying unit, the high-frequency magnetic field irradiating unit, and the receiving unit to control operations of two substances having different resonance frequencies
- a magnetic resonance imaging apparatus comprising: a measurement unit that obtains a nuclear magnetic resonance signal; an arithmetic unit that performs arithmetic processing on the nuclear magnetic resonance signal; and a display device that displays information after the arithmetic processing.
- a spectrum calculation unit for calculating the spectrum of each of the two magnetic resonance signals of the two substances having different resonance frequencies, and calculating temperature information in the subject based on each of the calculated spectrum peaks.
- a temperature information calculation unit that calculates temperature accuracy information that indicates the accuracy of the temperature information based on the calculated spectrum peaks, and the display based on the temperature information and the temperature accuracy information.
- a magnetic resonance imaging apparatus comprising: a display information generation unit that generates display information to be displayed on the apparatus.
- in-vivo temperature information and accuracy information thereof can be acquired in a short time and with a low burden on the subject.
- (A)-(c) is an external view of the magnetic resonance imaging apparatus of embodiment of this invention.
- 1 is a functional configuration diagram of a magnetic resonance imaging apparatus according to an embodiment of the present invention. It is a functional block diagram of the computer with which the magnetic resonance imaging apparatus of embodiment of this invention is provided. It is a flowchart for demonstrating the flow of the whole temperature information and temperature accuracy information calculation measurement of embodiment of this invention. It is explanatory drawing for demonstrating an example of the MRSI sequence of embodiment of this invention.
- (A)-(c) is explanatory drawing for demonstrating the area
- FIG. 1 is an external view of the MRI apparatus of this embodiment.
- FIG. 1A shows a horizontal magnetic field type MRI apparatus 100 using a tunnel magnet that generates a static magnetic field with a solenoid coil.
- FIG. 1B shows a hamburger type (open type) vertical magnetic field type MRI apparatus 120 in which magnets are separated into upper and lower sides in order to enhance the feeling of opening.
- FIG. 1C shows an MRI apparatus 130 that uses the same tunnel-type magnet as in FIG. 1A and has a feeling of openness by shortening the depth of the magnet and tilting it obliquely. In the present embodiment, any of these MRI apparatuses having these appearances can be used.
- the MRI apparatus of the present embodiment is not limited to these forms.
- various known MRI apparatuses can be used regardless of the form and type of the apparatus.
- the MRI apparatus 100 is representative.
- FIG. 2 is a functional configuration diagram of the MRI apparatus 100 of the present embodiment.
- the MRI apparatus 100 includes a static magnetic field generation unit including a static magnetic field coil 102 that generates a static magnetic field in a space where a subject 101 is placed, an x direction, and a y direction.
- a gradient magnetic field coil 103 gradient magnetic field application unit
- a shim coil 104 that adjusts the static magnetic field distribution
- a high frequency magnetic field for the measurement region of the subject 101 includes a static magnetic field generation unit including a static magnetic field coil 102 that generates a static magnetic field in a space where a subject 101 is placed, an x direction, and a y direction.
- a gradient magnetic field coil 103 gradient magnetic field application unit
- a shim coil 104 that adjusts the static magnetic field distribution
- a high frequency magnetic field for the measurement region of the subject 101 includes a high frequency magnetic field for the measurement region of the subject 101.
- Irradiating high-frequency magnetic field irradiation coil 105 (hereinafter simply referred to as a transmission coil; a high-frequency magnetic field irradiation unit) and a nuclear magnetic resonance signal receiving coil 106 for receiving a nuclear magnetic resonance signal generated from the subject 101 (hereinafter simply referred to as a reception coil); Receiver), transmitter 107, receiver 108, calculator 109, gradient magnetic field power supply 112, shim power supply 113, and sequence control device 1 It includes a 4, a.
- the gradient magnetic field coil 103 and the shim coil 104 are driven by a gradient magnetic field power supply unit 112 and a shim power supply unit 113, respectively.
- a case where separate transmission coils 105 and reception coils 106 are used will be described as an example.
- the transmission coil 105 and the reception coil 106 are configured as a single coil. May be.
- the high-frequency magnetic field irradiated by the transmission coil 105 is generated by the transmitter 107.
- the nuclear magnetic resonance signal detected by the receiving coil 106 is sent to the computer 109 through the receiver 108.
- the sequence control device 114 is configured such that the gradient magnetic field power supply unit 112 that is a power supply for driving the gradient magnetic field coil 103, the shim power supply unit 113 that is the power supply for driving the shim coil 104, the transmitter 107, and the receiver 108 in accordance with instructions from the computer 109. And the timing of application of the gradient magnetic field and high frequency magnetic field and reception of the nuclear magnetic resonance signal are controlled.
- the control time chart is called a pulse sequence, is preset according to measurement, and is stored in a storage device or the like included in the computer 109 described later.
- the computer 109 performs various arithmetic processes on the received nuclear magnetic resonance signal to generate image information, spectrum information, temperature information, and temperature accuracy information, and gives an instruction to the sequence control device 114, and the entire MRI apparatus 100 To control the operation.
- the computer 109 is an information processing device that includes a CPU, a memory, a storage device, and the like, and a display device 110 such as a display, an external storage device 111, an input device 115, and the like are connected to the computer 109.
- the display device 110 is an interface for displaying results obtained by the arithmetic processing to the operator.
- the input device 115 is an interface for an operator to input conditions, parameters, and the like necessary for the arithmetic processing performed in the present embodiment.
- the external storage device 111 holds, together with the storage device, data used for various arithmetic processes executed by the computer 109, data obtained by the arithmetic processes, input conditions, parameters, and the like.
- the MRI apparatus 100 calculates temperature information and an index indicating its accuracy (reliability) in one measurement.
- the function of the computer 109 of this embodiment that realizes this will be described.
- FIG. 3 is a functional block diagram of the computer 109 of this embodiment.
- the computer 109 of this embodiment includes a measurement unit 210 and a calculation unit 220.
- the measurement unit 210 operates the sequence control device 114 according to the pulse sequence and controls each unit to perform measurement to obtain a nuclear magnetic resonance signal.
- the calculation unit 220 performs various calculation processes on the nuclear magnetic resonance signal obtained by the measurement, and generates image information, spectrum information, temperature information, temperature accuracy information, and the like.
- the calculation unit 220 converts the nuclear magnetic resonance signal obtained by measurement into a spectrum, calculates temperature information and temperature accuracy information, and generates display information to be displayed on the display device 110.
- the calculation unit 220 of the present embodiment includes a spectrum calculation unit 230 that converts a nuclear magnetic resonance signal obtained by measurement into a spectrum, and a temperature information calculation unit that calculates temperature information inside the subject 101 from the spectrum. 240, a temperature accuracy information calculation unit 250 that calculates the accuracy of the calculated temperature information, and a display information generation unit 260.
- the temperature accuracy information calculation unit 250 also calculates temperature accuracy information using a spectrum from which the temperature information calculation unit 240 calculates temperature information.
- the various functions realized by the computer 109 are realized by the CPU loading a program held in the storage device into the memory and executing it.
- at least one of the various functions realized by the computer 109 is an information processing apparatus independent of the MRI apparatus 100 and is realized by an information processing apparatus capable of transmitting and receiving data to and from the MRI apparatus 100. It may be.
- FIG. 4 is a processing flow of the overall flow of temperature information and temperature accuracy information calculation measurement of the present embodiment.
- spectrum information of two substances having different resonance frequencies is used.
- NAA a metabolite
- the measurement unit 210 performs non-water suppression measurement (step S1101) to obtain a nuclear magnetic resonance signal of water. Thereafter, water suppression measurement is performed (step S1102), and a nuclear magnetic resonance signal of the metabolite is obtained.
- Both non-water suppression measurement and water suppression measurement are realized by controlling the sequence controller 114 according to a predetermined pulse sequence. An example of the predetermined pulse sequence will be described later.
- the nuclear magnetic resonance signals of water and NAA are obtained by separate measurements, but the water and NAA are measured by measuring while leaving a certain amount of signal without completely suppressing the water signal. These nuclear magnetic resonance signals may be obtained simultaneously.
- the spectrum calculation unit 230 performs Fourier transform on the obtained water and NAA nuclear magnetic resonance signals to calculate water and NAA spectra (step S1103).
- the temperature information calculation part 240 calculates the temperature (temperature information) in a subject from the spectrum of water and NAA (step S1104).
- the temperature accuracy information calculation unit 250 calculates an index (temperature accuracy information) indicating the temperature accuracy from the spectrum of water and NAA (step S1105).
- the display information generation unit 260 generates display information from the calculated temperature information and temperature accuracy information and displays the display information on the display device 110 (step S1106).
- a pulse sequence of non-water suppression measurement and water suppression measurement executed by the measurement unit 210 in steps S1101 and S1102 will be described.
- a pulse sequence (hereinafter referred to as MRSI sequence) of region selective magnetic resonance spectroscopic imaging for imaging metabolites will be described as an example.
- FIG. 5 is an example of the MRSI sequence 300.
- RF indicates the application timing of the high-frequency magnetic field pulse.
- Gx, Gy, and Gz indicate application timings of gradient magnetic field pulses in the x, y, and z directions, respectively.
- a / D indicates a signal measurement period.
- the MRSI sequence 300 shown in FIG. 5 is the same as the known MRSI sequence, and selectively excites a predetermined region of interest (voxel) using one excitation pulse RF1 and two inversion pulses RF2 and RF3.
- An FID signal (free induction attenuation) FID1 is obtained from this region of interest (voxel).
- FIGS. 6A and 6B are positioning scout images obtained by measurement performed prior to the main measurement.
- FIG. 6A shows a transformer image 410
- FIG. 6B shows a sagittal image 420
- FIG. Coronal image 430 the relationship between the operation of each unit and the excited region will be described with reference to FIGS. 5 and 6A, 6B, and 6C.
- a high-frequency magnetic field RF1 and gradient magnetic field pulses Gs1 and Gs1 'in the z direction are applied to excite the cross section 401 in the z direction.
- a high frequency magnetic field RF2 and a gradient magnetic field pulse Gs2 in the y direction are applied after TE / 4 (where TE is an echo time).
- TE is an echo time
- the high frequency magnetic field RF3 and the gradient magnetic field pulse Gs3 in the x direction are applied after TE / 2 from the application of the high frequency magnetic field RF2.
- the gradient magnetic field pulses Gd1 to Gd3 and Gd1 ′ to Gd3 ′ in each direction rephase the phase of nuclear magnetization excited by the high-frequency magnetic field RF1, and dephase the phase of nuclear magnetization excited by RF2 and RF3. It is a gradient magnetic field. Further, the phase encode gradient magnetic fields Gp1 and Gp2 are applied after the high-frequency magnetic field RF3. Thus, a nuclear magnetic resonance signal of the region of interest (voxel) 404 is obtained.
- the spectrum calculation unit 230 calculates the nuclear magnetic resonance signals of water and NAA of each region of interest (voxel) 404 measured in the MRSI sequence 300 in step S1103 in the time direction. Fourier transform is performed to calculate the spectrum of water and NAA in each region of interest (voxel) 404.
- the temperature information calculation unit 240 of the present embodiment calculates the resonance frequency of water and NAA, and converts the difference between the two (resonance frequency difference) into temperature, thereby obtaining temperature information of each region of interest (voxel) 404.
- the resonance frequencies of water and NAA are obtained by fitting water and NAA spectral peaks with a predetermined function.
- FIG. 7 is a process flow of the temperature information calculation process of the present embodiment.
- the temperature information calculation unit 240 calculates the resonance frequencies of water and NAA, respectively (step S4101).
- the obtained water and NAA spectral peaks are fitted using a Lorentzian function or the like shown in Formula (1) to obtain the resonance frequency ⁇ W of water and the resonance frequency ⁇ NAA of NAA , respectively.
- the temperature information calculation unit 240 calculates the resonance frequency difference ⁇ by taking the difference between the resonance frequency ⁇ W of water and the resonance frequency ⁇ NAA of NAA (step S4102).
- the temperature information calculation unit 240 calculates temperature information by converting the frequency difference to temperature using a temperature conversion formula for converting the frequency difference to temperature (step S4103).
- the temperature conversion formula is created in advance and held in a storage device or the like.
- An example of the temperature conversion formula used in this embodiment is shown in Formula (2).
- T p ⁇ ⁇ + q (2)
- T temperature
- p a coefficient having a temperature / frequency dimension
- q is a constant term.
- the temperature accuracy information calculation unit 250 of the present embodiment determines a model function that represents each of the spectral peaks of water and NAA, and obtains a plurality of virtual temperature information by adding random noise to the model function multiple times.
- the temperature accuracy information of the temperature information of each region of interest (voxel) 404 is obtained by statistically processing it.
- the model function adds and changes a plurality of different noises equivalent to the noise obtained from the spectrum obtained by measurement. In the present embodiment, for example, the standard deviation of noise obtained from the spectrum is calculated, noise groups having the same standard deviation are generated randomly, and added to the model function.
- FIG. 8 is a processing flow of temperature accuracy information calculation processing of the present embodiment.
- the temperature accuracy information calculation unit 250 calculates water peak information and a standard deviation of water noise (step S5101).
- the water peak information of this embodiment is each parameter of a function used for fitting.
- the resonance frequency ⁇ W of water, the half-value width a W of the spectrum peak, the height I W of the spectrum peak, the phase ⁇ W , and the constant term c W It is.
- the water peak measured using the Lorentz type function described in Formula (1) is fitted, and each parameter of the used function is determined. Specifically, the resonance frequency ⁇ W of water, the half-value width a W of the spectrum peak, the height I W of the spectrum peak, the phase ⁇ W , and the constant term c W are obtained.
- the standard deviation ⁇ W of water noise is calculated using a plurality of signal values N W in a noise region (region where there is no water or metabolite signal) of the obtained water spectrum.
- the signal-to-noise ratio (SNR W ) may be further calculated.
- the water signal-to-noise ratio SNR W is calculated by I W / ⁇ W.
- the temperature accuracy information calculation unit 250 calculates metabolite (NAA) peak information and the standard deviation of NAA noise (step S5102).
- the NAA peak information of this embodiment is the same as the water peak information. Therefore, similarly to the above-described method for calculating the water peak information, the NAA peak measured using the Lorentz type function described in Equation (1) is fitted to the fitting, and each parameter of the used function is determined. Specifically, the resonance frequency ⁇ NAA of the NAA peak, the half-value width a NAA of the spectrum peak, the height I NAA of the spectrum peak, the phase ⁇ NAA , and the constant term c NAA are calculated.
- the NAA noise standard deviation ⁇ NAA is calculated by using a plurality of signal values N NAA in the noise region of the obtained NAA spectrum. At this time, the signal-to-noise cost (SNR NAA ) may also be calculated for NAA .
- SNR NAA signal-to-noise cost
- the NAA signal-to-noise ratio SNR NAA is calculated as I NAA / ⁇ NAA .
- the temperature accuracy information calculation unit 250 calculates a model function related to the frequency ⁇ using the water peak information and the NAA peak information (step S5103).
- the Lorentz function shown in Expression (1) used for fitting is used as the model function.
- the water model information L W ( ⁇ ) is calculated using the water peak information
- the NAA model function L NAA ( ⁇ ) is calculated using the NAA peak information.
- the temperature accuracy information calculation unit 250 calculates a virtual temperature from each of a plurality of virtual spectra created by adding virtual noise to the model function.
- M M is a natural number
- the temperature accuracy information calculation unit 250 sets the counter m (m is an integer from 1 to M) to 1 (step S5104).
- the temperature accuracy information calculation unit 250 calculates virtual noises N W ( ⁇ ) and N NAA ( ⁇ ) having the same standard deviation as the noise standard deviations ⁇ W and ⁇ NAA , for example, A random number is used for random determination (step S5105).
- the determined mth virtual noises N W ( ⁇ , m) and N NAA ( ⁇ , m) satisfy the following equations (3) and (4).
- var (N W ( ⁇ , m)) ⁇ W (3)
- var (N NAA ( ⁇ , m )) ⁇ NAA ⁇ (4)
- var () is an operator that calculates the standard deviation in ().
- the temperature accuracy information calculation unit 250 converts the determined noises N W ( ⁇ , m) and N NAA ( ⁇ , m) into a model function L W ( ⁇ ) and L NAA ( ⁇ ) are added to generate model functions (virtual model functions) L W ( ⁇ , m) and L NAA ( ⁇ , m) after addition of the m-th virtual noise (step S5106).
- L W ( ⁇ , m) m L W ( ⁇ ) + N W ( ⁇ , m)
- L NAA ( ⁇ , m) m L NAA ( ⁇ ) + N NAA ( ⁇ , m) (6)
- the temperature accuracy information calculation unit 250 has a virtual model function L W ( ⁇ , m) in which the water spectrum and NAA spectrum calculated by the spectrum calculation unit 230 are equal to the noise standard deviations ⁇ W and ⁇ NAA. ) And L NAA ( ⁇ , m).
- the temperature accuracy information calculation unit 250 calculates a virtual temperature from the virtual model functions L W ( ⁇ , m) and L NAA ( ⁇ , m) (step S5107).
- the calculation of the virtual temperature is the same as the temperature information calculation process by the temperature information calculation unit 240.
- the virtual model functions L W ( ⁇ , m) and L NAA ( ⁇ , m) are fitted with the Lorentzian function of Equation (1), and the respective resonance frequencies ⁇ W (m) and ⁇ NAA (m ) Is calculated. Thereafter, a resonance frequency difference ⁇ (m) is calculated, and a virtual temperature in the mth subject 101 is calculated using a predetermined conversion formula (for example, formula (2)).
- the temperature accuracy information calculation unit 250 repeats the above steps S5104 to S5106 M times (steps S5108 and S5109), and calculates M virtual temperatures.
- step S5105 virtual noises N W ( ⁇ , m) and N NAA ( ⁇ , m) having the same standard deviation as the noise standard deviations ⁇ W and ⁇ NAA are determined randomly. .
- temperature accuracy information calculating unit 250 calculates a standard deviation sigma T of M virtual temperature (step S5110).
- the temperature accuracy information calculation unit 250 of this embodiment performs temperature accuracy information (which is an index of the temperature accuracy of the temperature calculated from the water and NAA spectra obtained by the measurements in Steps S1101 and S1102).
- the standard deviation is calculated.
- the standard deviation ⁇ T of M virtual temperatures is calculated as the temperature accuracy information, but the present invention is not limited to this.
- the obtained M virtual temperatures can be statistically processed, and various values serving as temperature accuracy indexes can be calculated as temperature accuracy information.
- the temperature accuracy information may be, for example, variance or standard error.
- FIG. 9 is an explanatory diagram for explaining display information displayed on the display screen 116 of the display device 110 according to the present embodiment.
- the display information generated by the display information generation unit 260 of the present embodiment is, for example, a temperature in which the temperature information of each region of interest (voxel) 404 calculated by the temperature information calculation unit 240 is associated with a matrix measured by the MRSI sequence 300.
- the temperature accuracy information of each region of interest (voxel) 404 calculated by the temperature image 502 using the values of the table 501 and the temperature table 501 as pixel values and the temperature accuracy information calculation unit 250 is associated with the matrix measured by the MRSI sequence 300.
- the temperature accuracy table 511 and the temperature accuracy image 512 using the values of the temperature accuracy table 511 as pixel values.
- a temperature difference image may be calculated as display information.
- the temperature difference image is obtained by subtracting the temperature of the reference voxel from the temperature of a voxel other than the reference voxel with an arbitrary voxel as a reference.
- the display information displayed on the display device 110 is not limited to the above.
- the display information may include various calculation results such as water peak information and NAA peak information obtained during the calculation of temperature information and temperature accuracy information. Further, temperature information, temperature accuracy information, water peak information, and NAA peak information may be displayed for each voxel of the matrix measured by the MRSI sequence 300. Further, the temperature image and the temperature difference image may be displayed by being superimposed on various images such as an MR image, a CT image, a PET image, and a SPECT image.
- the MRI apparatus 100 includes a static magnetic field generation unit (static magnetic field coil 102) that generates a static magnetic field in a space where a subject is placed, and a high-frequency magnetic field that irradiates the subject with a high-frequency magnetic field.
- An irradiation unit transmission coil 105
- a gradient magnetic field application unit gradient magnetic field coil 103
- a reception unit reception coil 106
- a measuring unit 210 that controls operations of the gradient magnetic field applying unit, the high-frequency magnetic field irradiating unit, and the receiving unit to obtain nuclear magnetic resonance signals of two substances having different resonance frequencies, and for the nuclear magnetic resonance signal
- the magnetic resonance imaging apparatus includes a calculation unit 220 that performs calculation processing and a display device 110 that displays information after the calculation processing, and the calculation unit 220 has different resonance frequencies.
- a temperature accuracy information calculation unit 250 that calculates temperature accuracy information indicating the accuracy of the temperature information based on each spectrum peak, and generates display information to be displayed on the display device 110 based on the temperature information and the temperature accuracy information.
- a display information generation unit 260 A display information generation unit 260.
- the temperature information calculation unit 240 determines a function representing each of the calculated spectrum peaks, and based on the function, calculates a resonance frequency of each of the two substances, and the two determined resonances A conversion unit that converts the frequency difference into temperature information.
- the temperature accuracy information calculation unit 240 adds a plurality of virtual noises to each of the determined model functions, a model function calculation unit that calculates a model function representing each calculated spectrum peak, and the 2 A virtual model function generation unit that generates a plurality of virtual model functions for each substance, and a virtual temperature calculation unit that calculates a virtual temperature from each of the virtual model functions, statistically processing the plurality of virtual temperatures, Obtain temperature accuracy information.
- each of a plurality of substances having different resonance frequencies is measured once by MRS / MRSI, thereby calculating the temperature in the subject and the accuracy information of the temperature.
- MRS / MRSI MRS / MRSI
- the temperature accuracy information calculation unit 250 obtains a plurality of virtual temperatures and calculates the temperature accuracy information, but the calculation procedure of the temperature accuracy information is not limited to this.
- the temperature accuracy information may be obtained numerically using, for example, an error propagation method based on a model function. A temperature calculation method using the error propagation method will be described below.
- FIG. 10 is a processing flow of temperature accuracy information calculation processing using the error propagation method of this embodiment.
- the temperature accuracy information calculation unit 250 calculates water peak information (step S5201). This is the same as step S5101. Further, metabolite (NAA) peak information is calculated (step S5202). This is also the same as step S5102. That is, each peak information is obtained by fitting each spectrum peak with a predetermined model function and determining the coefficient of the model function.
- NAA metabolite
- the frequency accuracy information is calculated as an index of the accuracy of the peak frequency (step S5203).
- the standard deviation of the peak frequency is calculated as the frequency accuracy information.
- the error propagation method of this modification when the measured spectral data is fitted with, for example, a Lorentz type function described in Equation (1), an error occurring in each parameter of the fitting function (hereinafter referred to as a fitting parameter). Estimate the amount and calculate the standard deviation.
- the deviation included in the measurement data of the signal intensity L (f k ) is denoted by ⁇ L.
- ⁇ L the deviation included in the measurement data of the signal intensity
- F1, F2 represent those obtained by substituting the end points of the frequency domain
- f k is the frequency point
- L (B, f) is each element b i of the B model function used for the fitting to be used for fitting.
- ⁇ B follows the distribution given by equation (8) below.
- D (B, F, k) i is a noise component N defined by the equation (8).
- This represents the coefficient of (f k ).
- Each element of ⁇ B than assume that the [Delta] b i variance var ([Delta] b i) is given by the following equation (9). Therefore, the standard deviation ⁇ w1 of the resonance frequency (water peak frequency) ⁇ w of the water peak of the present embodiment can be calculated by taking i 2 as the square root of var ( ⁇ b 1 ).
- the standard deviation ⁇ w1 of the water peak frequency ⁇ w and the standard deviation ⁇ NAA1 of the NAA peak frequency ⁇ NAA can be calculated.
- step S5204 frequency difference accuracy information (standard deviation) is used as an indicator of the accuracy of the peak frequency difference as a propagation error.
- the standard deviation ⁇ df of the peak frequency difference between the water peak and the NAA peak is calculated from the standard deviation ⁇ w1 of the resonance frequency ⁇ w of the water peak and the standard deviation ⁇ NAA1 of the resonance frequency ⁇ NAA of the NAA peak by the equation (10). It can be calculated.
- temperature accuracy information is calculated from the calculated standard deviation ⁇ df of the peak frequency difference (step S5205).
- the standard deviation ⁇ T of the calculated temperature can be calculated, for example, by converting ⁇ df to temperature using the temperature conversion formula described in Non-Patent Document 2. Specifically, the absolute value
- the temperature accuracy information (in this case, standard deviation) of the temperature calculated from the spectrum information of water obtained by measurement and the spectrum information of NAA can be calculated analytically.
- the temperature accuracy information calculation unit 250 calculates the model function representing the two calculated spectrum peaks, and the resonance frequency accuracy information of each of the two substances based on the model function.
- a frequency accuracy calculation unit that calculates the accuracy information of the resonance frequency difference between the two substances from the accuracy information of the two resonance frequencies, and a temperature difference from the accuracy information of the resonance frequency difference. The accuracy information may be calculated.
- the calculated temperature accuracy information may be other statistical values such as variance and standard error.
- the display information generation unit 260 generates a temperature table, a temperature image, and a temperature difference image as display information.
- the display information to be generated is not limited to these.
- a temperature image or a temperature difference image having an arbitrary spatial resolution may be generated as display information by interpolation or the like.
- the display information generation unit 260 may generate a high-resolution temperature image as the display information by interpolating the temperature information of each voxel. For example, a temperature image with higher resolution than the original temperature image is created by interpolating temperature information of adjacent target voxels to obtain temperature information of a new intermediate voxel.
- a high-resolution temperature difference image (high-resolution temperature difference image) is created by the following procedure from the original temperature difference image.
- the flow of processing in which the display information generation unit 260 generates a temperature difference image (high resolution temperature difference image) with higher resolution than the original resolution by interpolation will be described with reference to the flowchart of FIG.
- Step S6101 the selection of the reference voxel from the temperature image (hereinafter referred to as a low-resolution temperature image) matrix (each voxel) obtained from the result measured by the MRSI sequence 300 by the above-described method is accepted and set. (Step S6101). Two or more voxels may be selected.
- the temperature of the reference voxel is calculated (step S6102).
- the temperature of the voxel is set as the reference voxel temperature.
- an average value of the temperatures of the selected voxels is calculated.
- the sum of the water peak and NAA peak of the reference voxel is calculated, and then the resonance frequency of water and NAA is calculated by fitting, and the above equation (2), etc.
- the temperature may be calculated by a conversion formula of
- a target voxel for calculating a temperature difference image
- a target voxel There may be at least one target voxel, or the entire measurement region.
- step S6104 the difference between the temperature of each voxel included in the target voxel and the temperature of the reference voxel is calculated, and a temperature difference image is calculated (step S6104).
- a high-resolution temperature difference image is generated by interpolation from the temperature difference image (step S6105).
- a known method such as cubic interpolation or spline interpolation may be used in addition to linear interpolation.
- step S6105 the high-resolution temperature difference image calculated in step S6105 is displayed on the display device 110 (step S6106).
- the display information generation unit 260 receives a reference voxel receiving unit that receives a reference voxel as a reference from the voxels, and calculates a difference between the temperature information of each voxel and the temperature information of the reference voxel.
- a temperature difference information calculation unit that calculates information, and interpolates the temperature difference information of each voxel to generate a temperature difference image with a desired resolution as the display information.
- the method for generating the high-resolution temperature difference image is not limited to this.
- the flow of processing by the display information generation unit 260 for generating a high-resolution temperature difference image according to another procedure will be described with reference to the flowchart of FIG.
- a high-resolution temperature image is generated by interpolating a low-resolution temperature image obtained from the result measured by the MRSI sequence 300 (step S6201).
- a known method such as cubic interpolation or spline interpolation may be used in addition to linear interpolation.
- selection of a region of interest (ROI: Region Of Interest) as a reference on the interpolated high-resolution temperature image is received and set (step S6202).
- ROI Region Of Interest
- the number of voxels included in the selected reference ROI may be two or more.
- the shape of the reference ROI may be an arbitrary shape other than a circle, an ellipse, or a rectangle.
- the temperature (reference temperature) within the reference ROI is calculated (step S6203).
- the temperature information of each voxel included in the region selected as the reference ROI is set as the reference temperature.
- an average value of the temperatures of the selected voxels is calculated as a reference temperature.
- the selection of the ROI (herein referred to as the target ROI) for calculating the temperature difference image on the high resolution temperature image is received and set (step S6204).
- the target ROI There may be at least one target ROI. Moreover, the whole measurement area
- step S6205 a difference between the temperature of each voxel included in the target ROI and the reference temperature is calculated, and a high-resolution temperature difference image having the result as a pixel value is calculated.
- step S6205 the high-resolution temperature difference image calculated in step S6205 is displayed on the display device 110 (step S6206).
- the display information generation unit 260 interpolates the temperature information of each voxel to generate a temperature image with a desired resolution, and a reference region of interest on the temperature image with the desired resolution.
- a reference region-of-interest receiving unit that receives a selection; a reference temperature information calculating unit that calculates temperature information of the voxel in the reference region of interest as reference temperature information; and temperature information of each voxel of the temperature image of the desired resolution
- a temperature difference information calculation unit that calculates a difference between the reference temperature information and the reference temperature information as temperature difference information, and generates a temperature difference image using the temperature difference information of each voxel as a pixel value as the display information.
- high-resolution water and NAA spectrum data are calculated by spatially interpolating the water and NAA spectrum data obtained from the results measured by the MRSI sequence 300 (step S6301).
- a known method such as cubic interpolation or spline interpolation may be used in addition to linear interpolation.
- temperature information is calculated from the high-resolution water, NAA, and spectrum data by the method of the present embodiment, and a high-resolution temperature image (high-resolution temperature image) is generated (step S6302).
- the temperature information calculation process may be configured such that the display information generation unit 260 causes the temperature information calculation unit 240 to calculate.
- a region of interest ROI serving as a reference on the high resolution temperature image is received and set (step S6303).
- the number of voxels included in the reference ROI may be two or more.
- the shape of the reference ROI may be an arbitrary shape other than a circle, an ellipse, or a rectangle.
- the temperature (reference temperature) within the reference ROI is calculated (step S6304).
- the temperature information of each voxel included in the region selected as the reference ROI is set as the reference temperature.
- an average value of the temperatures of the selected voxels is calculated as a reference temperature.
- a target ROI selection and setting of an ROI for calculating a temperature difference image on the high resolution temperature image (referred to here as a target ROI) is received and set (step S6305).
- a target ROI There may be at least one target ROI. Moreover, the whole measurement area
- step S6306 a difference between the temperature of each voxel included in the target ROI and the reference temperature is calculated, and a high-resolution temperature difference image having the result as a pixel value is calculated.
- step S6306 the high-resolution temperature difference image calculated in step S6306 is displayed on the display device 110 (step S6307).
- the display information generation unit 260 spatially interpolates the spectrum of each voxel measured in the sequence to obtain spectrum data with a desired resolution, and the interpolated voxel from the interpolated spectrum data.
- a temperature image generation unit that generates a temperature image with the resolution using the calculated temperature information as pixel values, and a reference region of interest that receives a selection of a region of interest serving as a reference on the temperature image with the resolution
- a receiving unit a reference temperature information calculating unit that calculates the temperature information of each interpolated voxel in the reference region of interest as reference temperature information, and a voxel of each interpolated voxel in the temperature image of the desired resolution.
- a temperature difference information calculation unit that calculates a difference between temperature information and the reference temperature information as temperature difference information, and the temperature difference information of the voxel after each interpolation The generating a temperature difference image to the pixel value as the display information.
- an ROI for calculating a temperature image is selected on the high-resolution temperature image generated in step S6302, and each of the included ROIs included in the target ROI is selected.
- the voxel temperature be the pixel value of each pixel of the high-resolution temperature image.
- the display information generation unit 260 includes an interpolation unit that spatially interpolates the spectrum for each voxel measured in the sequence to obtain spectrum data with a desired resolution, and each of the interpolated spectrum data is subjected to interpolation. Voxel temperature information is calculated, and a temperature image having the resolution with the calculated temperature information as a pixel value is generated as the display information.
- the display information displayed on the display device 110 may be provided with a switch 520 that switches execution / non-execution of the interpolation process.
- the switch 520 is pressed, the display information generation unit 260 performs an interpolation process.
- the temperature accuracy information obtained by the method of the present embodiment may be used for determining measurement parameters. For example, a plurality of MR images are acquired by changing measurement parameters for one subject 101 in advance. Then, a plurality of pieces of temperature accuracy information are obtained by the above method using the plurality of acquired MR images. At this time, the obtained temperature accuracy information is stored in the storage device in association with the measurement parameter. At the time of measurement, measurement parameters that can obtain the best temperature accuracy under the constraint conditions such as measurement time are employed.
- the temperature accuracy information may be calculated for each VOI of various sizes at various positions, and the temperature accuracy information may be held in association with the positions.
- a histogram of pixel values of each voxel constituting the VOI for which the temperature accuracy information is calculated, instead of the measurement parameter and the measurement position, and the temperature accuracy information may be stored in association with each other.
- the pulse sequence used in the non-water suppression measurement and the water suppression measurement has been described as an example.
- the pulse sequence used in the non-water suppression measurement and the water suppression measurement is described here. Not limited. Any pulse sequence that can obtain the spectrum of the substance to be measured of each voxel may be used.
- a sequence for measuring a single region called a magnetic resonance spectrum copy pulse sequence (MRS sequence), a high-speed MRSI sequence using an oscillating gradient magnetic field called an echo planar spectroscopic image sequence (EPSI sequence), etc. It may be.
- MRS sequence magnetic resonance spectrum copy pulse sequence
- EPSI sequence echo planar spectroscopic image sequence
- DESCRIPTION OF SYMBOLS 100 MRI apparatus, 101: Subject, 102: Static magnetic field coil, 103: Gradient magnetic field coil, 104: Shim coil, 105: Transmission coil, 106: Reception coil, 107: Transmitter, 108: Receiver, 109: Calculator, 110: Display device, 111: External storage device, 112: Power supply unit for gradient magnetic field, 113: Power supply unit for shim, 114: Sequence control device, 115: Input device, 116: Display screen, 120: MRI device, 130: MRI Apparatus: 14: claim, 210: measurement unit, 220: calculation unit, 230: spectrum calculation unit, 240: temperature information calculation unit, 250: temperature accuracy information calculation unit, 260: display information generation unit, 300: MRSI sequence, 401: Cross section, 402: Cross section, 403: Cross section, 404: Voxel 410: Trans image, 420: Sagittal image, 430: Ronaru image, 501: temperature table, 50
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Radiology & Medical Imaging (AREA)
- Signal Processing (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
Abstract
Description
T=p×Δν+q (2)
ここで、Tは温度、pは温度/周波数の次元を持つ係数、qは定数項である。式(2)のpおよびqは、文献に記された公知の値や、実験的に求めた値を用いる。
var(NW(ν,m))=σW・・・(3)
var(NNAA(ν,m))=σNAA・・・(4)
なお、var()は()内の標準偏差を算出する演算子である。
LW(ν,m)m=LW(ν)+NW(ν,m)・・・(5)
LNAA(ν,m)m=LNAA(ν)+NNAA(ν,m)・・・(6)
Claims (14)
- 被検体が置かれる空間に静磁場を発生させる静磁場発生部と、前記被検体に高周波磁場を照射する高周波磁場照射部と、前記被検体に傾斜磁場を印加する傾斜磁場印加部と、前記被検体から発生する核磁気共鳴信号を受信する受信部と、前記傾斜磁場印加部と前記高周波磁場照射部と前記受信部との動作を制御して共鳴周波数が異なる2つの物質の核磁気共鳴信号を得る計測部と、前記核磁気共鳴信号に対する演算処理を行う演算部と、前記演算処理後の情報を表示する表示装置と、を備える磁気共鳴撮影装置であって、
前記演算部は、
前記共鳴周波数が異なる2つの物質それぞれの核磁気共鳴信号のスペクトルを各々算出するスペクトル算出部と、
前記算出した各々のスペクトルピークに基づいて前記被検体内の温度情報を算出する温度情報算出部と、
前記算出した各々のスペクトルピークに基づいて前記温度情報の精度を示す温度精度情報を算出する温度精度情報算出部と、
前記温度情報および前記温度精度情報に基づいて前記表示装置に表示する表示情報を生成する表示情報生成部と、を備えること
を特徴とする磁気共鳴撮影装置。 - 請求項1記載の磁気共鳴撮影装置であって、
前記温度情報算出部は、
前記算出した各々のスペクトルピークを表す関数を決定し、当該関数に基づき、前記2つの物質それぞれの共鳴周波数を算出する共鳴周波数算出部と、
前記決定した2つの共鳴周波数の差を温度情報に換算する換算部と、を備える
を特徴とする磁気共鳴撮影装置。 - 請求項1記載の磁気共鳴撮影装置であって、
前記温度精度情報算出部は、
前記算出した各々のスペクトルピークを表すモデル関数を算出するモデル関数算出部と、
前記モデル関数各々に複数の仮想的なノイズを付加し、前記2つの物質毎に複数の仮想モデル関数を生成する仮想モデル関数生成部と、
前記仮想モデル関数各々から仮想温度を算出する仮想温度算出部と、を備え、
前記複数の仮想温度を統計処理し、前記温度精度情報を得ること
を特徴とする磁気共鳴撮影装置。 - 請求項1記載の磁気共鳴撮影装置であって、
前記温度精度情報算出部は、
前記算出した各々のスペクトルピークを表すモデル関数を算出するモデル関数算出部と、
前記モデル関数に基づき、前記2つの物質それぞれの共鳴周波数の精度情報を算出する周波数精度算出部と、
前記2つの共鳴周波数の精度情報から2つの物質の共鳴周波数差の精度情報を算出する周波数差精度算出部と、を備え、
前記共鳴周波数差の精度情報から温度精度情報を算出すること
を特徴とする磁気共鳴撮影装置。 - 請求項1記載の磁気共鳴撮影装置であって、
前記共鳴周波数が異なる2つの物質は、水と代謝物であること
を特徴とする磁気共鳴撮影装置。 - 請求項1から5いずれか1項記載の磁気共鳴撮影装置であって、
前記計測部は、磁気共鳴スペクトロスコピーシーケンスに従って、前記各物質の磁気共鳴信号を得ること
を特徴とする磁気共鳴撮影装置。 - 請求項1から5いずれか1項記載の磁気共鳴撮影装置であって、
前記計測部は、磁気共鳴スペクトロスコピックイメージングシーケンスおよびエコープラナースペクトロスコピックイメージングシーケンスのいずれか一方のシーケンスに従って、前記各物質の磁気共鳴信号を得ること
を特徴とする磁気共鳴撮影装置。 - 請求項7記載の磁気共鳴撮影装置であって、
前記温度情報算出部は、前記温度情報を、前記シーケンスで計測したボクセル毎に算出し、
前記表示情報生成部は、前記各ボクセルの温度情報を補間することにより高分解能の温度画像を前記表示情報として生成すること
を特徴とする磁気共鳴撮像装置。 - 請求項7記載の磁気共鳴撮影装置であって、
前記温度情報算出部は、前記温度情報を、前記シーケンスで計測したボクセル毎に算出し、
前記表示情報生成部は、
前記各ボクセルの中から基準とする基準ボクセルの選択を受け付ける基準ボクセル受付部と、
前記各ボクセルの温度情報と前記基準ボクセルの温度情報との差分を温度差情報として算出する温度差情報算出部と、を備え、
前記各ボクセルの温度差情報を補間することにより所望の分解能の温度差画像を前記表示情報として生成すること
を特徴とする磁気共鳴撮影装置。 - 請求項7記載の磁気共鳴撮影装置であって、
前記温度情報算出部は、前記温度情報を、前記シーケンスで計測したボクセル毎に算出し、
前記表示情報生成部は、
前記各ボクセルの温度情報を補間することにより所望の分解能の温度画像を生成する補間部と、
前記所望の分解能の温度画像上で、基準となる関心領域の選択を受け付ける基準関心領域受付部と、
前記基準となる関心領域内の前記ボクセルの温度情報を基準温度情報として算出する基準温度情報算出部と、
前記所望の分解能の温度画像の各ボクセルの温度情報と前記基準温度情報との差分を温度差情報として算出する温度差情報算出部と、を備え、
前記各ボクセルの温度差情報を画素値とする温度差画像を、前記表示情報として生成すること
を特徴とする磁気共鳴撮像装置。 - 請求項7記載の磁気共鳴撮影装置であって、
前記表示情報生成部は、
前記シーケンスで計測したボクセル毎のスペクトルを空間的に補間して所望の分解能のスペクトルデータを得る補間部を備え、
前記補間後のスペクトルデータからそれぞれ補間後のボクセルの温度情報を算出し、算出した各温度情報を画素値とする前記分解能の温度画像を、前記表示情報として生成すること
を特徴とする磁気共鳴撮像装置。 - 請求項7記載の磁気共鳴撮影装置であって、
前記表示情報生成部は、
前記シーケンスで計測したボクセル毎のスペクトルを空間的に補間して所望の分解能のスペクトルデータを得る補間部と、
前記補間後のスペクトルデータからそれぞれ補間後のボクセルの温度情報を算出し、算出した各温度情報を画素値とする前記分解能の温度画像を生成する温度画像生成部と、
前記分解能の温度画像上で基準となる関心領域の選択を受け付ける基準関心領域受付部と、
前記基準となる関心領域内の各補間後のボクセルの前記温度情報を基準温度情報として算出する基準温度情報算出部と、
前記所望の分解能の温度画像内の各補間後のボクセルの温度情報と前記基準温度情報との差分を温度差情報として算出する温度差情報算出部と、を備え、
前記各補間後のボクセルの温度差情報を画素値とする温度差画像を前記表示情報として生成すること
を特徴とする磁気共鳴撮像装置。 - 請求項2記載の磁気共鳴撮影装置であって、
前記関数は、前記共鳴周波数を変数に持つ、ローレンツ型関数であること
を特徴とする磁気共鳴撮影装置。 - 請求項3記載の磁気共鳴撮影装置であって、
前記温度精度情報は、前記仮想温度の標準偏差であること
を特徴とする磁気共鳴撮影装置。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013556292A JP5797782B2 (ja) | 2012-01-30 | 2013-01-10 | 磁気共鳴撮影装置 |
US14/375,488 US9851424B2 (en) | 2012-01-30 | 2013-01-10 | Magnetic resonance imaging apparatus |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012016781 | 2012-01-30 | ||
JP2012-016781 | 2012-01-30 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2013114927A1 true WO2013114927A1 (ja) | 2013-08-08 |
Family
ID=48904971
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2013/050290 WO2013114927A1 (ja) | 2012-01-30 | 2013-01-10 | 磁気共鳴撮影装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9851424B2 (ja) |
JP (1) | JP5797782B2 (ja) |
WO (1) | WO2013114927A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017064175A (ja) * | 2015-09-30 | 2017-04-06 | 株式会社日立製作所 | 磁気共鳴イメージング装置、および、画像処理装置 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102013207390B4 (de) * | 2013-04-24 | 2015-03-05 | Siemens Aktiengesellschaft | Ermittlung einer Magnetresonanz-Ansteuersequenz mit einer automatischen Anpassung einer ersten und einer zweiten Selektionsrichtung |
DE102015200695B4 (de) * | 2015-01-19 | 2016-08-18 | Siemens Healthcare Gmbh | Generieren von Steuerinformationen für eine Magnetresonanz-Bildgebung unter Verwendung mehrerer Frequenzspektren von verschiedenen Spulenelementen |
US11176717B2 (en) * | 2019-09-26 | 2021-11-16 | Siemens Healthcare Gmbh | Guiding protocol development for magnetic resonance thermometry |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11253420A (ja) * | 1998-03-11 | 1999-09-21 | Technol Res Assoc Of Medical & Welfare Apparatus | 磁気共鳴診断装置 |
JP2004527351A (ja) * | 2001-05-30 | 2004-09-09 | インサイテック−テキソニックス リミテッド | 誤差補正を用いるmriベースの温度マッピング |
JP2005304558A (ja) * | 2004-04-16 | 2005-11-04 | Toshiba Corp | 磁気共鳴映像診断システム |
JP2008086525A (ja) * | 2006-09-29 | 2008-04-17 | Shiga Univ Of Medical Science | 小型マイクロ波温熱治療システム |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2619603A1 (en) * | 2010-09-21 | 2013-07-31 | Insightec Ltd. | Magnetic resonance thermometry using prf spectroscopy |
US20160192859A1 (en) * | 2013-09-09 | 2016-07-07 | Hitachi Medical Corporation | Magnetic resonance imaging apparatus and temperature information measurement method |
-
2013
- 2013-01-10 WO PCT/JP2013/050290 patent/WO2013114927A1/ja active Application Filing
- 2013-01-10 JP JP2013556292A patent/JP5797782B2/ja active Active
- 2013-01-10 US US14/375,488 patent/US9851424B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11253420A (ja) * | 1998-03-11 | 1999-09-21 | Technol Res Assoc Of Medical & Welfare Apparatus | 磁気共鳴診断装置 |
JP2004527351A (ja) * | 2001-05-30 | 2004-09-09 | インサイテック−テキソニックス リミテッド | 誤差補正を用いるmriベースの温度マッピング |
JP2005304558A (ja) * | 2004-04-16 | 2005-11-04 | Toshiba Corp | 磁気共鳴映像診断システム |
JP2008086525A (ja) * | 2006-09-29 | 2008-04-17 | Shiga Univ Of Medical Science | 小型マイクロ波温熱治療システム |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017064175A (ja) * | 2015-09-30 | 2017-04-06 | 株式会社日立製作所 | 磁気共鳴イメージング装置、および、画像処理装置 |
WO2017056996A1 (ja) * | 2015-09-30 | 2017-04-06 | 株式会社日立製作所 | 磁気共鳴イメージング装置、および、画像処理装置 |
Also Published As
Publication number | Publication date |
---|---|
US9851424B2 (en) | 2017-12-26 |
JP5797782B2 (ja) | 2015-10-21 |
JPWO2013114927A1 (ja) | 2015-05-11 |
US20150008925A1 (en) | 2015-01-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5449903B2 (ja) | 磁気共鳴イメージング装置 | |
JP6014770B2 (ja) | 磁気共鳴撮影装置および温度情報計測方法 | |
JP5221570B2 (ja) | 磁気共鳴イメージング装置及びマルチコントラスト画像取得方法 | |
US8754645B2 (en) | Method for spatially resolved determination of an MR parameter | |
US10215828B2 (en) | Magnetic resonance imaging apparatus | |
US20150338492A1 (en) | Magnetic resonance imaging device and quantitative susceptibility mapping method | |
CN107072586B (zh) | 磁共振成像装置 | |
JP6371554B2 (ja) | 磁気共鳴イメージング装置 | |
Hutter et al. | Slice-level diffusion encoding for motion and distortion correction | |
JP2022041978A (ja) | 磁気共鳴イメージング装置及び医用データ処理装置 | |
US10203387B2 (en) | MR imaging with enhanced susceptibility contrast | |
JP5797782B2 (ja) | 磁気共鳴撮影装置 | |
JP7115889B2 (ja) | 医用画像診断支援装置、および、磁気共鳴イメージング装置 | |
US10908235B2 (en) | Method of fast imaging of NMR parameters with variably-accelerated sensitivity encoding | |
JP5278914B2 (ja) | 磁気共鳴イメージング装置及び磁化率強調画像撮影方法 | |
EP2757386A1 (en) | Synthetic parameter-weigthed images based on relaxometry (SPARE) | |
JP7161377B2 (ja) | 画像処理装置、及び、それを含む磁気共鳴イメージング装置及び磁気共鳴イメージングシステム | |
JP6843706B2 (ja) | 磁気共鳴イメージング装置及び拡散強調画像の補正方法 | |
JP5636058B2 (ja) | 磁気共鳴撮影装置 | |
JP5684888B2 (ja) | 磁気共鳴イメージング装置 | |
JP6267684B2 (ja) | 磁気共鳴イメージング装置 | |
JP6487554B2 (ja) | 磁気共鳴イメージング装置 | |
JP6776217B2 (ja) | 画像処理装置、画像処理方法、画像処理プログラム及び磁気共鳴イメージング装置 | |
JP5881869B2 (ja) | 磁気共鳴イメージング装置 | |
JP6114846B2 (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: 13743643 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2013556292 Country of ref document: JP Kind code of ref document: A |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWE | Wipo information: entry into national phase |
Ref document number: 14375488 Country of ref document: US |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 13743643 Country of ref document: EP Kind code of ref document: A1 |