This application is a national phase application of PCT Application No. PCT/KR2020/006853, filed on 27 May 2020, which claims the benefit of and priority to Korean Patent Application No. 10-2020-0017998, filed on 13 Feb. 2020. The entire disclosures of the applications identified in this paragraph are incorporated herein by references.
The present invention relates to a magnetic resonance image processing apparatus and a magnetic resonance image processing method, and more particularly, to a magnetic resonance image processing apparatus for accelerating acquisition of a magnetic resonance image from a magnetic resonance signal using an artificial neural network, and a magnetic resonance image processing method.
In general, medical imaging apparatuses are apparatuses that acquire body information of a patient to provide an image. Medical imaging apparatuses include X-ray imaging apparatuses, ultrasound diagnosis apparatuses, computed tomography apparatuses, and magnetic resonance imaging (MRI) apparatuses.
Here, MRI enables a nuclear magnetic resonance phenomenon to occur in a hydrogen atomic nucleus in the human body by using a magnetic field and non-ionizing radiation which are harmless to the human body, thereby imaging density and physical/chemical properties of the atomic nucleus. The MRI apparatus holds an important position in a diagnosis field using medical images because image capturing conditions are not limited relatively, a variety of information about soft tissue is provided, and an image having excellent contrast is provided.
Specifically, an MRI apparatus is an image diagnostic apparatus that supplies energy having a constant frequency to an atomic nucleus in a state in which a constant magnetic field is applied to the atomic nucleus and converts energy emitted from the atomic nucleus into a signal to diagnose the interior of the human body. Since a proton constituting an atomic nucleus intrinsically has spin angular momentum and a magnetic dipole, when a magnetic field is applied on protons, the protons are arranged in a direction of the magnetic field, and the atomic nucleus performs a precession about the direction of the magnetic field. Images of the human body may be acquired through a nuclear magnetic resonance phenomenon, and the nuclear magnetic resonance phenomenon may be caused by such precession.
Meanwhile, imaging by an MRI apparatus may take 20 minutes to one hour or more according to an area to be imaged and a type of a magnetic resonance (MR) image. That is, there is a disadvantage in that an imaging time of the MRI apparatus is relatively long as compared with other medical imaging apparatuses. Such a disadvantage may put a burden on a patient, and in particular, the disadvantage makes it difficult to perform imaging on a patient with claustrophobia. Therefore, techniques for reducing an imaging time have been developed recently, and improvement is also required in terms of image quality.
The present invention is directed to facilitating a diagnosis of a lesion area of a patient by providing a magnetic resonance image processing apparatus and a magnetic resonance image processing method capable of acquiring a high-quality reconstructed image with shortening a magnetic resonance imaging time.
According to one embodiment of the present invention, a magnetic resonance image processing method performed by a magnetic resonance image processing apparatus includes acquiring a sub-sampled magnetic resonance signal, acquiring first k-space data from the sub-sampled magnetic resonance signal using a first parallel imaging technique, acquiring a first magnetic resonance image from the first k-space data using an inverse fast Fourier transform operation, preprocessing the first magnetic resonance image to generate first input image data, and acquiring a first output magnetic resonance image from the first input image data using a first artificial neural network model.
When the sub-sampled magnetic resonance signal is provided as a plurality of sub-sampled magnetic resonance signals, the first k-space data may be acquired as a plurality of pieces of first k-space data, and the first magnetic resonance image may be acquired as a plurality of first magnetic resonance images, and the preprocessing of the first magnetic resonance image to generate the first input image data may include combining the plurality of first magnetic resonance images to acquire a single second magnetic resonance image and preprocessing the second magnetic resonance image to generate the first input image data.
The combining of the plurality of first magnetic resonance images to acquire the single second magnetic resonance image may be performed using [Equations 1] below.
Ci=an ith-first magnetic resonance image.
Si=a sensitivity matrix of an ith coil.
m=the second magnetic resonance image according to complex value information.
mreal=the second magnetic resonance image according to the real value information.
The combining of the plurality of first magnetic resonance images to acquire the single second magnetic resonance image may be performed using [Equation 2] below.
m=the second magnetic resonance image according to complex value information.
Si*=a conjugate transposed matrix of a sensitivity matrix of an ith coil.
Ci=an itth-first magnetic resonance image.
The first magnetic resonance image may include a target image and an image aliased in a phase encoding direction, and the preprocessing of the first magnetic resonance image to generate the first input image data may include shifting the first magnetic resonance image such that the aliased image matches the target image to acquire the shifted first image, and stacking and grouping the first magnetic resonance image and the shifted first magnetic resonance image in a first direction perpendicular to each of the phase encoding direction and a readout direction to generate the first input image data.
The second magnetic resonance image may include a target image and an image aliased in a phase encoding direction, and the preprocessing of the second magnetic resonance image to generate the first input image data may include shifting the second magnetic resonance image such that the aliased image matches the target image to acquire the shifted second magnetic resonance image, and stacking and grouping the second magnetic resonance image and the shifted second magnetic resonance image in a first direction perpendicular to each of the phase encoding direction and a readout direction to generate the first input image data.
A noise pattern image generated based on the first magnetic resonance image may be input to the first artificial neural network model.
The first artificial neural network model may have an artificial neural network structure including a contracting path and an expanding path, the contracting path may include a plurality of convolution layers and pooling layers, the expanding path may include a plurality of convolution layers and unpooling layers, and the noise pattern image may be input to at least one of the plurality of convolution layers and the unpooling layers of the expanding path.
The magnetic resonance image processing method may include acquiring second k-space data from the first output magnetic resonance image using a fast Fourier transform operation, wherein the acquiring of the second k-space data includes decombining the first output magnetic resonance image to acquire a plurality of second output magnetic resonance images and acquiring the second k-space data from the second output magnetic resonance images using the fast Fourier transform operation.
The decombining of the first output magnetic resonance image to acquire the plurality of second output magnetic resonance images may be performed using [Equation 3] or [Equation 4] below.
{umlaut over (c)}i=Si·{umlaut over (m)} [Equation 3]
ζi=an ith second output magnetic resonance image.
Si=a sensitivity matrix of an ith coil.
{umlaut over (m)}=the first output to complex value information.
{umlaut over (c)}i=Si·(m{umlaut over (r)}eal·exp(i·∂)) [Equation 4]
ζi=the ith-second output magnetic resonance image.
Si=the sensitivity matrix of the ith coil.
m{umlaut over (r)}eal=the first output magnetic resonance image according to real value information.
The magnetic resonance image processing method may include acquiring third k-space data from the second k-space data using a second parallel imaging technique.
The magnetic resonance image processing method may include acquiring a third magnetic resonance image from the third k-space data using the inverse fast Fourier transform operation, preprocessing the third magnetic resonance image to generate second input image data, and acquiring a third output magnetic resonance image from the second input image data using a second artificial neural network model.
According to one embodiment of the present invention, a magnetic resonance image processing apparatus for performing a magnetic resonance image processing method includes a memory in which a magnetic resonance image processing program is stored, and a processor configured to execute the program, wherein the processor executes the program to acquire a sub-sampled magnetic resonance signal, acquire first k-space data from the sub-sampled magnetic resonance signal using a first parallel imaging technique, acquire a first magnetic resonance image from the first k-space data using an inverse fast Fourier transform operation, preprocess the first magnetic resonance image to generate first input image data, and acquire a first output magnetic resonance image from the first input image data using a first artificial neural network model.
According to a magnetic resonance image processing apparatus and a magnetic resonance image processing method according to the present invention, a diagnosis of a lesion area of a patient can be facilitated by providing the magnetic resonance image processing apparatus and the magnetic resonance image processing method capable of providing a high-quality reconstructed image even when a magnetic resonance imaging time is shortened.
In addition, a magnetic resonance image processing apparatus and a magnetic resonance image processing method according to the embodiments of the present invention acquire a plurality of magnetic resonance images and provide an output magnetic resonance image with high accuracy which is reconstructed from input data generated by preprocessing the plurality of magnetic resonance images by using a first artificial neural network model and a second artificial neural network model, thereby facilitating a diagnosis of a lesion area of a patient.
Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings so that those of ordinary skill in the art to which the present invention pertains can easily carry out the embodiments. However, the present invention may be embodied in different forms and is not limited to the embodiments set forth herein. In order to clarify the present invention, parts which are not related with the description will be omitted from the drawings, and like reference numbers will be used to refer to like parts throughout the drawings.
The present specification does not describe all elements of embodiments, and common knowledge in the technical field to which the present invention pertains or the same descriptions of the embodiments will be omitted. The term “part” or “portion” used in the specification may be implemented using hardware or software, and according to embodiments, one “part” or “portion” may be implemented as a single unit or element or include a plurality of units or elements.
In the present specification, an “image” may refer to multi-dimensional data composed of discrete image elements (for example, pixels in a two-dimensional (2D) image and voxels in a three-dimensional (3D) image). For example, the image may include a medical image of an object captured by a magnetic resonance imaging (MRI) apparatus, a computed tomography (CT) apparatus, an ultrasound imaging apparatus, or an X-ray imaging apparatus.
In the present specification, an “object” may be a target to be imaged and include a human, an animal, or a part of a human or animal. For example, the object may include a body part (an organ) or a phantom. The phantom refers to a material having a volume that is very close to a density and effective atomic number of an organism. The phantom may include a spherical phantom having similar properties to the human body.
An MRI system is a system which acquires a sectional image of a part of an object by expressing a strength of a magnetic resonance (MR) signal with respect to a radio frequency (RF) signal generated in a magnetic field having a specific strength through a contrast.
Specifically, in the MRI system, a main magnet generates a static magnetic field to align a magnetic dipole moment direction of a specific atomic nucleus of an object positioned in the static magnetic field in a direction of the static magnetic field direction. A gradient coil may apply a gradient signal to a static magnetic field to generate a gradient magnetic field, thereby inducing resonance frequencies differently according to parts of an object. An RF coil may emit an RF signal so as to match a resonance frequency of a part of an object whose image is to be acquired. In addition, when a gradient magnetic field is generated, the RF coil may receive MR signals having different resonance frequencies emitted from a plurality of parts of the object. The MRI system acquires an image by applying an image reconstruction technique to the MR signals received through such operations. In addition, the MRI system may perform serial or parallel signal processing on a plurality of MR signals received by a multi-channel RF coil to reconstruct the plurality of MR signals into image data.
Hereinafter, an MR image processing apparatus 100 according to one embodiment of the present invention will be described.
The MR image processing apparatus 100 according to one embodiment of the present invention relates to an MR image processing apparatus 100 which is applied directly or indirectly to the above-described MRI system to accelerate acquisition of an MR image from an MR signal.
The MR image processing apparatus 100 may include not only an MRI system capable of acquiring an MR image by detecting an MR signal by itself, but also an image processing apparatus for processing an externally acquired image, and a smartphone, a tablet personnel computer (PC), a PC, a smart television (TV), a micro-server, other home appliances, and other mobile or non-mobile computing devices which have a function of processing an MR image. In addition, the MR image processing apparatus 100 may include wearable devices, which have a communication function and a data processing function, such as a watch, glasses, a hair band, and a ring.
Referring to
The communication module 110 may provide a communication interface in the MR image processing apparatus 100 by interworking with a communication network and may serve to transmit and receive data to and from an MRI apparatus, a user terminal, and a management server. Here, the communication module may be a device including hardware and software necessary for transmitting and receiving signals such as control signals or data signals through a wired/wireless connection with other network devices.
The memory 120 may be a storage medium on which an MR image processing program is recorded. In addition, the memory 120 may perform a function of temporarily or permanently storing data processed by the processor. Here, the memory 120 may include a volatile storage medium or a non-volatile storage medium, but the scope of the present invention is not limited thereto.
The processor 130 may control an overall process performed by the MR image processing program in the MR image processing apparatus 100. Each operation of the process performed by the processor 130 will be described below with reference to
Here, the processor 130 may include any type of device capable of processing data, such as a processor. Here, the “processor” may refer to, for example, a data processing device that is embedded in hardware and includes a circuit physically structured to perform a function expressed in a code or command included in a program. Examples of the data processing device embedded in hardware may include a microprocessor, a central processing unit (CPU), a processor core, a multiprocessor, an application-specific integrated circuit (ASIC), a field programmable gate array (FPGA), and the like, but the scope of the present invention is not limited thereto.
The DB 140 may store a variety of data that is necessary for the MR image processing apparatus 100 to execute a program. For example, the DB may store learning data of an artificial neural network, data about an MR signal 310, k-space data, MR image data, and the like which will be described below.
Meanwhile, in the present invention, a “terminal” may include a wireless communication device with ensured portability and mobility, for example, any type of handheld-based wireless communication device such as a smartphone, a tablet PC, or a notebook computer. In addition, the “terminal” may include wearable devices, which have a communication function and a data processing function, such as a watch, glasses, a hair band, and a ring. Furthermore, the “terminal” may include a wired communication device such as a PC that may be connected to another terminal or server through a network.
In addition, the network refers to a connection structure that enables information exchange between nodes such as devices and servers and includes a local area network (LAN), a wide area network (WAN), Internet (WWW: World Wide Web), a wired or wireless data communication network, a telecommunication network, a wired or wireless television network, and the like. Examples of the wireless data communication network may include 3G, 4G, 5G, 3rd generation partnership project (3GPP), Long Term Evolution (LTE), World Interoperability for Microwave Access (WiMAX), Wi-Fi, Bluetooth communication, infrared communication, ultrasonic communication, visible light communication (VLC), LiFi, and the like, but the scope of the present invention is not limited thereto.
Hereinafter, an MR image processing method performed by the MR image processing apparatus 100 according to one embodiment of the present invention will be described in detail.
Referring to
Referring to
For example, the number of lines of a full-sampled MR signal 310 may be n, and the number of lines of the sub-sampled MR signal 310 may be n/2. Here, when a degree by which the number of sampled lines is decreased is a multiple of ½, it may be considered that an acceleration index of MR imaging is 2. When degrees by which the number of sampled lines is decreased are a multiple of ⅓ and a multiple of ¼, it may be considered that acceleration indexes are 3 and 4.
A sub-sampling method may be divided into a uniform sub-sampling method and a nonuniform sub-sampling method. The uniform sub-sampling method may mean that sampling is performed by constantly maintaining an interval between sampled lines. On the other hand, the nonuniform sub-sampling method may generally mean that sampling is performed more by gradually decreasing an interval between sampled lines in a direction toward a center of sampling data and sampling is performed less by gradually increasing the interval between the sampled lines in a direction away from the center of the sampling data.
Meanwhile, the MR image processing apparatus 100 may be included in an MRI system and may acquire input data corresponding to the sub-sampled MR signal 310 based on an MR signal received from an RF coil. In addition, the MR image processing apparatus 100 may acquire input data corresponding to the sub-sampled MR signal 310 from at least one of an external MR imaging apparatus, an external server, and a DB.
Referring again to
Here, a parallel imaging technique is a type of image reconstruction technique for acquiring k-space data and/or an MR image with high accuracy such as full-sampled k-space data and/or a full-sampled MR image from the sub-sampled MR signal 310 and/or sub-sampled k-space data.
In performing image reconstruction according to the parallel imaging technique, well-known techniques such as sensitivity profiles from an array of coils for encoding and reconstruction in parallel (SPACE RIP), simultaneous acquisition of spatial harmonics (SMASH), partially parallel imaging with localized sensitivity (PILS), generalized autocalibrating partially parallel acquisitions (GRAPPA), and iterative self-consistent parallel imaging reconstruction (SPIRiT) may be applied without limitation as long as the well-known techniques may be applied to the parallel imaging technique.
Referring to
That is, the first parallel imaging technique may be a GRAPPA technique. The GRAPPA may correct data about the sub-sampled MR signal 310 into high accuracy k-space data similar to full-sampled k-space data. In other words, when the number of lines of full-sampled image data is n and the number of lines of sub-sampled image data is n/2, the GRAPPA may estimate and generate the remaining unsampled n/2 lines from n/2 lines of the sub-sampled image data.
For example, in the sub-sampled MR signal 310 acquired from one channel of the RF coil, when a first line, a third line, and a fourth line are sampled and a second line is not sampled, the second line may be estimated by linearly combining the first line and the third line which are lines closest to the second line. As described above, the GRAPPA may estimate the remaining unsampled lines from a linear combination of adjacent lines, thereby correcting image data similarly to full-sampled image data. That is, it is possible to generate image data in a high accuracy k-space form for each channel.
Unlike such GRAPPA, the SPIRiT may estimate unsampled pixel data by linearly combining data of a pixel adjacent to a line including a pixel to be estimated and data of sampled lines adjacent to a corresponding line.
After operation S220, operation S230 of acquiring a first MR image from the first k-space data 320 using an inverse fast Fourier transform (IFFT) operation may be performed.
The GRAPPA may generate the first MR image, which is a reconstructed image, by performing an IFFT on image data in a complete k-space form. Here, the first MR image may include a target image and an image aliased in a phase encoding direction Ky.
Referring to
Meanwhile, when the sub-sampled MR signal 310 is provided as a plurality of sub-sampled MR signals 310, the first k-space data 320 may be acquired as a plurality of first k-space data 320 and the first MR image may be acquired as a plurality of first MR images. That is, when the RF coil for receiving the MR signal 310 is provided as a plurality of RF coils, a plurality of reconstructed images generated to correspond to the plurality of MR signals 310 received through respective channels become the first MR images.
After operation S230, operation S240 of preprocessing the first MR image to generate first input image data 340 may be performed.
Referring to
Referring to
Here, shifting the first MR image 720 may be performed by shifting the aliased images A1 and A2 so as to match the target image T. For example, the aliased image A2 of the object included in the first MR image 710 shifted leftward may be shifted to be disposed at the same position as a position of the target image T of the object included in the original first MR image 720 based on a plane formed by the phase encoding direction Ky and the readout direction Kx. In addition, the aliased image A1 of the object included in the first MR image 730 shifted rightward may be shifted to be disposed at the same position as the position of the target image T of the object included in the original first MR image 720 based on the plane formed by the phase encoding direction Ky and the readout direction Kx.
Referring to
Here, the first direction Kz is a direction parallel to a direction Kz of
Referring to
Specifically, operation S241′ of combining the plurality of first MR images to acquire the single second MR image 330 may be performed using [Equations 1] below.
Ci=an ith-first MR image.
Si=a sensitivity matrix of an ith coil.
m=a second MR image according to complex value information.
mreal=a second MR image according to real value information.
Here, the second MR image 330 may be acquired based on first MR image information and sensitivity information of the RF coil. The first MR image information may be obtained by multiplying information about the second MR image 330 by sensitivity information of each coil as a weight. Meanwhile, as a method of acquiring the second MR image 330, there is a method of acquiring the second MR image 330 by dividing the first MR image by sensitivity information of each corresponding coil, which may cause a problem of noise amplification. Accordingly, in order to prevent noise amplification and to simplify and accelerate a calculation, it may be easy to acquire the second MR image 330 according to real value information. In order to be able to acquire the second MR image 330 according to the real value information, the sum of sizes of sensitivity matrices of respective coils should satisfy 1.
Meanwhile, operation S241′ of combining the plurality of first MR images to acquire the single second MR image 330 may be performed using [Equation 2] below.
m=a second MR image according to complex value information.
Si*=a conjugate transposed matrix of a sensitivity matrix of an ith coil.
Ci is an ith-first MR image.
[Equation 2] may be used to acquire the second MR image 330 based on complex value information. In this case, a calculation may be more complicated as compared with a case in which [Equations 1] are applied to acquire the second MR image 330, but there is an advantage in that the more accurate second MR image 330 may be acquired.
Operation S240′ of preprocessing the first MR image to generate the first input image data 340 may include operation S242′ of preprocessing the combined second MR image 330 using [Equations 1] or [Equation 2] above to generate the first input image data 340.
The second MR image 330 in which the plurality of first MR images are combined may include a target image and an image aliased in the phase encoding direction Ky. In operation S242′ of preprocessing the second MR image 330 to generate the first input image data 340, the same process may be performed as in the above-described operation S240 of preprocessing the first MR image to generate the first input image data 340.
That is, operation S242′ of preprocessing the second MR image 330 to generate the first input image data 340 may include an operation (not shown) of shifting the second MR image 330 such that the aliased image matches the target image to acquire a shifted second MR image 330 and an operation (not shown) of stacking and grouping the second MR image 330 and the shifted second MR image 330 in the first direction Kz perpendicular to the phase encoding direction Ky and the readout direction Kx to generate the first input image data 340.
In summary, through a process of combining the plurality of first MR images to acquire the second MR image 330 and generating the first input image data 340 therefrom, the first input image data 340 with high accuracy may be generated based on a larger amount of data about the MR signals 31 and input to an artificial neural network to be described below.
After operation S240, operation S250 of acquiring a first output MR image 350 from the first input image data 340 using the first artificial neural network model 500 may be performed.
Referring to
For example, the first artificial neural network model 500 may be a model constructed by learning the correlation between at least one sub-sampled MR image and at least one full-sampled MR image in a unit of pixels of one or more sampled lines which are stacked in the phase encoding direction Ky using the neural network.
In addition, the first artificial neural network model 500 may be constructed using a variety of additional data in addition to the sub-sampled MR image and the full-sampled MR image. For example, as the additional data, at least one of k-space data corresponding to the first MR image, real image data, imaginary image data, magnitude image data, phase image data, sensitivity data of a multi-channel RF coil, and data about a noise pattern image NP may be used.
Referring to
The contracting path 510 may include a plurality of convolution layers Cony and pooling layers. For example, the contracting path 510 may include a plurality of 3×3 convolution layers as input layers and a plurality of combinations of pooling layers to which a wavelet transform is applied.
The expanding path 520 may include a plurality of convolution layers and un-pooling layers. For example, the expanding path 520 may include a plurality of combinations of a plurality of 3×3 convolution layers and a plurality of combinations of unpooling layers to which an inverse wavelet transform is applied and may further include a 1×1 convolution layer as an output layer. Meanwhile, the noise pattern image NP to be described below is input to each portion of the expanding path 520 to train the first artificial neural network.
Referring to
Specifically, referring to
Referring to
In addition, referring to
It can be seen that a portion of an image of an object is enlarged below each image. Here, it can be seen that, when an aliased image overlaps a target image, an enlarged portion h2 of the sub-sampled MR image 1220 is different from an enlarged portion h1 of the full-sampled MR image. Such a phenomenon is referred to as a hallucination phenomenon, and it can be confirmed that an image is different even in an enlarged portion h3 of the reconstructed image due to an influence of the sub-sampled MR image. In contrast, in an enlarged portion h4 of the reconstructed image 1240 using the noise pattern image NP, it can be confirmed that a hallucination phenomenon is mitigated similarly to the enlarged portion h1 of the full-sampled MR image 1210.
Referring to an effect of applying the noise pattern image NP, as an exemplary embodiment, the noise pattern image NP generated based on the first MR image may be input to the first artificial neural network model 500. That is, the first artificial neural network model 500 may perform a process of outputting the first output MR image 350 using the first input image data 340 and the noise pattern image NP as input data.
For example, the noise pattern image NP may be input to at least one of the plurality of convolution layers and the unpooling layers of the expanding path 520 of the first artificial neural network model 500. Specifically, the noise pattern image NP may be input to the convolution layer that is an output layer of the expanding path 520 and may be input by being pooled according to an image size of each unpooling layer.
In one embodiment, the MR image processing apparatus 100 may include the self-constructed first artificial neural network model 500. In another embodiment, the MR image processing apparatus 100 may obtain the first artificial neural network model 500 constructed in an external server or an external device from the external server or the external device. The MR image processing apparatus 100 may achieve the purpose of accelerating image acquisition rate and improving image quality by acquiring the first output MR image 350 that is a reconstructed image of a sub-sampled MR based on the first artificial neural network model 500 using the neural network.
Here, the acquired first output MR image 350 may be provided to a user terminal and utilized for diagnosing a lesion area of a patient. Meanwhile, an additional operation of processing an MR image may be performed to increase accuracy of the first output MR image 350.
Referring to
Referring to
Specifically, the operation of decombining the first output MR image 350 to acquire the plurality of second output MR images 400 may be performed using [Equation 3] or [Equation 4] below.
{umlaut over (c)}i=Si·{umlaut over (m)}
{umlaut over (c)}i=ith-second output MR image.
Si=a sensitivity matrix of an ith coil.
{umlaut over (m)}=a first output MR image according to complex value information.
[Equation 3] is an equation applicable when the first output MR image 350 is acquired according to complex value information.
{umlaut over (c)}i=Si·(m{umlaut over (r)}eal·exp(i·∂)) [Equation 4]
Si=the sensitivity matrix of the ith coil.
m{umlaut over (r)}eal=the first output magnetic resonance image according to real value information.
exp(i·φ)=phase information of the first output magnetic resonance image.
[Equation 4] is an equation applicable when the first output MR image 350 is acquired according to real value information. Here, phase information of the first output MR image 350 may be acquired by extracting complex value phase information in the above-described combining process. Accordingly, the first output MR images 350 may be decombined based on sensitivity information of each RF coil to generate a plurality of MR images.
Operation S260 of acquiring the second k-space data 360 may include operation S262 of acquiring the second k-space data 360 from the first output MR images 350 using a fast Fourier transform operation. Accordingly, a plurality of pieces of second k-space data 360 may be acquired by applying a fast Fourier transform to the plurality of MR images generated from the first output MR image 350.
After operation S260, operation S270 of acquiring third k-space data 370 from the second k-space data 360 using a second parallel imaging technique may be performed.
Here, the second parallel imaging technique may be a parallel imaging technique capable of acquiring the corrected full-sampled third k-space data 370 from the full-sampled second k-space data 360. As the second parallel imaging technique, the above-described SPIRiT technique may be used.
Here, since the GRAPPA is used as the first parallel imaging technique and the SPIRiT is used as the second parallel imaging technique, a reconstructed image with higher accuracy is acquired in a process of acquiring a reconstructed image.
After operation S270, operation S280 of acquiring a third MR image from the third k-space data 370 using an IFFT operation may be performed.
After operation S280, operation S290 of preprocessing the third MR image to generate second input image data 390 may be performed.
Here, in a method of preprocessing the third MR image to generate the second input image data 390, the same process may be performed as in the above-described operation S240 of preprocessing the first MR image to generate the first input image data 340.
Meanwhile, in a modified embodiment, operation S290 of preprocessing the third MR image to generate the second input image data 390 may include combining a plurality of third MR images to generate a single fourth MR image 380.
Here, the same method as the above-described method of combining the first MR images to acquire the single second MR image 330 may be used as a method of combining the plurality of third MR images to acquire the single fourth MR image 380.
In addition, operation S290 of preprocessing the third MR image to generate the second input image data 390 may include preprocessing the fourth MR image 380 to generate the second input image data 390.
Here, the same method as the above-described method of preprocessing the second MR image 330 to generate the first input image data 340 may be used as a method of preprocessing the fourth MR image 380 to generate the second input image data 390.
After operation S290, operation S300 of acquiring a third output MR image from the second input image data 390 using a second artificial neural network model 600 may be performed.
Here, the same artificial neural network model as the first artificial neural network model 500 may be used as the second artificial neural network model 600. Alternatively, the first artificial neural network model 500 additionally learning the first input image data 340 may be used, or an artificial neural network model different from the first artificial neural network model 500 may be used.
After operation S300, the second output MR image 400 is Fourier-transformed and/or decombined again and then is Fourier-transformed to acquire a fourth k-space data 410, and based on the second output MR image 400, processes from operation S260 to operation S300 may be repeatedly performed as one set.
In addition, the second output MR image 400 may be provided to a user terminal and used as an image for diagnosing a lesion area of a patient.
Hereinafter, Experimental Examples performed according to the MR image processing apparatus 100 and the MR image processing method according to one embodiment of the present invention will be described.
Referring to
Referring to
In the MR image processing apparatus 100 and the MR image processing method according to the embodiments of the present invention described above, even when a time taken for an MRI apparatus to capture an MR image is shortened, an MR image with high accuracy can be generated through image reconstruction.
In addition, the MR image processing apparatus 100 and the MR image processing method according to the embodiments of the present invention acquire a plurality of MR images, preprocess the plurality of MR images to generate input data, and provide an output MR image with high accuracy reconstructed using the first artificial neural network model 500 and the second artificial neural network model 600, thereby facilitating a diagnosis of a lesion area of a patient.
Meanwhile, the MR image processing method according to one embodiment of the present invention can be embodied in the form of a storage medium including commands executable by a computer such as a program module executed by the computer. A computer-readable medium may be any usable medium which can be accessed by the computer and includes all volatile/non-volatile and removable/non-removable media. In addition, the computer-readable media may include a computer storage medium. The computer storage medium includes all volatile/non-volatile and removable/non-removable media embodied by a certain method or technology for storing information such as a computer-readable command, a data structure, a program module, or other data. The method and system of the present invention have been described in relation to specific embodiments, but components or some or all of operations thereof can be embodied using a computer system having a general-purpose hardware architecture.
The above description is merely illustrative of the technical idea of the present invention, and those skilled in the art to which the present invention pertains may make various modifications and changes without departing from the essential characteristics of the present invention. Therefore, the embodiments disclosed in the present invention are not intended to limit the technical idea of the present invention but to describe the present invention, and the scope of the technical idea of the present invention is not limited by these embodiments. The scope of protection of the present invention should be interpreted by the following claims, and all technical ideas within the scope equivalent thereto should be construed as being included in the scope of the present invention.
Number | Date | Country | Kind |
---|---|---|---|
10-2020-0017998 | Feb 2020 | KR | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/KR2020/006853 | 5/27/2020 | WO |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2021/162176 | 8/19/2021 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
9196062 | Hwang et al. | Nov 2015 | B2 |
10198799 | Park et al. | Feb 2019 | B2 |
20170061620 | Park et al. | Mar 2017 | A1 |
20190369190 | Ye | Dec 2019 | A1 |
20200011951 | Shi et al. | Jan 2020 | A1 |
20200034998 | Schlemper | Jan 2020 | A1 |
20200072933 | Ye | Mar 2020 | A1 |
20200103483 | Hardy | Apr 2020 | A1 |
20210035337 | Kim | Feb 2021 | A1 |
20210166445 | Wang | Jun 2021 | A1 |
20210264645 | Polak | Aug 2021 | A1 |
20220075017 | Sabuncu | Mar 2022 | A1 |
Number | Date | Country |
---|---|---|
2015-533531 | Nov 2015 | JP |
2019-111322 | Jul 2019 | JP |
10-1330638 | Nov 2013 | KR |
10-1659578 | Sep 2016 | KR |
10-2018-0049166 | May 2018 | KR |
10-1923184 | Nov 2018 | KR |
10-2019-0124994 | Nov 2019 | KR |
10-2019-0138107 | Dec 2019 | KR |
WO-2019039708 | Feb 2019 | WO |
Entry |
---|
International Search Report from corresponding PCT Application No. PCT/KR2020/006853, dated Nov. 9, 2020. |
Sriram, A., et al.; “Grappanet:Combining Parallel Imaging With Deep Learning for Multi-Coil MRI Reconstruction”, Electrical Engineering and Systems Science, Dec. 18, 2019, pp. 1-18. |
Number | Date | Country | |
---|---|---|---|
20220309719 A1 | Sep 2022 | US |