CN114972145A - 基于深度学习的扩散张量重建方法、系统和电子设备 - Google Patents

基于深度学习的扩散张量重建方法、系统和电子设备 Download PDF

Info

Publication number
CN114972145A
CN114972145A CN202210581140.7A CN202210581140A CN114972145A CN 114972145 A CN114972145 A CN 114972145A CN 202210581140 A CN202210581140 A CN 202210581140A CN 114972145 A CN114972145 A CN 114972145A
Authority
CN
China
Prior art keywords
layer
images
image
network input
input image
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.)
Granted
Application number
CN202210581140.7A
Other languages
English (en)
Other versions
CN114972145B (zh
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.)
Xiamen University
Original Assignee
Xiamen 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 Xiamen University filed Critical Xiamen University
Priority to CN202210581140.7A priority Critical patent/CN114972145B/zh
Publication of CN114972145A publication Critical patent/CN114972145A/zh
Application granted granted Critical
Publication of CN114972145B publication Critical patent/CN114972145B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/04—Architecture, e.g. interconnection topology
    • G06N3/045—Combinations of networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
    • G06V10/774—Generating sets of training patterns; Bootstrap methods, e.g. bagging or boosting
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/82—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using neural networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10088—Magnetic resonance imaging [MRI]
    • G06T2207/10092—Diffusion tensor magnetic resonance imaging [DTI]

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Software Systems (AREA)
  • Computing Systems (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Linguistics (AREA)
  • Mathematical Physics (AREA)
  • Biophysics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Databases & Information Systems (AREA)
  • Medical Informatics (AREA)
  • Multimedia (AREA)
  • Image Analysis (AREA)

Abstract

本发明提供一种基于深度学习的扩散张量重建方法、系统和电子设备：获取T1加权图像、非扩散加权图像、扩散加权图像，分别简称为T1图像、b0图像、DW图像；拟合b0图像和DW图像，得到扩散张量；选用T1图像、部分b0图像和DW图像，以相邻三层为一组，获得原始网络输入图像组；对原始网络输入图像组进行排序，获得三层排序图像组和重建层排序图像组；构建并训练深度神经网络；对待成像对象进行图像采集，将采集的图像排序得到待重建的三层排序图像组和待重建的重建层排序图像组并输入训练好的深度神经网络重建得到扩散张量；由扩散张量计算得到DTI定量图像。本发明可减少扩散张量及其定量图重建所需DW图像数量，加快成像速度。

Description

基于深度学习的扩散张量重建方法、系统和电子设备
技术领域
本发明涉及磁共振成像领域，特别是指基于深度学习的扩散张量重建方法、系统和电子设备。
背景技术
扩散磁共振成像(DiffusionMagnetic Resonance Imaging，dMRI)利用水分子沿扩散梯度方向扩散的信息，可以提供组织的微结构信息。扩散张量成像(Diffusion TensorImaging，DTI)是临床常用的一种dMRI技术，在脑部疾病诊断和神经科学研究上具有重要应用价值。DTI由扩散加权(Diffusion-weighted，DW)图像拟合而来。由于DTI模型需要拟合6个参数，因此重建扩散张量至少需要6张DW图像。另一方面，DTI模型对噪声的敏感度很高，而DW图像通常信噪比较低，因此传统方法上高质量的DTI需要多于30张的DW图像来进行最小二乘拟合，导致DTI需要长的采样时间，而且在长时间的采样过程中，病人难免产生运动，从而导致图像伪影等问题。
发明内容
本发明的主要目的在于克服现有技术中的上述缺陷，提出一种基于深度学习的扩散张量重建方法、系统和电子设备，可大大减少扩散张量及其定量图重建所需的DW图像的数量，从而减少扫描获取DW图像的时间，加快获取高质量DTI的速度。
本发明采用如下技术方案：
一种基于深度学习的扩散张量重建方法，其特征在于，包括如下步骤：
S1获取图像，图像数量为P层，每层包括1张T1图像、L张b0图像和N张扩散方向不同的DW图像，T1图像为T1加权图像，b0图像为非扩散加权图像，DW图像为扩散加权图像，N个扩散方向由设计好的扩散梯度方向表确定；
S2对所有b0图像取平均，并将其连同所有DW图像，按照DTI拟合公式，采用最小二乘算法逐像素进行信号拟合，得到扩散张量标签；
S3沿层的方向抽取图像得到所有层循环DW图像，沿层的方向选取图像组成原始网络输入图像组；
S4对每个原始网络输入图像组进行排序，获得对应的三层排序网络输入图像组，并从三层排序网络输入图像组中抽取重建层图像，获得对应的重建层排序网络输入图像组；
S5构建能够适应三层排序网络输入图像组和重建层排序网络输入图像组同时输入的深度神经网络用于扩散张量重建，使用每个三层排序网络输入图像组和对应的重建层排序网络输入图像组作为网络输入，以扩散张量标签作为网络标签对深度神经网络进行训练，得到对应的训练好的深度神经网络；
S6以每层1张T1图像、1张b0图像和所有层循环DW图像的采集方案对待成像对象进行图像采集，将所采集的图像组成若干个待重建的三层排序网络输入图像组和对应的待重建的重建层排序网络输入图像组，将每个待重建的三层排序网络输入图像组和对应的待重建的重建层排序网络输入图像组输入训练好的深度神经网络，重建出每个重建层的扩散张量；
S7利用重建出的扩散张量进行特征值分解和公式计算，得到DTI定量图像。
优选的，步骤S3中沿层的方向抽取图像得到所有层循环DW图像，具体为：沿层的方向从任一端开始，以三层为一循环，按扩散梯度方向表的前M个扩散方向顺序每层抽取M/3张DW图像，获得所有层循环DW图像。
优选的，步骤S3中，沿层的方向选取图像组成原始网络输入图像组，具体为：
沿层的方向以所述抽取图像开始的同一端的第二层作为重建层，加上与其相邻的两层为一组，选取此三层的图像组成一个原始网络输入图像组；
不断向另一端移动重建层，并以新的重建层及其相邻两层为一组，组成新的原始网络输入图像组，直至重建层移动到另一端倒数第二层，获得P-2个所述原始网络输入图像组。
优选的，步骤S3中，以重建层和其相邻两层为一组，选取各层1张T1图像、第1张b0图像、循环DW图像，组成一个原始网络输入图像组，如此，一个原始网络输入图像组包含三层共3张T1图像、3张b0图像和M张DW图像；将图像按照3张T1图像-3张b0图像-M张DW图像的顺序排列，相同类型图像按照采集时层的顺序排列。
优选的，步骤S4具体包括如下：
以扩散梯度方向表的前M个扩散方向顺序作为DW图像的排列顺序，对每个原始网络输入图像组内的DW图像进行重新排序，从而形成新的三层排序；
以新的三层排序对原始网络输入图像组中的T1图像和b0图像进行排序，获得对应的三层排序网络输入图像组；
从三层排序网络输入图像组中选取重建层的T1图像、b0图像和DW图像组成重建层排序网络输入图像组。
一种基于深度学习的扩散张量重建系统，其特征在于，包括：
图像获取单元：获取图像，图像数量为P层，每层包括1张T1图像、L张b0图像和N张扩散方向不同的DW图像，T1图像为T1加权图像，b0图像为非扩散加权图像，DW图像为扩散加权图像，N个扩散方向由设计好的扩散梯度方向表确定；
标签获取单元：对所有b0图像取平均，并将其连同所有DW图像，按照DTI拟合公式，采用最小二乘算法逐像素进行信号拟合，得到扩散张量标签；
图像生成单元：沿层的方向抽取图像得到所有层循环DW图像，沿层的方向选取图像组成原始网络输入图像组；
图像处理单元：对每个原始网络输入图像组进行排序，获得对应的三层排序网络输入图像组，并从三层排序网络输入图像组中抽取重建层图像，获得对应的重建层排序网络输入图像组；
网络构建和训练单元：构建能够适应三层排序网络输入图像组和重建层排序网络输入图像组同时输入的深度神经网络用于扩散张量重建，使用每个三层排序网络输入图像组和对应的重建层排序网络输入图像组作为网络输入，以扩散张量标签作为网络标签对深度神经网络进行训练，得到对应的训练好的深度神经网络；
重建层扩散张量重建单元：以每层1张T1图像、1张b0图像和所有层循环DW图像的采集方案对待成像对象进行图像采集，将所采集的图像组成若干个待重建的三层排序网络输入图像组和对应的待重建的重建层排序网络输入图像组，将每个待重建的三层排序网络输入图像组和对应的待重建的重建层排序网络输入图像组输入训练好的深度神经网络，重建出每个重建层的扩散张量；
DTI定量图像生成单元：利用重建出的扩散张量进行特征值分解和公式计算，得到DTI定量图像。
优选的，所述图像生成单元还包括：
图像选取子单元：沿层的方向从任一端开始，以三层为一循环，按扩散梯度方向表的前M个扩散方向顺序每层抽取M/3张DW图像，获得所有层循环DW图像；沿层的方向以同一端的第二层作为重建层，加上与其相邻的两层为一组，选取此组各层1张T1图像、第1张b0图像、循环DW图像，组成一个原始网络输入图像组，如此，一个原始网络输入图像组包含三层共3张T1图像、3张b0图像和M张DW图像，不断向另一端移动重建层，并以新的重建层及其相邻两层为一组，组成新的原始网络输入图像组，直至重建层移动到另一端倒数第二层，获得P-2个所述原始网络输入图像组；
图像排序子单元：将原始网络输入图像组的图像按照3张T1图像-3张b0图像-M张DW图像的顺序排列，相同类型图像按照采集时层的顺序排列。
优选的，所述图像处理单元包括：
DW图像排序子单元：以扩散梯度方向表的前M个扩散方向顺序作为DW图像的排列顺序，对每个原始网络输入图像组内的DW图像进行重新排序，从而形成新的三层排序；
T1图像与b0图像排序子单元：以新的三层排序对原始网络输入图像组中的T1图像和b0图像进行排序，获得对应的三层排序网络输入图像组；
重建层选取子单元：从三层排序网络输入图像组中选取重建层的T1图像、b0图像和DW图像组成重建层排序网络输入图像组。
一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于：处理器执行计算机程序时实现所述的基于深度学习的扩散张量重建方法。
由上述对本发明的描述可知，与现有技术相比，本发明具有如下有益效果：
本发明包括：获取T1加权图像、非扩散加权图像和扩散加权图像；拟合b0图像和DW图像，得到扩散张量标签；选用T1图像、部分b0图像和DW图像，以相邻三层为一组，获得原始网络输入图像组；对原始网络输入图像组进行排序，获得三层排序图像组和重建层排序图像组；构建并训练深度神经网络；对待成像对象进行图像采集，将采集的图像排序得到待重建的三层排序图像组和待重建的重建层排序图像组并输入训练好的深度神经网络重建得到扩散张量；由扩散张量计算得到DTI定量图像。
本发明实现了快速重建扩散张量及其定量图，有效改善了因DW图像信噪比差造成利用少量DW图像重建DTI扩散张量及其定量图质量差的问题，提供了高质量的扩散张量及其定量图。
附图说明
图1为本发明方法流程图；
图2为本发明实施例每层使用3张DW图像得到的扩散张量及其定量图：(a)扩散张量的对角元素之一Dxx；(b)扩散张量的非对角元素之一Dxy；(c)MD图；(d)FA图；
图3为本发明实施例每层使用6张DW图像采用传统逐点最小二乘拟合方法得到的扩散张量及其定量图：(a)扩散张量的对角元素之一Dxx；(b)扩散张量的非对角元素之一Dxy；(c)MD图；(d)FA图；
图4为本发明实施例一种基于深度学习的扩散张量重建系统架构图；
图5为本发明实施例一种电子设备实施例示意图。
以下结合附图和具体实施例对本发明作进一步详述。
具体实施方式
以下通过具体实施方式对本发明作进一步的描述。
参见图1，一种基于深度学习的扩散张量重建方法，可大大减少扩散张量及其定量图重建时所需的DW图像的数量，从而减少扫描获取DW图像的时间，加快获取高质量DTI的速度，具体包括如下步骤：
S1获取图像，图像数量为P层，每层包括1张T1图像、L张b0图像和N张扩散方向不同的DW图像，T1图像为T1加权图像，b0图像为非扩散加权图像，DW图像为扩散加权图像。该步骤中，L不小于2，N一般不小于30，N个扩散方向由设计好的扩散梯度方向表确定，P不小于3。
在实际应用中，所述S1中，获取的T1图像、b0图像和DW图像均来源于人类连接组计划(Human Connectome Project，HCP)数据集。从数据集中每个志愿者抽取18张b0图像(即L＝18)、90张b＝1000s/mm2的DW图像(即N＝90)和1张T1图像。90张DW图像来自90个均匀分布的扩散加权方向，扩散梯度方向表由数据集给出。其中90个扩散加权方向由数据集的预处理进行了优化，抽取任意前9个方向(即M＝9)，这9个方向仍然保持均匀分布。数据集的预处理也包含DW图像的涡流校正和运动校正等。其中所有图像都为全脑图像，图像成像矩阵为174×145，每个志愿者共145层(即P＝145)，136个志愿者数据用于后续的训练，46个志愿者数据用于后续验证。
S2对所有b0图像取平均，并将其连同所有DW图像，按照DTI拟合公式，采用最小二乘算法逐像素进行信号拟合，得到扩散张量标签。具体如下：
按照DTI拟合公式，采用最小二乘算法，使用每层所有的b0图像和DW图像，逐像素进行信号拟合，获得扩散张量标签。DTI拟合公式为：
Figure BDA0003663774890000053
其中Si为DW图像的信号强度，S0为b0图像的信号强度，gi＝(gix,giy,giz)T为扩散梯度的单位方向向量，上标T表示向量或矩阵的转置，b为扩散梯度因子，D为扩散张量，可以表示为：
Figure BDA0003663774890000051
因此拟合公式可以改写为：
Figure BDA0003663774890000052
其中αi＝(gix 2,gix 2,giz 2,2gixgiy,2gixgiz,2giygiz)T，βi＝ln(S0/Si)/b，Dvec＝(Dxx,Dyy,Dzz,Dxy,Dxz,Dyz)T。在此基础上，拟合公式可以再进一步改写为：
Dvec＝A-1B
其中A＝(α1,α2,α3,…,αn)T，B＝(β1,β2,β3,…,βn)T；
最终通过最小二乘算法拟合出Dvec，即重建了扩散张量。
在实际应用中，Si是S1中从HCP数据集抽取的第i张DW图像的信号强度，S0是S1中抽取的18张b0图像的平均信号强度，gi是与第i张图像对应的扩散梯度的单位方向向量，b值取1000s/mm2。
S3沿层的方向抽取图像得到所有层循环DW图像，沿层的方向选取图像组成原始网络输入图像组。
其中，沿层的方向抽取图像得到所有层循环DW图像，具体为：沿层的方向从任一端开始，以三层为一循环，按扩散梯度方向表的前M个扩散方向顺序每层抽取M/3张DW图像，获得所有层循环DW图像，M/3远小于30且一般不大于6。实际应用中，可选取扩散梯度方向表的前9个扩散方向进行循环，按照顺序每层抽取3张DW图像。
沿层的方向选取图像组成原始网络输入图像组，具体为：
S3.1沿层的方向以抽取所有层循环DW图像开始的同一端的第二层作为重建层，加上与其相邻的两层为一组，选取此三层的图像组成一个原始网络输入图像组；具体的，以重建层和其相邻两层为一组，选取各层1张T1图像、第1张b0图像、循环DW图像，组成一个原始网络输入图像组，如此，一个原始网络输入图像组包含三层共3张T1图像、3张b0图像和M张DW图像。实际应用中，获得的原始网络输入图像组尺寸为174×145，共15个通道。
S3.2将图像按照3张T1图像-3张b0图像-M张DW图像的顺序排列，相同类型图像按照采集时层的顺序排列。
S3.3不断向另一端移动重建层，并以新的重建层及其相邻两层为一组，组成新的原始网络输入图像组，直至重建层移动到另一端倒数第二层，获得P-2个所述原始网络输入图像组。
S4对每个原始网络输入图像组进行排序，获得对应的三层排序网络输入图像组，并从三层排序网络输入图像组中抽取重建层图像，获得对应的重建层排序网络输入图像组。具体包括如下：
S4.1以扩散梯度方向表的前M个扩散方向顺序作为DW图像的排列顺序，对每个原始网络输入图像组内的DW图像进行重新排序，从而形成新的三层排序；
S4.2以新的三层排序对原始网络输入图像组中的T1图像和b0图像进行排序，获得对应的三层排序网络输入图像组。实际应用中，获得的三层排序网络输入图像组尺寸可为174×145，共15个通道。
S4.3从三层排序网络输入图像组中选取重建层的T1图像、b0图像和DW图像组成重建层排序网络输入图像组。实际应用中，获得的重建层排序网络输入图像组尺寸可为174×145，共5个通道。
S5构建能够适应三层排序网络输入图像组和重建层排序网络输入图像组同时输入的深度神经网络用于扩散张量重建，使用每个三层排序网络输入图像组和对应的重建层排序网络输入图像组作为网络输入，以扩散张量标签作为网络标签对深度神经网络进行训练，得到对应的训练好的深度神经网络。
实际应用中，深度神经网络使用常见的深度学习框架PyTorch搭建。神经网络模型基于五层的U-Net神经网络模型改进。在U-Net一个编码器一个解码器的网络结构基础上，将编码器拆分成两个编码器，从而适应同时两个输入。网络的每一层中均包含有相同的卷积、池化和采样结构，相同分辨率层面采用跳跃连接结构沟通。重建结果采用L1范数进行约束，训练过程不分块，批次数为16，总共训练100轮，其中每轮会完整地对所有训练数据和验证数据进行训练，并且更新网络权重。
对深度神经网络进行训练时，共有19720组样本用于训练，6670组样本用于验证。
S6以每层1张T1图像、1张b0图像和所有层循环DW图像的采集方案对待成像对象进行图像采集，将所采集的图像组成若干个待重建的三层排序网络输入图像组和对应的待重建的重建层排序网络输入图像组，将每个待重建的三层排序网络输入图像组和对应的待重建的重建层排序网络输入图像组输入训练好的深度神经网络，重建出每个重建层的扩散张量。
实际应用中，以每层1张T1图像、1张b0图像和所有层循环DW图像的采集方案对待成像对象进行图像采集指的是HCP数据集内，排除已经用于训练和验证的志愿者后，在额外的21个志愿者数据上以生成训练样本的相同方法生成的3045个样本。
S7利用重建出的扩散张量进行特征值分解和公式计算，得到DTI定量图像，包括但不限于MD、FA等。
该步骤中，重建的扩散张量进行特征值分解指的是对每个像素点的扩散张量单独进行特征值分解得到该像素点的扩散张量的特征值，进行进一步的公式计算。
本实施例每层使用3张DW图像得到的扩散张量及其定量图，分别展示在图2：(a)扩散张量的对角元素之一Dxx；(b)扩散张量的非对角元素之一Dxx；(c)MD图；(d)FA图。
为了评价本发明提出的方法，本发明提供了每层使用6张DW图像采用传统逐点最小二乘拟合方法得到的扩散张量及其定量图，分别展示在图3：(a)扩散张量的对角元素之一Dxx；(b)扩散张量的非对角元素之一Dxx；(c)MD图；(d)FA图。从图2和图3的对比可以看出，本发明可以较好地还原成像对象的细节，而且重建的扩散张量和定量图几乎没有噪点。同时，本发明每层使用的DW图像仅为3张，采集时间更少；本发明每次重建只需几秒时间，远快于传统算法每次需要3-4分钟的重建时间。
参见图4，本发明还提出一种基于深度学习的扩散张量重建系统，该系统采用上述的一种基于深度学习的扩散张量重建方法，包括：
图像获取单元10：获取图像，图像数量为P层，每层包括1张T1图像、L张b0图像和N张扩散方向不同的DW图像，T1图像为T1加权图像，b0图像为非扩散加权图像，DW图像为扩散加权图像。其中，L不小于2，N一般不小于30，N个扩散方向由设计好的扩散梯度方向表确定，P不小于3。
标签获取单元20：对所有b0图像取平均，并将其连同所有DW图像，按照DTI拟合公式，采用最小二乘算法逐像素进行信号拟合，得到扩散张量标签。具体为：按照DTI拟合公式，采用最小二乘算法，使用每层所有的b0图像和DW图像，逐像素进行信号拟合，获得扩散张量标签。DTI拟合公式为：
Figure BDA0003663774890000081
其中Si为DW图像的信号强度，S0为b0图像的信号强度，gi＝(gix,giy,giz)T为扩散梯度的单位方向向量，上标T表示向量或矩阵的转置，b为扩散梯度因子，D为扩散张量，可以表示为：
Figure BDA0003663774890000082
因此拟合公式可以改写为：
Figure BDA0003663774890000083
其中αi＝(gix 2,gix 2,giz 2,2gixgiy,2gixgiz,2giygiz)T，βi＝ln(S0/Si)/b，Dvec＝(Dxx,Dyy,Dzz,Dxy,Dxz,Dyz)T。在此基础上，拟合公式可以再进一步改写为：
Dvec＝A-1B
其中A＝(α1,α2,α3,…,αn)T，B＝(β1,β2,β3,…,βn)T；
最终通过最小二乘算法拟合出Dvec，即重建了扩散张量。
图像生成单元30：沿层的方向抽取图像得到所有层循环DW图像，沿层的方向选取图像组成原始网络输入图像组。具体的，其包括：
图像选取子单元31：沿层的方向从任一端开始，以三层为一循环，按扩散梯度方向表的前M个扩散方向顺序每层抽取M/3张DW图像，获得所有层循环DW图像；沿层的方向以同一端的第二层作为重建层，加上与其相邻的两层为一组，选取此组各层1张T1图像、第1张b0图像、循环DW图像，组成一个原始网络输入图像组，如此，一个原始网络输入图像组包含三层共3张T1图像、3张b0图像和M张DW图像，不断向另一端移动重建层，并以新的重建层及其相邻两层为一组，组成新的原始网络输入图像组，直至重建层移动到另一端倒数第二层，获得P-2个所述原始网络输入图像组。
图像排序子单元32：将原始网络输入图像组的图像按照3张T1图像-3张b0图像-M张DW图像的顺序排列，相同类型图像按照采集时层的顺序排列。
图像处理单元40：对每个原始网络输入图像组进行排序，获得对应的三层排序网络输入图像组，并从三层排序网络输入图像组中抽取重建层图像，获得对应的重建层排序网络输入图像组。具体的，其还包括：
DW图像排序子单元41：以扩散梯度方向表的前M个扩散方向顺序作为DW图像的排列顺序，对每个原始网络输入图像组内的DW图像进行重新排序，从而形成新的三层排序；
T1图像与b0图像排序子单元42：以新的三层排序对原始网络输入图像组中的T1图像和b0图像进行排序，获得对应的三层排序网络输入图像组；
重建层选取子单元43：从三层排序网络输入图像组中选取重建层的T1图像、b0图像和DW图像组成重建层排序网络输入图像组。
网络构建和训练单元50：构建能够适应三层排序网络输入图像组和重建层排序网络输入图像组同时输入的深度神经网络用于扩散张量重建，使用每个三层排序网络输入图像组和对应的重建层排序网络输入图像组作为网络输入，以扩散张量标签作为网络标签对深度神经网络进行训练，得到对应的训练好的深度神经网络。
重建层扩散张量重建单元60：以每层1张T1图像、1张b0图像和所有层循环DW图像的采集方案对待成像对象进行图像采集，将所采集的图像组成若干个待重建的三层排序网络输入图像组和对应的待重建的重建层排序网络输入图像组，将每个待重建的三层排序网络输入图像组和对应的待重建的重建层排序网络输入图像组输入训练好的深度神经网络，重建出每个重建层的扩散张量。
DTI定量图像生成单元70：利用重建出的扩散张量进行特征值分解和公式计算，得到DTI定量图像，包括但不限于MD、FA等。
本发明的目的是提供一种基于深度学习的快速高质量扩散张量及其定量图的重建方法及系统，可大大减少扩散张量及其定量图重建所需的DW图像的数量，从而减少扫描获取DW图像的时间，加快获取高质量DTI的速度。
参见图5，本发明实施例提供了一种电子设备500，包括存储器510、处理器520及存储在存储器510上并可在处理器520上运行的计算机程序511，处理器520执行计算机程序511时实现本发明实施例提供的一种基于深度学习的扩散张量重建方法。
上述仅为本发明的具体实施方式，但本发明的设计构思并不局限于此，凡利用此构思对本发明进行非实质性的改动，均应属于侵犯本发明保护范围的行为。

Claims (9)

1.一种基于深度学习的扩散张量重建方法，其特征在于，包括如下步骤：
S1获取图像，图像数量为P层，每层包括1张T1图像、L张b0图像和N张扩散方向不同的DW图像，T1图像为T1加权图像，b0图像为非扩散加权图像，DW图像为扩散加权图像，N个扩散方向由设计好的扩散梯度方向表确定；
S2对所有b0图像取平均，并将其连同所有DW图像，按照DTI拟合公式，采用最小二乘算法逐像素进行信号拟合，得到扩散张量标签；
S3沿层的方向抽取图像得到所有层循环DW图像，沿层的方向选取图像组成原始网络输入图像组；
S4对每个原始网络输入图像组进行排序，获得对应的三层排序网络输入图像组，并从三层排序网络输入图像组中抽取重建层图像，获得对应的重建层排序网络输入图像组；
S5构建能够适应三层排序网络输入图像组和重建层排序网络输入图像组同时输入的深度神经网络用于扩散张量重建，使用每个三层排序网络输入图像组和对应的重建层排序网络输入图像组作为网络输入，以扩散张量标签作为网络标签对深度神经网络进行训练，得到对应的训练好的深度神经网络；
S6以每层1张T1图像、1张b0图像和所有层循环DW图像的采集方案对待成像对象进行图像采集，将所采集的图像组成若干个待重建的三层排序网络输入图像组和对应的待重建的重建层排序网络输入图像组，将每个待重建的三层排序网络输入图像组和对应的待重建的重建层排序网络输入图像组输入训练好的深度神经网络，重建出每个重建层的扩散张量；
S7利用重建出的扩散张量进行特征值分解和公式计算，得到DTI定量图像。
2.如权利要求1所述的一种基于深度学习的扩散张量重建方法，其特征在于，步骤S3中沿层的方向抽取图像得到所有层循环DW图像，具体为：沿层的方向从任一端开始，以三层为一循环，按扩散梯度方向表的前M个扩散方向顺序每层抽取M/3张DW图像，获得所有层循环DW图像。
3.如权利要求2所述的一种基于深度学习的扩散张量重建方法，其特征在于，步骤S3中，沿层的方向选取图像组成原始网络输入图像组，具体为：
沿层的方向以所述抽取图像开始的同一端的第二层作为重建层，加上与其相邻的两层为一组，选取此三层的图像组成一个原始网络输入图像组；
不断向另一端移动重建层，并以新的重建层及其相邻两层为一组，组成新的原始网络输入图像组，直至重建层移动到另一端倒数第二层，获得P-2个所述原始网络输入图像组。
4.如权利要求3所述的一种基于深度学习的扩散张量重建方法，其特征在于，步骤S3中，以重建层和其相邻两层为一组，选取各层1张T1图像、第1张b0图像、循环DW图像，组成一个原始网络输入图像组，如此，一个原始网络输入图像组包含三层共3张T1图像、3张b0图像和M张DW图像；将图像按照3张T1图像-3张b0图像-M张DW图像的顺序排列,相同类型图像按照采集时层的顺序排列。
5.如权利要求1所述的一种基于深度学习的扩散张量重建方法，其特征在于，步骤S4具体包括如下：
以扩散梯度方向表的前M个扩散方向顺序作为DW图像的排列顺序，对每个原始网络输入图像组内的DW图像进行重新排序，从而形成新的三层排序；
以新的三层排序对原始网络输入图像组中的T1图像和b0图像进行排序，获得对应的三层排序网络输入图像组；
从三层排序网络输入图像组中选取重建层的T1图像、b0图像和DW图像组成重建层排序网络输入图像组。
6.一种基于深度学习的扩散张量重建系统，其特征在于，包括：
图像获取单元：获取图像，图像数量为P层，每层包括1张T1图像、L张b0图像和N张扩散方向不同的DW图像，T1图像为T1加权图像，b0图像为非扩散加权图像，DW图像为扩散加权图像，N个扩散方向由设计好的扩散梯度方向表确定；
标签获取单元：对所有b0图像取平均，并将其连同所有DW图像，按照DTI拟合公式，采用最小二乘算法逐像素进行信号拟合，得到扩散张量标签；
图像生成单元：沿层的方向抽取图像得到所有层循环DW图像，沿层的方向选取图像组成原始网络输入图像组；
图像处理单元：对每个原始网络输入图像组进行排序，获得对应的三层排序网络输入图像组，并从三层排序网络输入图像组中抽取重建层图像，获得对应的重建层排序网络输入图像组；
网络构建和训练单元：构建能够适应三层排序网络输入图像组和重建层排序网络输入图像组同时输入的深度神经网络用于扩散张量重建，使用每个三层排序网络输入图像组和对应的重建层排序网络输入图像组作为网络输入，以扩散张量标签作为网络标签对深度神经网络进行训练，得到对应的训练好的深度神经网络；
重建层扩散张量重建单元：以每层1张T1图像、1张b0图像和所有层循环DW图像的采集方案对待成像对象进行图像采集，将所采集的图像组成若干个待重建的三层排序网络输入图像组和对应的待重建的重建层排序网络输入图像组，将每个待重建的三层排序网络输入图像组和对应的待重建的重建层排序网络输入图像组输入训练好的深度神经网络，重建出每个重建层的扩散张量；
DTI定量图像生成单元：利用重建出的扩散张量进行特征值分解和公式计算，得到DTI定量图像。
7.如权利要求6所述的一种基于深度学习的扩散张量重建系统，其特征在于，所述图像生成单元还包括：
图像选取子单元：沿层的方向从任一端开始，以三层为一循环，按扩散梯度方向表的前M个扩散方向顺序每层抽取M/3张DW图像，获得所有层循环DW图像；沿层的方向以同一端的第二层作为重建层，加上与其相邻的两层为一组，选取此组各层1张T1图像、第1张b0图像、循环DW图像，组成一个原始网络输入图像组，如此，一个原始网络输入图像组包含三层共3张T1图像、3张b0图像和M张DW图像，不断向另一端移动重建层，并以新的重建层及其相邻两层为一组，组成新的原始网络输入图像组，直至重建层移动到另一端倒数第二层，获得P-2个所述原始网络输入图像组；
图像排序子单元：将原始网络输入图像组的图像按照3张T1图像-3张b0图像-M张DW图像的顺序排列，相同类型图像按照采集时层的顺序排列。
8.如权利要求6所述的一种基于深度学习的扩散张量的重建系统，其特征在于，所述图像处理单元包括：
DW图像排序子单元：以扩散梯度方向表的前M个扩散方向顺序作为DW图像的排列顺序，对每个原始网络输入图像组内的DW图像进行重新排序，从而形成新的三层排序；
T1图像与b0图像排序子单元：以新的三层排序对原始网络输入图像组中的T1图像和b0图像进行排序，获得对应的三层排序网络输入图像组；
重建层选取子单元：从三层排序网络输入图像组中选取重建层的T1图像、b0图像和DW图像组成重建层排序网络输入图像组。
9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于：处理器执行计算机程序时实现权利要求1至5中任一项所述的基于深度学习的扩散张量重建方法。
CN202210581140.7A 2022-05-26 2022-05-26 基于深度学习的扩散张量重建方法、系统和电子设备 Active CN114972145B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210581140.7A CN114972145B (zh) 2022-05-26 2022-05-26 基于深度学习的扩散张量重建方法、系统和电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210581140.7A CN114972145B (zh) 2022-05-26 2022-05-26 基于深度学习的扩散张量重建方法、系统和电子设备

Publications (2)

Publication Number Publication Date
CN114972145A true CN114972145A (zh) 2022-08-30
CN114972145B CN114972145B (zh) 2024-06-28

Family

ID=82955105

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210581140.7A Active CN114972145B (zh) 2022-05-26 2022-05-26 基于深度学习的扩散张量重建方法、系统和电子设备

Country Status (1)

Country Link
CN (1) CN114972145B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106780643A (zh) * 2016-11-21 2017-05-31 清华大学 磁共振多次激发扩散成像运动矫正方法
US20190333199A1 (en) * 2018-04-26 2019-10-31 The Regents Of The University Of California Systems and methods for deep learning microscopy
US20200096592A1 (en) * 2018-09-25 2020-03-26 Siemens Healthineers Ltd. Magnetic resonance diffusion tensor imaging method and device, and fiber tracking method and device
CN111738999A (zh) * 2020-06-15 2020-10-02 西安电子科技大学 基于弥散张量成像的大脑张量模板构建方法
CN112419203A (zh) * 2020-12-07 2021-02-26 贵州大学 基于对抗网络的扩散加权图像压缩感知恢复方法及装置
CN112763958A (zh) * 2020-12-10 2021-05-07 复旦大学 一种基于神经网络的多次激发平面回波磁共振成像方法
CN113490859A (zh) * 2019-02-25 2021-10-08 皇家飞利浦有限公司 使用多幅扩散加权mr图像来计算b0图像

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106780643A (zh) * 2016-11-21 2017-05-31 清华大学 磁共振多次激发扩散成像运动矫正方法
US20190333199A1 (en) * 2018-04-26 2019-10-31 The Regents Of The University Of California Systems and methods for deep learning microscopy
US20200096592A1 (en) * 2018-09-25 2020-03-26 Siemens Healthineers Ltd. Magnetic resonance diffusion tensor imaging method and device, and fiber tracking method and device
CN113490859A (zh) * 2019-02-25 2021-10-08 皇家飞利浦有限公司 使用多幅扩散加权mr图像来计算b0图像
CN111738999A (zh) * 2020-06-15 2020-10-02 西安电子科技大学 基于弥散张量成像的大脑张量模板构建方法
CN112419203A (zh) * 2020-12-07 2021-02-26 贵州大学 基于对抗网络的扩散加权图像压缩感知恢复方法及装置
CN112763958A (zh) * 2020-12-10 2021-05-07 复旦大学 一种基于神经网络的多次激发平面回波磁共振成像方法

Also Published As

Publication number Publication date
CN114972145B (zh) 2024-06-28

Similar Documents

Publication Publication Date Title
CN110461228B (zh) 使用多对比度和深度学习改进医学图像的质量
US10852379B2 (en) Artifact reduction by image-to-image network in magnetic resonance imaging
Ozcelik et al. Reconstruction of perceived images from fmri patterns and semantic brain exploration using instance-conditioned gans
CN111462020B (zh) 心脏图像的运动伪影校正方法、系统、存储介质和设备
CN110444277B (zh) 一种基于多生成多对抗的多模态脑部mri图像双向转换方法
CN106981051A (zh) 用于运动解析的mri的系统和方法
CN111598964B (zh) 一种基于空间自适应网络的定量磁化率图像重建方法
CN117036162B (zh) 轻量级胸部ct图像超分辨率的残差特征注意力融合方法
CN117223028A (zh) 用于带去噪磁共振图像重建的系统和方法
CN110853048A (zh) 基于粗、精训练的mri图像分割方法、装置和存储介质
WO2020252307A1 (en) Enhancements to quantitative magnetic resonance imaging techniques
CN113835059A (zh) 一种用于扩散加权磁共振成像的数据处理方法和系统
CN115861464A (zh) 基于多模mri同步生成的伪ct合成方法
CN116309524A (zh) 基于深度学习的心脏磁共振电影成像伪影抑制方法和系统
Wang et al. One for multiple: Physics-informed synthetic data boosts generalizable deep learning for fast MRI reconstruction
CN116912252B (zh) 面向模态缺失的影像数据的脑图谱个体化方法及系统
CN103985111B (zh) 一种基于双字典学习的4d‑mri超分辨率重构方法
CN111161370B (zh) 一种基于ai的人体多核dwi联合重建方法
CN114972145B (zh) 基于深度学习的扩散张量重建方法、系统和电子设备
CN116309754A (zh) 一种基于局部-全局信息协作的大脑医学图像配准方法及系统
CN107209242A (zh) 用于生成在其环境中的移动的对象的mr图像的方法和系统
KR20190117234A (ko) 인공신경망을 이용한 자기 공명 영상의 영상 프로토콜 선택 장치와 방법 및 프로그램이 기록된 컴퓨터 판독 가능한 기록매체
Qi et al. Multi-task MR imaging with iterative teacher forcing and re-weighted deep learning
WO2022120731A1 (zh) 基于循环生成对抗网络的mri-pet图像模态转换方法及系统
CN116705251A (zh) 一种基于深度学习的脑部mr影像多方位规范化系统

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