Magnetic resonance imaging device, image processing device, and imaging processing method

Abstract

To execute anatomical normalization for matching a functional image such as a blood flow image acquired by using an MRI device with a standard form accurately.SOLUTION: A form image in which a form is shown and a function image in which a function is shown are captured in the same imaging region of a subject. An arithmetic processing part executes processing for deforming the form image using deformation parameters, and moves a position of one or more structures in the form image to a position of a structure of a predetermined standard form respectively, and then deforms the function image using a value of the deformation parameters used when deforming the form image, thereby causing a position of a region in the function image to coincide with a position of a corresponding region of the standard form. Alternatively, the arithmetic processing part deforms the standard form in a reverse direction using the value of the deformation parameters, and causes the position of the region of the structure in the standard form to coincide with the position of the corresponding region in the function image.SELECTED DRAWING: Figure 5

Description

[Technical field]

[0001] The present invention relates to a magnetic resonance imaging (hereinafter referred to as "MRI") device or an image display device, and in particular to a technique for analyzing cerebral blood flow dynamics from a blood flow image. [Background technology]

[0002] An MRI device is a device that measures NMR signals generated by the nuclear spins that compose the tissues of a subject, particularly a human body, and creates two-dimensional or three-dimensional images of the morphology and function of the head, abdomen, limbs, etc. During imaging, the NMR signals are given different phase encodings by a gradient magnetic field and are frequency encoded, and are measured as time-series data. The measured NMR signals are reconstructed into an image by two-dimensional or three-dimensional Fourier transformation.

[0003] By evaluating the dynamics of cerebral blood flow, imaging diagnosis is performed for brain diseases in which abnormalities in cerebral blood flow are found, such as dementia, cerebral infarction, cerebrovascular disorders such as vascular stenosis, and epilepsy.

[0004] As an imaging method of cerebral blood flow using MRI, there is an arterial spin labeling (ASL) method (see Non-Patent Document 1). In the ASL method, the spin of protons in the blood passing through the label surface is inverted (labeling) by irradiating an RF pulse to a label surface set upstream of the blood flow from the imaging region. When this spin-inverted blood reaches the imaging region by the blood flow, it is exchanged with the spin of protons in the tissue of the imaging region in the capillaries, and the T1 relaxation time of the tissue changes. A blood flow image can be obtained by calculating the difference between a label image obtained by imaging the imaging region in this state and a control image obtained without inverting the blood protons. This method can label blood with an RF pulse, so that a blood flow image can be generated non-invasively.

[0005] Meanwhile, a method called anatomical standardization (normalization) is widely used to convert subjects' brain images obtained by PET (Positron Emission Tomography) or MRI, etc., into a standard brain coordinate system and align them to a standard brain. This anatomical standardization converts brain images with individual morphological differences into a standard brain, making it possible to compare blood flow and metabolic images pixel by pixel between subjects or between subject groups. Known methods of anatomical standardization include diffeomorphic anatomical registration through exponentiated Lie algebra (DARTEL) processing, which performs nonlinear transformation using multiple parameters.

[0006] Patent Document 1 also discloses a technique for obtaining a cerebral blood flow image by SPECT (Single Photon Emission Computed Tomography), which administers a drug labeled with a radioisotope to a subject, detects the emitted radiation to obtain projection data, and performs image reconstruction. In the technique of Patent Document 1, a standard brain image is obtained by transforming the reconstructed image to match a standard brain (anatomical normalization), and various quantitative value images are obtained by performing quantification of the standard brain image. It is also proposed to perform pathological analysis by quantitative evaluation by evaluating blood flow dynamics using the obtained quantitative value image. [Prior art documents] [Non-patent literature]

[0007] [Non-Patent Document 1] Kimura, Kabasawa, Yonekura, et al, Cerebral perfusion measurements using continuous arterial spin labeling: accuracy and limits of a quantitative approach, International Congress Series, 1256: 236-247, 2004 [Patent documents]

[0008] [Patent Document 1] JP 2006-119022 A Summary of the Invention [Problem to be solved by the invention]

[0009] The ASL method, which uses an MRI device to capture cerebral blood flow images, is advantageous for subjects in that it does not require the use of drugs and can be non-invasive compared to SPECT. In addition, if the cerebral blood flow images obtained by the ASL method can be converted to the standard brain coordinate system, it will be possible to compare the blood flow of subjects and healthy people in the same coordinate system, which will be useful for diagnosis, etc.

[0010] However, since the cerebral blood flow images taken by the ASL method do not contain information about the shape of the brain, it is difficult to convert the cerebral blood flow images accurately into the standard brain coordinate system.

[0011] Patent Document 1 discloses a method for determining deformation parameters for transforming a cerebral blood flow image obtained by SPECT to match a standard brain. However, because the cerebral blood flow image obtained by SPECT does not contain shape information, it is assumed that it is not easy to perform accurate coordinate conversion to match a standard brain.

[0012] Furthermore, the technique of obtaining a cerebral blood flow image by SPECT described in Patent Document 1 requires a longer examination time and is more expensive than imaging by an MRI device as in Non-Patent Document 1.

[0013] An object of the present invention is to perform anatomical normalization to accurately match functional images, such as blood flow images, acquired using an MRI apparatus to a standard form. [Means for solving the problem]

[0014] In order to achieve the above object, the present invention includes a static magnetic field generating unit that applies a static magnetic field to an imaging space in which a subject is placed, a gradient magnetic field generating unit that applies a gradient magnetic field to the imaging space, an irradiation coil that irradiates a high-frequency magnetic field to the subject in the imaging space, a receiving coil that receives a nuclear magnetic resonance signal from the subject, a measurement control unit that controls the gradient magnetic field generating unit, the irradiation coil, and the receiving coil to execute an imaging sequence and capture an image, and a calculation processing unit. The measurement control unit captures a morphological image showing a morphology and a functional image showing a function for the same imaging region of the subject. The calculation processing unit performs a process of transforming the morphological image using a deformation parameter and moving the positions of one or more structures in the morphological image to the positions of structures of a predetermined standard form, and then transforms the functional image using the value of the deformation parameter used when transforming the morphological image to match the position of the region in the functional image to the position of the corresponding region in the standard form, or transforms the standard form in the opposite direction using the value of the deformation parameter to match the position of the region of the structure of the standard form to the position of the corresponding region in the functional image. Effect of the Invention

[0015] According to the present invention, it is possible to perform anatomical normalization to accurately match functional images, such as blood flow images, acquired using an MRI device to a standard form, so that functional images can be easily evaluated by comparing them with functional images of healthy individuals. [Brief description of the drawings]

[0016] [Figure 1] 1 is a block diagram showing an example of the overall configuration of an MRI apparatus according to a first embodiment of the present invention. [Diagram 2] FIG. 2 is a diagram for explaining an example of an imaging pulse sequence in the ASL method. [Diagram 3] 4 is a flowchart showing the processing flow of the calculation processing unit 8 of the MRI apparatus of the embodiment 1. [Figure 4] FIG. 3 is a diagram showing an example of imaging conditions for the imaging pulse sequence of FIG. 2. [Diagram 5](a-1) to (a-3) are explanatory diagrams showing the procedure for anatomical normalization of 3DT1-weighted images, and (b-1) to (b-2) are explanatory diagrams showing the procedure for anatomical normalization of CBF images and ATT images. [Figure 6] FIG. 2 is a diagram showing an example of a GUI showing images and analysis results generated by a calculation processing unit 8 of the MRI apparatus of the first embodiment. [Figure 7] 10 is a flowchart showing the process flow of the calculation processing unit 8 of the MRI apparatus of the second embodiment. [Figure 8] 11 is a flowchart showing the process flow of the calculation processing unit 8 of the MRI apparatus of embodiment 3. [Figure 9] 10 is a flowchart showing the process flow of the calculation processing unit 8 of the MRI apparatus of embodiment 4. [Figure 10] 13 is a flowchart showing the process flow of the calculation processing unit 8 of the MRI apparatus of the fifth embodiment. [Figure 11] 13 is a flowchart showing the process flow of the calculation processing unit 8 of the MRI apparatus of the sixth embodiment. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0017] Hereinafter, preferred embodiments of the MRI apparatus of the present invention will be described in detail with reference to the accompanying drawings. In all the drawings for explaining the embodiments of the present invention, the same reference numerals are used for the components having the same functions, and the repeated description will be omitted.

[0018] <Embodiment 1> First, an overview of the MRI apparatus according to the first embodiment will be described.

[0019] The MRI apparatus of this embodiment 1 has a structure shown in FIG. 1 as an example, and includes a static magnetic field generating unit 2 that applies a static magnetic field to an imaging space in which a subject 1 is placed, a gradient magnetic field generating unit 3 that applies a gradient magnetic field to the imaging space, an irradiation coil 14a that irradiates a high-frequency magnetic field to the subject 1 in the imaging space, a receiving coil 14b that receives a nuclear magnetic resonance signal from the subject 1, a measurement control unit (sequencer) 4 that controls the gradient magnetic field generating unit 3, the irradiation coil 14a, and the receiving coil 14b to execute an imaging sequence and capture an image, and an arithmetic processing unit 8.

[0020] The measurement control unit 4 captures a morphological image showing the morphology and a functional image showing the function for the same imaging region of the subject 1.

[0021] The arithmetic processing unit 8 performs a process (anatomical normalization process) of deforming the morphological image using the transformation parameters and moving the positions of one or more structures in the morphological image to the positions of structures in a predetermined standard form using the transformation parameters. The arithmetic processing unit 8 deforms the functional image using the values ​​of the transformation parameters used when transforming the morphological image, thereby matching the positions of areas in the functional image to the positions of corresponding areas in the standard form.

[0022] Alternatively, the calculation processing unit 8 uses the values ​​of the above-mentioned transformation parameters to transform the standard form in the inverse direction, thereby making the position of the region of the structure in the standard form coincide with the position of the corresponding region in the functional image.

[0023] This allows anatomical normalization to be performed by accurately matching the coordinates of functional images that do not show the brain morphology, such as blood flow images, to the coordinates of the standard morphology using an MRI device. Therefore, functional images such as blood flow images of the subject 1 and a healthy person can be compared in the same coordinate system, which is useful for diagnosis, etc.

[0024] The morphological image referred to here is an image showing the morphology of the subject, and a T1 weighted image is used here. Note that the image is not limited to a T1 weighted image, and may be any image showing the morphology of the subject, such as an absolute value image, a T2 weighted image, or a proton density weighted image.

[0025] A functional image is an image that represents the functions of a subject. In this embodiment, an example will be described in which the functional image is an arterial spin labeling (ASL) image, particularly an arterial transit time (ATT) image calculated from a multi-PLD ASL image, and / or a cerebral blood flow (CBF) image of cerebral blood flow images. Note that the images are not limited to these, and any image that represents the functions of the subject may be used. For example, images using various physical property values and quantitative values such as T2, T2*, diffusion coefficient, flow velocity, magnetic susceptibility, elastic modulus, and contrast agent concentration as pixel values, or a fluid attenuated IR (FLAIR) image with water signal suppressed can be used.

[0026] On the other hand, as the standard form, any form may be used as long as the coordinates of one or more structures are determined. For example, the Talairach standard brain can be used.

[0027] As described above, in this embodiment, instead of directly normalizing a functional image such as a cerebral blood flow image with respect to the standard form (standard brain, etc.), first, the morphological image captured for the same imaging region is deformed to match the coordinates of the structures in the standard form (standard brain, etc.) (normalized), and the functional image is deformed with the value of the deformation parameter used for the deformation. Thus, even when the functional image does not include morphological information, the functional image can be accurately normalized.

[0028] Note that the morphological image is preferably captured for the same imaging region continuously with the capture of the functional image, but as long as it is captured for the same imaging region of the same subject, it may be captured before or after.

[0029] <Configuration of MRI apparatus> Hereinafter, the MRI apparatus of this embodiment will be described in detail.

[0030] The structure of the MRI apparatus according to the present invention will be described in detail with reference to Fig. 1. Fig. 1 is a block diagram showing the overall configuration of one embodiment of the MRI apparatus according to the present invention. This MRI apparatus obtains a tomographic image of a subject by utilizing the NMR phenomenon, and as shown in Fig. 1, the MRI apparatus is configured to include a static magnetic field generating system (static magnetic field generating unit) 2, a gradient magnetic field generating system (gradient magnetic field generating unit) 3, a transmitting system 5, a receiving system 6, a signal processing system 7, a sequencer (measurement control unit) 4, and a central processing unit (CPU) (arithmetic processing unit) 8.

[0031] The static magnetic field generating system 2 applies a static magnetic field to the imaging space in which the subject is placed. If the static magnetic field generating system 2 is a vertical magnetic field type, it generates a uniform static magnetic field in a direction perpendicular to the body axis of the subject 1, and if it is a horizontal magnetic field type, it generates a uniform static magnetic field in the direction of the body axis of the subject 1. A static magnetic field generating source of the permanent magnet type, normal conducting type, or superconducting type is arranged around the subject 1.

[0032] The gradient magnetic field generating system 3 applies a gradient magnetic field to the imaging space. The gradient magnetic field generating system 3 includes a gradient magnetic field coil 9 that applies gradient magnetic fields in three axial directions of X, Y, and Z, which are the coordinate system (stationary coordinate system) of the MRI apparatus, and a gradient magnetic field power supply 10 that drives each gradient magnetic field coil. The gradient magnetic field power supply 10 of each coil is driven according to an instruction from a sequencer 4 described later, thereby applying gradient magnetic fields Gx, Gy, and Gz in three axial directions of X, Y, and Z. During imaging, a slice selection gradient magnetic field pulse (Gs) is applied in a direction perpendicular to the slice plane (imaging cross section) to set a slice plane for the subject 1, and a phase encoding gradient magnetic field pulse (Gp) and a frequency encoding gradient magnetic field pulse (Gf) are applied in the remaining two directions perpendicular to the slice plane and perpendicular to each other, to encode position information in each direction in an echo signal.

[0033] The sequencer 4 controls the gradient magnetic field generating system 2, the transmission system 5, and the receiving system 6 to execute an imaging sequence and capture an image. That is, the sequencer 4 is a control means that repeatedly applies radio frequency magnetic field pulses (hereinafter referred to as "RF pulses") and gradient magnetic field pulses in a certain predetermined pulse sequence, and operates under the control of the calculation processing unit 8 to send various commands required for data collection of tomographic images of the subject 1 to the transmission system 5, the gradient magnetic field generating system 3, and the receiving system 6.

[0034] The transmission system 5 irradiates the subject 1 with an RF pulse (high frequency magnetic field) to cause nuclear magnetic resonance in the nuclear spins of atoms constituting the biological tissue of the subject 1, and is composed of a high frequency oscillator 11, a modulator 12, a high frequency amplifier 13, and a transmission side high frequency coil (transmission coil) 14a. The RF pulse output from the high frequency oscillator 11 is amplitude modulated by the modulator 12 at a timing according to a command from the sequencer 4, and this amplitude modulated RF pulse is amplified by the high frequency amplifier 13 and then supplied to the high frequency coil 14a arranged close to the subject 1, so that the RF pulse is irradiated to the subject 1.

[0035] The receiving system 6 receives an echo signal (NMR signal) emitted by the nuclear magnetic resonance of the atomic spins constituting the biological tissue of the subject 1, and is composed of a receiving radio frequency coil (receiving coil) 14b, a signal amplifier 15, a quadrature phase detector 16, and an A / D converter 17. The NMR signal of the response of the subject 1 induced by the electromagnetic wave irradiated from the transmitting radio frequency coil 14a is detected by the radio frequency coil 14b arranged close to the subject 1, amplified by the signal amplifier 15, and then split into two orthogonal systems of signals by the quadrature phase detector 16 at a timing according to a command from the sequencer 4. Each of the signals is converted into a digital quantity by the A / D converter 17 and sent to the signal processing system 7. The signal processing system 7 processes various data and displays and stores the results of the processing, etc. The signal processing system 7 has external storage devices such as an optical disk 19 and a magnetic disk 18, and a display 20 such as a CRT.

[0036] A processor (CPU) 8 receives the data from the receiving system 6, performs predetermined signal processing, and then executes image reconstruction processing to generate a tomographic image of the subject 1. The generated image is displayed on a display 20 and is also recorded in an external storage device such as a magnetic disk 18.

[0037] The operation unit 25, which is used to input various control information for the MRI apparatus and control information for processing performed by the signal processing system 7, is made up of a trackball or mouse 23 and a keyboard 24. This operation unit 25 is disposed close to the display 20, and an operator interactively controls various processes of the MRI apparatus through the operation unit 25 while looking at the display 20.

[0038] 1, the high frequency coil 14a on the transmitting side and the gradient magnetic field coil 9 are installed in the static magnetic field space of the static magnetic field generating system 2 into which the subject 1 is inserted, facing the subject 1 in the case of a vertical magnetic field system, or surrounding the subject 1 in the case of a horizontal magnetic field system. On the other hand, the high frequency coil 14b on the receiving side is installed facing or surrounding the subject 1.

[0039] The currently widely used MRI imaging system is the hydrogen nucleus (proton), which is the main constituent of the subject. By imaging the spatial distribution of proton density and the spatial distribution of the relaxation time of the excited state, the morphology and function of the human head, abdomen, limbs, etc. can be imaged in two or three dimensions.

[0040] <Capturing functional images> The measurement control unit (sequencer) 4 captures an image of cerebral blood flow as a functional image by the ASL (arterial spin labeling) method. An example of an imaging sequence is shown in FIG. 2. There are two types of ASL labeling methods: Pulsed ASL and pseudo-continuous ASL (pASL), and the sequence in FIG. 2 is an example of a sequence in the pASL method. The measurement control unit (sequencer) 4 executes the imaging sequence in FIG. 2 by controlling the gradient magnetic field generating unit 3, the irradiation coil 14a, the receiving coil 14b, etc. In the ASL method, by adjusting the RF pulse and the gradient magnetic field, the blood passing through the label surface set upstream of the blood flow from the imaging region is irradiated with a high-frequency pulse that inverts protons to label the blood, and then an image (label image) is captured at the timing when the delay time PLD (post labeling delay time: time elapsed after labeling) at which the blood reaches the imaging region has elapsed. Similarly, an image (control image) is captured for the same imaging region without inverting blood protons. Blood flow information is visualized by calculating the difference between the labeled image and the control image.

[0041] In the sequence of the pASL method shown in Figure 2, after labeling the blood passing through the label surface, imaging is performed after a delay time PLD (post labeling delay time) has elapsed in order to wait for the labeled blood to flow into the region of interest.

[0042] Specifically, the measurement control unit (sequencer) 4 executes a label imaging sequence for imaging a label image and a control imaging sequence for imaging a control image, as shown in FIG. 2. Both of these are sequences for sequentially irradiating and applying the following RF pulses 51 to 58 and gradient magnetic field pulses. First, a gradient magnetic field is applied to dephase the transverse magnetization while irradiating the imaging region with a pre-saturation pulse 51 that tilts the spins by 90 degrees. Next, an IR pulse 52, which is an inversion recovery pulse, is selectively or non-selectively irradiated to the imaging region. Next, in the label imaging sequence, an ASL pulse 53 that inverts blood protons by 180 degrees is selectively irradiated to the label surface (for example, a region 1-2 cm from the lower end of the cerebellum). In the control imaging sequence, a pulse 53 that does not invert blood protons is similarly irradiated to the label surface. Next, IR pulses 54 and 55, which are inversion recovery pulses, are selectively irradiated to the imaging region. Furthermore, the imaging region is irradiated with a fat suppression pulse 56, followed by non-selective irradiation with a blood flow signal suppression pulse 57. After the ASL pulse 53 is irradiated to the imaging region, a readout process 58 is executed at a timing when a predetermined delay time (PLD) has elapsed, and an NMR signal is acquired from the imaging region. Specifically, in the readout process 58, an RF pulse is irradiated to the imaging region, thereby exciting spins, and then an NMR signal is acquired by the receiving coil 14b while applying a readout gradient magnetic field and a phase encoding gradient magnetic field.

[0043] As a result, in the label imaging sequence, an NMR signal of the imaging region in a state where labeled blood has reached the imaging region is acquired, and the arithmetic processing unit 8 reconstructs a label image from the acquired NNR signal. Also, in the control imaging sequence, an NMR signal is acquired for the same imaging region without labeling the blood, and the arithmetic processing unit 8 reconstructs a control image from the acquired NNR signal.

[0044] Moreover, the measurement control unit (sequencer) 4 repeatedly executes the imaging sequence while changing the delay time (PLD) (the time from the end of application of the ASL pulse 53 to the start of the readout 58) (Multi-PLD ASL method). As a result, the calculation processing unit 8 reconstructs a plurality of label images with different delay times (PLD).

[0045] The calculation processing unit 8 calculates a cerebral blood flow (hereinafter referred to as CBF) image and an arterial transit time (hereinafter referred to as ATT) image from the obtained ASL image.

[0046] <Normalization of Functional Images> Next, the arithmetic processing unit 8 applies anatomical normalization to the ATT image and the CBF image. Furthermore, the arithmetic processing unit 8 performs ROI analysis and Voxel analysis on the anatomically normalized ATT image and the CBF image, and performs statistical analysis. As a result, the arithmetic processing unit 8 visualizes and displays abnormal blood flow based on the results of the statistical analysis.

[0047] <Control and arithmetic processing of the arithmetic processing unit 8> The control and calculation processing of the calculation processing unit 8 will be specifically described with reference to the flowchart of Fig. 3. Fig. 4 is a diagram showing an example of imaging conditions, Fig. 5 shows an image obtained in the steps of the flow of the first embodiment, and Fig. 6 shows an example of a screen displaying the imaging results of the first embodiment.

[0048] 3 of the arithmetic processing unit 8 is realized by software, in which a CPU in the arithmetic processing unit 8 reads and executes a program previously stored in a magnetic disk 18 or the like. It is also possible to realize a part or the whole of the arithmetic processing unit 8 by hardware. For example, a part or the whole of the arithmetic processing unit 8 may be configured using a custom IC such as an ASIC (Application Specific Integrated Circuit) or a programmable IC such as an FPGA (Field-Programmable Gate Array), and a circuit design may be performed to realize the function.

[0049] First, the operator places the subject 1 on a bed in the imaging area.

[0050] (Step S201) The calculation processing unit 8 accepts settings of imaging conditions for imaging sequences of Multi-PLD ASL images, 3D T1 weighted images, and proton density weighted images from the user via the operation unit 25. For example, the settings of the imaging conditions shown in FIG. 4 are accepted. The imaging region is the whole brain region in all of the imaging sequences of Multi-PLD ASL, TI weighted images, and proton density weighted images. In the imaging sequence of Multi-PLD ASL, seven types of delay times (PLD) are set here, ranging from 500 ms to 3000 ms.

[0051] (Step S202) The calculation processing unit 8 instructs the sequencer (measurement control unit 4) 4 to execute the Multi-PLD ASL imaging sequence shown in Fig. 2 under the imaging conditions set in step S201. As a result, the sequencer 4 controls the gradient magnetic field generating unit 3, the irradiation coil 14a, the receiving coil 14b, etc., and executes the imaging sequence. The calculation processing unit (CPU) 8 receives data from the receiving system 6 and reconstructs an ASL image for each delay time (PLD) of the Multi-PLD of the imaging region.

[0052] Moreover, the calculation processing unit 8 causes the sequencer (measurement control unit) 4 to execute, for example, a known gradient echo sequence before or after the Multi-PLD ASL imaging sequence. As a result, the calculation processing unit (CPU) 8 receives data from the receiving system 6 and reconstructs a 3D T1 weighted image of the imaging region. Furthermore, a proton density weighted image of the imaging region is captured using a known GRACE sequence or the like that is the same as the ASL imaging sequence. Note that the imaging region, FOV, and Matrix of the proton density weighted image are the same as those of the ASL. If they are different, registration of the ASL image and the proton density weighted image may be performed after reconstruction.

[0053] The arithmetic processing unit 8 stores each reconstructed image data on the magnetic disk 18. (Step S203) The calculation processing unit 8 receives a selection of an ASL image of a delay time (PLD) used for calculating CBF and ATT from the operator via the operation unit 25. Note that this selection of an ASL image may be performed automatically. (Step S204) The arithmetic processing unit 8 calculates the ATT and CBF for each pixel from the pixel values ​​of the Multi-PLD ASL image selected in step S203, and generates an ATT image and a CBF image. For example, the theoretical value of the ASL signal is calculated for each delay time (PLD) using a theoretical formula based on the two compartment model of formulas (1) and (2). The change in the theoretical value of the ASL signal due to the delay time (PLD) is fitted to the actual measured value of the ASL signal of each PLD (i.e., the pixel value of the ASL image captured for each PLD) by curve fitting based on the nonlinear least squares method or the like, thereby calculating the ATT and CBF for each pixel. From the ATT value and CBF value for each pixel, an ATT image displaying the value of ATT as a pixel value and a CBF image displaying the value of CBF can be generated. Note that the solution of the theoretical formulas (1) and (2) can be a known method described in the above-mentioned Non-Patent Document 1, so a detailed explanation will be omitted here.

number

number

[0054] In formula (1), t is time, T1m is the apparent T1 value inside the blood vessel, and T1e is the apparent T1 value outside the blood vessel. M0 is calculated by using the magnetization (pixel value) of the proton density weighted image, and M m0 is the value in the capillary of the proton density weighted image, and M e0 is the extravascular value of the proton density bridge image. Other variables and counts are as in Table 1 below. [Table 1]

[0055] Here, the ATT and CBF are calculated using the nonlinear least squares method based on the Two Compartment Model using equations (1) and (2), but they may be calculated using other known theoretical models or fitting techniques.

[0056] (Step S205) The calculation processing unit 8 executes a process (anatomical normalization) of deforming the 3D T1 weighted image captured in step S202 to match it with a standard brain, and obtains a deformation field for matching the 3D T1 weighted image with the standard brain.

[0057] Specifically, the arithmetic processor 8 first performs a process of moving the positions and shapes of structures (gray matter, white matter, etc.) in the 3D T1 weighted image to the coordinates of a predetermined standard brain structure using a known method, as shown in Figs. 5(a-1) to (a-3) (anatomical normalization (standardization)). For example, normalization is performed using diffeomorphic anatomical registration through exponentiated Lie algebra (DARTEL) processing, which performs nonlinear transformation using a large number of parameters. In the DARTEL processing, the target brain structure is deformed (anatomical normalization) according to a standard brain template created in advance.

[0058] Next, the calculation processing unit 8 saves the normalization parameters (a set of values ​​set as the parameters of the DARTEL transformation form) used when transforming the target brain as transformation parameters.

[0059] Here, the method of acquiring the anatomical normalization parameters is not limited to the DARTEL process, and may be any other known normalization method.

[0060] (Step S206) The calculation processing unit 8 registers the ATT image and the CBF image generated in step S204 with the 3D T1 weighted image captured in step S202, as shown in Fig. 5 (b-1) and (b-2). For example, the ATT image and the CBF image are registered with the 3D T1 weighted image by a known registration technique.

[0061] As a method for positioning, for example, the mutual information MI between the ATT image and the 3D T1 weighted image, and the CBF image and the 3D T1 weighted image are calculated using the following equations (3) to (6), and the positions are aligned (for example, by rotating and translating) so that the mutual information MI is maximized.

number

number

number

number

[0062] Here, ai is the pixel value (grayscale) of a pixel in the ATT image or CBF image, and bj is the pixel value (grayscale) of the corresponding pixel in the 3D T1 weighted image. h(ai,bj) is a two-dimensional histogram that counts the frequency of pixel value (grayscale) combinations for all corresponding pixel value combinations (ai,bj) and maps them. p(ai,bj) represents the probability that pixel values ​​ai and bj occur simultaneously (joint probability) and is calculated using formula (4). p(ai) and p(bj) represent the probability that pixel values ​​ai and bj occur, respectively (marginal probability), and are calculated using formulas (5) and (6), respectively.

[0063] (Step S207) The calculation processing unit 8 uses the values ​​of the deformation parameters stored in step S205 to deform each of the ATT image and the CBF image after the alignment in step S206.

[0064] This allows for anatomical normalization of ATT and CBF.

[0065] For the transformation process, it is preferable to use the same process as in step S205. Here, the DARTEL process is used, but other known transformation processes may be used.

[0066] (Step S208) The calculation processing unit 8 performs an analysis of cerebral blood flow dynamics using the ATT image and CBF image anatomically normalized in step S207. For example, from a healthy database (a database of anatomically normalized ATT images and CBF images of healthy people, for example, 50 people), it calculates the degree of deviation between the normalized ATT image and CBF image of the subject calculated in step 207. For example, it calculates a z-score for each voxel, and calculates the degree of deviation (deviation) of the normalized ATT image and CBF image of the subject from the healthy image as a numerical value (z-score).

[0067] (Step S209) The calculation processing unit 8 displays the deviation (Z score) of the blood flow rate (CBF image) and / or the blood flow arrival time (ATT image) calculated in step 208 on the display 20. As a display screen, for example, as shown in FIG. 6, each cross-sectional image of the CBF image and the ATT image and a 3D image are displayed, and the z-score is displayed in grayscale on top of them. The ATT image and / or the CBF image to be displayed may be a normalized image or an image before normalization. In this embodiment, the deviation (Z score) of the blood flow rate (CBF image) and / or the blood flow arrival time (ATT image) is displayed in grayscale in FIG. 6, but it is of course possible to display it in color.

[0068] The above is the description of the processing flow of the arithmetic processing unit 8 in the first embodiment.

[0069] As described above, in the first embodiment of the present invention, the calculation processing unit 8 can perform anatomical normalization after calculating functional images (ATT image and CBF image) that do not include shape information of the brain. This allows voxel analysis to be performed on the anatomically normalized functional images (ATT image and CBF image) and statistical analysis, etc. Therefore, from the results of the statistical analysis, abnormal blood flow that deviates from the image of a healthy person can be visualized and displayed, and the cerebral blood flow in dementia, cerebrovascular disorder, epilepsy, etc., where a decrease in cerebral blood flow is observed, can be evaluated using only the MRI device.

[0070] In the above-mentioned step 203, the selection of the images (Multi-PLD ASL images) to be used for calculating the functional images (ATT image and CBF image) is accepted from the user, but the calculation processing unit 8 may automatically select the images according to predetermined criteria.

[0071] It is also possible to store the selected image on the magnetic disk 18 or the like, and perform steps 204 and thereafter at a later date.

[0072] <<Embodiment 2>> An MRI apparatus according to a second embodiment of the present invention will be described.

[0073] The MRI apparatus of the second embodiment calculates an ATT image and a CBF image, and then applies anatomical normalization to the ATT image and the CBF image, as in the first embodiment. The MRI apparatus performs ROI analysis and statistical analysis on the anatomically normalized ATT image and CBF image, thereby visualizing and displaying abnormal blood flow.

[0074] Fig. 7 is a processing flow of the arithmetic processing unit 8 of embodiment 2. Hereinafter, the processing of the arithmetic processing unit 78 will be described with reference to Fig. 7. Note that in the flow of Fig. 7, the same steps as those in the flow of Fig. 3 of embodiment 1 are given the same step reference numerals, and the description thereof will be omitted.

[0075] (Steps S201 to S207) The calculation processing unit 8 performs steps S201 to S207 similar to steps S201 to S207 in the first embodiment to generate anatomically normalized ATT images and CBF images.

[0076] (Step S301) The arithmetic processor 8 sets one or more ROIs (regions of interest) for the ATT image and the CBF image anatomically normalized in step S207. The ROI may have the shape of a standard brain structure, or may have a desired shape drawn by the user. For example, an ROI from a known brain atlas such as the AAL (Automated Anatomical Labeling) atlas may be selected, or the user may draw an ROI on the image.

[0077] The calculation processing unit 8 performs an analysis of the cerebral blood flow dynamics in the ROI of the anatomically normalized ATT image and CBF image. For example, the degree of deviation of the ATT image and CBF image from the healthy database is calculated as a z-score for each ROI, and the degree of deviation of the blood flow arrival time and blood flow rate from healthy people is analyzed. Specifically, for example, a z-score is obtained that indicates how much the blood flow arrival time or blood flow rate of the ROI of the subject deviates from the average value of the blood flow arrival time or blood flow rate of the ROI of 50 people in the healthy database. This allows the value (z-score) indicating the degree of deviation to be calculated for each ROI when there are multiple regions (ROIs).

[0078] (Step S209) The calculation processing unit 8 displays the ATT image and / or CBF image and the z-score calculated for the ROI on the display 20, similar to step S209 in the first embodiment. For example, the z-score is displayed superimposed on the position of the ROI on the ATT image or CBF image.

[0079] As described above, in the second embodiment of the present invention, the z-score is calculated from the results of the statistical analysis and can be displayed together with the ATT image and / or the CBF image, so that the ROI with a large z-score, that is, deviating from the image of a healthy person, can be visualized and displayed. As a result, as in the first embodiment described above, it becomes possible to evaluate the cerebral blood flow in dementia, cerebrovascular disorder, epilepsy, etc., where a decrease in cerebral blood flow is observed, by using only MRI.

[0080] In addition, in the second embodiment, an ROI is set and the inside of the ROI is analyzed, so it is possible to determine easily and with a small amount of calculation whether the blood flow in the ROI is reduced compared to that of a healthy person.

[0081] In this embodiment, an example has been described in which a known ROI atlas such as an AAL atlas or an ROI drawn by a user is set as the ROI in step S301, but the ROI may be set by the calculation processing unit 8. For example, the calculation processing unit 8 may extract a region in which there is a significant difference in blood flow between the healthy database and the disease database in advance, define it as the ROI, and set it as the ROI in step S301.

[0082] Furthermore, the calculation processing unit 8 may compare the ATT image and / or CBF image generated in step S207 with images in the healthy subject database, extract a region having a significant difference, and set it as the ROI in step S301. <<Embodiment 3>> An MRI apparatus according to a third embodiment will be described.

[0083] The MRI apparatus of the third embodiment has a function of calculating the probability that a subject has a disease by using the ATT image and the CBF image.

[0084] Fig. 8 is a flowchart showing the process flow of the arithmetic processing unit 8 of embodiment 3. Each step in Fig. 8 will be described below.

[0085] (Steps S201 to S208) This is the same as in the first or second embodiment, and the calculation processing unit 8 generates anatomically normalized ATT images and CBF images.

[0086] (Step S401) The arithmetic processing unit 8 sets a predetermined ROI in advance on the ATT images and CBF images in the healthy database (a database of anatomically normalized ATT images and CBF images of healthy people) and the disease database (a database of anatomically normalized ATT images and CBF images of people with diseases) using, for example, a brain atlas, and obtains a feature amount for each ROI. For example, the average value of pixel values ​​(ATT value or CBF value) can be used as the feature amount. The arithmetic processing unit 8 obtains in advance a boundary (discrimination surface) between the distribution of feature amounts of the ATT images and CBF images of multiple healthy people and the distribution of feature amounts of multiple people with diseases. A known method (for example, a machine learning algorithm such as a Support Vector Machine) is used as a method for obtaining the discriminant surface.

[0087] The calculation processing unit 8 sets the above-mentioned predetermined ROI on the normalized ATT image and CBF image of the subject obtained in step 208, calculates the feature amount, and calculates whether the calculated feature amount is on the healthy person side or the diseased person side with respect to the discrimination plane, as well as the distance from the discrimination plane. Based on the calculated distance, the probability that the subject belongs to a person with a disease is calculated by a known method.

[0088] (Step S402) The calculation processing unit 8 displays the normalized ATT image and CBF image, and the result calculated in step S401 (probability that the subject belongs to a person with a disease). At this time, the distance calculated in step S401 may be displayed together with the probability.

[0089] The above is the description of the processing flow of the arithmetic processing unit 8 in the third embodiment.

[0090] According to the MRI apparatus of the third embodiment described above, the probability that the subject has a disease can be calculated from the ATT image and the CBF image, thereby assisting doctors in making diagnoses.

[0091] The method of calculating the discriminant surface (boundary) is not limited to machine learning algorithms such as Support Vector Machine, and any method that can distinguish the similarity of the ATT images and CBF images of a person with a disease and a healthy person may be used. For example, a deep learning method other than Support Vector Machine, clustering, statistical analysis, etc. may be used.

[0092] <<Embodiment 4>> An MRI apparatus according to a fourth embodiment will be described.

[0093] In the fourth embodiment, blood flow analysis is performed using a single PLD image (blood flow image) instead of a multi-PLD ASL. That is, a CBF image is calculated from a single PLD image, and anatomical normalization and statistical analysis are performed as in the first and second embodiments, and the results are displayed.

[0094] Fig. 9 is a flowchart showing the processing flow of the calculation processing unit 8 according to the fourth embodiment of the present invention. Each step in Fig. 9 will be described below. Note that, in the case where the same processing as that in the step in Fig. 3 is performed in the step in Fig. 9, the description will be omitted.

[0095] (Steps S901 to S903) Similar to steps 201 to 202 in FIG. 3, the calculation processing unit 8 accepts settings of imaging conditions from the user via the operation unit 25, and then executes the imaging sequence of the Multi-PLD ASL in FIG. 2 of the first embodiment in a Single PLD (one delay time PLD) to capture a label image and a control image.

[0096] Moreover, the calculation processing unit 8 captures a 3DT1 weighted image and a proton density weighted image.

[0097] The calculation processing unit 8 selects a label image and a control image captured with a single PLD (one delay time PLD) as images to be used in calculating the CBF.

[0098] (Step S904) The calculation processing unit 8 calculates a CBF image from the image of the single PLD and the proton density. For example, the calculation method uses equation (7) and calculates the CBF (f in equation (7)) by inputting the difference ΔM in pixel values ​​between the label image and the control image, the delay time PLD, and the pixel value of PD to equation (7). The calculation processing unit 8 generates a CBF image by calculating the CBF value for each pixel value.

[0099]

number

[0100] (Steps S205 to S209) The calculation processing unit 8 executes anatomical normalization to match the 3DT1 weighted image with a standard brain, similar to steps 205 to 209 in FIG. 3 of the first embodiment, and anatomically normalizes the CBF image using the transformation parameters at that time. Then, it calculates the degree of deviation (z-score) indicating how much the normalized CBF image deviates from the CBF image of the healthy subject database. The calculated degree of deviation is displayed together with the CBF image.

[0101] The above is the description of the processing flow of the arithmetic processing unit 8 in the fourth embodiment.

[0102] As described above, in the fourth embodiment of the present invention, the CBF of a single PLD can be used to analyze hemodynamics (CBF). Therefore, even when imaging of a multi-PLD ASL is difficult, the fourth embodiment of the present invention makes it possible to analyze hemodynamics using an ASL image of a single PLD.

[0103] <<Embodiment 5>> The MRI apparatus according to the fifth embodiment will be described.

[0104] In the fifth embodiment, the ROI shape of the brain atlas is transformed to match the brain morphology of an individual (subject), and then analysis is performed in the individual brain space. That is, unlike the first to fourth embodiments, in which the ATT image and the CBF image are transformed to match a standard brain using transformation parameters, the brain atlas is transformed to match the individual brain space of the subject by applying the transformation parameters in the reverse direction to the brain atlas, and the ATT image and the CBF image are analyzed using the ROI of the transformed brain atlas.

[0105] 10 is a flowchart showing the process flow of the arithmetic processing unit 8 according to the third embodiment of the present invention. Each step in FIG. 10 will be described in detail below.

[0106] (Steps S201 to S205) The calculation processing unit 8 performs steps S201 to S205 in the same manner as in the first embodiment.

[0107] (Step S601) The calculation processing unit 8 uses the deformation parameters stored in step S205 to deform the brain atlas in the reverse direction, whereby the brain atlas is deformed into a shape that matches the 3DT1 weighted image of the subject before normalization.

[0108] Furthermore, the calculation processing unit 8 registers the inversely transformed brain atlas to the ATT image and the CBF image. The transformation method used is, for example, the DARTEL method.

[0109] (Step S602) The calculation processing unit 8 calculates the degree of discrepancy from a healthy subject, in the same manner as in step S208, using the brain atlas that matches the brain space of the subject calculated in step S601, the ATT image, and the CBF image.

[0110] (Step S209) The calculation processing unit 8 displays the degree of deviation and the ATT image and CBF image, similarly to step S209 in the first embodiment.

[0111] The above is the description of the processing flow of the arithmetic processing unit 8 in the fifth embodiment.

[0112] As described above, in the fifth embodiment of the present invention, the atlas is inversely transformed to analyze the dynamics of cerebral blood flow. This makes it possible to avoid changes in values ​​of CBF and ATT due to the standard brain transformation, enabling more accurate blood flow evaluation.

[0113] <<Embodiment 6>> The MRI apparatus according to the sixth embodiment will be described.

[0114] In the sixth embodiment, in addition to the ATT images, CBF images, and 3D T1 weighted images, QSM (Quantitative Susceptibility Mapping) images are also used to perform statistical analysis.

[0115] 11 is a flowchart showing the process flow of the arithmetic processing unit 8 according to the third embodiment of the present invention. Each step in FIG. 11 will be described in detail below.

[0116] (Steps S1101 to S1102) The calculation processing unit 8 accepts settings of imaging conditions for a Multi-PLD ASL image, a 3D T1 weighted image, and a proton density weighted image, and performs imaging, similarly to steps S201 to S202 in the first embodiment; however, in the present embodiment, in addition to these, it also accepts settings of imaging conditions for a QSM image, and performs imaging.

[0117] (Steps S203 to S205) The calculation processing unit 8 performs steps S203 to S205 of the first embodiment to generate a CBF image and an ATT image, and obtains deformation parameters for matching the 3D T1 weighted image to the standard brain.

[0118] (Step S1106) The calculation processing unit 8 performs the QSM analysis. For example, a known method is used to create a global magnetic field map by removing phase wrapping from the collected phase images, and a local magnetic field map by removing background magnetic fields is obtained. In response to this, a magnetic susceptibility map (QSM) is obtained by estimating the dipole magnetic field. (Steps S1107 to S1109) Similar to steps S206 to S208 in embodiment 1, the calculation processing unit 8 aligns the CBF image, ATT image, and QSM image to the 3DT1 weighted image, and then anatomically normalizes them by deforming them using the deformation parameters saved in step S205.

[0119] Furthermore, the calculation processing unit 8 performs statistical analysis using the anatomically normalized ATT image, CBF image, QSM image, and 3DT1 weighted image. For example, similarly to the first embodiment, the degree of deviation from the image of a healthy subject is calculated.

[0120] The calculation processing unit 8 can perform analysis using multiple image types, such as ATT images and CBF images, as well as QSM images and 3DT1 weighted images, and can simultaneously calculate the discrepancy in blood flow, magnetic susceptibility, and brain volume in the same region.

[0121] (Step S209) The calculation processing unit 8 displays the analysis results as an image.

[0122] The above is the description of the processing flow of the arithmetic processing unit 8 in the sixth embodiment.

[0123] As described above, in the sixth embodiment of the present invention, by performing analysis using a plurality of image types, it becomes possible to simultaneously calculate the discrepancy of blood flow, magnetic susceptibility, and brain volume in the same region. By simultaneously analyzing blood flow, magnetic susceptibility, and brain volume, diagnosis can be facilitated.

[0124] Although the first to sixth embodiments of the present invention have been described above, the present invention is not limited to these embodiments. [Explanation of symbols]

[0125] 1: subject, 2: static magnetic field generating system, 3: gradient magnetic field generating system, 4: sequencer, 5: transmission system, 6: reception system, 7: signal processing system, 8: calculation processing unit (CPU), 9: gradient magnetic field coil, 10: gradient magnetic field power supply, 11: high frequency oscillator, 12: modulator, 13: high frequency amplifier, 14a: high frequency coil (transmission coil), 14b: high frequency coil (reception coil), 15: signal amplifier, 16: quadrature phase detector, 17: A / D converter, 18: magnetic disk, 19: optical disk, 20: display, 21: ROM, 22: RAM, 23: trackball or mouse, 24: keyboard

Claims

1. a static magnetic field generating unit that applies a static magnetic field to an imaging space in which a subject is placed, a gradient magnetic field generating unit that applies a gradient magnetic field to the imaging space, a transmitting coil that transmits a radio frequency magnetic field to the subject in the imaging space, a receiving coil that receives a nuclear magnetic resonance signal from the subject, a measurement control unit that controls the gradient magnetic field generating unit, the transmitting coil and the receiving coil to execute an imaging sequence and capture an image, and a calculation processing unit, the measurement control unit captures a morphological image showing a morphology and a functional image showing a function for the same imaging region of the subject; The magnetic resonance imaging apparatus is characterized in that the calculation processing unit deforms the morphological image using deformation parameters and performs a process of moving the positions of one or more structures in the morphological image to the positions of the structures in a predetermined standard form, and then deforms the functional image using the values ​​of the deformation parameters used when deforming the morphological image, thereby matching the positions of areas in the functional image to the positions of corresponding areas in the standard form, or deforms the standard form in the opposite direction using the values ​​of the deformation parameters, thereby matching the positions of the areas of the structures in the standard form to the positions of the corresponding areas in the functional image.

2. 2. A magnetic resonance imaging apparatus as described in claim 1, characterized in that the calculation processing unit performs processing to match the position of an area in the functional image with the position of a corresponding area in the morphological image before deforming the functional image.

3. 2. The magnetic resonance imaging apparatus according to claim 1, wherein the functional image is a cerebral blood flow image, A magnetic resonance imaging apparatus, wherein the standard form is a standard brain.

4. 2. The magnetic resonance imaging apparatus according to claim 1, wherein the morphological image is a T1 weighted image.

5. 4. The magnetic resonance imaging apparatus according to claim 3, wherein the imaging sequence is a sequence for implementing a Multi-PLD ASL (arterial spin labeling) method in which a plurality of types of time from blood labeling to imaging are set and imaging is performed for each of the plurality of types of time.

6. 6. The magnetic resonance imaging apparatus according to claim 5, wherein the functional image is an image of arterial transit time (ATT) and / or an image of cerebral blood flow (CBF).

7. 2. A magnetic resonance imaging apparatus as described in claim 1, characterized in that the calculation processing unit calculates the degree of deviation of the functional image of the subject after deformation from a healthy functional image obtained by deforming a functional image previously obtained for a healthy person so as to match the position of the corresponding area of ​​a standard morphology.

8. 8. The magnetic resonance imaging apparatus according to claim 7, wherein the calculation processing unit sets a region of interest in a functional image of the subject after deformation, calculates a degree of discrepancy between the functional image within the region of interest and a corresponding region of the healthy functional image, and displays the degree of discrepancy to a user.

9. 2. The magnetic resonance imaging apparatus according to claim 1, wherein the arithmetic processing unit sets a region of interest in the transformed functional image, calculates a predetermined feature amount for the functional image within the region of interest, and calculates on which side of a boundary between a distribution of the feature amounts of a plurality of healthy individuals and a distribution of the feature amounts of a plurality of individuals with a disease, which have been obtained in advance, the feature amount is located, as well as the distance from the boundary.

10. 2. The magnetic resonance imaging apparatus according to claim 1, wherein the imaging sequence is a sequence for realizing a Single PLD ASL (arterial spin labeling) method in which imaging is performed by setting a single time period from blood labeling to imaging, A magnetic resonance imaging apparatus, wherein the functional image is an image of cerebral blood flow (CBF).

11. 2. The magnetic resonance imaging apparatus according to claim 1, wherein the functional images include a cerebral blood flow image and a quantitative magnetic susceptibility mapping image, and the morphological images are T1-weighted images; The magnetic resonance imaging apparatus is characterized in that the calculation processing unit calculates the degree of deviation of the cerebral blood flow image, the quantitative susceptibility mapping image, and an image within the corresponding region of the T1-weighted image of the subject from the cerebral blood flow image, the quantitative susceptibility mapping image, and an image within the corresponding region of the T1-weighted image of a healthy subject.

12. a processor that receives a morphological image showing a morphology and a functional image showing a function, both of which are captured in the same imaging region of a subject; The image processing device is characterized in that the calculation processing unit performs a process of moving the position of one or more structures in the morphological image to the position of a structure in a predetermined standard form using deformation parameters, and then deforms the functional image using the values ​​of the deformation parameters to match the position of an area in the functional image to the position of a corresponding area in the standard form, or deforms the standard form in the opposite direction using the values ​​of the deformation parameters to match the position of the area of ​​the structure in the standard form to the position of the corresponding area in the functional image.

13. receiving a morphological image showing a morphology and a functional image showing a function, both of which are captured in the same imaging region of the subject; moving the positions of one or more structures in the morphological image to positions of structures of a predetermined standard shape using the deformation parameters; The functional image is deformed using the values ​​of the deformation parameters so that the position of the region in the functional image coincides with the position of the corresponding region in the standard form, or the standard form is deformed in the opposite direction using the values ​​of the deformation parameters so that the position of the region of the structure in the standard form coincides with the position of the corresponding region in the functional image.

13. An image processing method comprising:

Images