US20180028273A1 - Magnetic resonance imaging apparatus and quantitative-value computing program - Google Patents
Magnetic resonance imaging apparatus and quantitative-value computing program Download PDFInfo
- Publication number
- US20180028273A1 US20180028273A1 US15/614,782 US201715614782A US2018028273A1 US 20180028273 A1 US20180028273 A1 US 20180028273A1 US 201715614782 A US201715614782 A US 201715614782A US 2018028273 A1 US2018028273 A1 US 2018028273A1
- Authority
- US
- United States
- Prior art keywords
- value
- quantitative
- imaging
- pixel
- initial
- 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.)
- Granted
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B90/00—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
- A61B90/36—Image-producing devices or illumination devices not otherwise provided for
- A61B90/37—Surgical systems with images on a monitor during operation
-
- 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
-
- 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
-
- 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/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
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
- H04N19/10—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding
- H04N19/102—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding characterised by the element, parameter or selection affected or controlled by the adaptive coding
- H04N19/103—Selection of coding mode or of prediction mode
- H04N19/105—Selection of the reference unit for prediction within a chosen coding or prediction mode, e.g. adaptive choice of position and number of pixels used for prediction
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B90/00—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
- A61B90/36—Image-producing devices or illumination devices not otherwise provided for
- A61B90/37—Surgical systems with images on a monitor during operation
- A61B2090/374—NMR or MRI
-
- 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
Definitions
- the present invention relates to a technology that generates a quantitative-value image in response to a signal acquired from a magnetic resonance imaging apparatus.
- a magnetic resonance imaging (hereinafter, MRI) apparatus is a medical diagnostic imaging apparatus that mainly uses a nuclear magnetic resonance phenomenon of protons.
- the MRI apparatus is capable of non-invasively imaging any cross section and is capable of acquiring information related to vital functions such as bloodstream or metabolic functions, in addition to morphological information.
- a slice gradient magnetic field is applied to a subject positioned in a magnetostatic field and, simultaneously, a high-frequency magnetic field having a specific frequency is applied thereto, and thus nuclear magnetization is excited in a cross section that needs to be imaged.
- a phase-encoding gradient magnetic field and a lead-out gradient magnetic field are applied to the excited nuclear magnetization, thereby plane position information is assigned, and a nuclear magnetic resonance signal (echo) generated from the nuclear magnetization is measured.
- the measurement of the nuclear magnetic resonance signal is iterated until a measurement space referred to as a k space is filled with the signals.
- the signals, with which the k space is filled, are converted into an image through inverse Fourier transform.
- Pulses and gradient magnetic fields for generating an echo are applied, based on a preset pulse sequence.
- Various types of pulse sequences are known according to purposes.
- the phase-encoding gradient magnetic field is sequentially changed for every iteration time (hereinafter, TR) of the pulse sequence, and the number of nuclear magnetic resonance signals that is required for obtaining one tomogram is measured.
- TR iteration time
- a user selects and executes a pulse sequence, thereby acquiring a weighted image in which a relative difference in quantitative values of living tissue (referred to as subject parameters, for example, T1: longitudinal relaxation time, T2: transverse relaxation time, PD: proton density, or D: diffusion coefficient) is weighted.
- subject parameters for example, T1: longitudinal relaxation time, T2: transverse relaxation time, PD: proton density, or D: diffusion coefficient
- an imaging site that needs to be subjected to diagnosis using an image acquired by the imaging for positioning is set, the pulse sequence or the imaging parameter is changed, and a plurality of image types (for example, a T1 weighted image, a T2 weighted image, FLAIR: fluid attenuated inversion recovery, a diffusion weighted image, and MRA: magnetic resonance angiography) are acquired.
- the user adjusts a window level (WL) or a window width (WW) of an image acquired through a manual operation and performs clipping of a signal that interferes with a diagnosis as necessary, thereby generating an image for diagnosis.
- an object of the invention is to provide an MRI apparatus that is capable of reducing calculating costs in quantitative value estimating calculation, and to rapidly determine whether the imaging is successful or failed.
- an initial value of a quantitative value is set with reference to a pixel value obtained by performing calculation from a quantitative-value candidate group, and the quantitative value is computed through a localized optimization technique.
- the MRI apparatus of the invention includes: an imaging unit that acquires a nuclear magnetic resonance signal from a subject and performs imaging to acquire image data of the subject; a measurement control unit that controls the imaging unit such that the imaging unit performs imaging a plurality of times with different imaging parameter values in the same pulse sequence; a quantitative-value computing unit that computes a quantitative value which does not depend on the imaging parameter values, from a plurality of images having different pixel values obtained by performing the imaging the plurality of times; and a predicted-pixel-value computing unit that predicts a pixel value that is acquired from the plurality of imaging parameter values, for each of a plurality of predetermined quantitative-value candidate groups.
- the quantitative-value computing unit includes an initial-value selecting unit that selects an initial value of a quantitative value from the plurality of quantitative-value candidate groups with reference to the pixel value computed in the predicted-pixel-value computing unit, and an optimal-value computing unit that computes a quantitative value through a localized optimization technique using the initial value selected in the initial-value selecting unit.
- FIG. 1 is a functional block diagram of a calculator of a first embodiment.
- FIG. 2 is a diagram illustrating an example of flow of the calculator of the first embodiment.
- FIG. 3 is a diagram schematically illustrating a pixel-value predicting process of the first embodiment.
- FIG. 4 is a block diagram illustrating a common configuration of an MRI apparatus of a second embodiment.
- FIG. 5A illustrates an example of a functional block diagram of a calculator of the second embodiment.
- FIG. 5B illustrates another example of the functional block diagram of the calculator of the second embodiment.
- FIG. 6 illustrates an example of flow of the calculator of the second embodiment.
- FIG. 7 illustrates an example of a pulse sequence that is employed in the MRI apparatus of the embodiment.
- FIG. 8 illustrates effects of the second embodiment.
- FIG. 9 is a functional block diagram of a calculator of a third embodiment.
- FIG. 10 illustrates an example of flow of the calculator of the third embodiment.
- FIG. 11 is a functional block diagram of a calculator of a fourth embodiment.
- FIG. 12 illustrates an example of flow of the calculator of the fourth embodiment.
- FIG. 13 illustrates an example of presentation of an initial value in the fourth embodiment.
- FIG. 14 is a functional block diagram of a calculator of a fifth embodiment.
- FIG. 15 illustrates an example of flow of the calculator of the fifth embodiment.
- FIG. 16 illustrates an example of a user interface of the fifth embodiment.
- the embodiment is characterized by a calculator that performs quantitative value calculation using a plurality of images having different pixel values, which are obtained by performing imaging a plurality of times with different imaging parameter values in the same pulse sequence, the calculator including: a predicted-pixel-value computing unit that predicts a predicted value of a pixel value that is acquired from the imaging parameters in the imaging performed the plurality of times; an initial-value selecting unit that selects a quantitative-value candidate group which is used as initial values in quantitative value calculation from a plurality of predetermined quantitative-value candidate groups with reference to the pixel value predicted by the predicted-pixel-value computing unit; and an optimal-value computing unit that computes a quantitative value through a localized optimization technique using the initial value selected by the initial-value selecting unit.
- the embodiment is characterized by a program that is executed for performing the quantitative-value calculation and is installed in the calculator.
- the calculator may be incorporated in an MRI apparatus, or the calculator may be independently provided from the MRI apparatus.
- the “quantitative value” means a quantitative value that does not depend on a value of the imaging parameter, that is, a parameter depending on a subject which determines a pixel value of an MR image.
- the quantitative value includes a subject parameter such as a transverse relaxation time, a longitudinal relaxation time, spin density, a diffusion coefficient, or RF irradiation sensitivity, a proportionality coefficient of a signal function of a pulse sequence, or the like.
- a quantitative-value image means an image having, as a pixel value, any quantitative value of the quantitative values.
- a calculator 100 has a function of computing the quantitative value and includes: a predicted-pixel-value computing unit 220 that uses imaging parameters used in practical imaging and predicts a pixel value that is acquired from imaging parameters; an initial-value selecting unit 231 that uses a plurality of image data 310 acquired from the MRI apparatus and selects, as an initial value of the quantitative value computation, a predetermined candidate group from a plurality of predetermined quantitative-value candidate groups 320 ; and an optimal-value computing unit 232 that computes the quantitative value by using the initial value selected by the initial-value selecting unit 211 .
- functions of the initial-value selecting unit 231 and the optimal-value computing unit 232 are collectively referred to as a quantitative-value computing unit 230 .
- a function of the calculator 100 is mainly achieved by executing the program incorporated in a CPU.
- a part or the entirety of the quantitative-value computation may be performed by hardware such as an ASIC or an FPGA.
- the plurality of image data 310 and quantitative-value candidate groups 320 used by predicted-pixel-value computing unit 220 and the quantitative-value computing unit 230 are stored in a storage device 300 in advance.
- the storage device 300 may be an internal storage device in the MRI apparatus or an external storage device assembled to the MRI apparatus, or a storage device independently provided from the MRI apparatus.
- the image data 310 is data of an image that is imaged by the MRI apparatus in a pulse sequence and an imaging condition (imaging parameter) depending on a type of quantitative value as a purpose of the computation and is formed in response to an acquired signal, and contains a plurality of image data obtained by changing combinations of the imaging parameter values.
- the image data may be data obtained after image reconstruction is performed or data obtained before the image reconstruction is performed.
- the quantitative-value candidate groups 320 is a combination of predetermined quantitative values from a table in which a plurality of values are listed in associated with a plurality of types of quantitative values, respectively, and thus the plurality of values can be employed from values obtained in an empirical manner or literature.
- the table can have values different for each imaging site.
- Step S 101 Flow of a process of the calculator (quantitative-value computing unit 230 ) 100 of the embodiment is described with reference to FIG. 2 , based on such a configuration described above.
- the plurality of quantitative-value candidate groups 320 are acquired from the list of the quantitative values (Step S 101 ).
- the predicted-pixel-value computing unit 220 uses the imaging parameter values and the plurality of quantitative-value candidate groups 320 acquired in Step S 101 , and predicts a pixel value obtained in a case where imaging is performed with the imaging parameters (Step S 102 ).
- FIG. 3 schematically illustrates the prediction of the pixel value.
- FIG. 3 illustrates a case of computing a certain pixel value.
- the predicted-pixel-value computing unit 220 substitutes each of the L imaging parameters and the N quantitative-value candidate groups into a signal function, computes pixel values with respect to the L imaging parameters for each quantitative-value candidate group, and computes L pixel values.
- the signal function is a function generated for each pulse sequence, and a function that restores signal intensity of each pixel with the quantitative values and the imaging parameters as variables.
- the initial-value selecting unit 211 selects one quantitative-value candidate group from the plurality of quantitative-value candidate groups with reference to the predicted pixel value (Step S 103 ). Specifically, the predicted pixel value for each quantitative-value candidate group and a pixel value of the image data 310 that is practically imaged are matched, and the most matched quantitative-value candidate groups are selected.
- the optimal-value computing unit 232 starts the quantitative value computation using a localized optimization algorithm with the candidate group as initial values (Step S 104 ). Specifically, the pixel value of the image data 310 is fitted and the quantitative value is estimated as a variable of the signal function.
- the signal function is the same as the signal function used in Step S 102 of predicting the pixel value.
- the number of calculable quantitative values is smaller than or equal to the number of combinations (L sets in FIG. 3 ) having different imaging parameter values.
- a quantitative value as a value that does not depend on the imaging parameters is computed for each pixel.
- a quantitative value image (or, referred to as a subject parameter map) is obtained (Step S 105 ).
- the predicted pixel value is calculated for each imaging parameter, the initial value of the quantitative value is set with reference to the predicted pixel value, and thereby it is possible to compute the quantitative value with high accuracy by applying the localized optimization technique without resulting in a localized solution. In this manner, it is possible to significantly reduce calculation costs according to the quantitative-value computation.
- the MRI apparatus of the embodiment includes the imaging unit, a measurement control unit that controls the imaging unit and measures data that is required for the quantitative value computation, and the calculator that performs the quantitative value computation.
- the measurement control unit performs imaging a plurality of times with different imaging parameters in the same pulse sequence, and acquires a plurality of images (image data) having different pixel values.
- the calculator includes the quantitative-value computing unit that performs the quantitative value computation by using information of the plurality of image data acquired through imaging and the predetermined quantitative-value candidate groups.
- FIG. 4 is a block diagram illustrating a common configuration of an MRI apparatus 10 .
- the MRI apparatus 10 as the imaging unit, includes a magnet 101 that generates a magnetostatic field, a gradient coil 102 that generates a gradient magnetic field, an RF coil 107 that irradiates a subject (living body) 103 with a high-frequency magnetic pulse (hereinafter, RF pulse), an RF probe 108 that detects an echo signal that is generated from the subject 103 , and a bed (table) 115 on which the subject (for example, living body) 103 is mounted in a space of the magnetostatic field generated by the magnet 101 .
- RF pulse high-frequency magnetic pulse
- the MRI apparatus 10 (imaging unit) includes a gradient-magnetic-field power supply 105 that drives the gradient coil 102 , a high-frequency magnetic field generator 106 that drives the RF coil 107 , and a receiver 109 that receives an echo signal detected by the RF probe 108 .
- the MRI apparatus 10 may further include a shim coil 113 and a shim power supply 114 that drives the shim coil 113 .
- the shim coil 113 is formed to have a plurality of channels, and an additional magnetic field that corrects magnetostatic inhomogeneity with a current supplied from the shim power supply 114 .
- the currents flowing through the channels that configure the shim coil 113 during the adjusting of the magnetostatic homogeneity is controlled by a sequencer 104 which will be described below.
- the MRI apparatus 10 as a computation control system, includes the sequencer 104 that transmits commands to the gradient-magnetic-field power supply 105 and the high-frequency magnetic field generator 106 such that the gradient-magnetic-field power supply and the high-frequency magnetic field generator generate a gradient magnetic field and a high-frequency magnetic field, respectively, and sets a nuclear magnetic resonance frequency as a reference of detection to the receiver 109 , a calculator 110 that performs signal processing with respect to a detected signal, and performs control or the like of an entire operation of the MRI apparatus 10 , a display device 111 that displays a processing result from the calculator 110 , a storage device 112 that stores the processing result, and an input device 116 that receives instructions from a user.
- the sequencer 104 that transmits commands to the gradient-magnetic-field power supply 105 and the high-frequency magnetic field generator 106 such that the gradient-magnetic-field power supply and the high-frequency magnetic field generator generate a gradient magnetic field and a high-frequency magnetic field, respectively
- the display device 111 and the input device 116 are disposed to be close to each other or as an integrated console, and may function as a user interface.
- Various types of data that are required in the calculator 110 are stored in the storage device 112 .
- the storage device 112 not only the result of the signal processing, but also the detected signal, imaging conditions, or the like may be stored in the storage device 112 as necessary.
- control of the sequencer 104 causes the RF pulse to be applied to the subject 103 through the RF coil 107 , and the gradient magnetic pulse for assigning position information such as slice selection or phase encoding to the echo signal to be applied to the subject by the gradient coil 102 .
- a signal generated from the subject 103 is received by the RF probe 108 , the detected signal is transmitted to the calculator 110 , and the signal processing such as image reconstruction is performed on the signal.
- the calculator 110 performs not only the signal processing in which the received signal is processed, but also outputs an instruction to the sequencer 104 such that components operate at a timing and intensity which are programmed in advance, controls the operations of components that configure the imaging unit, and performs the measurement.
- an instruction to the sequencer 104 such that components operate at a timing and intensity which are programmed in advance, controls the operations of components that configure the imaging unit, and performs the measurement.
- a program in which the high-frequency magnetic field, the gradient magnetic field, the timing and intensity of signal reception are set, is referred to as the pulse sequence, and various pulse sequences are prepared in accordance with an imaging method or an imaging purpose. Measurement is performed depending on imaging parameters that are required for the pulse sequence and controlling the pulse sequence.
- the pulse sequence is generated in advance and is stored in the storage device 112 , the imaging parameters mainly include a flip angle (FA) of the RF pulse, iterating time (TR), echo time (TE), or the like, and values of the imaging parameters are input via a user interface.
- FA flip angle
- TR iterating time
- TE echo time
- control of a timing or intensity of the high-frequency magnetic field and the gradient magnetic field which is set in the pulse sequence enables any imaging cross section of an imaging target subject to be imaged.
- the pixel value of the MR image obtained in such a manner is determined depending on T1, T2, the diffusion coefficient, or the like of proton density or tissue of the subject 103 ; however, various changes in the values of the imaging parameters enable contrast of the tissue in an obtained image, a degree of T1 contrast, T2 contrast, or diffusion contrast to be changed.
- the calculator 110 of the embodiment controls the imaging unit such that the imaging unit performs the imaging a plurality of times in which the values of the imaging parameters are changes and the same pulse sequence is performed, and obtains a value that does not depend on the imaging parameters, that is, a quantitative value, through calculation from a plurality of images having different pixel values obtained in the imaging performed the plurality of times. In this manner, an image having the quantitative value as the pixel value, that is the quantitative value image, is obtained. Specifically, the pixel value of the image data obtained in the imaging is fitted to the signal function of the pulse sequence used in the imaging, and the quantitative value is obtained by estimating the quantitative value as a variable of the signal function.
- the quantitative-value candidate group selected through a predetermined technique is used as the initial value in the computation for estimating the quantitative value.
- FIG. 5 illustrates an example of a configuration of the calculator 110 for realizing the computation.
- the calculator 110 includes: the measurement control unit 210 that performs the imaging a plurality of times with different imaging parameters in the same pulse sequence and acquires a plurality of images with different pixel values; the predicted-pixel-value computing unit 220 that predicts and computes a pixel value that is acquired from the imaging in the quantitative-value candidate group; and the quantitative-value computing unit 230 that computes the quantitative value from the plurality of acquired images.
- the quantitative-value computing unit 230 of the embodiment includes: the initial-value selecting unit 231 that selects a predetermined quantitative-value candidate group from the plurality of quantitative-value candidate groups with reference to the pixel value computed in the predicted-pixel-value computing unit 220 ; and the optimal-value computing unit 232 that computes a quantitative value through a localized optimization technique using, as the initial value, the values of the quantitative-value candidate group selected in the initial-value selecting unit 231 .
- FIG. 5A illustrates only a portion related to the quantitative value computation; however, as illustrated in FIG.
- the calculator 110 includes: in addition to the units, an image reconstructing unit 240 that performs computation of image reconstruction, correction, or the like on measurement data from echo; and a display control unit 250 that causes the display device 111 to display the calculation result or the image from the calculator 110 ; or the like.
- the calculator 110 of the embodiment is configured to mainly include a CPU and a memory, and the functions of the calculator 110 are realized when the CPU loads and executes software (programs) stored in the storage device 112 in advance in the memory. Note that there is no need to realize all of the functions with software, and a part or the entirety thereof may be realized by hardware such as an application specific integrated circuit (ASIC). In addition, information that is required for executing processes realized by the functions, or information obtained during the processes or at the end of the processes is stored in the storage device 112 .
- ASIC application specific integrated circuit
- the measurement control unit 210 performs a predetermined pulse sequence while the imaging parameters change in response to an instruction of imaging start by the operator via the input device 116 , and acquires a plurality of images (Step S 601 ).
- the predicted-pixel-value computing unit 220 computes a predicted pixel value in the quantitative-value candidate group set in advance by using the used imaging parameters (Step S 602 ).
- the initial-value selecting unit 231 selects an initial quantitative value from the quantitative-value candidate group with reference to the predicted pixel value computed in the predicted-pixel-value computing unit 220 (Step S 603 ).
- the optimal-value computing unit 232 fits the pixel value of the image acquired in the measurement control unit 210 in the signal function through the localized optimization technique by using the initial value selected in the initial-value selecting unit 231 and computes the quantitative value (Step S 604 ).
- the measurement control unit 210 performs the predetermined pulse sequence by changing the values of the imaging parameters and performs the imaging the plurality of times.
- a case of using an RF-spoiled GRASS sequence as the predetermined pulse sequence is described as an example.
- the RF-spoiled GRASS sequence is a type of GrE system pulse sequence. As illustrated in FIG. 7 , first, a slice gradient magnetic pulse 701 is applied and a high-frequency magnetic (RF) pulse 702 is applied, thereby exciting predetermined slice magnetization in the subject 103 . Subsequently, a slice-encoding gradient magnetic pulse 703 and a phase-encoding gradient magnetic pulse 704 for adding position information of a slice direction and a phase-encoding direction to a phase of magnetization are applied. After a dephasing lead-out gradient magnetic field 705 is applied, while a lead-out gradient magnetic pulse 706 for adding position information of a lead-out direction is applied, a nuclear magnetic resonance signal (echo) 707 is measured in predetermined echo time (TE).
- TE predetermined echo time
- a rephasing slice-encoding gradient magnetic pulse 710 and a phase-encoding gradient magnetic pulse 709 are applied.
- the sequence is iterated for a predetermined iterating time (TR), and the number of echoes required for one image is measured.
- TR iterating time
- a phase ( ⁇ ) of the RF pulse is increased by a predetermined amount of increase during the iteration.
- the imaging parameters that are changeable in the RF-spoiled GRASS sequence are the flip angle (FA), the iterating time (TR), the echo time (TE), and an RF-phase image component value ( ⁇ ).
- FA flip angle
- TR iterating time
- TE echo time
- ⁇ RF-phase image component value
- the imaging is performed the plurality of times by changing the values of the imaging parameters, respectively.
- the combinations of imaging parameters of different values may be set by the operator for each imaging parameter; however, it is possible for the operator to select and set the imaging parameters from the combinations of the imaging parameters stored in the storage device 112 in advance.
- the imaging parameters used in the measurement control unit 210 are set by performing an operation of calling the parameters stored in the storage device 112 in advance, by the operator via the input device 116 . Note that it is possible to change the called imaging parameters via the input device 116 by the user before the start of the imaging or to store the changed parameters in the storage device 112 .
- the imaging parameter values used in the imaging is saved in the memory for being used in the predicted pixel value computation in the following Step S 602 or stored in the storage device 112 .
- the calculator 110 (image reconstructing unit) performs image reconstruction from measurement data obtained through the imaging the plurality of times with the different imaging parameter values described above, and obtains the plurality of image data.
- the plurality of image data is saved in the memory for the following computation or stored in the storage device 112 .
- the predicted-pixel-value computing unit 220 calculates a pixel value (predicted pixel value) obtained in a case where imaging is performed with the imaging parameter values used in the imaging, unlike the imaging described above.
- the calculation is performed to obtain an index used when an initial value (initial value of the quantitative value) used in computation in Step (Step S 604 ) of calculating the quantitative value is selected from the plurality of the quantitative-value candidate groups, and is performed using the signal function of the pulse sequence used in the imaging in Step S 601 , the imaging parameters used in the imaging, and the values of the plurality of predetermined quantitative-value candidate groups.
- the signal function is a function generated for each pulse sequence, and a function that restores the signal intensity of each pixel having, as variables, at least one of a parameter (subject parameter) as a physical property value of living tissue, which depends on the subject, and a parameter (apparatus parameter) that depends on the apparatus, and imaging conditions (imaging parameters) set when the user performs the pulse sequence.
- the subject parameters include the longitudinal relaxation time (T1), the transverse relaxation time (T2), and the spin density (PD).
- the apparatus parameters include magnetostatic field intensity (B0), irradiation intensity (B1) of a transmission coil, sensitivity (Sc) of a receiving coil, signal amplification by a pre-amplifier of the receiving coil, or a receiving coefficient (k) that is determined with a transmission loss, a receiving gain, or the like.
- this formula may be used as the signal function.
- a signal function is obtained as an interpolation function by numerical simulation in advance. The signal function obtained for each pulse sequence in advance is stored in the storage device 112 .
- a signal function fs of the RF-spoiled GRASS is expressed as in the following Expression (1) using the parameters described above.
- a represents a product of PD, Sc, and k, and is the proportionality coefficient of the signal function f.
- the signal function f is generated by interpolating a signal generated by exhaustively changing the imaging parameters (FA, TR, TE, and ⁇ ) with respect to any values of T1 and T2 of the subject parameters through the numerical simulation.
- B1 for example, 1
- B0 is the same as the magnetic field intensity (for example, 3T) of the apparatus used in the imaging.
- the numerical simulation is performed by using a subject model obtained by disposing spins on grid points, inputting the pulse sequence, the imaging parameters, and the apparatus parameters, and outputting an NMR signal by solving a Bloch equation which is a basic equation of the magnetic resonance phenomenon.
- the subject model is assigned with spatial distribution ( ⁇ , M0, T1, and T2) of the spins.
- ⁇ represents a gyromagnetic ratio
- M0 represents thermal equilibrium magnetization (spin density).
- the NMR signal is subjected to the image reconstruction, and thereby it is possible to obtain an image under assigned conditions.
- Bloch's equation is a first order linear ordinary differential equation, and is expressed in the following Expression (2).
- (x, y, z) represents a three-dimensional Cartesian coordinate system, and z represents an orientation of the magnetostatic field (intensity is B0).
- (Mx, My, Mz) represents spins
- H represents magnetic field intensity in a z direction
- H1 represents high-frequency magnetic-filed intensity
- Gx, Gy, and Gz represent gradient magnetic field intensity in subscript directions
- f0 represents a frequency in a rotating system of coordinates.
- the quantitative-value candidate groups used in Step S 602 of computing the predicted pixel value is a combination of the values of the quantitative values from the list of the quantitative values obtained in an empirical manner or literature.
- Table 2 shows an example of a list of the quantitative values.
- T1 [ms] 100, 250, 500, 700, 850, 1000, 1200, 1500, 1750, 2000, 2500, 3000, 3500, 4000, 5000 T2 [ms] 40, 50, 60, 70, 80, 100, 125, 150, 200, 400, 600, 1000, 1250, 1500, 2000, 2500 B1 [—] 0.5, 0.7, 0.9, 1.1, 1.3, 1.5 a [—] 0.1, 0.5, 1.0, 10, 100, 1000
- the number of quantitative-value candidate groups may be obtained as the same number as a number obtained by multiplying the numbers of listed values of the quantitative values ([the number of candidates of T1] ⁇ [the number of candidates of T2] ⁇ [the number of candidates of B1] ⁇ [the number of candidates of a]), all of the number thereof values may be the quantitative-value candidate groups, or the quantitative-value candidate groups may be narrowed to an appropriate number thereof.
- the predicted-pixel-value computing unit 220 computes a signal value I as the pixel value by substituting the values above and the values of the imaging parameters (FA, TR, and TE) used in the imaging in the signal function f, for each combination of the quantitative values (T1, T2, B1, and a) of the quantitative-value candidate group. In this manner, as illustrated in FIG. 3 , the predicted pixel value is computed for each quantitative-value candidate group.
- the computed values may be stored in the memory or the storage device 112 as the predicted pixel value data for each quantitative-value candidate group.
- the initial-value selecting unit 231 performs a matching process between the pixel value data obtained in calculation and the predicted pixel value data calculated in the predicted-pixel-value computing unit 220 in each pixel of the image acquired in the measurement control unit 210 .
- each quantitative value of the quantitative-value candidate group from which the predicted pixel value having the highest matching value is computed of the quantitative-value candidate groups is extracted as the initial value.
- P v represents a plurality of pixel value data obtained in the imaging performed the plurality of times
- a square error SSD i is calculated in the following expression.
- the optimal-value computing unit 232 computes the quantitative value through a localized optimization method using the imaging parameter values (TR, TE, FA, and ⁇ ) used in the imaging, and the initial value obtained in Step S 603 of selecting the initial value.
- TR j , TE j , FA j , and ⁇ j represent the imaging parameters in j-th imaging of imaging performed a total of L times
- P v,j represents obtained pixel value data
- T1, T2, B1, and a are computed when g expressed in the following expression has the smallest value.
- a fitting method of T1, T2, B1, and a is not limited thereto.
- any method may be employed.
- Steps S 601 to S 604 above causes the quantitative value to be obtained for each pixel.
- the quantitative value image (map) having the quantitative values as the pixel values is obtained.
- the quantitative value image may be displayed on the display device 111 as it is, or it is possible to estimate an image of imaging parameter values which is not imaged, using the computed quantitative value, and to generate an image having a different degree of contrast or intensity.
- the imaging parameters FA, TR, TE, and ⁇ are changed, the plurality of images are imaged by performing the RF-spoiled GRASS sequence, and a pixel value I for each pixel is fitted to the signal function f of Expression (1), thereby estimating the subject parameters T1 and T2, the apparatus parameter B1, and a as the product of the subject parameter and the apparatus parameter.
- the estimating quantitative values are four parameters, combination patterns of the imaging parameters need to be four or more patterns. The number of combinations of the imaging parameters changes depending on the estimating quantitative values.
- the computation technique of the quantitative values is not limited thereto.
- the spin echo sequence in which the signal intensity is formulated there is a method or the like in which the imaging is repeated while the imaging parameters TE are changed, and the quantitative value T2 and the proportionality coefficient a are computed by using a signal intensity function in which the signal intensity of the obtained image is formulated.
- the imaging parameters need to have two or more patterns.
- FIG. 8 illustrates an image 810 having a, B1, T1, and T2 which are computed in the embodiment, and an image 820 having a, B1, T1, and T2 which are computed by using a simulated annealing method with a global optimization algorithm as a fitting algorithm. Calculation time is as fast as about 15 times in the embodiment, and a difference between the calculation results is within 2%.
- the MRI apparatus of the embodiment includes: an imaging unit that acquires a nuclear magnetic resonance signal from a subject and acquires the image of the subject; a measurement control unit that controls the imaging unit such that the imaging unit performs imaging a plurality of times with different imaging parameter values in the same pulse sequence; a quantitative-value computing unit that computes a quantitative value which does not depend on the imaging parameter values, from a plurality of images having different pixel values obtained by performing the imaging the plurality of times; and further a predicted-pixel-value computing unit that predicts a pixel value that is acquired from the plurality of imaging parameter values, for each of a plurality of predetermined quantitative-value candidate groups.
- the quantitative-value computing unit includes an initial-value selecting unit that selects an initial value of a quantitative value from the plurality of quantitative-value candidate groups with reference to the pixel value computed in the predicted-pixel-value computing unit, and an optimal-value computing unit that computes a quantitative value through a localized optimization technique using the initial value selected in the initial-value selecting unit.
- the MRI apparatus of the embodiment further includes the storage device that stores the data required for the calculation in the quantitative-value computing unit.
- the storage device stores a quantitative-value list in which one or more candidates of each of quantitative values are listed in association with a plurality of types of quantitative values, and the quantitative-value computing unit uses the combinations of the candidates of the quantitative-value lists as the plurality of quantitative-value candidate groups.
- the localized optimization algorithm is advantageous in that calculation cost is lower than that of the global optimization algorithm, and is disadvantageous in that the calculation is likely to result in a localized solution.
- the initial value is set with high accuracy based on the predicted pixel value, and thereby advantages are achieved in that it is possible to achieve the calculation result as that obtained during the use of the global optimization algorithm without resulting in the localized solution, and it is possible to shorten quantitative-value computation time. In this manner, it is possible to shorten the time from the imaging to presenting a quantitative image to the operator, and it is possible for the operator to check the image without stress.
- the MRI apparatus of the embodiment is characterized in that the calculator is provided with separate means for computing a specific quantitative value, of the plurality of the quantitative-value candidates which are to be calculated, and thereby the computation speed is further improved.
- the MRI apparatus of the embodiment basically has the same configuration as the first embodiment. Hereinafter, a configuration different from that of the first embodiment will be described with a focus thereon. In addition, in the embodiment, a case where the separately computed specific quantitative value is the proportionality coefficient is described.
- a calculator 110 C of the embodiment includes the configuration illustrated in FIG. 5A or 5B , in which the quantitative-value computing unit 230 further has a proportionality-coefficient computing unit 233 that computes an initial value of the proportionality coefficient. Except for this, the components assigned with the same reference signs as in FIGS. 5A and 5B are the same as those in the second embodiment.
- a is not limited to PD ⁇ Sc ⁇ k.
- the functions of the components other than the proportionality-coefficient computing unit 233 are substantially the same as those in the first embodiment, and the repeated description thereof is omitted.
- the process flow of the quantitative value computation in the calculator 110 A of the embodiment is described with reference to FIG. 10 .
- the outline of the process is as follows.
- the measurement control unit 210 performs the predetermined pulse sequence while the imaging parameters change in response to the instruction of imaging start by the operator via the input device 116 , and acquires a plurality of images (Step S 611 ).
- the predicted-pixel-value computing unit 220 computes the predicted pixel value in the quantitative-value candidate group set in advance by using the used imaging parameters (Step S 612 ).
- the initial-value selecting unit 231 selects the initial quantitative value from the quantitative-value candidate group with reference to the predicted pixel value computed in the predicted-pixel-value computing unit 220 (Step S 613 ).
- the proportionality-coefficient computing unit 233 computes an initial proportionality coefficient for each pixel, based on the initial quantitative value selected in the initial-value selecting unit 231 and the image data acquired from the measurement control unit 210 . (Step S 614 ).
- the optimal-value computing unit 232 fits the pixel value of the image acquired in the measurement control unit 210 in the signal function through the localized optimization technique by using the initial value selected in the initial-value selecting unit 231 and the proportionality coefficient computed in the proportionality-coefficient computing unit 233 , and the quantitative value is computed (Step S 615 ).
- Step S 611 of imaging for acquiring the plurality of images is the same as Step S 601 in the second embodiment, the process proceeds to Step S 612 which is described.
- the predicted-pixel-value computing unit 220 predicts the pixel value which is obtained in the imaging parameters through the imaging performed every time.
- the prediction of the pixel value is performed for each of the plurality of quantitative-value candidate groups.
- Table 3 shows an example of the list of the quantitative values used in the embodiment.
- T1 [ms] 100, 250, 500, 700, 850, 1000, 1200, 1500, 1750, 2000, 2500, 3000, 3500, 4000, 5000 T2 [ms] 40, 50, 60, 70, 80, 100, 125, 150, 200, 400, 600, 1000, 1250, 1500, 2000, 2500 B1 [—] 0.5, 0.7, 0.9, 1.1, 1.3, 1.5 a [—] 1.0
- the number of the quantitative-value candidate groups is [the number of T1 values] ⁇ [the number of T2 values] ⁇ [the number of B1 values] at most.
- the N quantitative-value candidate groups are presented.
- the predicted pixel value is obtained for each quantitative-value candidate group.
- the initial-value selecting unit 231 performs matching of the N predicted pixel values and the pixel values of the image obtained through the actual imaging. Specifically, when focusing on a certain pixel, P v represents a plurality of pixel value data obtained in the imaging performed the plurality of times, P c,i (1, . . . N) represents predicted pixel value data obtained by the plurality of times of the imaging which is calculated for each of the N quantitative-value candidate groups, and a regularized cross correlation NCC i is expressed in the following expression. The regularized cross correlation is calculated, and P c,i is obtained when the maximum value of the correlation is obtained.
- TR j , TE j , FA j , and ⁇ j represent the imaging parameters in j-th imaging of imaging performed a total of L times
- P v,j represents obtained pixel value data
- P c,j represents predicted pixel value data that is calculated from the initial quantitative values T1 0 , T2 0 , B1 0 extracted by the initial-value selecting unit 231 and the imaging parameters TR j , TE j , FA j , and ⁇ j
- the proportionality coefficient a is calculated in the following Expression.
- P c,j is anyone (having the maximum correlation) of the plurality of predicted pixel values calculated in the predicted-pixel-value computing unit 220 , and thus there is no need to perform computation again in the step.
- the optimal-value computing unit 232 computes the quantitative value by using the localized optimization algorithm with the initial values T1 0 , T2 0 , and B1 0 extracted by the initial-value selecting unit 231 and the initial value a of the proportionality coefficient calculated by the proportionality-coefficient computing unit 233 .
- the proportionality coefficient is fixed in the quantitative-value candidate groups, the number of candidate groups is reduced, and thereby the calculation cost is reduced.
- the initial quantitative value that needs to be extracted and the proportionality coefficient for each pixel from the image data imaged are computed, and thereby it is possible to reset the initial value with high accuracy even with respect to the proportionality coefficient fixed in the candidate groups, and it is possible to compute without resulting in the localized solution even when the quantitative value is computed by using the localized optimization algorithm.
- the prediction of the pixel value and selection of the quantitative-value candidate group as the initial value of the quantitative value computation with reference to the result of the prediction are the same as those in the first to third embodiments; however, the MRI apparatus of the embodiment is characterized by having a function of presenting an initial value selected by the initial-value selecting unit. In this manner, operability in a quantitative-value diagnostic examination improves.
- the MRI apparatus of the embodiment basically has the same configuration as that of the second embodiment. Hereinafter, a configuration different from that of the second embodiment will be described with a focus thereon.
- a calculator 110 D of the embodiment includes the configuration illustrated in FIG. 5A or 5B , in which the quantitative-value computing unit 230 has an initial-value presenting unit 234 that causes an initial value image having, as the pixel value, the initial value selected by the initial-value selecting unit to be displayed on the display device 111 .
- the display control unit 250 FIG. 5B
- the measurement control unit 210 performs the predetermined pulse sequence while the imaging parameters change in response to the instruction of imaging start by the operator via the input device 116 , and acquires a plurality of images (Step S 621 ).
- the predicted-pixel-value computing unit 220 computes the predicted pixel value in the quantitative-value candidate group set in advance by using the used imaging parameters (Step S 622 ).
- the initial-value selecting unit 231 selects the initial quantitative value from the quantitative-value candidate group with reference to the predicted pixel value computed in the predicted-pixel-value computing unit 220 (Step S 623 ).
- the initial-value presenting unit 234 displays, as an initial-value image, the initial value selected by the initial-value selecting unit 231 on the display device 111 and presents the image to the operator (Step S 624 ).
- FIG. 13 illustrates an example of the initial-value image that is displayed on the display device 111 .
- the quantitative value is T1 and T2
- a T1 image 901 and a T2 image 902 are displayed.
- the operator checks the presented images, thereby it is possible to substantially check the structure, and it is possible to determine whether the imaging is successful or failed in the examination.
- This image is presented and the UI for inputting result of determination of whether the operator proceeds with the quantitative-value computation or performs imaging again, for example, instruction buttons 903 and 904 as illustrated in FIG. 13 , is displayed on the display device 111 .
- Step S 625 the process proceeds to the next step.
- the optimal-value computing unit 232 fits the pixel value of the image acquired in the measurement control unit 210 in the signal function through the localized optimization technique by using the initial value selected by the initial-value selecting unit 231 , and the quantitative value is computed (Step S 626 ).
- imaging is appropriately again performed by the instruction by the operator without performing S 626 of computation of the quantitative value.
- the MRI apparatus of the embodiment presents the initial value image, and thereby it is possible for the operator to determine whether or not the imaging needs to be performed without waiting for the result of the quantitative value calculation. In this manner, in the diagnosis by using the quantitative image, there is no need to wait for the determination of whether or not the imaging needs to be again performed until the end of the quantitative value calculation, and thus throughput improves.
- the calculator may be configured to separately compute a predetermined quantitative value, for example, the proportionality coefficient a, from the other quantitative values. In this case, it is possible to shorten time of correlation calculation for selecting the quantitative-value candidate groups by the initial-value selecting unit 231 .
- the initial-value image is presented and the operator is to determine whether or not the quantitative value calculation continues; however, for example, a value is presented when the value of the correlation (Expression (4) and Expression (7)) calculated in Step S 623 is the maximum, and thus whether or not the quantitative value calculation by the value continues is determined.
- a predetermined threshold value is set, and thereby whether or not the quantitative value calculation automatically continues on the apparatus side may be determined.
- the MRI apparatus of the embodiment further includes a region setting unit that receives designation of a range for which the quantitative value is calculated, and the optimal-value computing unit computes the quantitative value for the region that is received by the region setting unit.
- the optimal-value computing unit further includes a result presenting unit that presents a result obtained by calculating the quantitative value preferentially in association with the region which is received by the region setting unit.
- the MRI apparatus of the embodiment basically has the same configuration as that of the second embodiment. Hereinafter, a configuration different from that of the second embodiment will be described with a focus thereon.
- a calculator 110 E of the embodiment includes the configuration in FIG. 5A or 5B , and further includes a preferential-calculation-region designating unit 235 that designates a region for which the calculation is preferentially performed as illustrated in FIG. 14 , and a preferential-calculation-result presenting unit 236 that presents results of the preferential calculation.
- the preferential-calculation-region designating unit 235 generates the user interface (UI) such that the operator designates the preferential calculation region and displays the UI on the display device 111 .
- the preferential-calculation-result presenting unit 236 displays, on the display device 111 as the quantitative value image of the region, the result obtained after the quantitative computation is performed on the preferential calculation region.
- the functioning unit belongs to the quantitative-value computing unit 230 ; however, the functioning unit may belong to the function of the display control unit 250 .
- the measurement control unit 210 performs the predetermined pulse sequence while the imaging parameters change in response to the instruction of imaging start by the operator via the input device 116 , and acquires a plurality of images (Step S 631 ).
- the predicted-pixel-value computing unit 220 computes the predicted pixel value in the quantitative-value candidate group set in advance by using the used imaging parameters (Step S 632 ).
- the initial-value selecting unit 231 selects the initial quantitative value from the quantitative-value candidate group with reference to the predicted pixel value computed in the predicted-pixel-value computing unit 220 (Step S 633 ).
- the process flow described to here is the same as that of the second embodiment.
- the preferential-calculation-region designating unit 235 receives the instruction by the operator via the input device 116 , and designates the region on which the calculation needs to be preferentially performed (Step S 634 ).
- FIG. 16 illustrates an example of a User Interface 1200 of the preferential-calculation-region designating unit 235 .
- the preferential-calculation-region designating unit 235 displays a transverse section image 1201 , a coronal section image 1202 , a sagittal section image 1203 of the image which is first imaged by the measurement control unit 210 , and displays, on the display device 111 (display screen), a UI 1204 (for example, an enclosed shape in a dotted line) for designating the preferential calculation region.
- the operator operates a mouse or the like, thereby making it possible to perform operation of horizontal movement, rotation, contraction and expansion, or the like of the UI 1204 .
- types of UI 1204 illustrated in FIG. 16 are displayed on three sections, respectively and are associated with the coordinates. A result of one operation of the UI 1204 which overlaps the three sections is reflected to the other UI.
- the preferential-calculation-region designating unit 235 receives information of a UI operation and designates, as the preferential calculation region, a region surrounded by the UI 1204 .
- an “input completion button” 1205 for confirming the designation is displayed, the operator performs operation on the button 1205 , and thereby the designation of the preferential calculation region is completed.
- the optimal-value computing unit 232 fits the pixel value of the image acquired in the measurement control unit 210 in the signal function through the localized optimization technique by using the initial value selected by the initial-value selecting unit 231 and the quantitative value of the pixel is computed in the region designated by the preferential-calculation-region designating unit 235 (Step S 635 ).
- the preferential-calculation-result presenting unit 236 displays, on the display device 111 , the quantitative value in the region designated by the preferential-calculation-region designating unit 235 , and presents the value to the operator (Step S 636 ).
- a display mode of the quantitative value image there is not particular limitation on a display mode of the quantitative value image; however, for example, similar to an example ( FIG. 13 ) of displaying the initial value image in the third embodiment, one cross section of the preferential calculation region may be displayed, or, only the region designated as the preferential calculation region may be replaced with the quantitative value region on a three-section image illustrated in FIG. 16 .
- the UI on which the next process is selected, may be displayed on this screen or another screen.
- the processes after the process described above include a process of continuing the quantitative value computation for the region other than the preferential calculation region, image reconstruction, imaging end, or the like.
- the operator checks the calculation result of the preferential calculation region, and it is possible to perform the remaining quantitative value computation, or to perform the image reconstruction. In a case where the information that is wanted to be obtained is only the preferential calculation region, it is possible to end the imaging at that time.
- the optimal-value computing unit 232 computes the quantitative values for the remaining calculation region (Step S 637 ).
- the MRI apparatus of the embodiment it is possible to check image quality in the examination by the quantitative image diagnosis without waiting for the calculation of all of the regions, it is possible to determine the success or failure of the image, and the operability improves.
- the second embodiment to the fifth embodiment, in which the invention is applied to the MRI apparatus are described; however, the function of the calculator described in the embodiments can be performed not only the calculator incorporated in the MRI apparatus, but also by the calculator separately provided from the MRI apparatus.
- the storage unit of the data required for a part of the calculation or the calculation it is possible to use software or a recording medium that is built on the cloud or on the internet.
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- General Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- Signal Processing (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Animal Behavior & Ethology (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Multimedia (AREA)
- Gynecology & Obstetrics (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Description
- The present invention relates to a technology that generates a quantitative-value image in response to a signal acquired from a magnetic resonance imaging apparatus.
- A magnetic resonance imaging (hereinafter, MRI) apparatus is a medical diagnostic imaging apparatus that mainly uses a nuclear magnetic resonance phenomenon of protons. The MRI apparatus is capable of non-invasively imaging any cross section and is capable of acquiring information related to vital functions such as bloodstream or metabolic functions, in addition to morphological information. In general, a slice gradient magnetic field is applied to a subject positioned in a magnetostatic field and, simultaneously, a high-frequency magnetic field having a specific frequency is applied thereto, and thus nuclear magnetization is excited in a cross section that needs to be imaged. A phase-encoding gradient magnetic field and a lead-out gradient magnetic field are applied to the excited nuclear magnetization, thereby plane position information is assigned, and a nuclear magnetic resonance signal (echo) generated from the nuclear magnetization is measured. The measurement of the nuclear magnetic resonance signal is iterated until a measurement space referred to as a k space is filled with the signals. The signals, with which the k space is filled, are converted into an image through inverse Fourier transform.
- Pulses and gradient magnetic fields for generating an echo are applied, based on a preset pulse sequence. Various types of pulse sequences are known according to purposes. For example, in a gradient echo (GrE) type of fast imaging technique, the phase-encoding gradient magnetic field is sequentially changed for every iteration time (hereinafter, TR) of the pulse sequence, and the number of nuclear magnetic resonance signals that is required for obtaining one tomogram is measured.
- In general, in MR examination, a user selects and executes a pulse sequence, thereby acquiring a weighted image in which a relative difference in quantitative values of living tissue (referred to as subject parameters, for example, T1: longitudinal relaxation time, T2: transverse relaxation time, PD: proton density, or D: diffusion coefficient) is weighted. When a degree of weighting or a target quantitative value is changed, it is necessary to select another pulse sequence or to change an imaging parameter. In general MR examination, an imaging site that needs to be subjected to diagnosis using an image acquired by the imaging for positioning is set, the pulse sequence or the imaging parameter is changed, and a plurality of image types (for example, a T1 weighted image, a T2 weighted image, FLAIR: fluid attenuated inversion recovery, a diffusion weighted image, and MRA: magnetic resonance angiography) are acquired. The user adjusts a window level (WL) or a window width (WW) of an image acquired through a manual operation and performs clipping of a signal that interferes with a diagnosis as necessary, thereby generating an image for diagnosis.
- On the other hand, in recent years, with an early diagnosis of knee osteoarthritis, a diagnosis of a necrotic site of myocardium, or the like, as a target, clinical usefulness of using, as a diagnostic image, a quantitative value image having the quantitative values described above as pixel values is reported. In addition, when multicenter clinical study is performed by using the quantitative value image, there is no need to consider apparatus parameters depending on hardware, and thus it is easy to standardize an MR image, whereas it was difficult to standardize the image in the related art. For example, as a method of computing the quantitative value, a global optimization algorithm is used and a method of computing proton density is proposed (US-A-2009/0157350).
- In a method disclosed in US-A-2009/0157350, in order to avoid a case where an estimated quantitative value results in a localized solution, a quantitative value is computed by using the global optimization algorithm which is high in calculating costs. Compared to a case of the weighted image in which it is possible to check an image immediately after completion of the imaging, since the quantitative value starts to be calculated after the completion of the imaging, calculating time is standby time for determination of whether the imaging is successful or failed, and thus a heavy burden is imposed on an operator.
- In consideration of such problems described above, an object of the invention is to provide an MRI apparatus that is capable of reducing calculating costs in quantitative value estimating calculation, and to rapidly determine whether the imaging is successful or failed.
- According to the invention, in quantitative value computation, an initial value of a quantitative value is set with reference to a pixel value obtained by performing calculation from a quantitative-value candidate group, and the quantitative value is computed through a localized optimization technique.
- Specifically, the MRI apparatus of the invention includes: an imaging unit that acquires a nuclear magnetic resonance signal from a subject and performs imaging to acquire image data of the subject; a measurement control unit that controls the imaging unit such that the imaging unit performs imaging a plurality of times with different imaging parameter values in the same pulse sequence; a quantitative-value computing unit that computes a quantitative value which does not depend on the imaging parameter values, from a plurality of images having different pixel values obtained by performing the imaging the plurality of times; and a predicted-pixel-value computing unit that predicts a pixel value that is acquired from the plurality of imaging parameter values, for each of a plurality of predetermined quantitative-value candidate groups. The quantitative-value computing unit includes an initial-value selecting unit that selects an initial value of a quantitative value from the plurality of quantitative-value candidate groups with reference to the pixel value computed in the predicted-pixel-value computing unit, and an optimal-value computing unit that computes a quantitative value through a localized optimization technique using the initial value selected in the initial-value selecting unit.
- According to the invention, it is possible to shorten time taken to the quantitative value computation, and to expect improvement of operability by shortening the time taken from imaging to image checking.
-
FIG. 1 is a functional block diagram of a calculator of a first embodiment. -
FIG. 2 is a diagram illustrating an example of flow of the calculator of the first embodiment. -
FIG. 3 is a diagram schematically illustrating a pixel-value predicting process of the first embodiment. -
FIG. 4 is a block diagram illustrating a common configuration of an MRI apparatus of a second embodiment. -
FIG. 5A illustrates an example of a functional block diagram of a calculator of the second embodiment. -
FIG. 5B illustrates another example of the functional block diagram of the calculator of the second embodiment. -
FIG. 6 illustrates an example of flow of the calculator of the second embodiment. -
FIG. 7 illustrates an example of a pulse sequence that is employed in the MRI apparatus of the embodiment. -
FIG. 8 illustrates effects of the second embodiment. -
FIG. 9 is a functional block diagram of a calculator of a third embodiment. -
FIG. 10 illustrates an example of flow of the calculator of the third embodiment. -
FIG. 11 is a functional block diagram of a calculator of a fourth embodiment. -
FIG. 12 illustrates an example of flow of the calculator of the fourth embodiment. -
FIG. 13 illustrates an example of presentation of an initial value in the fourth embodiment. -
FIG. 14 is a functional block diagram of a calculator of a fifth embodiment. -
FIG. 15 illustrates an example of flow of the calculator of the fifth embodiment. -
FIG. 16 illustrates an example of a user interface of the fifth embodiment. - Next, embodiments of the invention will be described. In all of the figures with which the embodiments of the invention are described in reference, unless otherwise noted, the same reference signs are assigned to components having the same functions, and thus the repeated description thereof is omitted.
- The embodiment is characterized by a calculator that performs quantitative value calculation using a plurality of images having different pixel values, which are obtained by performing imaging a plurality of times with different imaging parameter values in the same pulse sequence, the calculator including: a predicted-pixel-value computing unit that predicts a predicted value of a pixel value that is acquired from the imaging parameters in the imaging performed the plurality of times; an initial-value selecting unit that selects a quantitative-value candidate group which is used as initial values in quantitative value calculation from a plurality of predetermined quantitative-value candidate groups with reference to the pixel value predicted by the predicted-pixel-value computing unit; and an optimal-value computing unit that computes a quantitative value through a localized optimization technique using the initial value selected by the initial-value selecting unit. In addition, the embodiment is characterized by a program that is executed for performing the quantitative-value calculation and is installed in the calculator. In the embodiment, the calculator may be incorporated in an MRI apparatus, or the calculator may be independently provided from the MRI apparatus.
- In the embodiment, the “quantitative value” means a quantitative value that does not depend on a value of the imaging parameter, that is, a parameter depending on a subject which determines a pixel value of an MR image. Specifically, the quantitative value includes a subject parameter such as a transverse relaxation time, a longitudinal relaxation time, spin density, a diffusion coefficient, or RF irradiation sensitivity, a proportionality coefficient of a signal function of a pulse sequence, or the like. A quantitative-value image means an image having, as a pixel value, any quantitative value of the quantitative values.
- Hereinafter, a configuration of the calculator of the embodiment and flow of the quantitative value computation will be described with reference to a functional block diagram in
FIG. 1 and flow inFIG. 2 . - As illustrated in
FIG. 1 , a calculator 100 has a function of computing the quantitative value and includes: a predicted-pixel-value computing unit 220 that uses imaging parameters used in practical imaging and predicts a pixel value that is acquired from imaging parameters; an initial-value selecting unit 231 that uses a plurality ofimage data 310 acquired from the MRI apparatus and selects, as an initial value of the quantitative value computation, a predetermined candidate group from a plurality of predetermined quantitative-value candidate groups 320; and an optimal-value computing unit 232 that computes the quantitative value by using the initial value selected by the initial-value selecting unit 211. Hereinafter, functions of the initial-value selecting unit 231 and the optimal-value computing unit 232 are collectively referred to as a quantitative-value computing unit 230. A function of the calculator 100 is mainly achieved by executing the program incorporated in a CPU. In addition, a part or the entirety of the quantitative-value computation may be performed by hardware such as an ASIC or an FPGA. - The plurality of
image data 310 and quantitative-value candidate groups 320 used by predicted-pixel-value computing unit 220 and the quantitative-value computing unit 230 are stored in astorage device 300 in advance. Thestorage device 300 may be an internal storage device in the MRI apparatus or an external storage device assembled to the MRI apparatus, or a storage device independently provided from the MRI apparatus. Theimage data 310 is data of an image that is imaged by the MRI apparatus in a pulse sequence and an imaging condition (imaging parameter) depending on a type of quantitative value as a purpose of the computation and is formed in response to an acquired signal, and contains a plurality of image data obtained by changing combinations of the imaging parameter values. Note that the image data may be data obtained after image reconstruction is performed or data obtained before the image reconstruction is performed. The quantitative-value candidate groups 320 is a combination of predetermined quantitative values from a table in which a plurality of values are listed in associated with a plurality of types of quantitative values, respectively, and thus the plurality of values can be employed from values obtained in an empirical manner or literature. In addition, the table can have values different for each imaging site. - Flow of a process of the calculator (quantitative-value computing unit 230) 100 of the embodiment is described with reference to
FIG. 2 , based on such a configuration described above. First, information of the plurality ofimage data 310 and imaging parameter values, which is acquired in combinations of different imaging parameter values, is acquired from thestorage device 300. In addition, the plurality of quantitative-value candidate groups 320 are acquired from the list of the quantitative values (Step S101). - The predicted-pixel-
value computing unit 220 uses the imaging parameter values and the plurality of quantitative-value candidate groups 320 acquired in Step S101, and predicts a pixel value obtained in a case where imaging is performed with the imaging parameters (Step S102).FIG. 3 schematically illustrates the prediction of the pixel value.FIG. 3 illustrates a case of computing a certain pixel value. As illustrated in the figure, when the image data having different L values with L imaging parameters is acquired and there are N quantitative-value candidate groups, the predicted-pixel-value computing unit 220 substitutes each of the L imaging parameters and the N quantitative-value candidate groups into a signal function, computes pixel values with respect to the L imaging parameters for each quantitative-value candidate group, and computes L pixel values. Note that the signal function is a function generated for each pulse sequence, and a function that restores signal intensity of each pixel with the quantitative values and the imaging parameters as variables. - Next, the initial-value selecting unit 211 selects one quantitative-value candidate group from the plurality of quantitative-value candidate groups with reference to the predicted pixel value (Step S103). Specifically, the predicted pixel value for each quantitative-value candidate group and a pixel value of the
image data 310 that is practically imaged are matched, and the most matched quantitative-value candidate groups are selected. - When the quantitative-value candidate groups are determined in Step S103, the optimal-
value computing unit 232 starts the quantitative value computation using a localized optimization algorithm with the candidate group as initial values (Step S104). Specifically, the pixel value of theimage data 310 is fitted and the quantitative value is estimated as a variable of the signal function. The signal function is the same as the signal function used in Step S102 of predicting the pixel value. The number of calculable quantitative values is smaller than or equal to the number of combinations (L sets inFIG. 3 ) having different imaging parameter values. - Finally, a quantitative value as a value that does not depend on the imaging parameters is computed for each pixel. In other words, a quantitative value image (or, referred to as a subject parameter map) is obtained (Step S105).
- According to the embodiment, the predicted pixel value is calculated for each imaging parameter, the initial value of the quantitative value is set with reference to the predicted pixel value, and thereby it is possible to compute the quantitative value with high accuracy by applying the localized optimization technique without resulting in a localized solution. In this manner, it is possible to significantly reduce calculation costs according to the quantitative-value computation.
- Based on the first embodiment described above, a more specific embodiment is described.
- The MRI apparatus of the embodiment includes the imaging unit, a measurement control unit that controls the imaging unit and measures data that is required for the quantitative value computation, and the calculator that performs the quantitative value computation. The measurement control unit performs imaging a plurality of times with different imaging parameters in the same pulse sequence, and acquires a plurality of images (image data) having different pixel values. The calculator includes the quantitative-value computing unit that performs the quantitative value computation by using information of the plurality of image data acquired through imaging and the predetermined quantitative-value candidate groups.
- First, an entire configuration of the MRI apparatus that is common in the embodiments of the invention including the embodiment is described with reference to
FIG. 4 .FIG. 4 is a block diagram illustrating a common configuration of anMRI apparatus 10. TheMRI apparatus 10, as the imaging unit, includes amagnet 101 that generates a magnetostatic field, agradient coil 102 that generates a gradient magnetic field, anRF coil 107 that irradiates a subject (living body) 103 with a high-frequency magnetic pulse (hereinafter, RF pulse), anRF probe 108 that detects an echo signal that is generated from the subject 103, and a bed (table) 115 on which the subject (for example, living body) 103 is mounted in a space of the magnetostatic field generated by themagnet 101. - Further, the MRI apparatus 10 (imaging unit) includes a gradient-magnetic-
field power supply 105 that drives thegradient coil 102, a high-frequencymagnetic field generator 106 that drives theRF coil 107, and areceiver 109 that receives an echo signal detected by theRF probe 108. When it is necessary to adjust magnetostatic homogeneity, theMRI apparatus 10 may further include ashim coil 113 and ashim power supply 114 that drives theshim coil 113. Theshim coil 113 is formed to have a plurality of channels, and an additional magnetic field that corrects magnetostatic inhomogeneity with a current supplied from theshim power supply 114. The currents flowing through the channels that configure theshim coil 113 during the adjusting of the magnetostatic homogeneity is controlled by asequencer 104 which will be described below. - In addition, the
MRI apparatus 10, as a computation control system, includes thesequencer 104 that transmits commands to the gradient-magnetic-field power supply 105 and the high-frequencymagnetic field generator 106 such that the gradient-magnetic-field power supply and the high-frequency magnetic field generator generate a gradient magnetic field and a high-frequency magnetic field, respectively, and sets a nuclear magnetic resonance frequency as a reference of detection to thereceiver 109, acalculator 110 that performs signal processing with respect to a detected signal, and performs control or the like of an entire operation of theMRI apparatus 10, adisplay device 111 that displays a processing result from thecalculator 110, astorage device 112 that stores the processing result, and aninput device 116 that receives instructions from a user. Thedisplay device 111 and theinput device 116 are disposed to be close to each other or as an integrated console, and may function as a user interface. Various types of data that are required in thecalculator 110 are stored in thestorage device 112. In addition, not only the result of the signal processing, but also the detected signal, imaging conditions, or the like may be stored in thestorage device 112 as necessary. - In the
MRI apparatus 10 having such a configuration, control of thesequencer 104 causes the RF pulse to be applied to the subject 103 through theRF coil 107, and the gradient magnetic pulse for assigning position information such as slice selection or phase encoding to the echo signal to be applied to the subject by thegradient coil 102. In addition, a signal generated from the subject 103 is received by theRF probe 108, the detected signal is transmitted to thecalculator 110, and the signal processing such as image reconstruction is performed on the signal. - The
calculator 110 performs not only the signal processing in which the received signal is processed, but also outputs an instruction to thesequencer 104 such that components operate at a timing and intensity which are programmed in advance, controls the operations of components that configure the imaging unit, and performs the measurement. Of programs described above, particularly, a program, in which the high-frequency magnetic field, the gradient magnetic field, the timing and intensity of signal reception are set, is referred to as the pulse sequence, and various pulse sequences are prepared in accordance with an imaging method or an imaging purpose. Measurement is performed depending on imaging parameters that are required for the pulse sequence and controlling the pulse sequence. The pulse sequence is generated in advance and is stored in thestorage device 112, the imaging parameters mainly include a flip angle (FA) of the RF pulse, iterating time (TR), echo time (TE), or the like, and values of the imaging parameters are input via a user interface. - In the
MRI apparatus 10, control of a timing or intensity of the high-frequency magnetic field and the gradient magnetic field which is set in the pulse sequence enables any imaging cross section of an imaging target subject to be imaged. Basically, the pixel value of the MR image obtained in such a manner is determined depending on T1, T2, the diffusion coefficient, or the like of proton density or tissue of the subject 103; however, various changes in the values of the imaging parameters enable contrast of the tissue in an obtained image, a degree of T1 contrast, T2 contrast, or diffusion contrast to be changed. In addition, it is possible to acquire information such as B1 as high-frequency magnetic distribution depending on the subject. - The
calculator 110 of the embodiment controls the imaging unit such that the imaging unit performs the imaging a plurality of times in which the values of the imaging parameters are changes and the same pulse sequence is performed, and obtains a value that does not depend on the imaging parameters, that is, a quantitative value, through calculation from a plurality of images having different pixel values obtained in the imaging performed the plurality of times. In this manner, an image having the quantitative value as the pixel value, that is the quantitative value image, is obtained. Specifically, the pixel value of the image data obtained in the imaging is fitted to the signal function of the pulse sequence used in the imaging, and the quantitative value is obtained by estimating the quantitative value as a variable of the signal function. The quantitative-value candidate group selected through a predetermined technique is used as the initial value in the computation for estimating the quantitative value. -
FIG. 5 illustrates an example of a configuration of thecalculator 110 for realizing the computation. As illustrated inFIG. 5A , thecalculator 110 includes: themeasurement control unit 210 that performs the imaging a plurality of times with different imaging parameters in the same pulse sequence and acquires a plurality of images with different pixel values; the predicted-pixel-value computing unit 220 that predicts and computes a pixel value that is acquired from the imaging in the quantitative-value candidate group; and the quantitative-value computing unit 230 that computes the quantitative value from the plurality of acquired images. - In addition, the quantitative-
value computing unit 230 of the embodiment includes: the initial-value selecting unit 231 that selects a predetermined quantitative-value candidate group from the plurality of quantitative-value candidate groups with reference to the pixel value computed in the predicted-pixel-value computing unit 220; and the optimal-value computing unit 232 that computes a quantitative value through a localized optimization technique using, as the initial value, the values of the quantitative-value candidate group selected in the initial-value selecting unit 231.FIG. 5A illustrates only a portion related to the quantitative value computation; however, as illustrated inFIG. 5B , thecalculator 110 includes: in addition to the units, animage reconstructing unit 240 that performs computation of image reconstruction, correction, or the like on measurement data from echo; and adisplay control unit 250 that causes thedisplay device 111 to display the calculation result or the image from thecalculator 110; or the like. - The
calculator 110 of the embodiment is configured to mainly include a CPU and a memory, and the functions of thecalculator 110 are realized when the CPU loads and executes software (programs) stored in thestorage device 112 in advance in the memory. Note that there is no need to realize all of the functions with software, and a part or the entirety thereof may be realized by hardware such as an application specific integrated circuit (ASIC). In addition, information that is required for executing processes realized by the functions, or information obtained during the processes or at the end of the processes is stored in thestorage device 112. - Hereinafter, flow of an operation in the embodiment will be described, based on the configuration of the
calculator 110 described above. - As illustrated in
FIG. 6 , as the outline of the operation, themeasurement control unit 210 performs a predetermined pulse sequence while the imaging parameters change in response to an instruction of imaging start by the operator via theinput device 116, and acquires a plurality of images (Step S601). Next, the predicted-pixel-value computing unit 220 computes a predicted pixel value in the quantitative-value candidate group set in advance by using the used imaging parameters (Step S602). Next, the initial-value selecting unit 231 selects an initial quantitative value from the quantitative-value candidate group with reference to the predicted pixel value computed in the predicted-pixel-value computing unit 220 (Step S603). Next, the optimal-value computing unit 232 fits the pixel value of the image acquired in themeasurement control unit 210 in the signal function through the localized optimization technique by using the initial value selected in the initial-value selecting unit 231 and computes the quantitative value (Step S604). - Hereinafter, processes in steps will be described in detail.
- Step S601
- The
measurement control unit 210 performs the predetermined pulse sequence by changing the values of the imaging parameters and performs the imaging the plurality of times. Here, a case of using an RF-spoiled GRASS sequence as the predetermined pulse sequence is described as an example. - The RF-spoiled GRASS sequence is a type of GrE system pulse sequence. As illustrated in
FIG. 7 , first, a slice gradientmagnetic pulse 701 is applied and a high-frequency magnetic (RF)pulse 702 is applied, thereby exciting predetermined slice magnetization in the subject 103. Subsequently, a slice-encoding gradientmagnetic pulse 703 and a phase-encoding gradientmagnetic pulse 704 for adding position information of a slice direction and a phase-encoding direction to a phase of magnetization are applied. After a dephasing lead-out gradientmagnetic field 705 is applied, while a lead-out gradientmagnetic pulse 706 for adding position information of a lead-out direction is applied, a nuclear magnetic resonance signal (echo) 707 is measured in predetermined echo time (TE). Finally, a rephasing slice-encoding gradientmagnetic pulse 710 and a phase-encoding gradientmagnetic pulse 709 are applied. The sequence is iterated for a predetermined iterating time (TR), and the number of echoes required for one image is measured. In the RF-spoiled GRASS sequence, a phase (θ) of the RF pulse is increased by a predetermined amount of increase during the iteration. - The imaging parameters that are changeable in the RF-spoiled GRASS sequence are the flip angle (FA), the iterating time (TR), the echo time (TE), and an RF-phase image component value (θ). Of the values, in general, θ is fixed to 117 degrees such that image contrast that less depends on T2 is obtained. When θ changes, dependency of the image contrast on T2 significantly changes.
- In the embodiment, the imaging is performed the plurality of times by changing the values of the imaging parameters, respectively. The combinations of imaging parameters of different values may be set by the operator for each imaging parameter; however, it is possible for the operator to select and set the imaging parameters from the combinations of the imaging parameters stored in the
storage device 112 in advance. In this case, the imaging parameters used in themeasurement control unit 210 are set by performing an operation of calling the parameters stored in thestorage device 112 in advance, by the operator via theinput device 116. Note that it is possible to change the called imaging parameters via theinput device 116 by the user before the start of the imaging or to store the changed parameters in thestorage device 112. -
- Specifically, as an example, Table 1 shows imaging parameters in a case where six images are acquired in the RF-spoiled GRASS sequence.
-
TABLE 1 No TR [ms] FA [deg.] θ [deg.] TE [ms] 1 10 30 8 5 2 30 30 22 5 3 40 10 2 5 4 40 30 2 5 5 40 10 8 5 6 40 10 5 5 - The imaging parameter values used in the imaging is saved in the memory for being used in the predicted pixel value computation in the following Step S602 or stored in the
storage device 112. - The calculator 110 (image reconstructing unit) performs image reconstruction from measurement data obtained through the imaging the plurality of times with the different imaging parameter values described above, and obtains the plurality of image data. The plurality of image data is saved in the memory for the following computation or stored in the
storage device 112. - Step S602
- The predicted-pixel-
value computing unit 220 calculates a pixel value (predicted pixel value) obtained in a case where imaging is performed with the imaging parameter values used in the imaging, unlike the imaging described above. The calculation is performed to obtain an index used when an initial value (initial value of the quantitative value) used in computation in Step (Step S604) of calculating the quantitative value is selected from the plurality of the quantitative-value candidate groups, and is performed using the signal function of the pulse sequence used in the imaging in Step S601, the imaging parameters used in the imaging, and the values of the plurality of predetermined quantitative-value candidate groups. - The signal function is a function generated for each pulse sequence, and a function that restores the signal intensity of each pixel having, as variables, at least one of a parameter (subject parameter) as a physical property value of living tissue, which depends on the subject, and a parameter (apparatus parameter) that depends on the apparatus, and imaging conditions (imaging parameters) set when the user performs the pulse sequence. The subject parameters include the longitudinal relaxation time (T1), the transverse relaxation time (T2), and the spin density (PD). The apparatus parameters include magnetostatic field intensity (B0), irradiation intensity (B1) of a transmission coil, sensitivity (Sc) of a receiving coil, signal amplification by a pre-amplifier of the receiving coil, or a receiving coefficient (k) that is determined with a transmission loss, a receiving gain, or the like.
- In a case of a pulse sequence in which the signal intensity is formulated, this formula may be used as the signal function. On the other hand, in a case of a pulse sequence in which the signal intensity is not formulated, a signal function is obtained as an interpolation function by numerical simulation in advance. The signal function obtained for each pulse sequence in advance is stored in the
storage device 112. - A signal function fs of the RF-spoiled GRASS is expressed as in the following Expression (1) using the parameters described above.
-
- Here, a represents a product of PD, Sc, and k, and is the proportionality coefficient of the signal function f. The signal function f is generated by interpolating a signal generated by exhaustively changing the imaging parameters (FA, TR, TE, and θ) with respect to any values of T1 and T2 of the subject parameters through the numerical simulation. At this time, since a change in B1 of the imaging target is the same as a change in FA as the imaging parameter, a constant is used as B1 (for example, 1). In addition, B0 is the same as the magnetic field intensity (for example, 3T) of the apparatus used in the imaging.
- The numerical simulation is performed by using a subject model obtained by disposing spins on grid points, inputting the pulse sequence, the imaging parameters, and the apparatus parameters, and outputting an NMR signal by solving a Bloch equation which is a basic equation of the magnetic resonance phenomenon.
- The subject model is assigned with spatial distribution (γ, M0, T1, and T2) of the spins. Here, γ represents a gyromagnetic ratio, and M0 represents thermal equilibrium magnetization (spin density). The NMR signal is subjected to the image reconstruction, and thereby it is possible to obtain an image under assigned conditions.
- Note that the Bloch's equation is a first order linear ordinary differential equation, and is expressed in the following Expression (2).
-
- In Expression (2), (x, y, z) represents a three-dimensional Cartesian coordinate system, and z represents an orientation of the magnetostatic field (intensity is B0). In addition, (Mx, My, Mz) represents spins, H represents magnetic field intensity in a z direction, H1 represents high-frequency magnetic-filed intensity, Gx, Gy, and Gz represent gradient magnetic field intensity in subscript directions, and f0 represents a frequency in a rotating system of coordinates.
- The quantitative-value candidate groups used in Step S602 of computing the predicted pixel value is a combination of the values of the quantitative values from the list of the quantitative values obtained in an empirical manner or literature. Table 2 shows an example of a list of the quantitative values.
-
TABLE 2 Value list T1 [ms] 100, 250, 500, 700, 850, 1000, 1200, 1500, 1750, 2000, 2500, 3000, 3500, 4000, 5000 T2 [ms] 40, 50, 60, 70, 80, 100, 125, 150, 200, 400, 600, 1000, 1250, 1500, 2000, 2500 B1 [—] 0.5, 0.7, 0.9, 1.1, 1.3, 1.5 a [—] 0.1, 0.5, 1.0, 10, 100, 1000 - In the example showed in Table 2, candidates of each of T1, T2, B1, and a are listed as a plurality of values. a represents the proportionality factor of the signal function f defined in Expression (3). The number of quantitative-value candidate groups may be obtained as the same number as a number obtained by multiplying the numbers of listed values of the quantitative values ([the number of candidates of T1]×[the number of candidates of T2]×[the number of candidates of B1]×[the number of candidates of a]), all of the number thereof values may be the quantitative-value candidate groups, or the quantitative-value candidate groups may be narrowed to an appropriate number thereof.
- The predicted-pixel-
value computing unit 220 computes a signal value I as the pixel value by substituting the values above and the values of the imaging parameters (FA, TR, and TE) used in the imaging in the signal function f, for each combination of the quantitative values (T1, T2, B1, and a) of the quantitative-value candidate group. In this manner, as illustrated inFIG. 3 , the predicted pixel value is computed for each quantitative-value candidate group. The computed values may be stored in the memory or thestorage device 112 as the predicted pixel value data for each quantitative-value candidate group. - Step S603
- In order to determine an initial value of iterative computation in the following Step S604 of calculating the quantitative value, the initial-
value selecting unit 231 performs a matching process between the pixel value data obtained in calculation and the predicted pixel value data calculated in the predicted-pixel-value computing unit 220 in each pixel of the image acquired in themeasurement control unit 210. As a result of the matching process, each quantitative value of the quantitative-value candidate group from which the predicted pixel value having the highest matching value is computed of the quantitative-value candidate groups is extracted as the initial value. There is no particular limitation on the match technique; however, for example, when focusing on a certain pixel, Pv represents a plurality of pixel value data obtained in the imaging performed the plurality of times, Pc,i (i=1, . . . N) represents predicted pixel value data obtained by the plurality of times of the imaging which is calculated for each of the N quantitative-value candidate groups, and a square error SSDi is calculated in the following expression. -
SDD i=Σ(P c,i −P v)2 (i=1, . . . ,N) (3) - When i=M and the minimum square error is obtained, an i=M-th quantitative-value candidate of the quantitative-value candidate groups is set as the initial value of the pixel.
- Step S604
- The optimal-
value computing unit 232 computes the quantitative value through a localized optimization method using the imaging parameter values (TR, TE, FA, and θ) used in the imaging, and the initial value obtained in Step S603 of selecting the initial value. - For example, when focusing on a certain pixel, TRj, TEj, FAj, and θj represent the imaging parameters in j-th imaging of imaging performed a total of L times, and Pv,j represents obtained pixel value data, T1, T2, B1, and a are computed when g expressed in the following expression has the smallest value.
-
- Specifically, with the quantitative values selected in the initial-
value selecting unit 231 as initial values T10, T20, B10, and a0, iterative calculation of the following expression is performed until g has a sufficiently small value. -
- Note that a fitting method of T1, T2, B1, and a is not limited thereto. In the embodiment, as long as the localized optimization algorithm represented by the Levenberg-Marquardt method is used, any method may be employed.
- The performance of Steps S601 to S604 above causes the quantitative value to be obtained for each pixel. In other words, the quantitative value image (map) having the quantitative values as the pixel values is obtained. The quantitative value image may be displayed on the
display device 111 as it is, or it is possible to estimate an image of imaging parameter values which is not imaged, using the computed quantitative value, and to generate an image having a different degree of contrast or intensity. - As described above, in the embodiment, the imaging parameters FA, TR, TE, and θ are changed, the plurality of images are imaged by performing the RF-spoiled GRASS sequence, and a pixel value I for each pixel is fitted to the signal function f of Expression (1), thereby estimating the subject parameters T1 and T2, the apparatus parameter B1, and a as the product of the subject parameter and the apparatus parameter. Here, since the estimating quantitative values are four parameters, combination patterns of the imaging parameters need to be four or more patterns. The number of combinations of the imaging parameters changes depending on the estimating quantitative values.
- In addition, the computation technique of the quantitative values is not limited thereto. For example, in the spin echo sequence in which the signal intensity is formulated, there is a method or the like in which the imaging is repeated while the imaging parameters TE are changed, and the quantitative value T2 and the proportionality coefficient a are computed by using a signal intensity function in which the signal intensity of the obtained image is formulated. In this case, since the estimating quantitative values are two parameters, the imaging parameters need to have two or more patterns.
-
FIG. 8 illustrates animage 810 having a, B1, T1, and T2 which are computed in the embodiment, and animage 820 having a, B1, T1, and T2 which are computed by using a simulated annealing method with a global optimization algorithm as a fitting algorithm. Calculation time is as fast as about 15 times in the embodiment, and a difference between the calculation results is within 2%. - The MRI apparatus of the embodiment includes: an imaging unit that acquires a nuclear magnetic resonance signal from a subject and acquires the image of the subject; a measurement control unit that controls the imaging unit such that the imaging unit performs imaging a plurality of times with different imaging parameter values in the same pulse sequence; a quantitative-value computing unit that computes a quantitative value which does not depend on the imaging parameter values, from a plurality of images having different pixel values obtained by performing the imaging the plurality of times; and further a predicted-pixel-value computing unit that predicts a pixel value that is acquired from the plurality of imaging parameter values, for each of a plurality of predetermined quantitative-value candidate groups. The quantitative-value computing unit includes an initial-value selecting unit that selects an initial value of a quantitative value from the plurality of quantitative-value candidate groups with reference to the pixel value computed in the predicted-pixel-value computing unit, and an optimal-value computing unit that computes a quantitative value through a localized optimization technique using the initial value selected in the initial-value selecting unit.
- In addition, the MRI apparatus of the embodiment further includes the storage device that stores the data required for the calculation in the quantitative-value computing unit. The storage device stores a quantitative-value list in which one or more candidates of each of quantitative values are listed in association with a plurality of types of quantitative values, and the quantitative-value computing unit uses the combinations of the candidates of the quantitative-value lists as the plurality of quantitative-value candidate groups.
- According to the embodiment, the following effects are to be achieved. In general, the localized optimization algorithm is advantageous in that calculation cost is lower than that of the global optimization algorithm, and is disadvantageous in that the calculation is likely to result in a localized solution. On the other hand, in the embodiment, the initial value is set with high accuracy based on the predicted pixel value, and thereby advantages are achieved in that it is possible to achieve the calculation result as that obtained during the use of the global optimization algorithm without resulting in the localized solution, and it is possible to shorten quantitative-value computation time. In this manner, it is possible to shorten the time from the imaging to presenting a quantitative image to the operator, and it is possible for the operator to check the image without stress.
- The MRI apparatus of the embodiment is characterized in that the calculator is provided with separate means for computing a specific quantitative value, of the plurality of the quantitative-value candidates which are to be calculated, and thereby the computation speed is further improved. The MRI apparatus of the embodiment basically has the same configuration as the first embodiment. Hereinafter, a configuration different from that of the first embodiment will be described with a focus thereon. In addition, in the embodiment, a case where the separately computed specific quantitative value is the proportionality coefficient is described.
- As illustrated in
FIG. 9 , a calculator 110C of the embodiment includes the configuration illustrated inFIG. 5A or 5B , in which the quantitative-value computing unit 230 further has a proportionality-coefficient computing unit 233 that computes an initial value of the proportionality coefficient. Except for this, the components assigned with the same reference signs as inFIGS. 5A and 5B are the same as those in the second embodiment. As described in the embodiment described above, the proportionality coefficient means “a” when “I=af” in the description as the signal function is expressed in Expression (1), and is a product of the spin density PD, the receiving coil sensitivity distribution Sc, and the receiving coefficient k that is determined by the receiving gain, or the like. Here, when “I=af” in the description, a is not limited to PD×Sc×k. The functions of the components other than the proportionality-coefficient computing unit 233 are substantially the same as those in the first embodiment, and the repeated description thereof is omitted. - The process flow of the quantitative value computation in the
calculator 110A of the embodiment is described with reference toFIG. 10 . The outline of the process is as follows. First, themeasurement control unit 210 performs the predetermined pulse sequence while the imaging parameters change in response to the instruction of imaging start by the operator via theinput device 116, and acquires a plurality of images (Step S611). Next, the predicted-pixel-value computing unit 220 computes the predicted pixel value in the quantitative-value candidate group set in advance by using the used imaging parameters (Step S612). Next, the initial-value selecting unit 231 selects the initial quantitative value from the quantitative-value candidate group with reference to the predicted pixel value computed in the predicted-pixel-value computing unit 220 (Step S613). Next, the proportionality-coefficient computing unit 233 computes an initial proportionality coefficient for each pixel, based on the initial quantitative value selected in the initial-value selecting unit 231 and the image data acquired from themeasurement control unit 210. (Step S614). Next, the optimal-value computing unit 232 fits the pixel value of the image acquired in themeasurement control unit 210 in the signal function through the localized optimization technique by using the initial value selected in the initial-value selecting unit 231 and the proportionality coefficient computed in the proportionality-coefficient computing unit 233, and the quantitative value is computed (Step S615). - Hereinafter, the processes of the embodiment will be described in detail. Since Step S611 of imaging for acquiring the plurality of images is the same as Step S601 in the second embodiment, the process proceeds to Step S612 which is described.
- Step S612
- The predicted-pixel-
value computing unit 220 predicts the pixel value which is obtained in the imaging parameters through the imaging performed every time. The prediction of the pixel value is performed for each of the plurality of quantitative-value candidate groups. In the quantitative-value candidate groups used in the embodiment, the value of the proportionality coefficient a is fixed, for example, as a=1, combinations of the other quantitative values (T1, T2, and B1) are generated from the list of quantitative values, and the quantitative-value candidate group is formed. Table 3 shows an example of the list of the quantitative values used in the embodiment. -
TABLE 3 Value list T1 [ms] 100, 250, 500, 700, 850, 1000, 1200, 1500, 1750, 2000, 2500, 3000, 3500, 4000, 5000 T2 [ms] 40, 50, 60, 70, 80, 100, 125, 150, 200, 400, 600, 1000, 1250, 1500, 2000, 2500 B1 [—] 0.5, 0.7, 0.9, 1.1, 1.3, 1.5 a [—] 1.0 - The number of the quantitative-value candidate groups is [the number of T1 values]×[the number of T2 values]×[the number of B1 values] at most. Here, the N quantitative-value candidate groups are presented. In this process, the predicted pixel value is obtained for each quantitative-value candidate group.
- Step S613
- The initial-
value selecting unit 231 performs matching of the N predicted pixel values and the pixel values of the image obtained through the actual imaging. Specifically, when focusing on a certain pixel, Pv represents a plurality of pixel value data obtained in the imaging performed the plurality of times, Pc,i (1, . . . N) represents predicted pixel value data obtained by the plurality of times of the imaging which is calculated for each of the N quantitative-value candidate groups, and a regularized cross correlation NCCi is expressed in the following expression. The regularized cross correlation is calculated, and Pc,i is obtained when the maximum value of the correlation is obtained. -
- When i=M and the maximum regularized cross correlation is obtained, an i=M-th quantitative-value candidate of the quantitative-value candidate groups is set as the initial quantitative value i the pixel. Since the regularized cross correlation is the matching process technique in which an influence of the proportionality coefficient is negligible, it is possible to fix the value of the proportionality coefficient in the quantitative value candidates. In this manner, since it is possible to reduce the number of the quantitative-value candidate groups, an advantage is achieved in that it is possible to improve the computation speed.
- Step S614
- In Step S611, the proportionality-
coefficient computing unit 233 computes the proportionality coefficient a fixed as a=1 by using the pixel value data obtained through imaging and the predicted pixel value data that is computed from the quantitative value selected by the initial-value selecting unit 231. For example, when focusing on a certain pixel, TRj, TEj, FAj, and θj represent the imaging parameters in j-th imaging of imaging performed a total of L times, Pv,j represents obtained pixel value data, and Pc,j represents predicted pixel value data that is calculated from the initial quantitative values T10, T20, B10 extracted by the initial-value selecting unit 231 and the imaging parameters TRj, TEj, FAj, and θj, the proportionality coefficient a is calculated in the following Expression. -
- Note that Pc,j is anyone (having the maximum correlation) of the plurality of predicted pixel values calculated in the predicted-pixel-
value computing unit 220, and thus there is no need to perform computation again in the step. - Step S615
- The optimal-
value computing unit 232 computes the quantitative value by using the localized optimization algorithm with the initial values T10, T20, and B10 extracted by the initial-value selecting unit 231 and the initial value a of the proportionality coefficient calculated by the proportionality-coefficient computing unit 233. - As described above, in the MRI apparatus of the embodiment, the proportionality coefficient is fixed in the quantitative-value candidate groups, the number of candidate groups is reduced, and thereby the calculation cost is reduced. In addition, the initial quantitative value that needs to be extracted and the proportionality coefficient for each pixel from the image data imaged are computed, and thereby it is possible to reset the initial value with high accuracy even with respect to the proportionality coefficient fixed in the candidate groups, and it is possible to compute without resulting in the localized solution even when the quantitative value is computed by using the localized optimization algorithm.
- Even in the embodiment, the prediction of the pixel value and selection of the quantitative-value candidate group as the initial value of the quantitative value computation with reference to the result of the prediction are the same as those in the first to third embodiments; however, the MRI apparatus of the embodiment is characterized by having a function of presenting an initial value selected by the initial-value selecting unit. In this manner, operability in a quantitative-value diagnostic examination improves. The MRI apparatus of the embodiment basically has the same configuration as that of the second embodiment. Hereinafter, a configuration different from that of the second embodiment will be described with a focus thereon.
- As illustrated in
FIG. 11 , acalculator 110D of the embodiment includes the configuration illustrated inFIG. 5A or 5B , in which the quantitative-value computing unit 230 has an initial-value presenting unit 234 that causes an initial value image having, as the pixel value, the initial value selected by the initial-value selecting unit to be displayed on thedisplay device 111. Although omitted in the figure, in a case where thecalculator 110D includes the display control unit 250 (FIG. 5B ), it is possible for the display control unit to function as the initial-value presenting unit 234. - The process flow of the quantitative value computation in the
calculator 110D of the embodiment is described with reference toFIG. 12 . - First, the
measurement control unit 210 performs the predetermined pulse sequence while the imaging parameters change in response to the instruction of imaging start by the operator via theinput device 116, and acquires a plurality of images (Step S621). Next, the predicted-pixel-value computing unit 220 computes the predicted pixel value in the quantitative-value candidate group set in advance by using the used imaging parameters (Step S622). Next, the initial-value selecting unit 231 selects the initial quantitative value from the quantitative-value candidate group with reference to the predicted pixel value computed in the predicted-pixel-value computing unit 220 (Step S623). - Next, the initial-
value presenting unit 234 displays, as an initial-value image, the initial value selected by the initial-value selecting unit 231 on thedisplay device 111 and presents the image to the operator (Step S624).FIG. 13 illustrates an example of the initial-value image that is displayed on thedisplay device 111. In this example, the quantitative value is T1 and T2, and aT1 image 901 and aT2 image 902 are displayed. The operator checks the presented images, thereby it is possible to substantially check the structure, and it is possible to determine whether the imaging is successful or failed in the examination. This image is presented and the UI for inputting result of determination of whether the operator proceeds with the quantitative-value computation or performs imaging again, for example,instruction buttons FIG. 13 , is displayed on thedisplay device 111. - In a case where the operator determines that the quantitative-value computation may proceed using the presented initial value of the quantitative value (Step S625), the process proceeds to the next step. Then, the optimal-
value computing unit 232 fits the pixel value of the image acquired in themeasurement control unit 210 in the signal function through the localized optimization technique by using the initial value selected by the initial-value selecting unit 231, and the quantitative value is computed (Step S626). On the other hand, in a case where the operator determines that the imaging needs to be again performed from the presented initial value image, imaging is appropriately again performed by the instruction by the operator without performing S626 of computation of the quantitative value. - The MRI apparatus of the embodiment presents the initial value image, and thereby it is possible for the operator to determine whether or not the imaging needs to be performed without waiting for the result of the quantitative value calculation. In this manner, in the diagnosis by using the quantitative image, there is no need to wait for the determination of whether or not the imaging needs to be again performed until the end of the quantitative value calculation, and thus throughput improves.
- In the description above, differences from the MRI apparatus of the second embodiment as the basis are described; however, similar to the MRI apparatus of the third embodiment, in the MRI apparatus of the embodiment, the calculator may be configured to separately compute a predetermined quantitative value, for example, the proportionality coefficient a, from the other quantitative values. In this case, it is possible to shorten time of correlation calculation for selecting the quantitative-value candidate groups by the initial-
value selecting unit 231. - In addition, in the description above, the initial-value image is presented and the operator is to determine whether or not the quantitative value calculation continues; however, for example, a value is presented when the value of the correlation (Expression (4) and Expression (7)) calculated in Step S623 is the maximum, and thus whether or not the quantitative value calculation by the value continues is determined. In this case, a predetermined threshold value is set, and thereby whether or not the quantitative value calculation automatically continues on the apparatus side may be determined.
- In the embodiment, a function of preferentially performing the quantitative value calculation is added according to a region designated by the operator, and thereby the operability in the quantitative value diagnostic examination is improved. In other words, the MRI apparatus of the embodiment further includes a region setting unit that receives designation of a range for which the quantitative value is calculated, and the optimal-value computing unit computes the quantitative value for the region that is received by the region setting unit. In addition, the optimal-value computing unit further includes a result presenting unit that presents a result obtained by calculating the quantitative value preferentially in association with the region which is received by the region setting unit. The MRI apparatus of the embodiment basically has the same configuration as that of the second embodiment. Hereinafter, a configuration different from that of the second embodiment will be described with a focus thereon.
- A
calculator 110E of the embodiment includes the configuration inFIG. 5A or 5B , and further includes a preferential-calculation-region designating unit 235 that designates a region for which the calculation is preferentially performed as illustrated inFIG. 14 , and a preferential-calculation-result presenting unit 236 that presents results of the preferential calculation. Specifically, the preferential-calculation-region designating unit 235 generates the user interface (UI) such that the operator designates the preferential calculation region and displays the UI on thedisplay device 111. The preferential-calculation-result presenting unit 236 displays, on thedisplay device 111 as the quantitative value image of the region, the result obtained after the quantitative computation is performed on the preferential calculation region. In the example illustrated in the figure, the functioning unit belongs to the quantitative-value computing unit 230; however, the functioning unit may belong to the function of thedisplay control unit 250. - The process flow of the quantitative value computation in the
calculator 110E of the embodiment is described with reference toFIG. 15 . - First, the
measurement control unit 210 performs the predetermined pulse sequence while the imaging parameters change in response to the instruction of imaging start by the operator via theinput device 116, and acquires a plurality of images (Step S631). Next, the predicted-pixel-value computing unit 220 computes the predicted pixel value in the quantitative-value candidate group set in advance by using the used imaging parameters (Step S632). Next, the initial-value selecting unit 231 selects the initial quantitative value from the quantitative-value candidate group with reference to the predicted pixel value computed in the predicted-pixel-value computing unit 220 (Step S633). The process flow described to here is the same as that of the second embodiment. - Next, the preferential-calculation-
region designating unit 235 receives the instruction by the operator via theinput device 116, and designates the region on which the calculation needs to be preferentially performed (Step S634). -
FIG. 16 illustrates an example of aUser Interface 1200 of the preferential-calculation-region designating unit 235. In the embodiment, an example of a head examination is described. The preferential-calculation-region designating unit 235 displays atransverse section image 1201, acoronal section image 1202, asagittal section image 1203 of the image which is first imaged by themeasurement control unit 210, and displays, on the display device 111 (display screen), a UI 1204 (for example, an enclosed shape in a dotted line) for designating the preferential calculation region. The operator operates a mouse or the like, thereby making it possible to perform operation of horizontal movement, rotation, contraction and expansion, or the like of theUI 1204. For example, types ofUI 1204 illustrated inFIG. 16 are displayed on three sections, respectively and are associated with the coordinates. A result of one operation of theUI 1204 which overlaps the three sections is reflected to the other UI. The preferential-calculation-region designating unit 235 receives information of a UI operation and designates, as the preferential calculation region, a region surrounded by theUI 1204. In addition, an “input completion button” 1205 for confirming the designation is displayed, the operator performs operation on thebutton 1205, and thereby the designation of the preferential calculation region is completed. - When the preferential calculation region is designated, the optimal-
value computing unit 232 fits the pixel value of the image acquired in themeasurement control unit 210 in the signal function through the localized optimization technique by using the initial value selected by the initial-value selecting unit 231 and the quantitative value of the pixel is computed in the region designated by the preferential-calculation-region designating unit 235 (Step S635). Next, the preferential-calculation-result presenting unit 236 displays, on thedisplay device 111, the quantitative value in the region designated by the preferential-calculation-region designating unit 235, and presents the value to the operator (Step S636). There is not particular limitation on a display mode of the quantitative value image; however, for example, similar to an example (FIG. 13 ) of displaying the initial value image in the third embodiment, one cross section of the preferential calculation region may be displayed, or, only the region designated as the preferential calculation region may be replaced with the quantitative value region on a three-section image illustrated inFIG. 16 . - Further, the UI, on which the next process is selected, may be displayed on this screen or another screen. For example, the processes after the process described above include a process of continuing the quantitative value computation for the region other than the preferential calculation region, image reconstruction, imaging end, or the like. The operator checks the calculation result of the preferential calculation region, and it is possible to perform the remaining quantitative value computation, or to perform the image reconstruction. In a case where the information that is wanted to be obtained is only the preferential calculation region, it is possible to end the imaging at that time.
- In a case where an instruction of continuation of the calculation is input, the optimal-
value computing unit 232 computes the quantitative values for the remaining calculation region (Step S637). - According to the MRI apparatus of the embodiment, it is possible to check image quality in the examination by the quantitative image diagnosis without waiting for the calculation of all of the regions, it is possible to determine the success or failure of the image, and the operability improves.
- As described above, the second embodiment to the fifth embodiment, in which the invention is applied to the MRI apparatus, are described; however, the function of the calculator described in the embodiments can be performed not only the calculator incorporated in the MRI apparatus, but also by the calculator separately provided from the MRI apparatus. In addition, it is possible to appropriately omit configurations other than the calculator of the components described in the embodiments or to add known components to the configurations, and this is included the invention. In addition, as the storage unit of the data required for a part of the calculation or the calculation, it is possible to use software or a recording medium that is built on the cloud or on the internet.
Claims (10)
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016-148607 | 2016-07-28 | ||
JP2016148607A JP6609226B2 (en) | 2016-07-28 | 2016-07-28 | Magnetic resonance imaging apparatus and quantitative value calculation program |
JPJP2016-148607 | 2016-07-28 |
Publications (2)
Publication Number | Publication Date |
---|---|
US20180028273A1 true US20180028273A1 (en) | 2018-02-01 |
US11187760B2 US11187760B2 (en) | 2021-11-30 |
Family
ID=61011854
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/614,782 Active 2039-01-06 US11187760B2 (en) | 2016-07-28 | 2017-06-06 | Magnetic resonance imaging apparatus and quantitative-value computing program |
Country Status (3)
Country | Link |
---|---|
US (1) | US11187760B2 (en) |
JP (1) | JP6609226B2 (en) |
CN (1) | CN107664754B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111803104A (en) * | 2020-07-20 | 2020-10-23 | 上海市第六人民医院 | Medical image display method, medium and electronic equipment |
US11582478B2 (en) * | 2020-09-08 | 2023-02-14 | Alibaba Group Holding Limited | Video encoding technique utilizing user guided information in cloud environment |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040066978A1 (en) * | 2001-04-19 | 2004-04-08 | Kyojiro Nanbu | Image processing method and image processing apparatus |
US20070280520A1 (en) * | 2005-12-29 | 2007-12-06 | Hiroshi Takai | Magnetic resonance imaging apparatus and image correction estimating method |
US7821266B2 (en) * | 2007-08-03 | 2010-10-26 | Siemens Aktiengesellschaft | Method and device for optimization of imaging parameters |
JP2011024926A (en) * | 2009-07-29 | 2011-02-10 | Hitachi Medical Corp | Magnetic resonance imaging apparatus |
US8355556B2 (en) * | 2009-11-16 | 2013-01-15 | University Health Network | Method of synthetic image generation for magnetic resonance imaging |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1004034A1 (en) * | 1997-08-13 | 2000-05-31 | Btg International Limited | Apparatus for and method of determining values of relaxation parameters |
CN101281241B (en) * | 2007-04-06 | 2012-06-20 | Ge医疗系统环球技术有限公司 | MRI equipment and control method thereof |
US20090157350A1 (en) | 2007-12-12 | 2009-06-18 | Chevron U.S.A. Inc. | Obtaining a proton density distribution from nuclear magnetic resonance data |
US9198599B2 (en) * | 2010-04-22 | 2015-12-01 | Hitachi Medical Corporation | Magnetic resonance imaging apparatus |
CN103064046B (en) * | 2012-12-25 | 2015-04-15 | 深圳先进技术研究院 | Image processing method based on sparse sampling magnetic resonance imaging |
CN104166109B (en) * | 2013-05-17 | 2018-08-03 | 上海联影医疗科技有限公司 | MR imaging method and device |
JP5946800B2 (en) * | 2013-07-22 | 2016-07-06 | 株式会社日立製作所 | Magnetic resonance imaging apparatus, image processing apparatus, image processing method, and image processing program |
JP6371554B2 (en) * | 2014-03-28 | 2018-08-08 | 株式会社日立製作所 | Magnetic resonance imaging system |
JP6250795B2 (en) * | 2014-04-22 | 2017-12-20 | 株式会社日立製作所 | Medical image diagnosis support apparatus, magnetic resonance imaging apparatus, and medical image diagnosis support method |
-
2016
- 2016-07-28 JP JP2016148607A patent/JP6609226B2/en active Active
-
2017
- 2017-06-06 US US15/614,782 patent/US11187760B2/en active Active
- 2017-06-07 CN CN201710423785.7A patent/CN107664754B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040066978A1 (en) * | 2001-04-19 | 2004-04-08 | Kyojiro Nanbu | Image processing method and image processing apparatus |
US20070280520A1 (en) * | 2005-12-29 | 2007-12-06 | Hiroshi Takai | Magnetic resonance imaging apparatus and image correction estimating method |
US7821266B2 (en) * | 2007-08-03 | 2010-10-26 | Siemens Aktiengesellschaft | Method and device for optimization of imaging parameters |
JP2011024926A (en) * | 2009-07-29 | 2011-02-10 | Hitachi Medical Corp | Magnetic resonance imaging apparatus |
US8355556B2 (en) * | 2009-11-16 | 2013-01-15 | University Health Network | Method of synthetic image generation for magnetic resonance imaging |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111803104A (en) * | 2020-07-20 | 2020-10-23 | 上海市第六人民医院 | Medical image display method, medium and electronic equipment |
US11582478B2 (en) * | 2020-09-08 | 2023-02-14 | Alibaba Group Holding Limited | Video encoding technique utilizing user guided information in cloud environment |
Also Published As
Publication number | Publication date |
---|---|
CN107664754B (en) | 2020-03-03 |
JP2018015307A (en) | 2018-02-01 |
JP6609226B2 (en) | 2019-11-20 |
US11187760B2 (en) | 2021-11-30 |
CN107664754A (en) | 2018-02-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6133537B2 (en) | MRI B1 mapping and B1 shimming | |
JP5735793B2 (en) | System for quantitatively separating seed signals in MR imaging | |
US8710840B2 (en) | Method and system for phase-sensitive magnetic resonance imaging | |
KR101955229B1 (en) | Method for setting an mri sequence | |
US8754645B2 (en) | Method for spatially resolved determination of an MR parameter | |
US10324153B2 (en) | System and method for flip angle determination in magnetic resonance imaging | |
US10132904B2 (en) | Generation of MR images with prospective motion correction and partial-volume-specific weighting of image information | |
US11122988B2 (en) | Magnetic resonance imaging apparatus | |
US7372268B2 (en) | Method and device for determination of coefficients of a magnetic resonance diffusion tensor | |
US20130307538A1 (en) | Magnetic resonance system and method to continuously correct phase errors of a magnetic resonance measurement sequence | |
US11393073B2 (en) | Magnetic resonance imaging apparatus and noise elimination method | |
US11187760B2 (en) | Magnetic resonance imaging apparatus and quantitative-value computing program | |
US11428769B2 (en) | Magnetic resonance imaging device, calculation device for generation of imaging parameter set, and imaging parameter set generation program | |
US11965945B2 (en) | Magnetic resonance system and shimming method and imaging method thereof | |
US20140253121A1 (en) | Magnetic resonance imaging apparatus and irradiation magnetic field distribution measurement method | |
JP7230149B2 (en) | Image processing method, image processing program, and computer obtained by magnetic resonance imaging | |
JP6721489B2 (en) | Magnetic resonance imaging apparatus and calculation image generation method | |
GB2574232A (en) | A method and apparatus for processing magnetic resonance images | |
US20210302522A1 (en) | Artificial intelligence based suppression of chemical species in magnetic resonance imaging | |
US10557909B2 (en) | Method and magnetic resonance tomography apparatus for diffusion imaging | |
JP6487554B2 (en) | Magnetic resonance imaging system | |
van Gorkum et al. | MRXCAT-CDTI: a numerical cardiac diffusion tensor imaging phantom | |
US10823803B2 (en) | Method and apparatus for reconstructing contrast levels from magnetic resonance acquisitions | |
US11454687B2 (en) | Multi-echo radial look-locker imaging for simultaneous T1 and fat fraction quantification | |
US11467238B2 (en) | Method for uniform reconstruction of multi-channel surface-coil magnetic resonance data without use of a reference scan |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: HITACHI, LTD., JAPAN Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:YOKOSAWA, SUGURU;TANIGUCHI, YO;OCHI, HISAAKI;AND OTHERS;REEL/FRAME:042612/0880 Effective date: 20170304 |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: FINAL REJECTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT VERIFIED |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
AS | Assignment |
Owner name: FUJIFILM HEALTHCARE CORPORATION, JAPAN Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:HITACHI, LTD.;REEL/FRAME:059115/0985 Effective date: 20220111 |