WO2023108484A1 - Magnetic resonance imaging method based on wave gradient coding field and deep learning model - Google Patents

Magnetic resonance imaging method based on wave gradient coding field and deep learning model Download PDF

Info

Publication number
WO2023108484A1
WO2023108484A1 PCT/CN2021/138359 CN2021138359W WO2023108484A1 WO 2023108484 A1 WO2023108484 A1 WO 2023108484A1 CN 2021138359 W CN2021138359 W CN 2021138359W WO 2023108484 A1 WO2023108484 A1 WO 2023108484A1
Authority
WO
WIPO (PCT)
Prior art keywords
magnetic resonance
wave
channel data
resonance imaging
model
Prior art date
Application number
PCT/CN2021/138359
Other languages
French (fr)
Chinese (zh)
Inventor
郑海荣
王海峰
梁栋
刘新
刘聪聪
崔卓须
Original Assignee
中国科学院深圳先进技术研究院
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 中国科学院深圳先进技术研究院 filed Critical 中国科学院深圳先进技术研究院
Priority to PCT/CN2021/138359 priority Critical patent/WO2023108484A1/en
Publication of WO2023108484A1 publication Critical patent/WO2023108484A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation

Definitions

  • the invention relates to the technical field of magnetic resonance imaging, in particular to a magnetic resonance imaging method, device, equipment and storage medium based on a wave gradient encoding field and a deep learning model.
  • Parallel imaging is an effective method of accelerating imaging.
  • Parallel imaging is a method that utilizes the spatial encoding capability of multi-channel phased array coils to reduce the number of encoding steps of the gradient magnetic field to achieve acceleration.
  • parallel imaging achieves the purpose of under-sampling through redundant information (prior information) during multi-channel k-space.
  • reconstruction methods can be divided into two categories, one is based on image domain anti-aliasing or solving methods, such as PILS, SENSE and ESPIRiT, etc.; the other is based on k-space data filling methods, such as SMASH, GRAPPA and SPIRiT et al.
  • the coil sensitivity information is estimated by using the acquired k-space reference line (ACS) at first, and the aliasing of the image is generated after the accelerated (under-sampled) k-space becomes the image domain, and then the image is aliased by using CSM performs anti-aliasing; the method based on k-space data filling is similar to the method in the image domain, and the convolution kernel (kernel) used to fill the original data of k-space is estimated through the collected ACS lines.
  • it is an acceleration method based on the image domain or k-space it is necessary to collect ACS lines for estimating CSM or kernel. Then collecting ACS takes a lot of time in practice, and after the ACS line is collected, it takes a lot of time to estimate CSM or kernel using commonly used algorithms (such as ESPIRiT, etc.), which is not friendly to practical applications.
  • a parallel imaging method using Wave-CAIPI (where CAIPI stands for Controlled Phase Shift or Misalignment Operation in Phase and Slice Direction) with controlled aliasing has a remarkable effect in 3D imaging by causing Aliasing, making full use of coil sensitivity changes in three directions, significantly high-magnification accelerated 3D imaging with reduced geometry factors and residual aliasing artifacts.
  • Wave-CAIPI combines Beam Phase Encoding (BPE) and Controlled Aliasing Parallel Imaging (CAIPIRINHA), which can significantly utilize the CSM information aliased in three directions in space, and use the extended SENSE model to solve, the system has a smaller condition number , so that the reconstructed geometric factor is smaller and the distribution is more uniform, so that the reconstructed image has a higher signal-to-noise ratio.
  • BPE Beam Phase Encoding
  • CAIPIRINHA Controlled Aliasing Parallel Imaging
  • the embodiment of the present application provides a magnetic resonance imaging method based on the Wave-CAIPI deep learning model, the method includes: introducing a deep generative network model into the magnetic resonance reconstruction model of the Wave-CAIPI gradient encoding field, and accelerating the magnetic resonance imaging method.
  • Resonance imaging make conjugate symmetry to the physical coil channel data through VCC to generate VCC channel data; merge the physical coil channel data and VCC channel data to reconstruct the geometric factor calculation model.
  • the method before introducing the deep generation network model into the magnetic resonance reconstruction model of the Wave-CAIPI gradient encoding field, the method further includes: using parallel imaging to acquire multiple images corresponding to different channels in accelerated magnetic resonance imaging Under-sampled k-space, based on different coil sensitivity information, the artifact-free image of the reconstructed under-sampled k-space is obtained.
  • the method further includes: introducing a deep generation model that does not require training.
  • the merging of the physical coil channel data and the VCC channel data to reconstruct the geometric factor calculation model includes: using the formula Get the geometric factor g-factor after VCC expansion, where, Indicates the coding matrix after VCC expansion.
  • the embodiment of the present application also provides a magnetic resonance imaging device based on the Wave-CAIPI gradient encoding field and the deep learning model
  • the device includes: an introduction unit for magnetic resonance reconstruction in the Wave-CAIPI gradient encoding field
  • the model introduces a deep generation network model to accelerate magnetic resonance imaging; the generation unit is used to perform conjugate symmetry on the physical coil channel data through VCC to generate VCC channel data; the reconstruction unit is used to merge the physical coil channel data and VCC channel data Rebuild the geometry factor calculation model.
  • the embodiment of the present application also provides a computer device, including a memory, a processor, and a computer program stored in the memory and operable on the processor.
  • the processor executes the program, it implements the The method described in any one of the descriptions of the examples.
  • the embodiment of the present application also provides a computer device, a computer-readable storage medium, on which a computer program is stored, and the computer program is used for: when the computer program is executed by a processor, the computer program according to the present application is implemented.
  • a computer device a computer-readable storage medium, on which a computer program is stored, and the computer program is used for: when the computer program is executed by a processor, the computer program according to the present application is implemented.
  • the magnetic resonance imaging method based on the Wave-CAIPI gradient coding field and the deep learning model provided by the present invention combines Wave-CAIPI, VCC and DGM, which not only utilizes the advantages of Wave-CAIPI and VCC to reduce the system condition number, but is applicable
  • the reconstructed g-factor is smaller and more uniform, so that the reconstructed graphics have a higher signal-to-noise ratio, and there is no need to collect ACS lines to estimate CSM or kernel during the solution process, which greatly shortens the traditional convolutional neural network.
  • the data collection time of training data training network parameters avoids the errors introduced by inaccurate data collection, and does not require time-consuming estimation of CSM or kernel in traditional reconstruction, which can be more suitable for clinical needs.
  • Fig. 1 shows the schematic flow chart of the magnetic resonance imaging method based on Wave-CAIPI gradient encoding field and deep learning model provided by the embodiment of the present application;
  • Fig. 2 shows an exemplary structural block diagram of a magnetic resonance imaging device 200 based on a Wave-CAIPI gradient encoding field and a deep learning model according to an embodiment of the present application;
  • FIG. 3 shows a schematic structural diagram of a computer system suitable for implementing a terminal device according to an embodiment of the present application
  • Fig. 4 shows the structural diagram of the depth generation model provided by the embodiment of the present application.
  • FIG. 5 shows a schematic diagram of the simulation results combined with Wave-CAIPI and VCC provided by the embodiment of the present application
  • Fig. 6 shows a schematic diagram of the reconstruction and comparison results of WV-DGM and WV-SENSE (Wave-VCC SENSE) provided in the embodiment of the present application on brain data.
  • first and second are used for descriptive purposes only, and cannot be interpreted as indicating or implying relative importance or implicitly specifying the quantity of indicated technical features.
  • the features defined as “first” and “second” may explicitly or implicitly include at least one of these features.
  • “plurality” means at least two, such as two, three, etc., unless otherwise specifically defined.
  • the first feature may be in direct contact with the first feature or the first and second feature may be in direct contact with the second feature through an intermediary. touch.
  • “above”, “above” and “above” the first feature on the second feature may mean that the first feature is directly above or obliquely above the second feature, or simply means that the first feature is higher in level than the second feature.
  • “Below”, “beneath” and “beneath” the first feature may mean that the first feature is directly below or obliquely below the second feature, or simply means that the first feature is less horizontally than the second feature.
  • FIG. 1 shows a schematic flowchart of a magnetic resonance imaging method based on a Wave-CAIPI gradient encoding field and a deep learning model provided by an embodiment of the present application.
  • the method includes:
  • Step 110 introducing a deep generative network model into the magnetic resonance reconstruction model of the Wave-CAIPI gradient encoding field to accelerate magnetic resonance imaging;
  • Step 120 performing conjugate symmetry on the physical coil channel data through VCC to generate VCC channel data
  • step 130 the physical coil channel data and the VCC channel data are merged to reconstruct a geometric factor calculation model.
  • the method before introducing the deep generative network model in the magnetic resonance reconstruction model of the Wave-CAIPI gradient encoding field in the present application, the method further includes: using parallel imaging to acquire multiple images corresponding to different channels in accelerated magnetic resonance imaging The undersampled k-space of , by using different coil sensitivity information, obtains the artifact-free image of the reconstructed undersampled k-space.
  • parallel imaging is an effective way to reduce the geometric factor of accelerated imaging, which simultaneously acquires multiple under-sampled k-spaces corresponding to different channels, and the coil sensitivity distribution of each channel is Different, so there is a difference in coil sensitivity encoding between each k-space, that is, there is complementary information, by using different coil sensitivity information, an image without artifacts can be reconstructed from multiple under-sampled k-spaces , in essence, parallel imaging mainly uses redundant information (prior information) between different coil k-spaces to achieve the purpose of undersampling and accelerated acquisition.
  • the signal model is as follows,
  • the effect of the Wave-CAIPI gradient field can be represented by Psf, which is a three-dimensional phase diagram with periodic sinusoidal changes, where t is linearly corresponding to the code k x of the k-space readout direction, Psf(t,y,z) can also be written in discrete form Psf(k x ,y,z), and Psf(k x ,y,z) is Fourier transformed in the y and z directions to become Psf(k x , k y , k z ), it can describe the offset of the Wave-CAIPI track relative to the Cartesian track, and perform an inverse Fourier transform of Psf(k x ,y, z) in the readout direction, which becomes Psf (x, y, z), which can describe the diffusion effect of the Wave-CAIPI gradient field in the readout direction.
  • Psf is a three-dimensional phase diagram with periodic sinusoidal changes, where t is linearly corresponding to the code
  • formulas (3) and (4) can be equivalently written as, (5), where Psf[x,y,z] is where Psf(k x ,y,z) is the inverse Fourier transform in the readout direction.
  • Wave-CAIPI spatial encoding gradient into the encoding matrix, combined with CAIPIRINHA, pixels can be moved (aliased) in 3D space, thereby reducing the condition number of the solution system and shrinking the geometric factor (g-factor).
  • VCC VCC
  • the specific method is to generate the data of the virtual conjugate coil channel by performing conjugate symmetry on the data of the physical coil channel, and merge the two together for reconstruction, which can further Reduce the number of conditions of the system, reduce the geometric factor, and improve the image SNR.
  • wave*(x,y,z) represents the data expanded by VCC, Indicates the Psf that expands the number of channels on the basis of the original Psf[x,y,z].
  • the combination of the physical coil channel data and the VCC channel data in the present application to reconstruct the geometric factor calculation model includes: through the formula Get the geometric factor g-factor after VCC expansion, where, Indicates the coding matrix after VCC expansion.
  • the g-factor after VCC extension is obtained through (5) as follows, (7), where, Represents the coding matrix after VCC extension, (7) represents the g-factor calculation model after VCC extension, the calculation method of the original g-factor is similar to (7), except that E is the model without VCC extension.
  • the method further includes: introducing a deep generation model that does not require training.
  • the solution method based on the image domain is dominant. While using the advantages of Wave-CAIPI, it still needs to estimate the CSM or kernel (k-space data completion method) , it takes a certain amount of time to collect ACS in clinical practice, and it is estimated that CSM or kernel will also take a lot of time.
  • a deep generative model (DGM) that does not require training is introduced, combined with (1), (3) and (4), the model proposed by the present invention is named as Wave-VCC-DGM,
  • G( ⁇ ) represents DGM
  • represents fixed random noise
  • the output result after G( ⁇ ) operation is the image to be solved.
  • the core idea is to use the network parameters of DGM to fit the image to be solved.
  • DGM It is a simple convolution generator (of course, in this invention, it is not limited to the DGM network architecture. DGM is only a model suitable for the image domain. If it is solved in k-space, the model needs to be changed to a network that conforms to k-space solution Structure), consisting of upsampling (Upsampling), class convolution (Convolution-like), BatchNormalization (BN), and activation function layers, the structure is shown in Figure 4.
  • the deep generative network model used in the present invention can make an effective approximation to the traditional CSM or k-space kernel based on the image domain, and contains all the previous assumptions.
  • the model itself is low-parameterized and does not involve Convolution, and has a simple structure.
  • DMG does not use convolution operations, its structure is closely related to convolutional neural networks. Specifically, the network combines pixels between different channels. Since the coupling between linearly arranged pixels is not provided, Therefore, DGM is not a convolution operation in the traditional sense, but similar to convolution operations in convolutional neural networks, weights are shared between spatial locations. The coupling between pixels in DMG comes from upsampling.
  • the present invention combines DGM, Wave-CAIPI and VCC to propose a new network model Wave-VCC-DGM that does not need to be trained, which can not only utilize the reduction system of Wave-CAIPI
  • the advantage of the condition number, and the advantage of VCC combined with the target background phase, further reduces the g-factor of the reconstructed image and improves the SNR of the image.
  • DGM can take advantage of the advantage of not requiring training, and directly use ADAM (not limited to ADAM) to solve directly to meet the requirements of final image fitting, which can greatly reduce the cost of collecting training data in clinical practice. Time, and does not need to consume a lot of time for the estimation of CSM, or the calculation of the kernel, which is very friendly for rapid implementation in the clinic.
  • the Wave-VCC-DGM model can make parallel imaging more flexible.
  • Wave-VCC has better effect
  • the result is shown in Figure 5 can see that compared with the traditional image-based domain (taking SENSE as an example)
  • the reconstruction results of , combined with Wave-CAIPI and VCC have better results.
  • the data is uniformly sampled from the template, achieving a 4x speedup.
  • the simulated Wave coding model is used for preliminary experiments, and the traditional reconstruction based on the image domain (taking SENSE as an example) is carried out respectively.
  • the traditional reconstruction based on the image domain taking SENSE as an example
  • the present invention has also carried out the reconstruction test of data encoded by real Wave space.
  • the Wave encoding amplitude is 1.2mT/m, and the number of cycles is 7. This is just an example, and Wave is not limited to the above parameters.
  • the result is shown in Figure 6 , the results show that WV-DGM has a better effect than traditional WV-SENSE.
  • FIG. 2 shows an exemplary structural block diagram of a magnetic resonance imaging apparatus 200 based on a Wave gradient encoding field and a deep learning model according to an embodiment of the present application.
  • the device includes:
  • the introduction unit 210 is used to introduce a deep generation network model into the magnetic resonance reconstruction model of the Wave-CAIPI gradient encoding field to accelerate magnetic resonance imaging;
  • the generating unit 220 is configured to perform conjugate symmetry on the physical coil channel data through the VCC to generate VCC channel data;
  • the reconstruction unit 230 is configured to combine the physical coil channel data and the VCC channel data to reconstruct the geometric factor calculation model.
  • the units or modules recorded in the device 200 correspond to the steps in the method described with reference to FIG. 1 . Therefore, the operations and features described above for the method are also applicable to the device 200 and the units contained therein, and will not be repeated here.
  • the apparatus 200 may be pre-implemented in the browser of the electronic device or other security applications, and may also be loaded into the browser of the electronic device or its security applications by downloading or other means.
  • the corresponding units in the apparatus 200 may cooperate with the units in the electronic device to implement the solutions of the embodiments of the present application.
  • FIG. 3 shows a schematic structural diagram of a computer system 300 suitable for implementing a terminal device or a server according to an embodiment of the present application.
  • a computer system 300 includes a central processing unit (CPU) 301 that can operate according to a program stored in a read-only memory (ROM) 302 or loaded from a storage section 308 into a random-access memory (RAM) 303 Instead, various appropriate actions and processes are performed.
  • ROM read-only memory
  • RAM random-access memory
  • various programs and data necessary for the operation of the system 300 are also stored.
  • the CPU 301 , ROM 302 , and RAM 303 are connected to each other via a bus 304 .
  • An input/output (I/O) interface 305 is also connected to the bus 304 .
  • the following components are connected to the I/O interface 305: an input section 306 including a keyboard, a mouse, etc.; an output section 307 including a cathode ray tube (CRT), a liquid crystal display (LCD), etc., and a speaker; a storage section 308 including a hard disk, etc. and a communication section 309 including a network interface card such as a LAN card, a modem, or the like.
  • the communication section 309 performs communication processing via a network such as the Internet.
  • a drive 310 is also connected to the I/O interface 305 as needed.
  • a removable medium 311, such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, etc., is mounted on the drive 310 as necessary so that a computer program read therefrom is installed into the storage section 308 as necessary.
  • an embodiment of the present disclosure includes a magnetic resonance imaging method based on a Wave-CAIPI gradient encoding field and a deep learning model, which includes a computer program tangibly embodied on a machine-readable medium, the computer program including a method for executing Program code for the method in Figure 1.
  • the computer program may be downloaded and installed from a network via communication portion 309 and/or installed from removable media 311 .
  • each block in a flowchart or block diagram may represent a module, program segment, or part of code that includes one or more logical functions for implementing specified executable instructions.
  • the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or they may sometimes be executed in the reverse order, depending upon the functionality involved.
  • each block of the block diagrams and/or flowchart illustrations, and combinations of blocks in the block diagrams and/or flowchart illustrations can be implemented by a dedicated hardware-based system that performs the specified functions or operations , or may be implemented by a combination of dedicated hardware and computer instructions.
  • the units or modules involved in the embodiments described in the present application may be implemented by means of software or by means of hardware.
  • the described units or modules can also be set in a processor, for example, it can be described as: a processor includes a first sub-region generating unit, a second sub-region generating unit, and a display area generating unit.
  • a processor includes a first sub-region generating unit, a second sub-region generating unit, and a display area generating unit.
  • the names of these units or modules do not constitute limitations on the units or modules themselves in some cases, for example, the display area generation unit can also be described as "used to generate The cell of the display area of the text".
  • the present application also provides a computer-readable storage medium, which may be the computer-readable storage medium contained in the aforementioned devices in the above-mentioned embodiments; computer-readable storage media stored in the device.
  • the computer-readable storage medium stores one or more programs, and the aforementioned programs are used by one or more processors to execute the text generation method applied to transparent window envelopes described in this application.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

The present application discloses a magnetic resonance imaging method and apparatus based on a wave (Wave-CAIPI) gradient coding field and a deep learning model, a device, and a storage medium. The method comprises: introducing a deep generative model (DGM) into a magnetic resonance reconstruction model of the Wave-CAIPI gradient coding field to accelerate magnetic resonance imaging; performing conjugate symmetry on physical coil channel data by means of a virtual coil concept (VCC) to generate VCC channel data; and merging the physical coil channel data and the VCC channel data to reconstruct a geometric factor calculation model. According to the solution provided by the present application, a VCC (Wave-CAIPI) technology in a Wave-CAIPI gradient coding occasion is combined with the DGM, the advantage of reducing the number of system conditions by the Wave-CAIPI and the VCC is utilized, and g-factor suitable for reconstruction is smaller and more uniform, such that a reconstructed graph has a higher signal-to-noise ratio, and the step of needing a large amount of training data to train network parameters in a conventional convolutional neural network is shortened.

Description

基于波浪梯度编码场和深度学习模型的磁共振成像方法Magnetic Resonance Imaging Method Based on Wave Gradient Encoded Field and Deep Learning Model 技术领域technical field
本发明涉及磁共振成像技术领域,具体涉及一种基于波浪梯度编码场和深度学习模型的磁共振成像方法、装置、设备及其存储介质。The invention relates to the technical field of magnetic resonance imaging, in particular to a magnetic resonance imaging method, device, equipment and storage medium based on a wave gradient encoding field and a deep learning model.
背景技术Background technique
在传统的磁共振成像(MRI)中,扫描时间长是固有的特点。现有的磁共振成像加速方法中,并行成像是一种有效的加速成像的方法,并行成像是一种利用多通道的相控阵列线圈的空间编码能力,减少梯度磁场的编码步数,达到加速的目的,实际中,并行成像是通过多通道k空间之期间的冗余信息(先验信息),达到欠采样的目的。在并行成像中,重建方法可以分为两类,一类是基于图像域解混叠或者求解方法,比如PILS、SENSE和ESPIRiT等;另外一类是基于k空间数据填充的方法,比如SMASH、GRAPPA和SPIRiT等。在图像域的重建方法中,首先使用采集的k空间参考线(ACS)估计线圈敏感度信息(CSM),加速后(欠采样)的k空间变成图像域后产生图像的混叠,然后利用CSM进行解混叠;基于k空间数据填充的方法,类似于图像域的方法,通过采集的ACS线,估计出用于填充k空间原始数据的卷积核(kernel)。无论是基于图像域还基于k空间的加速方法,都需要采集ACS线用于估计CSM或者kernel。然后采集ACS在实际中是需要消耗大量时间的,而且采集完成ACS线之后使用常用的算法(比如ESPIRiT等)估计CSM 或者kernel时需要消耗大量的时间,这对于实际应用是不友好的。Long scan times are inherent in conventional magnetic resonance imaging (MRI). Among the existing magnetic resonance imaging acceleration methods, parallel imaging is an effective method of accelerating imaging. Parallel imaging is a method that utilizes the spatial encoding capability of multi-channel phased array coils to reduce the number of encoding steps of the gradient magnetic field to achieve acceleration. In practice, parallel imaging achieves the purpose of under-sampling through redundant information (prior information) during multi-channel k-space. In parallel imaging, reconstruction methods can be divided into two categories, one is based on image domain anti-aliasing or solving methods, such as PILS, SENSE and ESPIRiT, etc.; the other is based on k-space data filling methods, such as SMASH, GRAPPA and SPIRiT et al. In the reconstruction method of the image domain, the coil sensitivity information (CSM) is estimated by using the acquired k-space reference line (ACS) at first, and the aliasing of the image is generated after the accelerated (under-sampled) k-space becomes the image domain, and then the image is aliased by using CSM performs anti-aliasing; the method based on k-space data filling is similar to the method in the image domain, and the convolution kernel (kernel) used to fill the original data of k-space is estimated through the collected ACS lines. Whether it is an acceleration method based on the image domain or k-space, it is necessary to collect ACS lines for estimating CSM or kernel. Then collecting ACS takes a lot of time in practice, and after the ACS line is collected, it takes a lot of time to estimate CSM or kernel using commonly used algorithms (such as ESPIRiT, etc.), which is not friendly to practical applications.
一种使用Wave-CAIPI(其中CAIPI表示控制相位在相位和选层方向进行偏移或错位操作)可控混叠的并行成像方法在3D成像中具有显著的效果,其通过在三个空间方向造成混叠,充分利用三个方向的线圈敏感度变化,显著高倍加速三维成像的降低几何因子和残留混叠伪影。Wave-CAIPI结合了波束相位编码(BPE)和可控混叠并行成像(CAIPIRINHA),能够显著利用空间上三个方向混叠的CSM信息,利用扩展的SENSE模型求解,系统具有更小的条件数,使得重建的几何因子更小以及分布更加均匀,从而重建图像具有更高的信噪比。尽管Wave-CAIPI能够兼顾更高信噪比的同时实现更高的加速倍数,但是重建过程仍然需要采集ACS估计CSM,这个过程仍然没有避免在基于图像域或者k空间中消耗时间的过程,这在临床的应用中是不友好的。A parallel imaging method using Wave-CAIPI (where CAIPI stands for Controlled Phase Shift or Misalignment Operation in Phase and Slice Direction) with controlled aliasing has a remarkable effect in 3D imaging by causing Aliasing, making full use of coil sensitivity changes in three directions, significantly high-magnification accelerated 3D imaging with reduced geometry factors and residual aliasing artifacts. Wave-CAIPI combines Beam Phase Encoding (BPE) and Controlled Aliasing Parallel Imaging (CAIPIRINHA), which can significantly utilize the CSM information aliased in three directions in space, and use the extended SENSE model to solve, the system has a smaller condition number , so that the reconstructed geometric factor is smaller and the distribution is more uniform, so that the reconstructed image has a higher signal-to-noise ratio. Although Wave-CAIPI can achieve a higher acceleration factor while taking into account a higher signal-to-noise ratio, the reconstruction process still needs to collect ACS to estimate CSM. This process still does not avoid the time-consuming process based on the image domain or k-space. It is unfriendly in clinical application.
发明内容Contents of the invention
鉴于现有技术中的上述缺陷或不足,期望提供一种基于Wave-CAIPI梯度编码场和深度学习模型的磁共振成像方法、装置、设备及其存储介质。In view of the above-mentioned defects or deficiencies in the prior art, it is desired to provide a magnetic resonance imaging method, device, equipment and storage medium based on the Wave-CAIPI gradient encoding field and deep learning model.
第一方面,本申请实施例提供了一种基于Wave-CAIPI深度学习模型的磁共振成像方法,该方法包括:在Wave-CAIPI梯度编码场的磁共振重建模型中引入深度生成网络模型,加速磁共振成像;通过VCC对物理线圈通道数据做共轭对称,生成VCC通道数据;将物理线圈通道数据和VCC通道数据合并重建几何因子计算模型。In the first aspect, the embodiment of the present application provides a magnetic resonance imaging method based on the Wave-CAIPI deep learning model, the method includes: introducing a deep generative network model into the magnetic resonance reconstruction model of the Wave-CAIPI gradient encoding field, and accelerating the magnetic resonance imaging method. Resonance imaging; make conjugate symmetry to the physical coil channel data through VCC to generate VCC channel data; merge the physical coil channel data and VCC channel data to reconstruct the geometric factor calculation model.
在其中一个实施例中,所述在Wave-CAIPI梯度编码场的磁共振重建模型中引入深度生成网络模型之前,该方法还包括:利用并行成像在磁共振加速成像采集多个对应于不同通道的欠采样k空间,通过基于不同的线圈敏感度信息,得到重建后的欠采样k空间的无伪影的图 像。In one of the embodiments, before introducing the deep generation network model into the magnetic resonance reconstruction model of the Wave-CAIPI gradient encoding field, the method further includes: using parallel imaging to acquire multiple images corresponding to different channels in accelerated magnetic resonance imaging Under-sampled k-space, based on different coil sensitivity information, the artifact-free image of the reconstructed under-sampled k-space is obtained.
在其中一个实施例中,所述在Wave-CAIPI梯度编码场的磁共振重建模型中引入深度生成网络模型,加速磁共振成像,包括:在磁共振成像建立Wave梯度场(Wave)的编码模型:wave(x,y,z)=Em(x,y,z),其中,E是编码矩阵,
Figure PCTCN2021138359-appb-000001
M表示在相位编码方向由可控混叠并行成像采样模式导致的偏移混叠,F x表示在x方向的傅里叶变换,S表示线圈敏感度信息,Psf(k x,y,z)为Wave-CAIPI梯度场的作用效果;通过wave(x,y,z)=Em(x,y,z)和
Figure PCTCN2021138359-appb-000002
Figure PCTCN2021138359-appb-000003
得到
Figure PCTCN2021138359-appb-000004
其中,Psf[x,y,z]为Psf(k x,y,z)在读出方向的反傅里叶变换。
In one of the embodiments, the introduction of the deep generation network model in the magnetic resonance reconstruction model of the Wave-CAIPI gradient encoding field accelerates the magnetic resonance imaging, including: setting up the encoding model of the Wave gradient field (Wave) in the magnetic resonance imaging: wave(x,y,z)=Em(x,y,z), where E is the encoding matrix,
Figure PCTCN2021138359-appb-000001
M represents the offset aliasing caused by the controlled aliasing parallel imaging sampling mode in the phase encoding direction, F x represents the Fourier transform in the x direction, S represents the coil sensitivity information, Psf(k x ,y,z) is the effect of the Wave-CAIPI gradient field; through wave(x,y,z)=Em(x,y,z) and
Figure PCTCN2021138359-appb-000002
Figure PCTCN2021138359-appb-000003
get
Figure PCTCN2021138359-appb-000004
Wherein, Psf[x,y,z] is the inverse Fourier transform of Psf(k x ,y,z) in the readout direction.
在其中一个实施例中,所述通过VCC对物理线圈通道数据做共轭对称,生成VCC通道数据,包括:通过wave(x,y,z)=Em(x,y,z)和
Figure PCTCN2021138359-appb-000005
Figure PCTCN2021138359-appb-000006
得到
Figure PCTCN2021138359-appb-000007
Figure PCTCN2021138359-appb-000008
其中,wave*(x,y,z)表示通过VCC扩展后的数据,
Figure PCTCN2021138359-appb-000009
表示在现有的Psf[x,y,z]基础上扩展通道数的Psf。
In one of the embodiments, the VCC is used to perform conjugate symmetry on the physical coil channel data to generate the VCC channel data, including: through wave(x, y, z)=Em(x, y, z) and
Figure PCTCN2021138359-appb-000005
Figure PCTCN2021138359-appb-000006
get
Figure PCTCN2021138359-appb-000007
Figure PCTCN2021138359-appb-000008
Among them, wave*(x,y,z) represents the data expanded by VCC,
Figure PCTCN2021138359-appb-000009
Indicates the Psf that expands the number of channels on the basis of the existing Psf[x,y,z].
在其中一个实施例中,所述生成VCC通道数据之后,该方法还包括:引入无需训练的深度生成模型。In one of the embodiments, after the VCC channel data is generated, the method further includes: introducing a deep generation model that does not require training.
在其中一个实施例中,所述将物理线圈通道数据和VCC通道数据合并重建几何因子计算模型,包括:通过公式
Figure PCTCN2021138359-appb-000010
Figure PCTCN2021138359-appb-000011
得到经过VCC扩展后的几何因子g-factor,其中,
Figure PCTCN2021138359-appb-000012
表示经过VCC扩展后的编码矩阵。
In one of the embodiments, the merging of the physical coil channel data and the VCC channel data to reconstruct the geometric factor calculation model includes: using the formula
Figure PCTCN2021138359-appb-000010
Figure PCTCN2021138359-appb-000011
Get the geometric factor g-factor after VCC expansion, where,
Figure PCTCN2021138359-appb-000012
Indicates the coding matrix after VCC expansion.
第二方面,本申请实施例还提供了一种基于Wave-CAIPI梯度编码场和深度学习模型的磁共振成像装置,该装置包括:引入单元,用于在Wave-CAIPI梯度编码场的磁共振重建模型中引入深度生成网络模 型,加速磁共振成像;生成单元,用于通过VCC对物理线圈通道数据做共轭对称,生成VCC通道数据;重建单元,用于将物理线圈通道数据和VCC通道数据合并重建几何因子计算模型。In the second aspect, the embodiment of the present application also provides a magnetic resonance imaging device based on the Wave-CAIPI gradient encoding field and the deep learning model, the device includes: an introduction unit for magnetic resonance reconstruction in the Wave-CAIPI gradient encoding field The model introduces a deep generation network model to accelerate magnetic resonance imaging; the generation unit is used to perform conjugate symmetry on the physical coil channel data through VCC to generate VCC channel data; the reconstruction unit is used to merge the physical coil channel data and VCC channel data Rebuild the geometry factor calculation model.
第三方面,本申请实施例还提供了一种计算机设备,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如本申请实施例描述中任一所述的方法。In the third aspect, the embodiment of the present application also provides a computer device, including a memory, a processor, and a computer program stored in the memory and operable on the processor. When the processor executes the program, it implements the The method described in any one of the descriptions of the examples.
第四方面,本申请实施例还提供了一种计算机设备一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序用于:所述计算机程序被处理器执行时实现如本申请实施例描述中任一所述的方法。In a fourth aspect, the embodiment of the present application also provides a computer device, a computer-readable storage medium, on which a computer program is stored, and the computer program is used for: when the computer program is executed by a processor, the computer program according to the present application is implemented. The method described in any one of the descriptions of the examples.
本发明的有益效果:Beneficial effects of the present invention:
本发明提供的基于Wave-CAIPI梯度编码场和深度学习模型的磁共振成像方法,将Wave-CAIPI、VCC和DGM相结合,其不仅仅利用了Wave-CAIPI和VCC降低系统条件数的优势,适用重建的g-factor更小,更加均匀,从而重建图形具有更高的信噪比,而且在求解的过程中不需要采集ACS线估计CSM或者kernel,进而大大缩短传统卷积神经网络中需要大量的训练数据训练网络参数的数据采集时间,并且规避了采集数据不准确而引入的误差,也不需要传统重建中需要消耗时间去估计CSM或者kernel,这能更加适应于临床需求。The magnetic resonance imaging method based on the Wave-CAIPI gradient coding field and the deep learning model provided by the present invention combines Wave-CAIPI, VCC and DGM, which not only utilizes the advantages of Wave-CAIPI and VCC to reduce the system condition number, but is applicable The reconstructed g-factor is smaller and more uniform, so that the reconstructed graphics have a higher signal-to-noise ratio, and there is no need to collect ACS lines to estimate CSM or kernel during the solution process, which greatly shortens the traditional convolutional neural network. The data collection time of training data training network parameters avoids the errors introduced by inaccurate data collection, and does not require time-consuming estimation of CSM or kernel in traditional reconstruction, which can be more suitable for clinical needs.
附图说明Description of drawings
通过阅读参照以下附图所作的对非限制性实施例所作的详细描述,本申请的其它特征、目的和优点将会变得更明显:Other characteristics, objects and advantages of the present application will become more apparent by reading the detailed description of non-limiting embodiments made with reference to the following drawings:
图1示出了本申请实施例提供的基于Wave-CAIPI梯度编码场和深度学习模型的磁共振成像方法的流程示意图;Fig. 1 shows the schematic flow chart of the magnetic resonance imaging method based on Wave-CAIPI gradient encoding field and deep learning model provided by the embodiment of the present application;
图2示出了根据本申请一个实施例的基于Wave-CAIPI梯度编码场 和深度学习模型的磁共振成像装置200的示例性结构框图;Fig. 2 shows an exemplary structural block diagram of a magnetic resonance imaging device 200 based on a Wave-CAIPI gradient encoding field and a deep learning model according to an embodiment of the present application;
图3示出了适于用来实现本申请实施例的终端设备的计算机系统的结构示意图;FIG. 3 shows a schematic structural diagram of a computer system suitable for implementing a terminal device according to an embodiment of the present application;
图4示出了本申请实施例提供的深度生成模型结构图;Fig. 4 shows the structural diagram of the depth generation model provided by the embodiment of the present application;
图5示出了本申请实施例提供的结合Wave-CAIPI和VCC的仿真结果示意图;FIG. 5 shows a schematic diagram of the simulation results combined with Wave-CAIPI and VCC provided by the embodiment of the present application;
图6示出了本申请实施例提供的WV-DGM和WV-SENSE(Wave-VCC SENSE)在大脑数据上的重建比较结果的示意图。Fig. 6 shows a schematic diagram of the reconstruction and comparison results of WV-DGM and WV-SENSE (Wave-VCC SENSE) provided in the embodiment of the present application on brain data.
具体实施方式Detailed ways
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图对本发明的具体实施方式做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本发明。但是本发明能够以很多不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似改进,因此本发明不受下面公开的具体实施例的限制。In order to make the above objects, features and advantages of the present invention more comprehensible, specific implementations of the present invention will be described in detail below in conjunction with the accompanying drawings. In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention. However, the present invention can be implemented in many other ways different from those described here, and those skilled in the art can make similar improvements without departing from the connotation of the present invention, so the present invention is not limited by the specific embodiments disclosed below.
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”、“顺时针”、“逆时针”、“轴向”、“径向”、“周向”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。In describing the present invention, it should be understood that the terms "center", "longitudinal", "transverse", "length", "width", "thickness", "upper", "lower", "front", " Back", "Left", "Right", "Vertical", "Horizontal", "Top", "Bottom", "Inner", "Outer", "Clockwise", "Counterclockwise", "Axial" , "radial", "circumferential" and other indicated orientations or positional relationships are based on the orientations or positional relationships shown in the drawings, which are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying the referred device or Elements must have certain orientations, be constructed and operate in certain orientations, and therefore should not be construed as limitations on the invention.
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此, 限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。In addition, the terms "first" and "second" are used for descriptive purposes only, and cannot be interpreted as indicating or implying relative importance or implicitly specifying the quantity of indicated technical features. Thus, the features defined as "first" and "second" may explicitly or implicitly include at least one of these features. In the description of the present invention, "plurality" means at least two, such as two, three, etc., unless otherwise specifically defined.
在本发明中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或成一体;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系,除非另有明确的限定。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。In the present invention, unless otherwise clearly specified and limited, terms such as "installation", "connection", "connection" and "fixation" should be understood in a broad sense, for example, it can be a fixed connection or a detachable connection , or integrated; it may be mechanically connected or electrically connected; it may be directly connected or indirectly connected through an intermediary, and it may be the internal communication of two components or the interaction relationship between two components, unless otherwise specified limit. Those of ordinary skill in the art can understand the specific meanings of the above terms in the present invention according to specific situations.
在本发明中,除非另有明确的规定和限定,第一特征在第二特征“上”或“下”可以是第一和第二特征直接接触,或第一和第二特征通过中间媒介间接接触。而且,第一特征在第二特征“之上”、“上方”和“上面”可是第一特征在第二特征正上方或斜上方,或仅仅表示第一特征水平高度高于第二特征。第一特征在第二特征“之下”、“下方”和“下面”可以是第一特征在第二特征正下方或斜下方,或仅仅表示第一特征水平高度小于第二特征。In the present invention, unless otherwise clearly specified and limited, the first feature may be in direct contact with the first feature or the first and second feature may be in direct contact with the second feature through an intermediary. touch. Moreover, "above", "above" and "above" the first feature on the second feature may mean that the first feature is directly above or obliquely above the second feature, or simply means that the first feature is higher in level than the second feature. "Below", "beneath" and "beneath" the first feature may mean that the first feature is directly below or obliquely below the second feature, or simply means that the first feature is less horizontally than the second feature.
需要说明的是,当元件被称为“固定于”或“设置于”另一个元件,它可以直接在另一个元件上或者也可以存在居中的元件。当一个元件被认为是“连接”另一个元件,它可以是直接连接到另一个元件或者可能同时存在居中元件。本文所使用的术语“垂直的”、“水平的”、“上”、“下”、“左”、“右”以及类似的表述只是为了说明的目的,并不表示是唯一的实施方式。It should be noted that when an element is referred to as being “fixed on” or “disposed on” another element, it may be directly on the other element or there may be an intervening element. When an element is referred to as being "connected to" another element, it can be directly connected to the other element or intervening elements may also be present. As used herein, the terms "vertical", "horizontal", "upper", "lower", "left", "right" and similar expressions are for the purpose of illustration only and are not intended to represent the only embodiments.
请参考图1,图1示出了本申请实施例提供的基于Wave-CAIPI梯度编码场和深度学习模型的磁共振成像方法的流程示意图。Please refer to FIG. 1 , which shows a schematic flowchart of a magnetic resonance imaging method based on a Wave-CAIPI gradient encoding field and a deep learning model provided by an embodiment of the present application.
如图1所示,该方法包括:As shown in Figure 1, the method includes:
步骤110,在Wave-CAIPI梯度编码场的磁共振重建模型中引入深度生成网络模型,加速磁共振成像; Step 110, introducing a deep generative network model into the magnetic resonance reconstruction model of the Wave-CAIPI gradient encoding field to accelerate magnetic resonance imaging;
步骤120,通过VCC对物理线圈通道数据做共轭对称,生成VCC通道数据; Step 120, performing conjugate symmetry on the physical coil channel data through VCC to generate VCC channel data;
步骤130,将物理线圈通道数据和VCC通道数据合并重建几何因子计算模型。In step 130, the physical coil channel data and the VCC channel data are merged to reconstruct a geometric factor calculation model.
采用上述技术方案,将Wave-CAIPI、VCC和DGM相结合,其不仅仅利用了Wave-CAIPI和VCC降低系统条件数的优势,适用重建的g-factor更小,更加均匀,从而重建图形具有更高的信噪比,而且在求解的过程中不需要采集ACS线估计CSM或者kernel,进而大大缩短传统卷积神经网络中需要大量的训练数据训练网络参数的数据采集时间,并且规避了采集数据不准确而引入的误差,也不需要传统重建中需要消耗时间去估计CSM或者kernel,这能更加适应于临床需求。Using the above technical solution, combining Wave-CAIPI, VCC and DGM, it not only takes advantage of the advantages of Wave-CAIPI and VCC to reduce the system condition number, but the g-factor applicable to reconstruction is smaller and more uniform, so that the reconstructed graphics have more It has a high signal-to-noise ratio, and does not need to collect ACS lines to estimate CSM or kernel during the solution process, thereby greatly shortening the data collection time that requires a large amount of training data to train network parameters in traditional convolutional neural networks, and avoiding the need for data collection. Accurate and introduced errors do not require time-consuming estimation of CSM or kernel in traditional reconstruction, which can be more suitable for clinical needs.
在一些实施例中,本申请中的在Wave-CAIPI梯度编码场的磁共振重建模型中引入深度生成网络模型之前,该方法还包括:利用并行成像在磁共振加速成像采集多个对应于不同通道的欠采样k空间,通过基于不同的线圈敏感度信息,得到重建后的欠采样k空间的无伪影的图像。In some embodiments, before introducing the deep generative network model in the magnetic resonance reconstruction model of the Wave-CAIPI gradient encoding field in the present application, the method further includes: using parallel imaging to acquire multiple images corresponding to different channels in accelerated magnetic resonance imaging The undersampled k-space of , by using different coil sensitivity information, obtains the artifact-free image of the reconstructed undersampled k-space.
具体地,在磁共振加速成像中,并行成像是一种有效的减低几何因子的加速成像的手段,其同时采集多个对应于不同通道的欠采样k空间,每个通道的线圈敏感度分布是不同的,所以每个k空间之间存在线圈敏感度编码的差异,即存在互补的信息,通过利用不同的线圈敏感度信息,可以从多个通的欠采样k空间重建得到无伪影的图像,实质上,并行成像主要是利用的是不同线圈k空间之间的冗余信息(先验信息),达到欠采样加速采集的目的,其信号模型如下,Specifically, in accelerated magnetic resonance imaging, parallel imaging is an effective way to reduce the geometric factor of accelerated imaging, which simultaneously acquires multiple under-sampled k-spaces corresponding to different channels, and the coil sensitivity distribution of each channel is Different, so there is a difference in coil sensitivity encoding between each k-space, that is, there is complementary information, by using different coil sensitivity information, an image without artifacts can be reconstructed from multiple under-sampled k-spaces , in essence, parallel imaging mainly uses redundant information (prior information) between different coil k-spaces to achieve the purpose of undersampling and accelerated acquisition. The signal model is as follows,
Figure PCTCN2021138359-appb-000013
Figure PCTCN2021138359-appb-000013
其中,
Figure PCTCN2021138359-appb-000014
(2)其中,g y(t)和g z(t)是在读出期间一对相位差π/2的Wave-CAIPI梯度场。
in,
Figure PCTCN2021138359-appb-000014
(2) where g y (t) and g z (t) are a pair of Wave-CAIPI gradient fields with a phase difference of π/2 during readout.
在一些实施例中,本申请中的在Wave-CAIPI梯度编码场的磁共振重建模型中引入深度生成网络模型,加速磁共振成像,包括:在磁共振成像建立Wave-CAIPI梯度场的编码模型:wave(x,y,z)=Em(x,y,z),其中,E是编码矩阵,
Figure PCTCN2021138359-appb-000015
M表示在相位编码方向由可控混叠并行成像采样模式导致的偏移混叠,f x表示在x方向的傅里叶变换,S表示线圈敏感度信息,Psf(k x,y,z)为Wave-CAIPI梯度场的作用效果;通过wave(x,y,z)=Em(x,y,z)和
Figure PCTCN2021138359-appb-000016
Figure PCTCN2021138359-appb-000017
得到
Figure PCTCN2021138359-appb-000018
其中,Psf[x,y,z]为Psf(k x,y,z)在读出方向的反傅里叶变换。
In some embodiments, the deep generation network model is introduced into the magnetic resonance reconstruction model of the Wave-CAIPI gradient encoding field in this application to accelerate magnetic resonance imaging, including: establishing a coding model of the Wave-CAIPI gradient field in magnetic resonance imaging: wave(x,y,z)=Em(x,y,z), where E is the encoding matrix,
Figure PCTCN2021138359-appb-000015
M represents the offset aliasing caused by the controllable aliasing parallel imaging sampling mode in the phase encoding direction, f x represents the Fourier transform in the x direction, S represents the coil sensitivity information, Psf(k x ,y,z) is the effect of the Wave-CAIPI gradient field; through wave(x,y,z)=Em(x,y,z) and
Figure PCTCN2021138359-appb-000016
Figure PCTCN2021138359-appb-000017
get
Figure PCTCN2021138359-appb-000018
Wherein, Psf[x,y,z] is the inverse Fourier transform of Psf(k x ,y,z) in the readout direction.
具体地,在MRI中,Wave-CAIPI梯度场的作用效果可以用Psf表示,它是一个周期正弦变化的三维相位图,其中t与k空间读出方向的编码k x是一一线性对应的,Psf(t,y,z)也可以写作离散形式Psf(k x,y,z),将Psf(k x,y,z)在y和z方向做傅里叶变换,变成Psf(k x,k y,k z),它可以描述Wave-CAIPI轨迹相对于笛卡尔轨迹的偏移量,将Psf(k x,y,z)在读出方向做反傅里叶变换,即变成Psf(x,y,z),其可以描述Wave-CAIPI梯度场在读出方向的扩散效应。Wave-CAIPI的编码模型可以简化为,wave(x,y,z)=Em(x,y,z)(3),其中E是编码矩阵,
Figure PCTCN2021138359-appb-000019
(4),其中S表示CSM,F x表示在x方向的傅里叶变换,M表示在相位编码方向由CAIPIRINHA采样模式导致的偏移混叠。实际上,公式(3)和(4)可以等价写成,
Figure PCTCN2021138359-appb-000020
Figure PCTCN2021138359-appb-000021
(5),其中Psf[x,y,z]是在Psf(k x,y,z)是在读出方向的反傅里叶变换。通过在编码矩阵中引入Wave-CAIPI空间编码梯度,并且结合CAIPIRINHA,可以使得像素在3D空间中产生移动(混叠),进而降低求解系统条件数,缩小几何因子(g-factor)。
Specifically, in MRI, the effect of the Wave-CAIPI gradient field can be represented by Psf, which is a three-dimensional phase diagram with periodic sinusoidal changes, where t is linearly corresponding to the code k x of the k-space readout direction, Psf(t,y,z) can also be written in discrete form Psf(k x ,y,z), and Psf(k x ,y,z) is Fourier transformed in the y and z directions to become Psf(k x , k y , k z ), it can describe the offset of the Wave-CAIPI track relative to the Cartesian track, and perform an inverse Fourier transform of Psf(k x ,y, z) in the readout direction, which becomes Psf (x, y, z), which can describe the diffusion effect of the Wave-CAIPI gradient field in the readout direction. The encoding model of Wave-CAIPI can be simplified as, wave(x,y,z)=Em(x,y,z)(3), where E is the encoding matrix,
Figure PCTCN2021138359-appb-000019
(4), where S denotes the CSM, F x denotes the Fourier transform in the x direction, and M denotes the offset aliasing caused by the CAIPIRINHA sampling pattern in the phase encoding direction. In fact, formulas (3) and (4) can be equivalently written as,
Figure PCTCN2021138359-appb-000020
Figure PCTCN2021138359-appb-000021
(5), where Psf[x,y,z] is where Psf(k x ,y,z) is the inverse Fourier transform in the readout direction. By introducing the Wave-CAIPI spatial encoding gradient into the encoding matrix, combined with CAIPIRINHA, pixels can be moved (aliased) in 3D space, thereby reducing the condition number of the solution system and shrinking the geometric factor (g-factor).
在一些实施例中,本申请中的通过VCC对物理线圈通道数据做共轭对称,生成VCC通道数据,包括:通过wave(x,y,z)=Em(x,y,z)和
Figure PCTCN2021138359-appb-000022
Figure PCTCN2021138359-appb-000023
得到
Figure PCTCN2021138359-appb-000024
Figure PCTCN2021138359-appb-000025
其中,wave*(x,y,z)表示通过VCC扩展后的数据,
Figure PCTCN2021138359-appb-000026
表示在现有的Psf[x,y,z]基础上扩展通道数的Psf。
In some embodiments, the conjugate symmetry of the physical coil channel data through VCC in this application is used to generate VCC channel data, including: through wave(x, y, z)=Em(x, y, z) and
Figure PCTCN2021138359-appb-000022
Figure PCTCN2021138359-appb-000023
get
Figure PCTCN2021138359-appb-000024
Figure PCTCN2021138359-appb-000025
Among them, wave*(x,y,z) represents the data expanded by VCC,
Figure PCTCN2021138359-appb-000026
Indicates the Psf that expands the number of channels on the basis of the existing Psf[x,y,z].
具体地,使用Wave-CIAPI加速MRI成像的过程中能够有效降低g-factor,提高SNR,VCC(VCC)技术是提高编码算子的系统条件数的另外一种并行成像方法,其中VCC主要利用目标的背景相位和线圈相位,提供额外的编码能力,具体做法是通过对物理线圈通道的数据做共轭对称,生成虚拟的共轭线圈通道的数据,将两者合并到一起用于重建,能够进一步降低系统的条件数量,降低几何因子,提高图像SNR,通过在Wave-CAIPI中引入VCC技术之后,并且结合(3)和(4)可以得到,
Figure PCTCN2021138359-appb-000027
(6)其中,wave*(x,y,z)表示通过VCC扩展后的数据,
Figure PCTCN2021138359-appb-000028
表示在原来的Psf[x,y,z]基础上扩展通道数的Psf。
Specifically, the use of Wave-CIAPI to accelerate MRI imaging can effectively reduce the g-factor and improve SNR. VCC (VCC) technology is another parallel imaging method to improve the system condition number of the encoding operator. VCC mainly uses the target The background phase and the coil phase provide additional encoding capabilities. The specific method is to generate the data of the virtual conjugate coil channel by performing conjugate symmetry on the data of the physical coil channel, and merge the two together for reconstruction, which can further Reduce the number of conditions of the system, reduce the geometric factor, and improve the image SNR. After introducing the VCC technology in Wave-CAIPI, and combining (3) and (4), it can be obtained,
Figure PCTCN2021138359-appb-000027
(6) Among them, wave*(x,y,z) represents the data expanded by VCC,
Figure PCTCN2021138359-appb-000028
Indicates the Psf that expands the number of channels on the basis of the original Psf[x,y,z].
在一些实施例中,本申请中的将物理线圈通道数据和VCC通道数据合并重建几何因子计算模型,包括:通过公式
Figure PCTCN2021138359-appb-000029
Figure PCTCN2021138359-appb-000030
得到经过VCC扩展后的几何因子g-factor,其中,
Figure PCTCN2021138359-appb-000031
表示经过VCC扩展后的编码矩阵。
In some embodiments, the combination of the physical coil channel data and the VCC channel data in the present application to reconstruct the geometric factor calculation model includes: through the formula
Figure PCTCN2021138359-appb-000029
Figure PCTCN2021138359-appb-000030
Get the geometric factor g-factor after VCC expansion, where,
Figure PCTCN2021138359-appb-000031
Indicates the coding matrix after VCC expansion.
具体地,通过(5)得到经过VCC扩展后的g-factor如下,
Figure PCTCN2021138359-appb-000032
Figure PCTCN2021138359-appb-000033
(7),其中,
Figure PCTCN2021138359-appb-000034
表示经过VCC扩展后的编码矩阵,(7)表示经过VCC扩展后的g-factor计算模型,原始g-factor的计算方式和(7)类似,只是E是没有经过VCC扩展后的模型。通过比较传统的g-factor和
Figure PCTCN2021138359-appb-000035
可以得到经过VCC扩展后 的模型会对g-factor有效的降低,提高图像重建的信噪比。
Specifically, the g-factor after VCC extension is obtained through (5) as follows,
Figure PCTCN2021138359-appb-000032
Figure PCTCN2021138359-appb-000033
(7), where,
Figure PCTCN2021138359-appb-000034
Represents the coding matrix after VCC extension, (7) represents the g-factor calculation model after VCC extension, the calculation method of the original g-factor is similar to (7), except that E is the model without VCC extension. By comparing traditional g-factor and
Figure PCTCN2021138359-appb-000035
It can be obtained that the model after VCC expansion will effectively reduce the g-factor and improve the signal-to-noise ratio of image reconstruction.
在一些实施例中,本申请中的所述生成VCC通道数据之后,该方法还包括:引入无需训练的深度生成模型。In some embodiments, after the VCC channel data is generated in the present application, the method further includes: introducing a deep generation model that does not require training.
具体地,在(6)中的求解模型中,使用基于图像域的求解方式占主要地位,这在使用Wave-CAIPI的优势的同时,仍然需要估计CSM或者kernel(k空间的数据补全方法),在临床中需要消耗一定的时间去进行采集ACS,估计CSM或者kernel也需要大量的时间。在本发明中,保留Wave-CAIPI的优势的同时引入一种不需要训练的深度生成模型(DGM),结合(1)、(3)和(4)式,本发明提出的模型为,命名为Wave-VCC-DGM,
Figure PCTCN2021138359-appb-000036
Specifically, in the solution model in (6), the solution method based on the image domain is dominant. While using the advantages of Wave-CAIPI, it still needs to estimate the CSM or kernel (k-space data completion method) , it takes a certain amount of time to collect ACS in clinical practice, and it is estimated that CSM or kernel will also take a lot of time. In the present invention, while retaining the advantages of Wave-CAIPI, a deep generative model (DGM) that does not require training is introduced, combined with (1), (3) and (4), the model proposed by the present invention is named as Wave-VCC-DGM,
Figure PCTCN2021138359-appb-000036
其中,G(ξ)表示DGM,ξ表示固定的随机噪声,经过G(ξ)操作之后输出的结果是待求解的图像,其中核心思想为用DGM的网络参数去拟合需要求解的图像,DGM是一个简单的卷积生成器(当然在本发明中不仅仅局限于DGM网络架构,DGM只是一个适用于图像域的模型,如果在k空间进行求解的话,模型需要更改为符合k空间求解的网络结构),由上采样(Upsampling)、类卷积(Convolution-like)、BatchNormalization(BN)、和激活函数层组成,结构如图4所示。Among them, G(ξ) represents DGM, ξ represents fixed random noise, and the output result after G(ξ) operation is the image to be solved. The core idea is to use the network parameters of DGM to fit the image to be solved. DGM It is a simple convolution generator (of course, in this invention, it is not limited to the DGM network architecture. DGM is only a model suitable for the image domain. If it is solved in k-space, the model needs to be changed to a network that conforms to k-space solution Structure), consisting of upsampling (Upsampling), class convolution (Convolution-like), BatchNormalization (BN), and activation function layers, the structure is shown in Figure 4.
在本发明中使用的深度生成网络模型,其本身能够对传统的基于图像域的CSM或者k空间的kernel进行一个有效的近似,并且包含之前所有的假设,模型本身是低参数化的,不涉及卷积,并且具有简单的结构。虽然DMG不使用卷积操作,但是其结构与卷积神经网络密切相关,具体来说,该网络将不同通道之间的像素组合,由于按照线性排列的像素之间的耦合性并没有被提供,因此DGM并不是传统意义上的卷积操作,但是和卷积神经网络中的卷积操作类似,权重在空间位置之间共享。DMG中像素之间的耦合来源于上采样。在(8)中可以看出,本发明将DGM、Wave-CAIPI以及VCC相结合,提出一种新的并不 需要进行训练的网络模型Wave-VCC-DGM,不仅可以利用Wave-CAIPI的降低系统条件数的优势,而且利用VCC结合目标背景相位的优势,进一步降低重建图像的g-factor,提高图像的SNR。The deep generative network model used in the present invention can make an effective approximation to the traditional CSM or k-space kernel based on the image domain, and contains all the previous assumptions. The model itself is low-parameterized and does not involve Convolution, and has a simple structure. Although DMG does not use convolution operations, its structure is closely related to convolutional neural networks. Specifically, the network combines pixels between different channels. Since the coupling between linearly arranged pixels is not provided, Therefore, DGM is not a convolution operation in the traditional sense, but similar to convolution operations in convolutional neural networks, weights are shared between spatial locations. The coupling between pixels in DMG comes from upsampling. In (8), it can be seen that the present invention combines DGM, Wave-CAIPI and VCC to propose a new network model Wave-VCC-DGM that does not need to be trained, which can not only utilize the reduction system of Wave-CAIPI The advantage of the condition number, and the advantage of VCC combined with the target background phase, further reduces the g-factor of the reconstructed image and improves the SNR of the image.
除此之外,DGM能够使用不需要训练的优势,直接使用ADAM(不局限于ADAM)这类求解算法直接进行求解,达到最终图像拟合的要求,这在临床中可以大大减低采集训练数据的时间,并且不需要消耗大量的时间用于CSM的估计,或者kernel的计算,这对于临床中快速实现是很友好的。Wave-VCC-DGM模型能够使得并行成像更加灵活,在传统的求解并行成像中的两类方法中(基于k空间和基于图像域),都会增加求解过程中的复杂度,比如基于k空间的计算kernel和基于图像域的估计CSM,在非笛卡尔坐标系(Non-Cartesian)中计算kernel是相对比较复杂的,但是本发明中提出的模型可以避免这个问题,因为本发明使用Calibration-free的方式进行求解,大大降低求解过程的复杂度。In addition, DGM can take advantage of the advantage of not requiring training, and directly use ADAM (not limited to ADAM) to solve directly to meet the requirements of final image fitting, which can greatly reduce the cost of collecting training data in clinical practice. Time, and does not need to consume a lot of time for the estimation of CSM, or the calculation of the kernel, which is very friendly for rapid implementation in the clinic. The Wave-VCC-DGM model can make parallel imaging more flexible. In the two traditional methods of solving parallel imaging (based on k-space and based on image domain), it will increase the complexity of the solution process, such as calculation based on k-space The kernel and the estimated CSM based on the image domain are relatively complicated to calculate the kernel in the non-Cartesian coordinate system (Non-Cartesian), but the model proposed in the present invention can avoid this problem, because the present invention uses the Calibration-free method Solving, greatly reducing the complexity of the solving process.
本发明经过实验证明,Wave-VCC具有更好的效果,首先使用仿真的Wave编码(Psf为仿真),结果如图5所示,可以看到相对于传统的基于图像域(以SENSE为例)的重建结果,结合Wave-CAIPI和VCC具有更好的结果。数据以均匀采样模板,实现4倍加速。The present invention proves through experiment, Wave-VCC has better effect, first uses the Wave coding of emulation (Psf is emulation), the result is shown in Figure 5, can see that compared with the traditional image-based domain (taking SENSE as an example) The reconstruction results of , combined with Wave-CAIPI and VCC have better results. The data is uniformly sampled from the template, achieving a 4x speedup.
在图5中使用仿真的Wave编码模型进行初步试验,分别进行了传统基于图像域的重建(以SENSE为例),使用4倍加速的均匀采样mask,可以看到Wave-CAIPI和VCC(Wave-VCC)结合的结果具有更好的重建信噪比。本发明还进行了经过真实Wave空间编码的数据的重建试验,Wave编码幅值为1.2mT/m,cycle的数量为7,这仅仅是示例,Wave不局限于以上参数,结果如图6所示,结果显示WV-DGM相对于传统的WV-SENSE具有更好的效果。In Figure 5, the simulated Wave coding model is used for preliminary experiments, and the traditional reconstruction based on the image domain (taking SENSE as an example) is carried out respectively. Using a 4-fold accelerated uniform sampling mask, you can see that Wave-CAIPI and VCC (Wave- VCC) combined results have a better reconstruction signal-to-noise ratio. The present invention has also carried out the reconstruction test of data encoded by real Wave space. The Wave encoding amplitude is 1.2mT/m, and the number of cycles is 7. This is just an example, and Wave is not limited to the above parameters. The result is shown in Figure 6 , the results show that WV-DGM has a better effect than traditional WV-SENSE.
进一步地,参考图2,图2示出了根据本申请一个实施例的基于 Wave梯度编码场和深度学习模型的磁共振成像装置200的示例性结构框图。Further, referring to FIG. 2 , FIG. 2 shows an exemplary structural block diagram of a magnetic resonance imaging apparatus 200 based on a Wave gradient encoding field and a deep learning model according to an embodiment of the present application.
如图2所示,该装置包括:As shown in Figure 2, the device includes:
引入单元210,用于在Wave-CAIPI梯度编码场的磁共振重建模型中引入深度生成网络模型,加速磁共振成像;The introduction unit 210 is used to introduce a deep generation network model into the magnetic resonance reconstruction model of the Wave-CAIPI gradient encoding field to accelerate magnetic resonance imaging;
生成单元220,用于通过VCC对物理线圈通道数据做共轭对称,生成VCC通道数据;The generating unit 220 is configured to perform conjugate symmetry on the physical coil channel data through the VCC to generate VCC channel data;
重建单元230,用于将物理线圈通道数据和VCC通道数据合并重建几何因子计算模型。The reconstruction unit 230 is configured to combine the physical coil channel data and the VCC channel data to reconstruct the geometric factor calculation model.
应当理解,装置200中记载的诸单元或模块与参考图1描述的方法中的各个步骤相对应。由此,上文针对方法描述的操作和特征同样适用于装置200及其中包含的单元,在此不再赘述。装置200可以预先实现在电子设备的浏览器或其他安全应用中,也可以通过下载等方式而加载到电子设备的浏览器或其安全应用中。装置200中的相应单元可以与电子设备中的单元相互配合以实现本申请实施例的方案。It should be understood that the units or modules recorded in the device 200 correspond to the steps in the method described with reference to FIG. 1 . Therefore, the operations and features described above for the method are also applicable to the device 200 and the units contained therein, and will not be repeated here. The apparatus 200 may be pre-implemented in the browser of the electronic device or other security applications, and may also be loaded into the browser of the electronic device or its security applications by downloading or other means. The corresponding units in the apparatus 200 may cooperate with the units in the electronic device to implement the solutions of the embodiments of the present application.
下面参考图3,其示出了适于用来实现本申请实施例的终端设备或服务器的计算机系统300的结构示意图。Referring now to FIG. 3 , it shows a schematic structural diagram of a computer system 300 suitable for implementing a terminal device or a server according to an embodiment of the present application.
如图3所示,计算机系统300包括中央处理单元(CPU)301,其可以根据存储在只读存储器(ROM)302中的程序或者从存储部分308加载到随机访问存储器(RAM)303中的程序而执行各种适当的动作和处理。在RAM303中,还存储有系统300操作所需的各种程序和数据。CPU301、ROM302以及RAM303通过总线304彼此相连。输入/输出(I/O)接口305也连接至总线304。As shown in FIG. 3 , a computer system 300 includes a central processing unit (CPU) 301 that can operate according to a program stored in a read-only memory (ROM) 302 or loaded from a storage section 308 into a random-access memory (RAM) 303 Instead, various appropriate actions and processes are performed. In the RAM 303, various programs and data necessary for the operation of the system 300 are also stored. The CPU 301 , ROM 302 , and RAM 303 are connected to each other via a bus 304 . An input/output (I/O) interface 305 is also connected to the bus 304 .
以下部件连接至I/O接口305:包括键盘、鼠标等的输入部分306;包括诸如阴极射线管(CRT)、液晶显示器(LCD)等以及扬声器等的输出部分307;包括硬盘等的存储部分308;以及包括诸如LAN卡、调制 解调器等的网络接口卡的通信部分309。通信部分309经由诸如因特网的网络执行通信处理。驱动器310也根据需要连接至I/O接口305。可拆卸介质311,诸如磁盘、光盘、磁光盘、半导体存储器等等,根据需要安装在驱动器310上,以便于从其上读出的计算机程序根据需要被安装入存储部分308。The following components are connected to the I/O interface 305: an input section 306 including a keyboard, a mouse, etc.; an output section 307 including a cathode ray tube (CRT), a liquid crystal display (LCD), etc., and a speaker; a storage section 308 including a hard disk, etc. and a communication section 309 including a network interface card such as a LAN card, a modem, or the like. The communication section 309 performs communication processing via a network such as the Internet. A drive 310 is also connected to the I/O interface 305 as needed. A removable medium 311, such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, etc., is mounted on the drive 310 as necessary so that a computer program read therefrom is installed into the storage section 308 as necessary.
特别地,根据本公开的实施例,上文参考图1描述的过程可以被实现为计算机软件程序。例如,本公开的实施例包括一种基于Wave-CAIPI梯度编码场和深度学习模型的磁共振成像方法,其包括有形地包含在机器可读介质上的计算机程序,所述计算机程序包含用于执行图1的方法的程序代码。在这样的实施例中,该计算机程序可以通过通信部分309从网络上被下载和安装,和/或从可拆卸介质311被安装。In particular, according to an embodiment of the present disclosure, the process described above with reference to FIG. 1 may be implemented as a computer software program. For example, an embodiment of the present disclosure includes a magnetic resonance imaging method based on a Wave-CAIPI gradient encoding field and a deep learning model, which includes a computer program tangibly embodied on a machine-readable medium, the computer program including a method for executing Program code for the method in Figure 1. In such an embodiment, the computer program may be downloaded and installed from a network via communication portion 309 and/or installed from removable media 311 .
附图中的流程图和框图,图示了按照本发明各种实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段、或代码的一部分,前述模块、程序段、或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个接连地表示的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或操作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。The flowchart and block diagrams in the Figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present invention. In this regard, each block in a flowchart or block diagram may represent a module, program segment, or part of code that includes one or more logical functions for implementing specified executable instructions. It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or they may sometimes be executed in the reverse order, depending upon the functionality involved. It should also be noted that each block of the block diagrams and/or flowchart illustrations, and combinations of blocks in the block diagrams and/or flowchart illustrations, can be implemented by a dedicated hardware-based system that performs the specified functions or operations , or may be implemented by a combination of dedicated hardware and computer instructions.
描述于本申请实施例中所涉及到的单元或模块可以通过软件的方式实现,也可以通过硬件的方式来实现。所描述的单元或模块也可以设置在处理器中,例如,可以描述为:一种处理器包括第一子区域生 成单元、第二子区域生成单元以及显示区域生成单元。其中,这些单元或模块的名称在某种情况下并不构成对该单元或模块本身的限定,例如,显示区域生成单元还可以被描述为“用于根据第一子区域和第二子区域生成文本的显示区域的单元”。The units or modules involved in the embodiments described in the present application may be implemented by means of software or by means of hardware. The described units or modules can also be set in a processor, for example, it can be described as: a processor includes a first sub-region generating unit, a second sub-region generating unit, and a display area generating unit. Wherein, the names of these units or modules do not constitute limitations on the units or modules themselves in some cases, for example, the display area generation unit can also be described as "used to generate The cell of the display area of the text".
作为另一方面,本申请还提供了一种计算机可读存储介质,该计算机可读存储介质可以是上述实施例中前述装置中所包含的计算机可读存储介质;也可以是单独存在,未装配入设备中的计算机可读存储介质。计算机可读存储介质存储有一个或者一个以上程序,前述程序被一个或者一个以上的处理器用来执行描述于本申请的应用于透明窗口信封的文本生成方法。As another aspect, the present application also provides a computer-readable storage medium, which may be the computer-readable storage medium contained in the aforementioned devices in the above-mentioned embodiments; computer-readable storage media stored in the device. The computer-readable storage medium stores one or more programs, and the aforementioned programs are used by one or more processors to execute the text generation method applied to transparent window envelopes described in this application.
以上描述仅为本申请的较佳实施例以及对所运用技术原理的说明。本领域技术人员应当理解,本申请中所涉及的发明范围,并不限于上述技术特征的特定组合而成的技术方案,同时也应涵盖在不脱离前述发明构思的情况下,由上述技术特征或其等同特征进行任意组合而形成的其它技术方案。例如上述特征与本申请中公开的(但不限于)具有类似功能的技术特征进行互相替换而形成的技术方案。The above description is only a preferred embodiment of the present application and an illustration of the applied technical principle. Those skilled in the art should understand that the scope of the invention involved in this application is not limited to the technical solution formed by the specific combination of the above-mentioned technical features, and should also cover the technical solutions formed by the above-mentioned technical features or without departing from the aforementioned inventive concept. Other technical solutions formed by any combination of equivalent features. For example, a technical solution formed by replacing the above-mentioned features with technical features with similar functions disclosed in (but not limited to) this application.

Claims (10)

  1. 一种基于波浪梯度编码场和深度学习模型的磁共振成像方法,其特征在于,该方法包括:A magnetic resonance imaging method based on a wave gradient encoding field and a deep learning model, characterized in that the method comprises:
    在波浪梯度编码场的磁共振重建模型中引入深度生成网络模型,加速磁共振成像;Introduce a deep generative network model into the magnetic resonance reconstruction model of the wave gradient encoding field to accelerate magnetic resonance imaging;
    通过虚拟共轭线圈对物理线圈通道数据做共轭对称,生成虚拟共轭线圈通道数据;Make conjugate symmetry to the physical coil channel data through the virtual conjugate coil, and generate the virtual conjugate coil channel data;
    将物理线圈通道数据和虚拟共轭线圈通道数据合并重建几何因子计算模型。The physical coil channel data and the virtual conjugate coil channel data are combined to reconstruct the geometric factor calculation model.
  2. 根据权利要求1所述基于波浪梯度编码场和深度学习模型的磁共振成像方法,其特征在于,所述在波浪梯度编码场的磁共振重建模型中引入深度生成网络模型之前,该方法还包括:According to the magnetic resonance imaging method based on the wave gradient encoding field and the deep learning model according to claim 1, it is characterized in that, before the introduction of the deep generation network model in the magnetic resonance reconstruction model of the wave gradient encoding field, the method also includes:
    利用并行成像在磁共振加速成像采集多个对应于不同通道的欠采样k空间,Using parallel imaging to acquire multiple undersampled k-spaces corresponding to different channels in Accelerated Magnetic Resonance Imaging,
    通过基于不同的线圈敏感度信息,得到重建后的欠采样k空间的无伪影的图像。Based on different coil sensitivity information, artifact-free images of reconstructed under-sampled k-space are obtained.
  3. 根据权利要求2所述基于波浪梯度编码场和深度学习模型的磁共振成像方法,其特征在于,所述在波浪梯度编码场的磁共振重建模型中引入深度生成网络模型,加速磁共振成像,包括:According to the magnetic resonance imaging method based on the wave gradient encoding field and the deep learning model according to claim 2, it is characterized in that, in the magnetic resonance reconstruction model of the wave gradient encoding field, the depth generation network model is introduced to accelerate the magnetic resonance imaging, including :
    在磁共振成像建立波浪梯度场的编码模型:wave(x,y,z)=Em(x,y,z),其中,E是编码矩阵,
    Figure PCTCN2021138359-appb-100001
    M表示在相位编码方向由可控混叠并行成像采样模式导致的偏移混叠,F x表示在x方向的傅里叶变换,S表示线圈敏感度信息,Psf(k x,y,z)为波浪梯度场的作用效果;
    Establish the encoding model of the wave gradient field in magnetic resonance imaging: wave(x,y,z)=Em(x,y,z), where E is the encoding matrix,
    Figure PCTCN2021138359-appb-100001
    M represents the offset aliasing caused by the controlled aliasing parallel imaging sampling mode in the phase encoding direction, F x represents the Fourier transform in the x direction, S represents the coil sensitivity information, Psf(k x ,y,z) is the effect of wave gradient field;
    通过wave(x,y,z)=Em(x,y,z)和
    Figure PCTCN2021138359-appb-100002
    得到
    Figure PCTCN2021138359-appb-100003
    其中,Psf[x,y,z]为Psf(k x,y,z)在读出方向的反傅里叶变换。
    By wave(x,y,z)=Em(x,y,z) and
    Figure PCTCN2021138359-appb-100002
    get
    Figure PCTCN2021138359-appb-100003
    Wherein, Psf[x,y,z] is the inverse Fourier transform of Psf(k x ,y,z) in the readout direction.
  4. 根据权利要求3所述基于波浪梯度编码场和深度学习模型的磁共振成像方法,其特征在于,所述通过虚拟共轭线圈对物理线圈通道数据做共轭对称,生成虚拟共轭线圈通道数据,包括:According to the magnetic resonance imaging method based on the wave gradient encoding field and the deep learning model according to claim 3, it is characterized in that, the virtual conjugate coil is used to perform conjugate symmetry on the physical coil channel data to generate virtual conjugate coil channel data, include:
    通过wave(x,y,z)=Em(x,y,z)和
    Figure PCTCN2021138359-appb-100004
    得到
    Figure PCTCN2021138359-appb-100005
    其中,wave*(x,y,z)表示通过虚拟共轭线圈扩展后的数据,
    Figure PCTCN2021138359-appb-100006
    表示在现有的Psf[x,y,z]基础上扩展通道数的Psf。
    By wave(x,y,z)=Em(x,y,z) and
    Figure PCTCN2021138359-appb-100004
    get
    Figure PCTCN2021138359-appb-100005
    Among them, wave*(x,y,z) represents the data expanded by the virtual conjugate coil,
    Figure PCTCN2021138359-appb-100006
    Indicates the Psf that expands the number of channels on the basis of the existing Psf[x,y,z].
  5. 根据权利要求4所述基于波浪梯度编码场和深度学习模型的磁共振成像方法,其特征在于,所述生成虚拟共轭线圈通道数据之后,该方法还包括:According to the magnetic resonance imaging method based on the wave gradient encoding field and the deep learning model according to claim 4, it is characterized in that, after the virtual conjugate coil channel data is generated, the method also includes:
    引入无需训练的深度生成模型。Introduce training-free deep generative models.
  6. 根据权利要求4所述基于波浪梯度编码场和深度学习模型的磁共振成像方法,其特征在于,所述将物理线圈通道数据和虚拟共轭线圈通道数据合并重建几何因子计算模型,包括:According to the magnetic resonance imaging method based on the wave gradient encoding field and the deep learning model according to claim 4, it is characterized in that, the physical coil channel data and the virtual conjugate coil channel data are merged to reconstruct the geometric factor calculation model, including:
    通过公式
    Figure PCTCN2021138359-appb-100007
    得到经过VCC扩展后的几何因子g-factor,其中,
    Figure PCTCN2021138359-appb-100008
    表示经过VCC扩展后的编码矩阵。
    by formula
    Figure PCTCN2021138359-appb-100007
    Get the geometric factor g-factor after VCC expansion, where,
    Figure PCTCN2021138359-appb-100008
    Indicates the coding matrix after VCC expansion.
  7. 一种基于基于波浪梯度编码场和深度学习模型的磁共振成像装置,其特征在于,该装置包括:A magnetic resonance imaging device based on a wave gradient encoding field and a deep learning model, characterized in that the device includes:
    引入单元,用于在波浪梯度编码场的磁共振重建模型中引入深度生成网络模型,加速磁共振成像;The introduction unit is used to introduce a deep generation network model into the magnetic resonance reconstruction model of the wave gradient encoding field to accelerate magnetic resonance imaging;
    生成单元,用于通过虚拟共轭线圈对物理线圈通道数据做共轭对称,生成虚拟共轭线圈通道数据;A generation unit is used to perform conjugate symmetry on the physical coil channel data through the virtual conjugate coil to generate virtual conjugate coil channel data;
    重建单元,用于将物理线圈通道数据和虚拟共轭线圈通道数据合并重建几何因子计算模型。The reconstruction unit is used to combine the physical coil channel data and the virtual conjugate coil channel data to reconstruct the geometric factor calculation model.
  8. 根据权利要求7所述基于基于波浪梯度编码场和深度学习模型的磁共振成像装置,其特征在于,所述在波浪梯度编码场的磁共振重建模型中引入深度生成网络模型之前,该装置还包括:According to claim 7, based on the magnetic resonance imaging device based on wave gradient coding field and deep learning model, it is characterized in that, before the introduction of deep generation network model in the magnetic resonance reconstruction model of wave gradient coding field, the device also includes :
    利用并行成像在磁共振加速成像采集多个对应于不同通道的欠采样k空间,Using parallel imaging to acquire multiple undersampled k-spaces corresponding to different channels in Accelerated Magnetic Resonance Imaging,
    通过基于不同的线圈敏感度信息,得到重建后的欠采样k空间的无伪影的图像。Based on different coil sensitivity information, artifact-free images of reconstructed under-sampled k-space are obtained.
  9. 一种计算机设备,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1-6中任一所述的方法。A computer device, comprising a memory, a processor, and a computer program stored on the memory and operable on the processor, characterized in that, when the processor executes the program, it implements any of claims 1-6 described method.
  10. 一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序用于:A computer-readable storage medium having stored thereon a computer program for:
    所述计算机程序被处理器执行时实现如权利要求1-6中任一所述的方法。When the computer program is executed by the processor, the method according to any one of claims 1-6 is implemented.
PCT/CN2021/138359 2021-12-15 2021-12-15 Magnetic resonance imaging method based on wave gradient coding field and deep learning model WO2023108484A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
PCT/CN2021/138359 WO2023108484A1 (en) 2021-12-15 2021-12-15 Magnetic resonance imaging method based on wave gradient coding field and deep learning model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2021/138359 WO2023108484A1 (en) 2021-12-15 2021-12-15 Magnetic resonance imaging method based on wave gradient coding field and deep learning model

Publications (1)

Publication Number Publication Date
WO2023108484A1 true WO2023108484A1 (en) 2023-06-22

Family

ID=86775017

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2021/138359 WO2023108484A1 (en) 2021-12-15 2021-12-15 Magnetic resonance imaging method based on wave gradient coding field and deep learning model

Country Status (1)

Country Link
WO (1) WO2023108484A1 (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108872899A (en) * 2018-06-28 2018-11-23 中国科学院深圳先进技术研究院 Reconstructing device, method, magnetic resonance system, equipment and the medium of magnetic resonance image
CN108957375A (en) * 2018-07-23 2018-12-07 深圳先进技术研究院 MR imaging method, device, equipment and storage medium
CN113298902A (en) * 2021-05-13 2021-08-24 中国科学院深圳先进技术研究院 Method for reconstructing magnetic resonance image in convoluted field of view, computer device and storage medium
CN113313778A (en) * 2021-05-13 2021-08-27 中国科学院深圳先进技术研究院 Magnetic resonance image reconstruction method, computer device, and storage medium

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108872899A (en) * 2018-06-28 2018-11-23 中国科学院深圳先进技术研究院 Reconstructing device, method, magnetic resonance system, equipment and the medium of magnetic resonance image
CN108957375A (en) * 2018-07-23 2018-12-07 深圳先进技术研究院 MR imaging method, device, equipment and storage medium
CN113298902A (en) * 2021-05-13 2021-08-24 中国科学院深圳先进技术研究院 Method for reconstructing magnetic resonance image in convoluted field of view, computer device and storage medium
CN113313778A (en) * 2021-05-13 2021-08-27 中国科学院深圳先进技术研究院 Magnetic resonance image reconstruction method, computer device, and storage medium

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
ZHILANG QIU; SEN JIA; SHI SU; YANJIE ZHU; XIN LIU; HAIRONG ZHENG; DONG LIANG; HAIFENG WANG: "Highly accelerated parallel MRI using wave encoding and virtual conjugate coils", MAGNETIC RESONANCE IN MEDICINE, WILEY-LISS, US, vol. 86, no. 3, 15 April 2021 (2021-04-15), US , pages 1345 - 1359, XP071674783, ISSN: 0740-3194, DOI: 10.1002/mrm.28803 *

Similar Documents

Publication Publication Date Title
Lin et al. Artificial intelligence for MR image reconstruction: an overview for clinicians
Plenge et al. Super‐resolution methods in MRI: can they improve the trade‐off between resolution, signal‐to‐noise ratio, and acquisition time?
US7840045B2 (en) Method and system for parallel reconstruction in the K-space domain for application in imaging systems
WO2020019412A1 (en) Magnetic resonance imaging method and apparatus, and device and storage medium
US8890522B2 (en) Accelerated pseudo-random data magnetic resonance imaging system and method
US11131738B2 (en) Multiband, multishot magnetic resonance elastography
WO2022236851A1 (en) Method for reconstructing wrap-around field-of-view magnetic resonance image, computer device, and storage medium
Chen et al. Accelerated two-dimensional cine DENSE cardiovascular magnetic resonance using compressed sensing and parallel imaging
Li et al. SuperMAP: Deep ultrafast MR relaxometry with joint spatiotemporal undersampling
WO2022213666A1 (en) Joint k-space and image-space reconstruction imaging method and device
Levac et al. Fse compensated motion correction for mri using data driven methods
CN109920017A (en) The parallel MR imaging reconstructing method of the full variation Lp pseudonorm of joint from consistency based on feature vector
Zaid Alkilani et al. FD‐Net: An unsupervised deep forward‐distortion model for susceptibility artifact correction in EPI
WO2023108484A1 (en) Magnetic resonance imaging method based on wave gradient coding field and deep learning model
CN117495992A (en) Reconstruction method of undersampled nuclear magnetic resonance image
KR102163337B1 (en) Method for accelerating multiple-acquisition magnetic resonance imaging by varying undersampling-dimension and device for the same
CN113298901B (en) Method for reconstructing magnetic resonance image in convoluted field of view, computer device and storage medium
US20180286088A1 (en) Iterative reconstruction of quantitative magnetic resonance images
CN113506258B (en) Under-sampling lung gas MRI reconstruction method for multitask complex value deep learning
WO2022027419A1 (en) Magnetic resonance imaging method and apparatus, and computer storage medium
JP6618786B2 (en) Magnetic resonance imaging apparatus and image processing apparatus
CN114565690A (en) Magnetic resonance imaging method based on wave gradient coding field and deep learning model
WO2024103414A1 (en) Method and apparatus for reconstructing magnetic resonance image
WO2023092319A1 (en) Magnetic resonance multi contrast-ratio parameter imaging method based on wave-type gradient
Lee et al. Accelerated 3D myelin water imaging using joint spatio‐temporal reconstruction

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 21967621

Country of ref document: EP

Kind code of ref document: A1