[go: up one dir, main page]

CN118096920B - Abdominal organ imaging method and device based on deep learning and coil sensitivity - Google Patents

Abdominal organ imaging method and device based on deep learning and coil sensitivity Download PDF

Info

Publication number
CN118096920B
CN118096920B CN202410260672.XA CN202410260672A CN118096920B CN 118096920 B CN118096920 B CN 118096920B CN 202410260672 A CN202410260672 A CN 202410260672A CN 118096920 B CN118096920 B CN 118096920B
Authority
CN
China
Prior art keywords
image
resolution
diffusion
deep learning
reconstruction
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202410260672.XA
Other languages
Chinese (zh)
Other versions
CN118096920A (en
Inventor
王鹤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fudan University
Original Assignee
Fudan University
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 Fudan University filed Critical Fudan University
Priority to CN202410260672.XA priority Critical patent/CN118096920B/en
Publication of CN118096920A publication Critical patent/CN118096920A/en
Application granted granted Critical
Publication of CN118096920B publication Critical patent/CN118096920B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/008Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • G06N3/0455Auto-encoder networks; Encoder-decoder networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/0464Convolutional networks [CNN, ConvNet]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/09Supervised learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10088Magnetic resonance imaging [MRI]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20081Training; Learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20084Artificial neural networks [ANN]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明公开了基于深度学习和线圈灵敏度的腹部器官成像方法及装置,属于磁共振成像领域,步骤包括:采集高分辨弥散加权数据;通过传统高分辨率弥散重建方法进行图像重建;选择质量高的图像;生成与被选择图像对应的被伪影污染的高分辨率图像;对生成的成对图像输入深度学习的网络进行泛化和训练;将基于传统高分辨率弥散重建的腹部器官多次激发弥散加权图像输入S5构建的网络模型,得到高信噪比、高分辨率及残留伪影少的腹部器官高分辨弥散加权图像;还公开了成像装置。本发明采用上述一种基于深度学习和线圈灵敏度的腹部器官成像方法及装置,提高了重建后的高分辨弥散加权图像质量,为临床诊断提供更精细的解剖信息。

The present invention discloses an abdominal organ imaging method and device based on deep learning and coil sensitivity, which belongs to the field of magnetic resonance imaging, and the steps include: collecting high-resolution diffusion-weighted data; reconstructing images by traditional high-resolution diffusion reconstruction methods; selecting high-quality images; generating high-resolution images contaminated by artifacts corresponding to the selected images; inputting the generated paired images into a deep learning network for generalization and training; inputting the abdominal organ multiple-excitation diffusion-weighted images based on traditional high-resolution diffusion reconstruction into a network model constructed by S5, and obtaining high-resolution diffusion-weighted images of abdominal organs with high signal-to-noise ratio, high resolution and few residual artifacts; and an imaging device is also disclosed. The present invention adopts the above-mentioned abdominal organ imaging method and device based on deep learning and coil sensitivity, improves the quality of the reconstructed high-resolution diffusion-weighted images, and provides more detailed anatomical information for clinical diagnosis.

Description

基于深度学习和线圈灵敏度的腹部器官成像方法及装置Abdominal organ imaging method and device based on deep learning and coil sensitivity

技术领域Technical Field

本发明涉及磁共振成像技术领域,尤其是涉及一种基于深度学习和线圈灵敏度的腹部器官成像方法及装置。The present invention relates to the technical field of magnetic resonance imaging, and in particular to an abdominal organ imaging method and device based on deep learning and coil sensitivity.

背景技术Background Art

临床上常规的弥散加权图像采集方法是通过单次激发平面回波成像(single-shot Echo Planar Imaging,SSH-EPI)实现,但其图像空间分辨率低,导致重建后图像几何失真严重,图像清晰度低。因此目前的改进策略是通过多次激发(沿相位编码方向或者频率编码方向两种)实现,如多次激发平面回波弥散加权成像(Multi-shot Echo PlanarImaging Diffusion Weighted Imaging,MSH-EPIDWI)技术,其通过完成较高的空间分辨率采集,从而重建出较少的几何变形和高分辨率图像,为临床诊断提供更精细的解剖信息。The conventional diffusion-weighted image acquisition method in clinical practice is achieved through single-shot echo planar imaging (SSH-EPI), but its image spatial resolution is low, resulting in severe geometric distortion and low image clarity after reconstruction. Therefore, the current improvement strategy is to achieve it through multiple excitations (along the phase encoding direction or the frequency encoding direction), such as the multi-shot echo planar imaging diffusion weighted imaging (MSH-EPIDWI) technology, which completes the acquisition with higher spatial resolution, thereby reconstructing less geometric deformation and high-resolution images, providing more detailed anatomical information for clinical diagnosis.

但是,对多次激发技术而言,不同次激发之间会存在线性和非线性相位变化,扩散梯度则放大了这种相位差异进而造成多次激发的图像重建中的混叠伪影,而腹部器官相比于头部和四肢等其他部位,在成像中不可避免的会由于呼吸、心跳等运动的存在加剧这种相位差异,从而引入更多的运动伪影,使得腹部器官高分辨弥散成像极具挑战。However, for multiple excitation technology, there will be linear and nonlinear phase changes between different excitations, and the diffusion gradient will amplify this phase difference and thus cause aliasing artifacts in the image reconstruction of multiple excitations. Compared with other parts such as the head and limbs, the abdominal organs will inevitably have this phase difference aggravated during imaging due to the presence of movements such as breathing and heartbeat, thereby introducing more motion artifacts, making high-resolution diffusion imaging of abdominal organs extremely challenging.

目前传统的高分辨率弥散重建方法,如MUSE及一系列的变种算法:POCS-MUSE、self-feeding MUSE和DL-MUSE都是针对脑部等运动很微小部位的高分辨率弥散成像的重建方法,但在运动剧烈且复杂的腹部器官的应用中,效果并不理想,仍然存在非常严重的图像伪影。At present, traditional high-resolution diffusion reconstruction methods, such as MUSE and a series of variant algorithms: POCS-MUSE, self-feeding MUSE and DL-MUSE are all reconstruction methods for high-resolution diffusion imaging of parts with very small movements such as the brain. However, in applications of abdominal organs with violent and complex movements, the effect is not ideal and there are still very serious image artifacts.

发明内容Summary of the invention

本发明的目的是提供基于深度学习和线圈灵敏度的腹部器官成像方法及装置,通过对成对的高质量伪影较少的图像和对应的的带伪影的图像进行深度学习网络训练和模型构建,提高重建后的高分辨弥散加权图像质量,为临床诊断提供更精细的解剖信息。The purpose of the present invention is to provide an abdominal organ imaging method and device based on deep learning and coil sensitivity. By performing deep learning network training and model construction on paired high-quality images with fewer artifacts and corresponding images with artifacts, the quality of reconstructed high-resolution diffusion-weighted images is improved, thereby providing more detailed anatomical information for clinical diagnosis.

为实现上述目的,本发明提供了一种基于深度学习和线圈灵敏度的腹部器官成像方法,步骤包括:To achieve the above object, the present invention provides an abdominal organ imaging method based on deep learning and coil sensitivity, the steps comprising:

S1、采集进行多次激发的腹部器官高分辨弥散加权数据;S1, collecting high-resolution diffusion-weighted data of abdominal organs with multiple excitations;

S2、通过传统高分辨率弥散重建方法对采集的数据进行图像重建;S2, reconstructing the collected data using a traditional high-resolution diffusion reconstruction method;

S3、选择重建后伪影少、质量高的图像;S3, select images with fewer artifacts and higher quality after reconstruction;

S4、使用步骤S3选择的图像和与选择图像对应的线圈灵敏度图,进行腹部器官运动伪影仿真和图像重建自身过程的混叠伪影仿真,生成对应的被伪影污染的高分辨率弥散图像;S4, using the image selected in step S3 and the coil sensitivity map corresponding to the selected image, performing motion artifact simulation of abdominal organs and aliasing artifact simulation of the image reconstruction process itself, and generating a corresponding high-resolution diffusion image contaminated by artifacts;

S5、将步骤S3和步骤S4生成的成对图像输入深度学习网络进行泛化和训练,构建高分辨弥散图像伪影矫正的网络模型;S5, inputting the paired images generated in step S3 and step S4 into the deep learning network for generalization and training, and constructing a network model for high-resolution diffuse image artifact correction;

S6、将基于传统高分辨率弥散重建的腹部器官多次激发弥散加权图像输入步骤S5构建的网络模型,得到高信噪比、高分辨率及残留伪影少的腹部器官高分辨弥散加权图像。S6. Inputting the multi-excitation diffusion-weighted image of the abdominal organs based on traditional high-resolution diffusion reconstruction into the network model constructed in step S5, to obtain a high-resolution diffusion-weighted image of the abdominal organs with high signal-to-noise ratio, high resolution and few residual artifacts.

优选的,所述步骤S3具体包括:Preferably, the step S3 specifically includes:

S31、对重建后的图像进行分割;S31, segmenting the reconstructed image;

S32、基于预训练的分割模型计算图像评价指标;S32, calculating image evaluation indicators based on the pre-trained segmentation model;

S33、按照图像评价指标的某个阈值进行图像筛选。S33, performing image screening according to a certain threshold of the image evaluation index.

优选的,所述图像筛选基于伪影信号比GSR、信噪比SNR、对比度噪声比CNR中一种或多种类似图像评价指标的组合作为图像质量筛选方法。Preferably, the image screening is based on a combination of one or more similar image evaluation indicators such as artifact signal ratio GSR, signal-to-noise ratio SNR, and contrast-to-noise ratio CNR as an image quality screening method.

优选的,所述步骤S6中基于传统高分辨率弥散重建的腹部器官多次激发弥散加权图像存在的问题包括但不限于存在严重伪影的、信噪比低、因运动而出现信号缺失和图像模糊。Preferably, problems existing in the abdominal organ multiple excitation diffusion weighted image based on traditional high-resolution diffusion reconstruction in step S6 include but are not limited to severe artifacts, low signal-to-noise ratio, signal loss due to motion, and image blur.

一种基于深度学习和线圈灵敏度的腹部器官成像装置,包括图像准备模块、运动仿真模块和图像矫正模块;An abdominal organ imaging device based on deep learning and coil sensitivity, comprising an image preparation module, a motion simulation module and an image correction module;

图像准备模块,用于对被测目标进行多次激发弥散数据的磁共振信号的激发、采集和重建;An image preparation module, used for exciting, collecting and reconstructing magnetic resonance signals of the target object for multiple excitation diffusion data;

运动仿真模块,用于通过对被测目标的各种运动情况进行仿真;Motion simulation module, used to simulate various motion conditions of the target under test;

图像矫正模块,用于根据图像准备模块和运动仿真模块提供的信息,使用矫正算法或者网络模型来对高分辨弥散加权图像进行图像伪影矫正。The image correction module is used to correct image artifacts of the high-resolution diffusion-weighted image using a correction algorithm or a network model according to the information provided by the image preparation module and the motion simulation module.

优选的,所述图像准备模块包括高分辨弥散数据重建单元和图像筛选单元;Preferably, the image preparation module includes a high-resolution diffusion data reconstruction unit and an image screening unit;

高分辨弥散数据重建单元,用于重建高分辨弥散加权图像;A high-resolution diffusion data reconstruction unit, used for reconstructing a high-resolution diffusion-weighted image;

图像筛选单元,用于筛选高质量的高分辨弥散加权图像。The image screening unit is used to screen high-quality high-resolution diffusion-weighted images.

优选的,所述运动仿真模块包括物体运动仿真单元和重建运动仿真单元;Preferably, the motion simulation module includes an object motion simulation unit and a reconstruction motion simulation unit;

物体运动仿真单元,用于利用仿真运动来直接模拟被测物体运动对采集信号的影响;An object motion simulation unit, used to directly simulate the influence of the motion of the measured object on the collected signal by using simulated motion;

重建运动仿真单元,用于利用仿真运动来模拟被测腹部器官的运动对高分辨率弥散重建方法本身的影响。The reconstruction motion simulation unit is used to use simulated motion to simulate the influence of the motion of the abdominal organ under test on the high-resolution diffusion reconstruction method itself.

优选的,所述图像矫正模块包括对图像本身幅值、伪影、相位进行计算矫正的模块。Preferably, the image correction module includes a module for calculating and correcting the amplitude, artifacts and phase of the image itself.

因此,本发明采用上述结构的基于深度学习和线圈灵敏度的腹部器官成像方法及装置,具有以下有益效果:Therefore, the present invention adopts the above-mentioned structure based on deep learning and coil sensitivity of the abdominal organ imaging method and device, which has the following beneficial effects:

(1)使用深度学习方式代替传统高分辨率弥散重建方法基于数学模型的近似求解方式,能够更加准确有效的消除挑战性更高的腹部器官的高分辨弥散加权图像重建中的伪影,且重建结果质量更高,同时速度更快;(1) Using deep learning to replace the approximate solution method based on mathematical models in traditional high-resolution diffusion reconstruction methods can more accurately and effectively eliminate artifacts in the high-resolution diffusion-weighted image reconstruction of more challenging abdominal organs, and the reconstruction results are of higher quality and faster speed;

(2)传统高分辨率弥散成像重建方法重建的图像质量和加速倍数会受到g-factor的制约,因此在多次弥散加权序列应用中,其激发次数有限,限制了分辨率的提高,本发明提出的基于深度学习的方法则不受这种模型计算的制约,可以对更高激发次数、更高分辨率的应用有效;(2) The image quality and acceleration factor of the traditional high-resolution diffusion imaging reconstruction method are restricted by the g-factor. Therefore, in the application of multiple diffusion-weighted sequences, the number of excitations is limited, which limits the improvement of resolution. The deep learning-based method proposed in the present invention is not restricted by this model calculation and can be effective for applications with higher excitation times and higher resolution.

(3)成熟稳健的网络模型将大大缩减腹部器官的重建时间并提高其重建精度。(3) A mature and robust network model will greatly reduce the reconstruction time of abdominal organs and improve their reconstruction accuracy.

下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。The technical solution of the present invention is further described in detail below through the accompanying drawings and embodiments.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

图1为本发明实施例的肝脏成像整体流程图和具体网络结构图;FIG1 is an overall flow chart and a specific network structure diagram of liver imaging according to an embodiment of the present invention;

图2为本发明实施例中一例脂质不良的血管平滑肌脂肪瘤患者的多个b值(0,600和800s/mm2)的不同重建方法的结果对比示意图;FIG2 is a schematic diagram showing the comparison of the results of different reconstruction methods at multiple b values (0, 600 and 800 s/mm 2 ) of an angiomyolipoma patient with dyslipidemia according to an embodiment of the present invention;

图3为本发明实施例中两个健康受试的多个b值(0,600和800s/mm2)的不同重建方法的代表性重建结果对比;FIG3 is a comparison of representative reconstruction results of different reconstruction methods at multiple b values (0, 600 and 800 s/mm 2 ) of two healthy subjects in an embodiment of the present invention;

图4为本发明实施例基于深度学习和线圈灵敏度的腹部器官成像装置结构示意图。FIG4 is a schematic diagram of the structure of an abdominal organ imaging device based on deep learning and coil sensitivity according to an embodiment of the present invention.

具体实施方式DETAILED DESCRIPTION

实施例Example

为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments.

参照图1-3,一种基于深度学习和线圈灵敏度的腹部器官成像方法,包括步骤:1-3, a method for abdominal organ imaging based on deep learning and coil sensitivity includes the steps of:

S1、采集50例健康受试进行多次激发的腹部器官高分辨弥散加权数据,其中每个受试都采集了单次和四次激发数据。S1. High-resolution diffusion-weighted data of abdominal organs were collected from 50 healthy subjects with multiple excitations, and single and four-shot data were collected for each subject.

S2、通过传统高分辨率弥散重建方法对采集的数据进行图像重建,重建后,图像出现很多的混叠或者运动伪影,无法用于临床的诊断。S2. The collected data is reconstructed using the traditional high-resolution diffusion reconstruction method. After reconstruction, the image has many aliasing or motion artifacts and cannot be used for clinical diagnosis.

S3、选择重建后伪影少、质量高的图像作为步骤S4的准备图像。具体为,每个受试的肝脏弥散加权图像数据,以包含有3个b值(以0,600和800s/mm2为例)和沿着头角方向覆盖全肝扫描20层为例,针对重建后的弥散加权图像数据,选取其b=0s/mm2无弥散加权图像为例进行图像分割获取肝脏掩模,如图1所示,使用的是预训练好的nnUNet模型为例进行的图形分割;同时,针对这批健康受试的所有不同b值加权、不同层的肝脏弥散加权数据进行混合,使用预训练好的nnUNet模型计算全部高分辨率弥散加权图像的GSR评价指标,并按照肝脏区域的GSR评价指标的20%阈值为例进行筛选。S3, select the image with few artifacts and high quality after reconstruction as the preparation image for step S4. Specifically, for each subject's liver diffusion-weighted image data, take the image data containing 3 b values (taking 0, 600 and 800s/ mm2 as examples) and covering 20 layers of the whole liver scan along the head angle direction as an example, for the reconstructed diffusion-weighted image data, select its b=0s/ mm2 non-diffusion-weighted image as an example to perform image segmentation to obtain the liver mask, as shown in Figure 1, using the pre-trained nnUNet model as an example for graphic segmentation; at the same time, for this batch of healthy subjects, all liver diffusion-weighted data with different b value weightings and different layers are mixed, and the pre-trained nnUNet model is used to calculate the GSR evaluation index of all high-resolution diffusion-weighted images, and the 20% threshold of the GSR evaluation index of the liver area is used as an example for screening.

S4、使用步骤S3中得到的图像和与选择图像对应的线圈灵敏度图,进行腹部器官运动仿真和图像重建过程的混叠伪影仿真,生成对应的被伪影污染的高分辨率图像,具体为:S4, using the image obtained in step S3 and the coil sensitivity map corresponding to the selected image, perform abdominal organ motion simulation and aliasing artifact simulation in the image reconstruction process to generate a corresponding high-resolution image contaminated by artifacts, specifically:

针对腹部器官运动仿真,如图1的肝脏运动仿真部分所示,考虑到数据采集期间腹部器官运动(呼吸和心跳等)的非均匀性,使用基于几何空间的各种刚体和非刚体运动仿真(包括但不限于平移、旋转、变形、加速等)。同时为了确保最终图像重建的高质量,考虑到在这种运动下引起的磁场对重建图像的不均匀性的影响,加入了与腹部图像对应的线圈灵敏度图(CSM)的运动情况。最终产生对运动影响充分仿真的被伪影污染的高分辨弥散加权图像。本实施例中,采用了三种运动仿真:分别为图像沿x或y轴的位移小于3个像素,图像沿x或y轴旋转小于3度,及具有弹性高斯特性的变形仿真。针对图像采集期间腹部呼吸和心脏运动对场不均匀性的影响,对线圈灵敏度图CSM也采用了三种运动仿真:分别为CSM沿x或y轴的位移小于3个像素,CSM沿x或y轴旋转小于3度,CSM的弹性高斯特性变形仿真。考虑到数据采集过程中,不同次激发间实际发生的运动不一致的情况,对四次激发的数据图像的重建,每次激发的图像和CSM的运动都随机选取了三种运动中的两种进行组合,以确保更真实地呈现生理动态特性。For the simulation of abdominal organ motion, as shown in the liver motion simulation part of FIG1, various rigid and non-rigid motion simulations based on geometric space (including but not limited to translation, rotation, deformation, acceleration, etc.) are used in consideration of the non-uniformity of abdominal organ motion (respiration and heartbeat, etc.) during data acquisition. At the same time, in order to ensure the high quality of the final image reconstruction, the motion of the coil sensitivity map (CSM) corresponding to the abdominal image is added in consideration of the influence of the magnetic field caused by such motion on the non-uniformity of the reconstructed image. Finally, a high-resolution diffusion-weighted image contaminated by artifacts that fully simulates the influence of motion is generated. In this embodiment, three types of motion simulation are used: the displacement of the image along the x or y axis is less than 3 pixels, the rotation of the image along the x or y axis is less than 3 degrees, and the deformation simulation with elastic Gaussian characteristics. In view of the influence of abdominal breathing and cardiac motion on field non-uniformity during image acquisition, three types of motion simulation are also used for the coil sensitivity map CSM: the displacement of the CSM along the x or y axis is less than 3 pixels, the rotation of the CSM along the x or y axis is less than 3 degrees, and the elastic Gaussian deformation simulation of the CSM. Taking into account the inconsistency of the actual motion between different excitations during data acquisition, for the reconstruction of the data images of the four excitations, two of the three motions were randomly selected for each excitation image and CSM motion to ensure a more realistic presentation of the physiological dynamic characteristics.

S5、将步骤S3和步骤S4生成的50例健康受试成对(高质量和运动污染)图像输入深度学习网络进行泛化和训练,构建高分辨弥散图像伪影矫正的网络模型。S5. Input the paired images (high quality and motion contaminated) of 50 healthy subjects generated in steps S3 and S4 into the deep learning network for generalization and training, and construct a network model for high-resolution diffusion image artifact correction.

具体的,伪影矫正网络整体结构可采用4级编码器-解码器监督网络通过Pytorch实现、训练在NVIDIA Tesla P100×4GPU(每个处理器有16GB内存)上执行为例,也可采用本领域技术人员熟悉的其他编程或训练平台实现。卷积层权重以高斯随机分布初始化为例,使用AdamW优化器去伪影网络为例,也可采用其他本领域技术人员熟悉的深度学习优化器。训练过程包括利用L1损失进行150000次迭代为例,从初始学习率3×10-4开始为例,通过余弦退火逐渐降低到10-6为例。所有的权值可使用凯明法初始化或其他本领域技术人员熟悉的初始化方法,并在高质量图像和运动模拟图像之间的L1损失监督下进行更新。整个训练时间约为46小时,根据损失度量选择模型。Specifically, the overall structure of the artifact correction network can be implemented by Pytorch using a 4-level encoder-decoder supervision network, and the training can be performed on an NVIDIA Tesla P100×4 GPU (each processor has 16GB of memory) as an example, or it can be implemented using other programming or training platforms familiar to those skilled in the art. The weights of the convolutional layer are initialized using a Gaussian random distribution as an example, and the AdamW optimizer is used to remove the artifact network as an example, or other deep learning optimizers familiar to those skilled in the art can be used. The training process includes using L1 loss for 150,000 iterations as an example, starting from an initial learning rate of 3× 10-4 as an example, and gradually reducing it to 10-6 through cosine annealing as an example. All weights can be initialized using the Kaiming method or other initialization methods familiar to those skilled in the art, and updated under the supervision of L1 loss between high-quality images and motion simulation images. The entire training time is about 46 hours, and the model is selected based on the loss metric.

S6、将基于传统高分辨率弥散重建的存在严重伪影的、信噪比低、因运动而出现信号缺失和图像模糊等问题的腹部器官多次激发弥散加权图像输入步骤S5中训练好的网络模型,并进行训练和泛化性能优化,得到高信噪比、高分辨率及残留伪影少的腹部器官高分辨弥散加权图像。S6. Input the multiple-excitation diffusion-weighted images of abdominal organs based on traditional high-resolution diffusion reconstruction, which have serious artifacts, low signal-to-noise ratio, signal loss and image blur due to movement, etc., into the network model trained in step S5, and perform training and generalization performance optimization to obtain high-resolution diffusion-weighted images of abdominal organs with high signal-to-noise ratio, high resolution and few residual artifacts.

具体的,使用进行四次激发的高分辨弥散加权图像数据对第一批10个病人和10个健康受试为例进行肝脏弥散加权数据的多通道数据的测试应用,使用传统高分辨肝脏弥散图像进行重建,图像会带有很多伪影无法用于临床诊断。将这新的一批肝脏的数据和第一批50个健康受试进行高质量图像筛选后剩余的80%数据作为深度学习网络模型的参数输入,进行泛化和训练,最终得到可用于临床诊断的、高分辨率、高精度的极少伪影的肝脏高分辨弥散加权图像。Specifically, the high-resolution diffusion-weighted image data of the first batch of 10 patients and 10 healthy subjects were used as an example to test the multi-channel data of liver diffusion-weighted data. The traditional high-resolution liver diffusion image was used for reconstruction, and the image would have many artifacts and could not be used for clinical diagnosis. The data of this new batch of livers and the remaining 80% of the data after high-quality image screening of the first batch of 50 healthy subjects were used as parameter inputs of the deep learning network model for generalization and training, and finally high-resolution, high-precision, and artifact-free liver high-resolution diffusion-weighted images that can be used for clinical diagnosis were obtained.

使用进行四次激发的高分辨弥散加权图像数据对第二批10个健康和10个病人受试和第一批筛选后的80%健康受试对应的数据进行传统高分辨率重建,经训练好的网络模型得到肝脏高分辨弥散加权图像。The high-resolution diffusion-weighted image data with four excitations were used to perform traditional high-resolution reconstruction on the data corresponding to the second batch of 10 healthy and 10 patient subjects and 80% of the healthy subjects after the first screening, and the high-resolution diffusion-weighted image of the liver was obtained through the trained network model.

将得到肝脏高分辨弥散加权图像和低分辨率图像进行定性和定量对比分析。通过观察,发现按本发明提供的方法得到的图像质量上有显著的提高和明显的改进(随着b值的增大,由于呼吸和心脏运动引起的伪影和图像模糊,特别是肝脏左叶等部位的信号缺失等情况的发生都明显大大的降低和消除)。The high-resolution diffusion-weighted image and low-resolution image of the liver were subjected to qualitative and quantitative comparative analysis. Through observation, it was found that the image quality obtained by the method provided by the present invention was significantly improved (with the increase of b value, the occurrence of artifacts and image blurring caused by breathing and cardiac movement, especially the signal loss in the left lobe of the liver, etc., was significantly reduced and eliminated).

图2和图3分别展示了在一个病人和两个健康受试身上的对比结果。图2为一例脂质不良的血管平滑肌脂肪瘤患者基于传统高分辨弥散重建方法(C-E)、临床常用的低分辨率重建方法(G-I)和使用本发明的运动伪影矫正的高分辨弥散重建方法(K-M)在多个b值(0,600和800s/mm2)的不同重建方法的结果对比示意图。其中参考图像为临床常规扫描的T2加权(A)、T1加权(B)和DCE图像(F为动脉期,J为门脉期,N为延迟期)。实心的红色箭头指向的是位于肝脏边缘的尺寸较小的病灶,空心的红色箭头指向的是明显的运动伪影,空心的蓝色箭头指向的是因存在运动伪影而引起的信号缺失和信号强度错误的地方。图3中,为使用运动伪影矫正网络的高分辨弥散重建方法(A的第一行和B的第一行)和临床常用的低分辨率弥散重建方法,空心的红色箭头指向的是传统高分辨率方法因存在腹部、特别是肝脏左叶靠近心脏区域的运动伪影而引起的信号缺失和信号强度错误的地方。通过重建结果对比可以看出,本发明提供的方法,为临床肝脏疾病的诊断提供卓越的图像质量,可以准确描绘小囊肿、恶性和良性肿瘤,即使是严重的铁超载和肝脏对比度差,其增强的诊断能力也扩展到位于肝脏左叶、血管附近或受心脏和血管脉动影响的小尺寸病变或感兴趣区域。Figures 2 and 3 show the comparison results on a patient and two healthy subjects, respectively. Figure 2 is a schematic diagram showing the comparison of the results of different reconstruction methods at multiple b values (0, 600 and 800 s/mm2) based on the traditional high-resolution diffusion reconstruction method (C-E), the clinically commonly used low-resolution reconstruction method (G-I) and the high-resolution diffusion reconstruction method (K-M) using the motion artifact correction of the present invention for a patient with lipid dyslipidemia. The reference images are T2-weighted (A), T1-weighted (B) and DCE images (F is the arterial phase, J is the portal phase, and N is the delayed phase) of routine clinical scans. The solid red arrow points to a smaller lesion at the edge of the liver, the hollow red arrow points to an obvious motion artifact, and the hollow blue arrow points to a place where the signal is missing and the signal intensity is wrong due to the presence of motion artifacts. In Figure 3, the high-resolution diffusion reconstruction method using the motion artifact correction network (the first row of A and the first row of B) and the low-resolution diffusion reconstruction method commonly used in clinical practice are shown. The hollow red arrows point to the areas where the traditional high-resolution method causes signal loss and signal intensity errors due to motion artifacts in the abdomen, especially in the left lobe of the liver near the heart. By comparing the reconstruction results, it can be seen that the method provided by the present invention provides excellent image quality for the diagnosis of clinical liver diseases, and can accurately depict small cysts, malignant and benign tumors, even in severe iron overload and poor liver contrast. Its enhanced diagnostic capability is also extended to small-sized lesions or regions of interest located in the left lobe of the liver, near blood vessels, or affected by heart and blood vessel pulsation.

参照图4,一种基于深度学习和线圈灵敏度的腹部器官成像装置,包括:Referring to FIG. 4 , an abdominal organ imaging device based on deep learning and coil sensitivity includes:

图像准备模块,用于对被测目标进行多次激发弥散数据的磁共振信号的激发、采集和重建,同时进行高质量的高分辨率弥散数据的准备。图像准备模块包括高分辨弥散数据重建单元,用于重建传统高分辨肝脏器官弥散加权图像;图像筛选单元,使用基于伪影信号比GSR的筛选量纲来进行20%阈值为例的高质量高分辨弥散图像为例的筛选。The image preparation module is used to excite, collect and reconstruct the magnetic resonance signals of the target object for multiple excitation diffusion data, and prepare high-quality high-resolution diffusion data at the same time. The image preparation module includes a high-resolution diffusion data reconstruction unit, which is used to reconstruct traditional high-resolution liver organ diffusion weighted images; an image screening unit, which uses a screening dimension based on the artifact signal ratio GSR to screen high-quality high-resolution diffusion images with a threshold of 20% as an example.

运动仿真模块,用于通过对被测目标的各种运动情况进行仿真,从而达到模拟磁共振信号采集过程中的被测物体运动对实际图像重建的影响。运动仿真模块包括物体运动仿真单元,用于利用图像本身包括但不限于平移、旋转、变形或加速等运动仿真来直接模拟肝脏受到的呼吸、心跳等运动对数据采集的信号的影响;重建运动仿真单元,用于利用线圈敏感度CSM包括但不限于平移、旋转、变形或加速等运动仿真形式来模拟上述刚体或非刚体运动对重建过程本身的影响。The motion simulation module is used to simulate various motion conditions of the target to be measured, so as to simulate the influence of the motion of the measured object during the magnetic resonance signal acquisition process on the actual image reconstruction. The motion simulation module includes an object motion simulation unit, which is used to directly simulate the influence of the liver's breathing, heartbeat and other movements on the data acquisition signal by using the image itself, including but not limited to translation, rotation, deformation or acceleration and other motion simulation forms; a reconstruction motion simulation unit, which is used to simulate the influence of the above rigid or non-rigid body motion on the reconstruction process itself by using the coil sensitivity CSM, including but not limited to translation, rotation, deformation or acceleration and other motion simulation forms.

图像矫正模块,用于根据图像准备模块和运动仿真模块提供的信息,使用矫正算法或者网络模型来对高分辨弥散加权图像进行图像伪影矫正,以提高得到高质量图像。图像矫正模块包括对图像本身幅值、伪影、相位进行计算矫正的模块。The image correction module is used to correct image artifacts of high-resolution diffusion-weighted images using correction algorithms or network models based on the information provided by the image preparation module and the motion simulation module to improve the quality of the image. The image correction module includes a module for calculating and correcting the amplitude, artifacts, and phase of the image itself.

在本发明的描述中,需要理解的是,术语“x轴”、“y轴”、“平移”、“旋转”、“变形”、等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的方法、装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。In the description of the present invention, it is necessary to understand that the orientations or positional relationships indicated by terms such as "x-axis", "y-axis", "translation", "rotation", and "deformation" are based on the orientations or positional relationships shown in the accompanying drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the methods, devices, or elements referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore should not be understood as limiting the present invention.

最后应说明的是:以上实施例仅用以说明本发明的技术方案而非对其进行限制,尽管参照较佳实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对本发明的技术方案进行修改或者等同替换,而这些修改或者等同替换亦不能使修改后的技术方案脱离本发明技术方案的精神和范围。Finally, it should be noted that the above embodiments are only used to illustrate the technical solution of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that they can still modify or replace the technical solution of the present invention with equivalents, and these modifications or equivalent replacements cannot cause the modified technical solution to deviate from the spirit and scope of the technical solution of the present invention.

Claims (7)

1.基于深度学习和线圈灵敏度的腹部器官成像方法,其特征在于,步骤包括:1. An abdominal organ imaging method based on deep learning and coil sensitivity, characterized in that the steps include: S1、采集进行多次激发的腹部器官高分辨弥散加权数据;S1, collecting high-resolution diffusion-weighted data of abdominal organs with multiple excitations; S2、通过传统高分辨率弥散重建方法对采集的数据进行图像重建;S2, reconstructing the collected data using a traditional high-resolution diffusion reconstruction method; S3、选择重建后伪影少、质量高的图像;S3, select images with fewer artifacts and higher quality after reconstruction; S4、使用步骤S3选择的图像和与选择图像对应的线圈灵敏度图,进行腹部器官运动伪影仿真和图像重建自身过程的混叠伪影仿真,生成对应的被伪影污染的高分辨率弥散图像;S4, using the image selected in step S3 and the coil sensitivity map corresponding to the selected image, performing motion artifact simulation of abdominal organs and aliasing artifact simulation of the image reconstruction process itself, and generating a corresponding high-resolution diffusion image contaminated by artifacts; S5、将步骤S3和步骤S4生成的成对图像输入深度学习网络进行泛化和训练,构建高分辨弥散图像伪影矫正的网络模型;S5, inputting the paired images generated in step S3 and step S4 into the deep learning network for generalization and training, and constructing a network model for high-resolution diffuse image artifact correction; S6、将基于传统高分辨率弥散重建的腹部器官多次激发弥散加权图像输入步骤S5构建的网络模型,得到高信噪比、高分辨率及残留伪影少的腹部器官高分辨弥散加权图像。S6. Inputting the multi-excitation diffusion-weighted image of the abdominal organs based on traditional high-resolution diffusion reconstruction into the network model constructed in step S5, to obtain a high-resolution diffusion-weighted image of the abdominal organs with high signal-to-noise ratio, high resolution and few residual artifacts. 2.根据权利要求1所述的基于深度学习和线圈灵敏度的腹部器官成像方法,其特征在于,所述步骤S3具体包括:2. The abdominal organ imaging method based on deep learning and coil sensitivity according to claim 1, characterized in that step S3 specifically comprises: S31、对重建后的图像进行分割;S31, segmenting the reconstructed image; S32、基于预训练的分割模型计算图像评价指标;S32, calculating image evaluation indicators based on the pre-trained segmentation model; S33、按照图像评价指标的某个阈值进行图像筛选。S33, performing image screening according to a certain threshold of the image evaluation index. 3.根据权利要求2所述的基于深度学习和线圈灵敏度的腹部器官成像方法,其特征在于:所述图像筛选基于伪影信号比GSR、信噪比SNR、对比度噪声比CNR中一种或多种组合作为图像质量筛选方法。3. The abdominal organ imaging method based on deep learning and coil sensitivity according to claim 2 is characterized in that: the image screening is based on one or more combinations of artifact signal ratio GSR, signal-to-noise ratio SNR, and contrast-to-noise ratio CNR as an image quality screening method. 4.基于深度学习和线圈灵敏度的腹部器官成像装置,应用上述权利要求1-3任一项所述的基于深度学习和线圈灵敏度的腹部器官成像方法,其特征在于:包括图像准备模块、运动仿真模块和图像矫正模块;4. An abdominal organ imaging device based on deep learning and coil sensitivity, applying the abdominal organ imaging method based on deep learning and coil sensitivity as described in any one of claims 1 to 3, characterized in that it includes an image preparation module, a motion simulation module and an image correction module; 图像准备模块,用于对被测目标进行多次激发弥散数据的磁共振信号的激发、采集和重建;An image preparation module, used for exciting, collecting and reconstructing magnetic resonance signals of the target object for multiple excitation diffusion data; 运动仿真模块,用于通过对被测目标的各种运动情况进行仿真;Motion simulation module, used to simulate various motion conditions of the target under test; 图像矫正模块,用于根据图像准备模块和运动仿真模块提供的信息,使用矫正算法或者网络模型来对高分辨弥散加权图像进行图像伪影矫正。The image correction module is used to correct image artifacts of the high-resolution diffusion-weighted image using a correction algorithm or a network model according to the information provided by the image preparation module and the motion simulation module. 5.根据权利要求4所述的基于深度学习和线圈灵敏度的腹部器官成像装置,其特征在于:所述图像准备模块包括高分辨弥散数据重建单元和图像筛选单元;5. The abdominal organ imaging device based on deep learning and coil sensitivity according to claim 4, characterized in that: the image preparation module includes a high-resolution diffusion data reconstruction unit and an image screening unit; 高分辨弥散数据重建单元,用于重建高分辨弥散加权图像;A high-resolution diffusion data reconstruction unit, used for reconstructing a high-resolution diffusion-weighted image; 图像筛选单元,用于筛选高质量的高分辨弥散加权图像。The image screening unit is used to screen high-quality high-resolution diffusion-weighted images. 6.根据权利要求4所述的基于深度学习和线圈灵敏度的腹部器官成像装置,其特征在于:所述运动仿真模块包括物体运动仿真单元和重建运动仿真单元;6. The abdominal organ imaging device based on deep learning and coil sensitivity according to claim 4, characterized in that: the motion simulation module includes an object motion simulation unit and a reconstruction motion simulation unit; 物体运动仿真单元,用于利用仿真运动来直接模拟被测物体运动对采集信号的影响;An object motion simulation unit, used to directly simulate the influence of the motion of the measured object on the collected signal by using simulated motion; 重建运动仿真单元,用于利用仿真运动来模拟被测腹部器官的运动对高分辨率弥散重建方法本身的影响。The reconstruction motion simulation unit is used to use simulated motion to simulate the influence of the motion of the abdominal organ under test on the high-resolution diffusion reconstruction method itself. 7.根据权利要求4所述的基于深度学习和线圈灵敏度的腹部器官成像装置,其特征在于:所述图像矫正模块包括对图像本身幅值、伪影、相位进行计算矫正的模块。7. The abdominal organ imaging device based on deep learning and coil sensitivity according to claim 4 is characterized in that: the image correction module includes a module for calculating and correcting the amplitude, artifacts, and phase of the image itself.
CN202410260672.XA 2024-03-07 2024-03-07 Abdominal organ imaging method and device based on deep learning and coil sensitivity Active CN118096920B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410260672.XA CN118096920B (en) 2024-03-07 2024-03-07 Abdominal organ imaging method and device based on deep learning and coil sensitivity

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410260672.XA CN118096920B (en) 2024-03-07 2024-03-07 Abdominal organ imaging method and device based on deep learning and coil sensitivity

Publications (2)

Publication Number Publication Date
CN118096920A CN118096920A (en) 2024-05-28
CN118096920B true CN118096920B (en) 2024-09-06

Family

ID=91154663

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410260672.XA Active CN118096920B (en) 2024-03-07 2024-03-07 Abdominal organ imaging method and device based on deep learning and coil sensitivity

Country Status (1)

Country Link
CN (1) CN118096920B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117011409A (en) * 2023-08-10 2023-11-07 厦门大学 Multi-position physical intelligent high-definition diffusion magnetic resonance data generation method

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2697994C2 (en) * 2014-07-03 2019-08-21 Конинклейке Филипс Н.В. Multi-frame magnetic resonance (mr) tomography system and method of operation thereof
US11100621B2 (en) * 2017-10-20 2021-08-24 Imaging Biometrics, Llc Simulated post-contrast T1-weighted magnetic resonance imaging
CN108594146B (en) * 2018-04-24 2020-07-28 朱高杰 Magnetic resonance diffusion weighted imaging method based on deep learning and convex set projection
US11308661B2 (en) * 2018-05-21 2022-04-19 The Board Of Trustees Of The Leland Stanford Junior University Motion robust reconstruction of multi-shot diffusion-weighted images without phase estimation via locally low-rank regularization
CN110895320B (en) * 2019-10-31 2021-12-24 清华大学 Deformation-free single-shot planar echo imaging method and device based on deep learning
CN112051531B (en) * 2020-09-14 2022-10-28 首都医科大学附属北京天坛医院 Multi-excitation navigation-free magnetic resonance diffusion imaging method and device
CN112763958B (en) * 2020-12-10 2022-06-21 复旦大学 A Multi-Excitation Echo Plane Magnetic Resonance Imaging Method Based on Neural Network

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117011409A (en) * 2023-08-10 2023-11-07 厦门大学 Multi-position physical intelligent high-definition diffusion magnetic resonance data generation method

Also Published As

Publication number Publication date
CN118096920A (en) 2024-05-28

Similar Documents

Publication Publication Date Title
Chen et al. Deep learning for image enhancement and correction in magnetic resonance imaging—state-of-the-art and challenges
Lyu et al. Cine cardiac MRI motion artifact reduction using a recurrent neural network
Hoch et al. Advanced MRI of the optic nerve
US11415655B2 (en) Reduced field-of-view perfusion imaging with high spatiotemporal resolution
CN109696647A (en) The K space acquisition method and method for reconstructing of three-dimensional repeatedly excitation Diffusion-Weighted MR Imaging
Jiang et al. Respiratory motion correction in abdominal MRI using a densely connected U-Net with GAN-guided training
US20230380714A1 (en) Method and system for low-field mri denoising with a deep complex-valued convolutional neural network
Rodríguez‐Soto et al. Correction of artifacts induced by B0 inhomogeneities in breast MRI using reduced‐field‐of‐view echo‐planar imaging and enhanced reversed polarity gradient method
CN115359144B (en) Magnetic resonance plane echo imaging artifact simulation method and system
CN118096920B (en) Abdominal organ imaging method and device based on deep learning and coil sensitivity
JP7529892B2 (en) Method and apparatus for acquiring and reconstructing a sequence of volumetric diffusion-weighted magnetic resonance images - Patents.com
CN112489150B (en) Multi-scale sequential training method of deep neural network for rapid MRI
CN118261822A (en) Self-supervised MRI diffusion weighted image denoising method based on voxel replacement
Chen et al. Assessment of rat sciatic nerve using diffusion-tensor imaging with readout-segmented echo planar imaging
CN116430287A (en) A method of lung magnetic resonance dynamic imaging
Baker et al. Rapid 2D 23Na MRI of the calf using a denoising convolutional neural network
CN113341357A (en) vBM3 d-based magnetic resonance image reconstruction method
Tabari et al. Ultrafast brain imaging with deep learning multi-shot EPI: preliminary clinical evaluation
CN114998470B (en) Image imaging method, system and storage medium based on motion correction
CN113066145B (en) Deep learning-based rapid whole-body diffusion weighted imaging method and related equipment
CN117036520A (en) Motion artifact correction method for magnetic resonance image
CN118330529A (en) A fast magnetic resonance imaging and image reconstruction method
Liu Plug-and-Play Methods for Magnetic Resonance Image Reconstruction
Liu et al. Feasibility of the application of deep learning-reconstructed ultra-fast respiratory-triggered T2-weighted imaging at 3 T in liver imaging
Yu et al. Using Diffusion-Tensor Imaging With Readout-Segmented Echo Planar

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant