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

Abstract

To obtain a volume rendering image from a 3D MRI image by automatically performing appropriate opacity setting.SOLUTION: A three-dimensional image of a subject is received to obtain a distribution of pixel values of the three-dimensional image, a pixel value of a feature amount predetermined according to the distribution of pixel values is calculated, and opacity is set for each pixel value included in the three-dimensional image on the basis of the pixel value of the feature amount. Thus, opacity can be set automatically. Opacity is used to generate a volume rendering image of the three-dimensional image.SELECTED DRAWING: Figure 6

Description

[Technical field]

[0001] The present invention relates to a technique for creating a volume rendering image of a desired tissue from a three-dimensional (3D) image captured by a magnetic resonance imaging (hereinafter, referred to as "MRI") device or the like. [Background technology]

[0002] An MRI device is a device that measures nuclear magnetic resonance (NMR) signals generated by atomic spins that constitute 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. In 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] For image diagnosis of cerebral vascular lesions such as cerebral aneurysms, a method is used in which images of cerebral blood vessels with high signal intensity are obtained without using contrast agents by imaging using a 3D TOF (time of flight) sequence. By projecting the obtained 3D TOF image at a desired angle using the MIP (maximum intensity projection) method, blood vessel images can be viewed from a desired angle. Note that, before projection using the MIP method, signals from unnecessary tissues other than the blood vessels to be observed, such as subcutaneous fat and superficial blood vessels, are removed (clipping process), and further, opacity is set to make the pixels of the brain parenchyma transparent so that they are not displayed, thereby obtaining an image in which the morphology of the cerebral blood vessels can be observed more clearly. However, the clipping process and opacity setting are performed manually, which is a cumbersome and time-consuming task.

[0004] Therefore, in order to perform the above clipping process automatically, Patent Document 1 discloses a technique in which a 3D MRI image is processed to generate a brain extraction mask image for extracting the brain region and a blood vessel mask image for extracting the blood vessel region, and an image of the brain region and blood vessels in the brain is extracted from the 3D image using a mask that combines both mask images. With this method, blood vessels in the brain region and blood vessels at the base of the skull, which are necessary for observing cerebral blood vessels, can be extracted (clipped) from one 3D MRI image.

[0005] On the other hand, Patent Document 2 discloses a technique for identifying the brain contour in an X-ray CT image, acquiring the MR signal intensity of pixels located on the brain contour in the corresponding MRI image, and setting an opacity curve based on this signal intensity so that the pixels on the brain contour are displayed opaquely and the other areas are displayed transparently. This opacity curve indicates the relationship between the MR signal value and the opacity of the display pixel value, and specifically, the range of MRI values ​​indicated by the brain contour pixels is set to an opacity of 1, and the other MRI value ranges are set to an opacity of 0. This makes it possible to display only the brain contour as a pseudo 3D image. [Prior art documents] [Patent documents]

[0006] [Patent Document 1] JP 2020-39507 A [Patent Document 2] JP 2012-100955 A Summary of the Invention [Problem to be solved by the invention]

[0007] The technology of Patent Document 1 is a technology for automatically performing clipping processing, and when creating a volume rendering image, opacity setting must be performed manually.

[0008] On the other hand, the method of Patent Document 2 requires two types of images, an X-ray CT image and an MRI image. In addition, the method of Patent Document 2 discloses an opacity setting method for making only the MRI signal indicating the region of the brain contour opaque, but does not disclose any method for setting opacity so as to make only the contour of the blood vessels in the brain opaque. Since the blood vessels in the brain are more complicated and intricate than the contour shape of the brain, it is presumed that it is not easy to set opacity using the technology of Patent Document 2.

[0009] Furthermore, since the method of Patent Document 2 does not perform clipping processing, if the method of Patent Document 2 is applied to setting the opacity of cerebral blood vessels, signals from unnecessary tissues other than the target blood vessels, such as subcutaneous fat and superficial blood vessels, will remain in the MRI image, requiring cumbersome clipping.

[0010] In addition, assuming that the clipping process of Patent Document 1 and the opacity process of Patent Document 2 are combined, the method of Patent Document 2 sets the opacity based on the histogram of MR signal intensity of the contour pixels of the region to be extracted, whereas the technology of Patent Document 1 uses an integrated mask to extract the brain and cerebral blood vessel images from the 3D image and performs clipping to remove unnecessary tissues, so the histogram of MR signal intensity changes from the histogram before clipping. Therefore, there is a possibility that an appropriate opacity for evaluating blood vessels cannot be set from the histogram of the changed MRI signal intensity.

[0011] The objective of this application is to obtain a volume rendering image from a 3D MRI image by automatically setting appropriate opacity. [Means for solving the problem]

[0012] In order to achieve the above object, the magnetic resonance imaging apparatus of the present invention includes a static magnetic field generating unit that generates a static magnetic field in a space in which an imaging region of a subject is disposed, a measurement control unit that applies a gradient magnetic field and a high frequency magnetic field to the subject and detects a nuclear magnetic resonance signal generated from the imaging region, and a calculation processing unit that reconstructs a three-dimensional image of the subject using the detected nuclear magnetic resonance signal. The calculation processing unit obtains a distribution of pixel values ​​of the three-dimensional image, calculates pixel values ​​of a predetermined feature amount from the distribution of pixel values, and sets an opacity for each pixel value included in the three-dimensional image based on the pixel values ​​of the feature amount. Effect of the Invention

[0013] According to the present invention, since an appropriate opacity can be automatically set from a 3D MRI image, an MRI volume rendering image can be automatically created. This reduces the amount of complicated work for the operator and improves efficiency. [Brief description of the drawings]

[0014] [Figure 1] 1 is a block diagram showing an example of the overall configuration of an MRI apparatus according to the present invention. [Diagram 2] An example of a 3D TOF imaging pulse sequence. [Diagram 3] 4 is a flowchart showing the processing of a calculation processing unit according to the first embodiment of the present invention. [Figure 4] 5 is a flowchart showing an opacity setting process of the arithmetic processing unit according to the first embodiment of the present invention. [Diagram 5] 1A to 1C are diagrams showing images before, during, and after processing according to the first embodiment of the present invention. [Figure 6] 4 is a graph showing a histogram and an opacity curve generated for opacity setting by the calculation processing unit according to the first embodiment of the present invention. [Figure 7] 10 is a flowchart showing the processing of a calculation processing unit according to the second embodiment of the present invention. [Figure 8] 11 is a flowchart showing the processing of a calculation processing unit according to a third embodiment of the present invention. [Figure 9]10 is a flowchart showing the processing of a calculation processing unit according to a fourth embodiment of the present invention. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0015] 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.

[0016] The MRI apparatus of this embodiment captures 3D TOF images, automatically performs clipping processing and transparency setting, and creates an optimal volume rendering image.

[0017] First, an overall outline of an example of an MRI apparatus according to the present invention will be described with reference to Fig. 1. Fig. 1 is a block diagram showing the overall configuration of an embodiment of an MRI apparatus according to the present invention. This MRI apparatus obtains a tomographic image of a subject by utilizing the NMR phenomenon, and as shown in Fig. 1, the MRI apparatus is configured with a static magnetic field generating unit 2, a gradient magnetic field generating unit 3, a transmitting unit 5, a receiving unit 6, a signal processing unit 7, a sequencer 4, and a calculation processing unit (CPU) 8. The gradient magnetic field generating unit 3, the transmitting unit 5, the receiving unit 6, and the sequencer 4 are collectively referred to as a measurement control unit 100.

[0018] The static magnetic field generating unit 2 generates a static magnetic field in a space in which the imaging region of the subject 1 is located. The static magnetic field generating unit 2 is a static magnetic field generating source of a permanent magnet type, a resistive conductive type, or a superconductive type, which is arranged around the subject 1. The static magnetic field generating source generates a uniform static magnetic field in a direction perpendicular to the body axis of the subject 1 if it is a vertical magnetic field type, or in the direction of the body axis of the subject 1 if it is a horizontal magnetic field type.

[0019] The gradient magnetic field generating unit 3 is configured with a gradient magnetic field coil 9 that applies gradient magnetic fields in three axial directions of X, Y, and Z, which are the coordinate part (static coordinate part) of the MRI apparatus, and a gradient magnetic field power supply 10 that drives each gradient magnetic field coil, and applies gradient magnetic fields Gx, Gy, and Gz in three axial directions of X, Y, and Z by driving the gradient magnetic field power supply 10 of each coil according to an instruction from a sequencer 4 described later. 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.

[0020] The transmitter 5 irradiates the subject 1 with a radio frequency magnetic field pulse (hereinafter referred to as "RF pulse") in order to cause nuclear magnetic resonance in the nuclear spins of atoms constituting the biological tissue of the subject 1, and is configured with a radio frequency oscillator 11, a modulator 12, a radio frequency amplifier 13, and a radio frequency coil (transmitting coil) 14a on the transmitting side. The RF pulse output from the radio frequency oscillator 11 is amplitude-modulated by the modulator 12 at a timing according to a command from the sequencer 4 described later, and this amplitude-modulated RF pulse is amplified by the radio frequency amplifier 13 and then supplied to the radio frequency coil 14a arranged close to the subject 1. As a result, the RF pulse is irradiated from the radio frequency coil 14a to the subject 1.

[0021] The receiving unit 6 detects echo signals (hereinafter referred to as NMR signals) emitted by nuclear magnetic resonance (NMR) of atomic nuclear spins constituting the biological tissue of the subject 1, and is configured with a receiving side high-frequency coil (receiving coil) 14b, a signal amplifier 15, a quadrature phase detector 16, and an A / D converter 17.

[0022] The NMR signal of the response of the subject 1 induced by the electromagnetic waves irradiated from the transmitting coil 14a is detected by the receiving coil 14b arranged close to the subject 1, amplified by the signal amplifier 15, and then split into two orthogonal signals by the quadrature phase detector 16 at a timing specified by the 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 unit 7.

[0023] The sequencer 4 is a control means for repeatedly applying RF pulses and gradient magnetic field pulses according to a certain pulse sequence, and receiving the NMR signals generated thereby at a certain timing. The sequencer 4 operates under the control of the calculation processing unit 8, and sends various commands to the transmitting unit 5, the gradient magnetic field generating unit 3, and the receiving unit 6 according to the pulse sequence, and collects NMR signal data required for generating a tomographic image of the subject 1.

[0024] The signal processing unit 7 performs various data processing and displays and stores the processing results, and has an external storage device such as an optical disk 19 or a magnetic disk 18, and a display 20 such as a CRT.

[0025] When the NMR signal data is input from the receiving unit 6 to a processor (CPU) 8, the processor 8 executes signal processing, image reconstruction, and other processes to reconstruct a tomographic image of the subject 1. The processor 8 displays the tomographic image of the subject 1 on a display 20, and records it on an external storage device such as a magnetic disk 18.

[0026] The operation unit 25 is used to input various control information for the MRI apparatus and control information for processing performed by the signal processing unit 7, and is composed 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.

[0027] 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 unit 2 into which the subject 1 is inserted, facing the subject 1 in the case of a vertical magnetic field type, or surrounding the subject 1 in the case of a horizontal magnetic field type. On the other hand, the high frequency coil 14b on the receiving side is installed facing or surrounding the subject 1. As described above, the gradient magnetic field generating unit 3, the transmitting unit 5, the receiving unit 6, and the sequencer 4 are collectively referred to as the measurement control unit 100.

[0028] 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.

[0029] As a method for imaging cerebral blood vessels using MRI, there is a pulse sequence of the 3D Time of Flight (TOF) method. FIG. 2 is an example of a sequence diagram of the 3D TOF method. In FIG. 2, RF indicates the timing of application of an RF pulse 401 and acquisition of an NMR signal 402, and Gs, Gp, and Gr indicate the application timing and magnitude of gradient magnetic field pulses 403 to 408 in the slice direction, phase encoding direction, and readout direction. The number after the hyphen after the three-digit number indicates the number of repetitions. For example, the "1" in 402 "-1" indicates the pulse applied and the signal acquired in the first repetition, and the "2" in "-2" indicates the pulse applied and the signal acquired in the second repetition. 2, a sequence of irradiating an RF pulse 401 while applying a gradient magnetic field pulse 403 in a slice direction Gs, applying a gradient magnetic field pulse 404 in a slice direction Gs, a gradient magnetic field pulse 405 in a phase encoding direction Gp, and a gradient magnetic field pulse 408 in a readout direction Gr, and then receiving an NMR signal 402 while applying a gradient magnetic field pulse 407 in the readout direction is repeated a predetermined number of times. Each time the sequence is repeated, the magnitude of the gradient magnetic field pulse 404 in the slice direction Gs is changed, and the magnitudes of the gradient magnetic field pulses 405 and 406 in the phase encoding direction Gp are changed, thereby changing the slice position and phase encoding to obtain an NMR signal required for generating a 3D TOF image.

[0030] In the MRI apparatus of this embodiment, the arithmetic processing unit 8 automatically performs clipping processing and opacity setting from the 3D TOF image to create a volume rendering image will be described below. The processing performed by the arithmetic processing unit 8 in each embodiment is image processing, and the arithmetic processing unit 8 functions as an image processing device.

[0031] <<Embodiment 1>> The clipping process and the opacity setting process of the MRI apparatus of the embodiment 1 will be described in detail with reference to the flowcharts of Fig. 3 and Fig. 4. Also, Fig. 5(a) to (c) show images during processing, and Fig. 6 shows a histogram of pixel values.

[0032] In the first embodiment, the arithmetic processing unit 8 obtains a distribution of pixel values ​​of a three-dimensional image (3D TOF image), calculates a pixel value of a predetermined feature from the distribution of pixel values, and sets an opacity for each pixel value included in the three-dimensional image based on the pixel value of the feature. When the imaging site is the brain, the arithmetic processing unit 8 sets an opacity for each pixel value such that an opacity of 0 is set for a pixel value range corresponding to an image of the brain parenchyma in the three-dimensional image, and an opacity of 1 is set for a pixel value range corresponding to an image of blood vessels.

[0033] Specifically, the pixel value of the feature is the average value of the distribution of pixel values ​​of the three-dimensional image. The calculation processing unit 8 sets an opacity curve in which the opacity changes from 0 to 1 as the pixel value increases, with a pixel value that is a predetermined value greater than the pixel value (average value) of the feature as the boundary. For example, the predetermined value is set based on the variance (σ) calculated from the distribution of pixel values.

[0034] More specifically, the calculation processing unit 8 sets an opacity curve in which the opacity changes linearly or nonlinearly from 0 to 1, with a pixel value that is a predetermined value (3σ) greater than the pixel value of the feature amount as the boundary.

[0035] The processing performed by the arithmetic processing unit 8 will now be described in further detail.

[0036] The processing flow of the arithmetic processing unit 8 shown in Figures 3 and 4 is stored in advance as a program on the magnetic disk 18, and is realized by software when the arithmetic processing unit (CPU) 8 reads the program from the magnetic disk 1 and executes it. 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.

[0037] (Step S201) The operator sets up the imaging region (here, the brain) of the subject 1 in the static magnetic field space generated by the static magnetic field generating unit 2. The calculation processing unit 8 receives, via the operation unit 25, an input from the operator of settings of imaging conditions for a sequence for capturing a 3D TOF image.

[0038] (Step S202) The arithmetic processing unit 8 instructs the sequencer 4 to execute the 3D TOF pulse sequence of FIG. 2 under the imaging conditions set in step S201. As a result, imaging of the subject 1 is performed using the 3D TOF pulse sequence. The arithmetic processing unit 8 processes the acquired NMR signals to reconstruct a 3D TOF image. The generated image data is stored in the magnetic disk 18. Here, an image of the brain is captured as an example.

[0039] (Step S203) The arithmetic processing unit 8 receives a selection of a 3D TOF image to be processed from among the images stored in the magnetic disk 18 from the operator via the operation unit 25 (FIG. 5(a)). The arithmetic processing unit 8 starts a process of creating a volume rendering image for the selected 3D TOF image. Note that, without performing this step S203, if image capture is completed in step S202, the arithmetic processing unit 8 may select the captured image as the processing target.

[0040] (Step S204) The processor 8 automatically performs clipping processing on the 3D TOF image acquired in step S202 using a known method (FIG. 5(b)). In this embodiment, for example, the method of Patent Document 1 is used to remove subcutaneous fat regions, skull regions, and the like, and an image including the brain parenchyma region including blood vessels and the blood vessels at the base of the skull is obtained.

[0041] (Step S205) The calculation processing unit 8 calculates the opacity of the 3D TOF image after the clipping process is applied, which is acquired in step S204. The process of step S205 will be described in detail with reference to the flow of FIG.

[0042] (Step S41) The calculation processing unit 8 normalizes the 3D TOF image after clipping processing has been applied. For example, normalization is performed so that the signal values ​​(pixel values) fall within the range of 0 to 255.

[0043] (Step S42) The calculation processing unit 8 uses a discriminant analysis method to calculate a threshold value 61 (see FIG. 6 ) for removing background regions so as to maximize the degree of separation of pixel values ​​of the normalized 3D TOF image, and replaces pixel values ​​below the threshold value 61 with 0 to remove the background regions of the 3D TOF image.

[0044] (Step S43) The arithmetic processing unit 8 generates a histogram (distribution of pixel values) of the pixel values ​​and pixel counts of the 3D TOF image after the background region has been removed. An example of the histogram is shown in FIG. 6. The 3D TOF image after the background region has been removed is composed almost entirely of images of the brain parenchyma region and the blood vessel region. The pixels of the blood vessel region have a large signal value but a small number of pixels, while the pixels of the brain parenchyma region have a smaller signal value than the blood vessel region but a much larger number of pixels than the blood vessel region. Therefore, as shown in FIG. 6, the histogram of pixel values ​​of the 3D TOF image after the background region has been removed shows that the pixels of the brain parenchyma have a peaked distribution, and the pixels of the blood vessel region are distributed in a region with a larger pixel value than the distribution of the pixels of the brain parenchyma.

[0045] (Step S44) In order to set opacity 0 (transparent) for the pixel value range of the brain parenchyma and opacity 1 (opaque) for the pixels of the blood vessels, the calculation processing unit 8 calculates the average value (Mean) and variance (σ) of the pixel values ​​(signal values) as the feature amount of the histogram, and sets the opacity 0 (transparent) for the pixel value range of the brain parenchyma from the minimum pixel value 0 to the pixel value of Mean (average value) + 3σ (= pixel value C1), and sets the opacity to a value greater than 0 for the pixel value range of the blood vessels. Here, as shown in FIG. 6, an opacity curve A is set in which the pixel value range of Mean (average value) + 3σ (= pixel value C1) or less is opacity 0, the opacity changes linearly or nonlinearly from 0 to 1 in the range of Mean + 3σ (= pixel value C1) to Mean + 5σ (= pixel value C2), and the opacity is 1 (opaque) for Mean + 5σ (= pixel value C2) or more. FIG. 6 shows, as an example, an opacity curve A in which the opacity value increases linearly (proportional to the pixel value) in the range from Mean+3σ (=pixel value C1) to Mean+5σ (=pixel value C2).

[0046] (Step S206) The arithmetic processing unit 8 creates a volume rendering image using the 3D TOF image to which clipping has been applied, obtained in step S204. Here, the value of the opacity curve set in step S205 is used when creating the volume rendering image. That is, for each pixel of the 3D TOF image after clipping obtained in step S204, if the pixel value is between 0 and Mean+3σ, opacity 0 is set, and if the pixel value is Mean (average value)+5σ or more, opacity 1 is set, and the transparency of the linear or nonlinear opacity curve A is set for the pixels between them, and volume rendering is performed. As a result, as shown in FIG. 5(c), a volume rendering image of the blood vessel image can be obtained by making the pixels of the brain parenchyma transparent.

[0047] (Step S207) The calculation processing unit 8 displays the volume rendering image, for which the clipping and opacity calculated in step S206 have been set, on the display 20, and stores the image.

[0048] As described above, in the first embodiment, the calculation processing unit 8 can automatically perform clipping processing and automatic opacity setting from 3D TOF images, and can automate the creation of MRI volume rendering images. This reduces the amount of work for the operator and improves efficiency.

[0049] In the above embodiment, the opacity is set in step 205 as shown in Figure 6 by setting an opacity curve A that changes linearly from opacity 0 to 1 between Mean (average value) + 3σ and Mean (average value) + 5σ, but it may also be a sigmoid curve like curve B in Figure 6.

[0050] It is also possible to set an opacity curve in which the opacity value changes stepwise from 0 to 1, with a pixel value that is a predetermined value greater than the pixel value of the feature (for example, Mean (average value) + 3σ) as the boundary.

[0051] In this embodiment, the average value (Mean) of the pixel values ​​of the 3D image after clipping is set as the feature value, and the pixel values ​​C1 and C2 of the boundary for changing the opacity are set based on this feature value, but the feature value does not necessarily have to be the average value. For example, the peak value of the histogram in Fig. 6 may be set as the feature value, and the pixel values ​​C1 and C2 may be pixel values ​​that are a predetermined value away from the peak value.

[0052] <<Embodiment 2>> The clipping process and the opacity setting process of the MRI apparatus of the second embodiment will be described in detail with reference to the flowchart of FIG.

[0053] The second embodiment is similar to the first embodiment in that clipping processing and opacity setting are automatically performed on 3D TOF images, and the creation of MRI volume rendering images is automated. Furthermore, in the second embodiment, step S301 is performed between steps S205 and S206 in the flow of FIG. 3 of the first embodiment, so that pixel values ​​C1 and C2, which are the boundaries for changing the opacity from 0 to 1, are searched for to be more optimal pixel values.

[0054] That is, the calculation processing unit 8 sets multiple opacity curves while shifting the boundary pixel values ​​C1 and C2, performs calculations on the distribution of pixel values ​​of a three-dimensional image using the multiple opacity curves, and selects one opacity curve based on the calculation results.

[0055] A specific description will be given below with reference to the flow in Fig. 7. In Fig. 7, the same steps as those in the first embodiment are given the same step reference numerals, and the description will be omitted.

[0056] (Steps S201 to S205) Since this is the same as in the first embodiment, the description will be omitted. In step S205, based on the histogram of pixel values ​​of the 3D TOF image, an opacity curve A is set in which the opacity changes linearly from 0 to 1 between the average pixel value Mean+3σ (=pixel value C1) and the average pixel value Mean+5σ (=pixel value C2).

[0057] (Step S301) The arithmetic processing unit 8 further exploratively determines pixel values ​​C1 and C2 that change the opacity from 0 to 1. For example, the pixel values ​​C1 and C2 determined in step S205 are set as initial values, and the pixel values ​​C1 and C2 (C1=Mean+3σ, C2=Mean+5σ) that change the opacity from 0 to 1 are changed (shifted) to set an opacity curve A in which the opacity changes, for example, linearly from 0 to 1. For example, to search for C1, step S206 is performed while changing the value of C1 to obtain a volume rendering image as shown in FIG. 5(c), the volume of the area remaining in the volume rendering image is calculated, and the amount of change is used as a feature to perform the search. Specifically, the pixel value C1 is set as an initial value, and the amount of change in the volume remaining in the obtained volume rendering image is obtained by performing step S206 each time while changing the pixel value, and if the amount of change is smaller than a predetermined threshold, the pixel value at that time is determined as the value of C1 of the opacity curve A. This makes it possible to determine a more appropriate pixel value C1 in which no image of the brain parenchyma remains in the volume rendering image, and only blood vessels remain.

[0058] Similarly, while changing the pixel value C2 in the same manner, step S206 is performed to obtain a volume rendering image, and the volume of the area remaining in the volume rendering image is calculated.If the amount of change is smaller than a predetermined threshold, the pixel value at that time is determined as the value of C2 of the opacity curve A.

[0059] This makes it possible to set appropriate pixel values ​​C1 and C2 that allow only blood vessels to remain in the volume rendering image.

[0060] In addition, in step S301, instead of the volume of the region remaining in the volume rendering image, it is also possible to use the area of ​​a specified cross section of the volume rendering image or the area of ​​a two-dimensional image obtained by projecting the volume rendering image onto a specified plane, and search for pixel values ​​C1 and C2 using the amount of change as a feature.

[0061] Furthermore, the search method is not limited to these methods, and may involve multiplying the set opacity curve A and a histogram, and searching for pixel values ​​C1 and C2 where the area of ​​the histogram curve after the multiplication falls within a predetermined area range of only the major blood vessels.

[0062] (Steps S206~S207) Volume rendering is created using the pixel values ​​calculated in step S301 in the same manner as in the first embodiment.

[0063] As described above, in the second embodiment of the present invention, the calculation processing unit 8 automatically performs clipping processing and searches for pixel values ​​for opacity setting from the 3D TOF image, sets a more optimal opacity, and automates the creation of an MRI volume rendering image. This makes it possible to obtain a volume rendering image with a more optimal opacity setting.

[0064] In the second embodiment, instead of the linear opacity curve A, the curve B such as the sigmoid curve described in the first embodiment may be used.

[0065] <<Embodiment 3>> The processing of 3D TOF images by the MRI apparatus of the third embodiment will be described in detail with reference to the flowchart of FIG.

[0066] In the third embodiment, similar to the first embodiment, clipping processing and automatic opacity setting are performed on 3D TOF images to automate the creation of MRI volume rendering images. However, by performing step S401 between steps S206 and S207 in the flow of FIG. 3 in the first embodiment, regions unnecessary for evaluation that remain even after volume rendering are deleted.

[0067] That is, the calculation processing unit 8 performs a process of extracting an unnecessary region by performing a predetermined image processing on the volume rendering image, and then performs a process of deleting the extracted unnecessary region. For example, the predetermined image processing is a process of extracting a discontinuous region.

[0068] A specific description will be given with reference to the flow in Fig. 8. In Fig. 8, the same steps as those in the first embodiment are given the same step reference numerals, and the description will be omitted.

[0069] (Steps S201 to S206) Since this is the same as in the first embodiment, a description will be omitted, but after clipping and opacity setting are automatically performed on the 3D TOF image, a volume rendering image is generated.

[0070] (Step S401) The calculation processing unit 8 performs a process of deleting regions that are unnecessary for diagnosis, such as subcutaneous fat regions, remaining in the volume rendering image. For example, the calculation processing unit 8 labels all remaining regions, extracts discontinuous regions, and regards the extracted regions as regions other than blood vessels and deletes them.

[0071] Alternatively, a method may be used in which a feature amount for each region of a volume rendering image is calculated, and regions other than blood vessels are identified and removed based on the feature amount.

[0072] Alternatively, a method may be used in which a learning model that has been trained to distinguish between blood vessels and unnecessary regions through machine learning or the like using a large number of volume rendering images is used to remove regions other than blood vessels based on image recognition.

[0073] (Step S207) The volume rendering image after the unnecessary regions have been deleted is displayed and saved in the same manner as in the first embodiment.

[0074] As described above, in the third embodiment of the present invention, the calculation processing unit 8 automatically performs clipping processing and opacity setting on the 3D TOF image, and further automatically deletes the remaining areas unnecessary for diagnosis, thereby making it possible to create a volume rendering image more suitable for diagnosis without the operator having to perform cumbersome operations.

[0075] <<Embodiment 4>> The processing of 3D TOF images by the MRI apparatus of the fourth embodiment will be described in detail with reference to the flowchart of FIG.

[0076] In the fourth embodiment, similar to the first embodiment, clipping processing and automatic opacity setting are automatically performed on 3D TOF images to automate the creation of MRI volume rendering images, but step S501 is performed between steps S206 and S207 in the flow of the first embodiment in Fig. 3 to further extract a predetermined specific blood vessel image from the blood vessel image in the volume rendering image. The specific blood vessel image is a predetermined main blood vessel or a blood vessel on the left or right side of the subject.

[0077] A specific description will be given with reference to the flow in Fig. 9. In Fig. 9, the same steps as those in the first embodiment are given the same step reference numerals, and the description will be omitted.

[0078] (Steps S201 to S206) Since this is the same as in the first embodiment, a description will be omitted, but after clipping and opacity setting are automatically performed on the 3D TOF image, a volume rendering image is generated.

[0079] (Step S501) The arithmetic processing unit 8 creates a volume rendering image by extracting an image of a predetermined major blood vessel and / or a blood vessel on the left side and / or the right side of the imaging site (here, the brain) of the subject 1 from the volume rendering image generated in step 206. Specifically, the initial points of the major blood vessels or the left and right blood vessels are extracted from the volume rendering image, and the blood vessels are extracted by tracing the blood vessel image from the initial points using the continuity of the blood vessels.

[0080] For example, in the case of left and right separately, the approximate left and right center positions of the body (brain) are obtained from the image position of the DICOM information attached to the MRI image and the patient position, the left and right regions are set as volume rendering images, and the pixels with the highest signal values ​​on the blood vessels in the left and right regions are obtained, respectively, and these are set as the initial points. From these initial points, consecutive pixels on the image are traced using a region growing method or the like to extract consecutive blood vessels in the left region and consecutive blood vessels in the right region separately.

[0081] In addition, when extracting major blood vessels, first, the initial point is set to the pixel with the largest pixel value, and the blood vessel is extracted from the initial point by a region growing method or the like using the continuity of the blood vessel, and then the initial point is set to the pixel with the largest pixel value in the remaining image, and the second blood vessel is extracted. By repeating this process, multiple blood vessels are extracted in order. After labeling the multiple extracted blood vessels, the largest continuous area is regarded as a major blood vessel such as the middle cerebral artery, and the next largest continuous area is regarded as a blood vessel such as the vertebral artery, using the position information and continuity, and each is extracted. These blood vessels may be extracted by labeling by deep learning using a trained learning model, or the like.

[0082] (Step S207) One or more volume rendering images from which desired blood vessels have been extracted are displayed and stored in the same manner as in the first embodiment.

[0083] As described above, in embodiment 4 of the present invention, the calculation processing unit 8 automatically performs clipping processing and opacity setting from 3D TOF images, and further automates the creation of volume rendering images for the left and right and for major blood vessels, thereby making it possible to create volume rendering images more suitable for diagnosis without the operator having to perform complicated tasks.

[0084] Although the embodiments of the present invention have been described above, it goes without saying that the present invention is not limited to these. [Explanation of symbols]

[0085] 1: subject, 2: static magnetic field generator, 3: gradient magnetic field generator, 4: sequencer, 5: transmitter, 6: receiver, 7: signal processor, 8: arithmetic processor (image processor, 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 (transmitting coil), 14b: high frequency coil (receiving 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 generates a static magnetic field in a space in which an imaging region of a subject is disposed; a measurement control unit that applies a gradient magnetic field and a high frequency magnetic field to the subject and detects a nuclear magnetic resonance signal generated from the imaging region; and a calculation processing unit that generates a three-dimensional image of the subject using the detected nuclear magnetic resonance signal, The arithmetic processing unit is generating a 3D TOF image (three-dimensional time of flight image) of a brain parenchyma region including blood vessels as the three-dimensional image; determining a histogram of pixel values ​​of the 3D TOF image; removing pixels of a background region of the 3D TOF image having a pixel value equal to or less than a predetermined threshold value from the histogram to obtain a histogram in which the pixels of the brain parenchyma show a distribution having a peak and the pixels of the blood vessels are distributed in a region having a pixel value larger than that of the distribution having the peak of the pixels of the brain parenchyma; Calculate the average value (Mean) and variance (σ) of the pixel values ​​of the histogram, and set an opacity value of zero as the pixel value range of the brain parenchyma, the pixel value range being equal to or less than a pixel value C1 (=Mean+3σ), which is the sum of the average value (Mean) and three times the variance (σ); A magnetic resonance imaging apparatus characterized in that, in order to set the opacity to a value greater than 0, the range of pixel values ​​greater than pixel value C1 (= Mean + 3σ) is defined as the range of pixel values ​​of the blood vessels, the opacity value is changed from zero to 1 between the pixel value C1 (= Mean + 3σ) and a pixel value C2 (= Mean + 5σ), which is the sum of the average value (Mean) and 5 times the variance (σ), and an opacity curve is set in which the opacity value is 1 for pixel values ​​greater than or equal to pixel value C2.

2. 2. The magnetic resonance imaging apparatus according to claim 1, wherein the calculation processing unit generates a volume rendering image of the blood vessels by volume rendering the three-dimensional image using the opacity value.

3. A static magnetic field generating unit that generates a static magnetic field in a space in which an imaging portion of a subject is disposed, a measurement control unit that applies a gradient magnetic field and a high frequency magnetic field to the subject and detects a nuclear magnetic resonance signal generated from the imaging portion, and a calculation processing unit that generates a three-dimensional image of the subject using the detected nuclear magnetic resonance signal; The arithmetic processing unit is generating a 3D TOF image (three-dimensional time of flight image) of a brain parenchyma region including blood vessels as the three-dimensional image; A histogram of pixel values ​​of the 3D TOF image is obtained, a mean value (Mean) and a variance (σ) of the pixel values ​​of the histogram are calculated, and an opacity value of zero is set in a range of pixel values ​​equal to or less than a pixel value C1 (=Mean+3σ) which is the sum of the mean value (Mean) and three times the variance (σ); an opacity curve is set in which the opacity value is changed from zero to one between the pixel value C1 (= Mean + 3σ) and a pixel value C2 (= Mean + 5σ) which is the sum of the average value (Mean) and five times the variance (σ), and the opacity value is one for pixel values ​​equal to or greater than the pixel value C2; The magnetic resonance imaging apparatus is characterized in that the calculation processing unit generates a volume rendering image of the three-dimensional image using the opacity of the opacity curve while shifting at least one of the pixel values ​​C1 and C2 of the opacity curve, calculates the volume of the image remaining in the volume rendering image and / or the area of ​​a two-dimensional image generated from the volume rendering image, and determines the values ​​of the pixel values ​​C1 and / or C2 to be used in the opacity curve based on the volume and / or area.

4. 3. The magnetic resonance imaging apparatus according to claim 2, wherein the arithmetic processing unit extracts unnecessary regions by performing predetermined image processing on the volume rendering image, and performs processing to delete the extracted unnecessary regions.

5. a static magnetic field generating unit that generates a static magnetic field in a space in which an imaging region of a subject is disposed; a measurement control unit that applies a gradient magnetic field and a high frequency magnetic field to the subject and detects a nuclear magnetic resonance signal generated from the imaging region; and a calculation processing unit that reconstructs a three-dimensional image of the subject using the detected nuclear magnetic resonance signal, The arithmetic processing unit is determining a distribution of pixel values ​​of the three-dimensional image, calculating a pixel value of a predetermined feature amount from the distribution of pixel values, and setting an opacity curve in which an opacity changes from 0 to 1 as the pixel value increases, with a pixel value that is a predetermined value larger than the pixel value of the feature amount as a boundary; setting a plurality of opacity curves while shifting the pixel value that is the boundary, performing a calculation on the distribution of pixel values ​​of the three-dimensional image using the plurality of opacity curves, and selecting one opacity curve based on the calculation result; volume-rendering the three-dimensional image using the selected opacity curve, extracting an unnecessary region by performing predetermined image processing on the obtained volume-rendered image, and deleting the extracted unnecessary region; 4. A magnetic resonance imaging apparatus, comprising: a first image processing unit that performs image processing for extracting discontinuous regions;

6. 3. The magnetic resonance imaging apparatus according to claim 2, wherein the arithmetic processing unit further performs a process of extracting an image of a predetermined specific blood vessel from the volume rendering image of the blood vessel.

7. A processor receives a 3D TOF image (three-dimensional time of flight image) of the subject's head and processes it; The arithmetic processing unit is generating a 3D TOF image of a brain parenchyma region including blood vessels from the 3D TOF image of the head; A histogram of pixel values ​​of a 3D TOF image of the brain parenchyma region including the blood vessels is obtained; removing pixels of a background region of the 3D TOF image having a pixel value equal to or less than a predetermined threshold value from the histogram to obtain a histogram in which the pixels of the brain parenchyma show a distribution having a peak and the pixels of the blood vessels are distributed in a region having a pixel value larger than that of the distribution having the peak of the pixels of the brain parenchyma; Calculate the average value (Mean) and variance (σ) of the pixel values ​​of the histogram, and set an opacity value of zero by assuming that the pixel value range of the brain parenchyma is equal to or less than a pixel value C1 (=Mean+3σ), which is the sum of the average value (Mean) and three times the variance (σ); A range of pixel values ​​greater than a pixel value C1 (=Mean+3σ) is considered to be the range of pixel values ​​of the blood vessels, and in order to set the opacity to a value greater than 0, the opacity value is changed from zero to one between the pixel value C1 (=Mean+3σ) and a pixel value C2 (=Mean+5σ), which is the sum of the mean value (Mean) and five times the variance (σ), and an opacity curve is set in which the opacity value is one for pixel values ​​equal to or greater than the pixel value C2.

13. An image processing device comprising:

8. Receive a 3D TOF image (three-dimensional time of flight image) of the subject's head, generating a 3D TOF image of a brain parenchyma region including blood vessels from the 3D TOF image of the head; A histogram of pixel values ​​of a 3D TOF image of the brain parenchyma region including the blood vessels is obtained; removing pixels of a background region of the 3D TOF image having a pixel value equal to or less than a predetermined threshold value from the histogram to obtain a histogram in which the pixels of the brain parenchyma show a distribution having a peak and the pixels of the blood vessels are distributed in a region having a pixel value larger than that of the distribution having the peak of the pixels of the brain parenchyma; Calculate the average value (Mean) and variance (σ) of the pixel values ​​of the histogram, and set an opacity value of zero by assuming that the pixel value range of the brain parenchyma is equal to or less than a pixel value C1 (=Mean+3σ), which is the sum of the average value (Mean) and three times the variance (σ); A range of pixel values ​​greater than a pixel value C1 (=Mean+3σ) is considered to be the range of pixel values ​​of the blood vessels, and in order to set the opacity to a value greater than 0, the opacity value is changed from zero to one between the pixel value C1 (=Mean+3σ) and a pixel value C2 (=Mean+5σ), which is the sum of the average value (Mean) and five times the variance (σ), and an opacity curve is set in which the opacity value is one for pixel values ​​equal to or greater than the pixel value C2; An image processing method, comprising: generating a volume rendering image of a 3D TOF image of a brain parenchymal region including the blood vessels using the opacity value.

Images