CN110827331B - 图像配准模型的训练方法、图像配准方法和计算机设备 - Google Patents

图像配准模型的训练方法、图像配准方法和计算机设备 Download PDF

Info

Publication number
CN110827331B
CN110827331B CN201911063804.5A CN201911063804A CN110827331B CN 110827331 B CN110827331 B CN 110827331B CN 201911063804 A CN201911063804 A CN 201911063804A CN 110827331 B CN110827331 B CN 110827331B
Authority
CN
China
Prior art keywords
image
image registration
loss
deformation field
network
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
CN201911063804.5A
Other languages
English (en)
Other versions
CN110827331A (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.)
Shanghai United Imaging Intelligent Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Intelligent Healthcare Co Ltd
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 Shanghai United Imaging Intelligent Healthcare Co Ltd filed Critical Shanghai United Imaging Intelligent Healthcare Co Ltd
Priority to CN201911063804.5A priority Critical patent/CN110827331B/zh
Publication of CN110827331A publication Critical patent/CN110827331A/zh
Application granted granted Critical
Publication of CN110827331B publication Critical patent/CN110827331B/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
    • G06T7/00—Image analysis
    • G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
    • 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/10081—Computed x-ray tomography [CT]
    • 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/20—Special algorithmic details
    • G06T2207/20081—Training; Learning

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)
  • Medical Treatment And Welfare Office Work (AREA)

Abstract

本申请涉及一种图像配准模型的训练方法、图像配准方法和计算机设备。图像配准模型的训练方法包括：获取n个训练样本图像；n≥3；将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场；根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失，以及根据n个变形场计算信息约束损失；根据图像配准损失和信息约束损失计算网络损失，利用网络损失对初始图像配准网络进行训练；当网络损失达到满足预设条件时，初始图像配准网络训练完成，得到图像配准模型。该方法中，除考虑图像配准时的配准损失，还考虑了各变形场之间的形变联系，由此提高得到的变形场的准确性，也使得变形场具有形变一致性。

Description

图像配准模型的训练方法、图像配准方法和计算机设备
技术领域
本申请涉及图像处理技术领域，特别是涉及一种图像配准模型的训练方法、图像配准方法和计算机设备。
背景技术
图像配准是将两幅图像进行空间匹配的过程，若要将图像A配准到图像B，则是将B作为参考图像，A作为浮动图像，得到图像A配准到图像B的变形场。医学影像配准技术是医学影像处理的基础，在图像信息融合、辅助诊断、手术规划以及医学基础理论研究等领域发挥着十分重要的作用。早期的医学影像配准方法主要是基于灰度的配准方法，根据光流场理论将待配准影像作为连续运动的影像序列的两帧，然后计算变形场的大小，从而实现影像的配准。但是此方法在医学影像质量较差、且形变较大的情况下，配准精度明显下降。
随着人工智能技术的不断发展，基于深度学习的医学影像配准方法也应运而生，其需要大量的训练样本通过不断迭代的方式学习数据的特征。传统技术通常是将变形场的一阶梯度作为正则损失函数，指导图像配准网络训练过程中网络参数的更新，以训练得到收敛的图像配准网络。
但是，传统技术只考虑了变形场局部的平滑性，得到的图像配准网络精度较低，从而使用图像配准网络进行图像配准的结果准确度也较低，且不具备良好的形变一致性。
发明内容
基于此，有必要针对统技术图像配准的结果准确度较低，且不具备良好的形变一致性的问题，提供一种图像配准模型的训练方法、图像配准方法和计算机设备。
第一方面，本申请实施例提供一种图像配准模型的训练方法，包括：
获取n个训练样本图像；n≥3；
将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场；
根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失，以及根据n个变形场计算信息约束损失；
根据图像配准损失和信息约束损失计算网络损失，利用网络损失对初始图像配准网络进行训练；当网络损失达到满足预设条件时，初始图像配准网络训练完成，得到图像配准模型。
在其中一个实施例中，上述方法还包括：
根据第n个变形场计算变形场平滑损失；
根据图像配准损失、信息约束损失和变形场平滑损失计算网络损失。
在其中一个实施例中，将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场，包括：
将第i个训练样本图像和第i+1个训练样本图像输入初始图像配准网络进行图像配准，得到n-1个变形场；其中，1≤i≤n-1；
将第1个训练样本图像和第n个训练样本图像输入初始图像配准网络进行图像配准，得到第n个变形场。
在其中一个实施例中，根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失，包括：
利用第n个变形场对第n个变形场对应的第一训练样本图像进行图像空间变换，得到配准图像；
计算配准图像与第n个变形场对应的第二训练样本图像之间的差异值，将差异值作为图像配准损失。
在其中一个实施例中，根据n个变形场计算信息约束损失，包括：
将第1个变形场至第n-1个变形场进行变形场叠加，得到叠加后的变形场；
根据叠加后的变形场和第n个变形场计算信息约束损失。
在其中一个实施例中，根据叠加后的变形场和第n个变形场计算信息约束损失，包括：
将叠加后的变形场与第n个变形场中对应位置的像素点值进行求差，得到求差结果，将求差结果作为信息约束损失。
在其中一个实施例中，当n＝3时，将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场，包括：
将第1个训练样本图像作为参考图像，第2个训练样本图像作为浮动图像输入初始图像配准网络进行图像配准，得到第一个变形场；
将第2个训练样本图像作为参考图像，第3个训练样本图像作为浮动图像输入初始图像配准网络进行图像配准，得到第二个变形场；
将第1个训练样本图像作为参考图像，第3个训练样本图像作为浮动图像输入初始图像配准网络进行图像配准，得到第三个变形场。
在其中一个实施例中，根据第n个变形场计算变形场平滑损失，包括：
计算第n个变形场一阶梯度的所有像素点值的平方和，将平方和作为变形场平滑损失。
第二方面，本申请实施例提供一种图像配准方法，包括：
获取待配准图像和参考图像；
将待配准图像和参考图像输入图像配准模型进行图像配准，得到待配准图像的配准结果；其中，图像配准模型的训练方式包括：
获取n个训练样本图像；n≥3；
将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场；
根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失，以及根据n个变形场计算信息约束损失；
根据图像配准损失和信息约束损失计算网络损失，利用网络损失对初始图像配准网络进行训练；当网络损失达到满足预设条件时，初始图像配准网络训练完成，得到图像配准模型。
第三方面，本申请实施例提供一种图像配准模型的训练装置，包括：
第一获取模块，用于获取n个训练样本图像；n≥3；
第一图像配准模块，用于将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场；
计算模块，用于根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失，以及根据n个变形场计算信息约束损失；
训练模块，用于根据图像配准损失和信息约束损失计算网络损失，利用网络损失对初始图像配准网络进行训练；当网络损失达到满足预设条件时，初始图像配准网络训练完成，得到图像配准模型。
第四方面，本申请实施例提供一种图像配准装置，包括：
第二获取模块，用于获取待配准图像和参考图像；
第二图像配准模块，用于将待配准图像和参考图像输入图像配准模型进行图像配准，得到待配准图像的配准结果。
第五方面，本申请实施例提供一种计算机设备，包括存储器和处理器，存储器存储有计算机程序，处理器执行计算机程序时实现以下步骤：
获取n个训练样本图像；n≥3；
将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场；
根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失，以及根据n个变形场计算信息约束损失；
根据图像配准损失和信息约束损失计算网络损失，利用网络损失对初始图像配准网络进行训练；当网络损失达到满足预设条件时，初始图像配准网络训练完成，得到图像配准模型。
第六方面，本申请实施例还提供一种计算机设备，包括存储器和处理器，存储器存储有计算机程序，处理器执行计算机程序时实现以下步骤：
获取待配准图像和参考图像；
将待配准图像和参考图像输入图像配准模型进行图像配准，得到待配准图像的配准结果；其中，图像配准模型的训练方式包括：
获取n个训练样本图像；n≥3；
将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场；
根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失，以及根据n个变形场计算信息约束损失；
根据图像配准损失和信息约束损失计算网络损失，利用网络损失对初始图像配准网络进行训练；当网络损失达到满足预设条件时，初始图像配准网络训练完成，得到图像配准模型。
第七方面，本申请实施例提供一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
获取n个训练样本图像；n≥3；
将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场；
根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失，以及根据n个变形场计算信息约束损失；
根据图像配准损失和信息约束损失计算网络损失，利用网络损失对初始图像配准网络进行训练；当网络损失达到满足预设条件时，初始图像配准网络训练完成，得到图像配准模型。
第八方面，本申请实施例还提供一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
获取待配准图像和参考图像；
将待配准图像和参考图像输入图像配准模型进行图像配准，得到待配准图像的配准结果；其中，图像配准模型的训练方式包括：
获取n个训练样本图像；n≥3；
将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场；
根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失，以及根据n个变形场计算信息约束损失；
根据图像配准损失和信息约束损失计算网络损失，利用网络损失对初始图像配准网络进行训练；当网络损失达到满足预设条件时，初始图像配准网络训练完成，得到图像配准模型。
上述图像配准模型的训练方法、图像配准方法、装置、计算机设备和可读存储介质，首先获取n个训练样本图像；n≥3；将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场；根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失，以及根据n个变形场计算信息约束损失；根据图像配准损失和信息约束损失计算网络损失，利用网络损失对初始图像配准网络进行训练；当网络损失达到满足预设条件时，初始图像配准网络训练完成，得到图像配准模型。该方法中，除考虑图像配准时的配准损失，保证了图像配准模型的精度，还考虑了各变形场之间的形变联系，即前n-1个变形场的叠加与第n个变形场相等的性质，由此进一步提高根据图像配准模型得到的变形场的准确性，也使得变形场具有形变一致性。
附图说明
图1为一个实施例提供的图像配准模型的训练方法的流程示意图；
图2为另一个实施例提供的图像配准模型的训练方法的流程示意图；
图3为又一个实施例提供的图像配准模型的训练方法的流程示意图；
图3a为一个实施例提供的变形场叠加过程的示意图；
图3b为一个实施例提供的n＝3时图像配准模型的训练过程的示意图；
图4为一个实施例提供的图像配准方法的流程示意图；
图5为一个实施例提供的图像配准模型的训练装置的结构示意图；
图6为一个实施例提供的图像配准装置的结构示意图；
图7为一个实施例提供的计算机设备的内部结构示意图。
具体实施方式
本申请实施例提供的图像配准模型的训练方法，可以适用于进行图像配准的网络模型的训练过程，此图像可以为医学图像，也可以为其他领域中的图像。图像配准主要分为刚性和非刚性两类，刚性配准主要解决简单的图像整体移动，如平移、旋转等问题；非刚性配准主要解决的是图像的弹性变换的问题，它容许变换过程中任意两个像素点之间的对应关系发生变动。医学图像配准技术是手术计划定制、影像引导放射治疗以及医学影像融合处理等方面的基础和前提，具有重要的临床应用价值；同时，人体的大部分组织和器官具有非刚性特征，因而针对非刚性医学影像的研究具有重要意义。
医学图像配准是指对于一幅医学图像寻求一种或者一系列的空间变换，使它与另外一幅医学图像或者多幅图像上的对应点达到空间上的一致。这种一致是建立两幅图像中相同解剖结构的匹配对应关系。通常配准算法的策略是，假定两幅待配准的图像分别为M和F，定义一个目标函数(也称能量函数)E＝E(M,T(F))，其中T代表对图像F作空间变换，一般认为E值的大小代表了图像的配准程度，对于特定的配准任务，图像M和F是确定的，因此E是空间变换T的函数，即E＝E(T)。这样配准的问题就转化为寻求一定的空间变换T使E＝E(T)取最优值得问题。根据E的定义不同，最优值对应为最大值或最小值，在定义了E之后，还必须解决一个问题，即如何找到空间变换T，使得E取得最优值。在非刚性配准的过程中，如果输入的图像具有一定的关系，那么得到的变形场也应该具有相对应的关系。例如将图像B配到图像A的变形场为fBA，将图像C配到图像B的变形场为fCB，将图像C配到图像A的变形场为fCA，那么变形场fCA＝fCB+T(fBA),其中T代表一种空间变换，对于非刚性配准任务来讲，这种空间变换关系在非刚性配准模型得到的变形场中得以体现，即变形场应该具有形变一致性。
近年来，随着计算机计算能力的大幅增强，具备了处理海量数据的能力，人工智能神经网络技术在实际应用中发挥了巨大的作用，同时深度学习也成为了人工智能领域中的核心话题。在医疗健康领域，深度学习在医学图像的应用主要表现为辅助诊断或者影像检测。基于深度学习的医学图像配准直接从大量的训练样本中通过不断迭代的方式学习数据的特征，通过组合低层特征形成更加抽象的高层特征，从而得到网络模型实现图像的配准。对于基于深度学习的医学图像配准而言，T为网络结构，E为损失函数，E的优化问题为根据损失函数的损失值通过反向传播算法迭代更新网络结构的相关参数。
基于深度学习的非刚性医学图像配准方法主要分为两大类：监督学习和非监督学习。然而监督的学习方法需要大量的带标注的数据，医学图像配准任务很难获取金标准数据。所以，基于深度学习的非监督方法成为医学图像非刚性配准任务的首选，其常用的网络结构为基于卷积神经网络的UNet网络结构，通过回归学习，实现待配准图像对到形变场的映射。另外，为了得到更加平滑的变形场，传统技术通常将变形场的一阶梯度作为正则损失函数，指导网络训练过程中的参数的更新，使得学习出的变形场更加平滑。但是这种控制变形场平滑的方法只考虑了变形场局部的平滑性，无法避免变形场出现交叉的情况，进行图像配准的结果准确度较低，也无法从根本上保证非刚性配准模型得到的变形场具有形变一致性。本申请提供的图像配准模型的训练方法、图像配准方法和计算机设备，旨在解决上述技术问题。
为了使本申请的目的、技术方案及优点更加清楚明白，通过下述实施例并结合附图，对本申请实施例中的技术方案进行进一步的详细说明。应当理解，此处所描述的具体实施例仅用以解释本申请，并不用于限定本申请。
需要说明的是，下述方法实施例的执行主体可以是图像配准模型的训练装置，该装置可以通过软件、硬件或者软硬件结合的方式实现成为计算机设备的部分或者全部。下述方法实施例以执行主体为计算机设备为例进行说明，该计算机设备可以是终端，也可以是服务器，可以为单独的计算设备，也可以集成于医学成像设备上，只要能完成图像配准模型的训练即可，本实施例对此不做限定。
图1为一个实施例提供的图像配准模型的训练方法的流程示意图。本实施例涉及的是计算机设备利用训练样本图像对初始图像配准网络进行训练，得到图像配准模型的具体过程。如图1所示，该方法包括：
S101，获取n个训练样本图像；n≥3。
具体的，计算机设备首先需要获取大量的训练样本图像，本实施例所指出的训练样本图像的个数大于等于3个，是基于变形场的形变一致性需至少3个图像才可满足而言，并不代表训练样本图像的总个数。其中，可以将大量的训练样本图像分为多组，每组包括n个训练样本图像，然后对每组训练样本图像都执行以下方法步骤。
可选的，训练样本图像可以为核磁共振成像(Nuclear Magnetic ResonanceImaging，MRI)、正电子发射型计算机断层成像(Positron Emission ComputedTomography，PET)及电子计算机断层扫描成像(Computed Tomography，CT)等，可以是脑部图像、胸部图像、腹部图像等，但值得说明的是，n个训练样本图像需为同一类型的图像。可选的，计算机设备获取训练样本图像的方式可以为从计算机设备的存储器中直接调取，也可以为从影像归档和通信系统(Picture Archiving and Communication Systems，PACS)中获取，本实施例对此不做限制。
S102，将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场。
具体的，计算机设备将上述n个训练样本图像输入初始图像配准网络进行循环两两图像配准，可以得到n个变形场。其中，图像配准过程为一个图像作为浮动图像、另一个图像作为参考图像输入初始图像配准网络，得到对应的变形场。
其中，循环两两图像配准是指：假设n个训练样本图像的编号分别为1、2、3...n，可以将1和2、2和3、3和4、4和5、...、n-1和n、n和1进行两两图像配准；但值得说明的是，本实施例并不限于是编号相邻的两个训练样本图像进行配准，只要前后两次图像配准时，前一次的浮动图像为下一次的参考图像，将n个训练样本图像进行不重复的两两图像配准后，再将最后一个训练样本图像和第一个训练样本图像进行图像配准即可。示例性的，假设有A、B、C、D四个训练样本图像，可以进行(AB、BC、CD、AD)的循环配准，得到四个对应的变形场；也可以进行(AC、CD、DB、AB)的循环配准，得到四个对应的变形场。
可选的，上述初始图像配准网络可以为神经网络模型，可以为卷积神经网络(Convolutional Neural Networks，CNN)、全卷积网络(Fully Convolutional Networks，FCN)，也可以为其他网络模型。可选的，该初始图像配准网络为UNet网络，整个网络结构分为压缩路径和非压缩路径两部分，即缩小和扩大图像的特征映射，并将缩小端的低层特征送入扩大端的相应位置，为模型训练提供更多的有效信息，加速训练过程中模型的收敛。
S103，根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失，以及根据n个变形场计算信息约束损失。
具体的，计算机设备可以根据第n个变形场和第n个变形场对应的两个训练样本图像计算图像配准损失，此处的第n个变形场为初始图像配准网络最后一次进行图像配准得到的变形场，依据的是输入初始图像配准网络的次序，并不特指是哪个变形场。对于上述(AB、BC、CD、AD)的示例，则第n个变形场为对图像A和D进行配准得到的fAD，那么根据fAD和训练样本图像A可得到配准之后的图像，然后将配准之后的图像与训练样本图像D进行比较，得到图像配准损失。可选的，图像配准损失的损失函数可以为均方误差损失(MSE Loss)或者归一化交叉相关损失(Normalization cross correlation，NCC Loss)。
然后，计算机设备还根据得到的n个变形场计算信息约束损失，假设n为3时，得到的变形场为fBA、fCB和fCA，计算机设备将fBA和fCB进行叠加后得到fCA’，理论上fCA’应与fCA相等，即fCA’＝fCA。但因此时的初始图像配准网络还未收敛，fCA’与fCA有一定的差异，利用该差异可得到对应的信息约束损失。
S104，根据图像配准损失和信息约束损失计算网络损失，利用网络损失对初始图像配准网络进行训练；当网络损失达到满足预设条件时，初始图像配准网络训练完成，得到图像配准模型。
具体的，计算机设备可以根据图像配准损失和信息约束损失计算网络损失，可选的，可以将图像配准损失和信息约束损失进行直接求和，或者平均求和，或者加权求和，得到网络损失。当为加权求和时，每个损失都有各自不同的权重，权重的大小代表该损失的重要程度；权重越大，该损失在训练过程中起的作用越大。然后计算机设备利用得到的网络损失调整初始图像配准网络的网络参数，可选的，可以采用反向梯度传播的方法调整网络参数。
以上述方式迭代训练初始图像配准网络，当某次迭代过程的网络损失满足预设条件时，表征初始图像配准网络训练完成，便得到了收敛的图像配准模型。可选的，预设条件可以为网络损失达到收敛，也可以为网络损失小于或者等于预设阈值，对此本实施例不做限制。
本实施例提供的图像配准模型的训练方法，计算机设备将获取的n个训练样本图像输入初始图像配准网络进行循环两两配准，得到n个变形场；然后根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失，以及根据n个变形场计算信息约束损失；最后根据图像配准损失和信息约束损失计算网络损失，利用网络损失对初始图像配准网络进行训练，直至得到收敛的图像配准模型。该方法中，除考虑图像配准时的配准损失，保证了图像配准模型的精度，还考虑了各变形场之间的形变联系，即前n-1个变形场的叠加与第n个变形场相等的性质，由此进一步提高根据图像配准模型得到的变形场的准确性，也使得变形场具有形变一致性。
可选的，在其中一些实施例中，上述方法还包括：根据第n个变形场计算变形场平滑损失；根据图像配准损失、信息约束损失和变形场平滑损失计算网络损失。
可选的，计算机设备可以计算第n个变形场一阶梯度的所有像素点的平方和，将该平方和作为变形场平滑损失。可选的，计算机设备可以将图像配准损失、信息约束损失和变形场平滑损失进行直接求和，或者平均求和，或者加权求和，得到上述网络损失。由此本申请还考虑了变形场的局部平滑性，进一步提高了得到的图像配准模型的精度。
可选的，在其中一些实施例中，将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场，包括：将第i个训练样本图像和第i+1个训练样本图像输入初始图像配准网络进行图像配准，得到n-1个变形场；其中，1≤i≤n-1；将第1个训练样本图像和第n个训练样本图像输入初始图像配准网络进行图像配准，得到第n个变形场。
示例性的，当n＝3时，计算机设备可以将第1个训练样本图像作为参考图像，第2个训练样本图像作为浮动图像输入初始图像配准网络进行图像配准，得到第一个变形场；将第2个训练样本图像作为参考图像，第3个训练样本图像作为浮动图像输入初始图像配准网络进行图像配准，得到第二个变形场；将第1个训练样本图像作为参考图像，第3个训练样本图像作为浮动图像输入初始图像配准网络进行图像配准，得到第三个变形场。
图2为另一个实施例提供的图像配准模型的训练方法的流程示意图。本实施例涉及的是计算机设备根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失的具体过程。在上述实施例的基础上，可选的，如图2所示，S103可以包括：
S201，利用第n个变形场对第n个变形场对应的第一训练样本图像进行图像空间变换，得到配准图像。
具体的，计算机设备利用第n个变形场对其对应的第一训练样本图像进行图像空间变换，即根据第n个变形场提供的坐标位置对应关系，利用插值算法对第一训练样本图像进行计算，得到配准图像；需要说明的是，该第一训练样本图像为得到第n个变形场时的浮动图像。
S202，计算配准图像与第n个变形场对应的第二训练样本图像之间的差异值，将差异值作为图像配准损失。
具体的，计算机设备计算配准图像与第二训练样本图像(该第二训练样本图像为得到第n个变形场时的参考图像)之间的差异值，可选的，可以计算配准图像与第二训练样本图像对应坐标位置的像素差异值，将该差异值作为图像配准损失。对于上述示例，第n个变形场为第三个变形场，第一训练样本图像为第3个训练样本图像，第二训练样本图像为第1个训练样本图像。
本实施例提供的图像配准模型的训练方法，计算机设备利用第n个变形场对第一训练样本图像进行图像空间变换，得到配准图像，然后计算配准图像与第二训练样本图像之间的差异值，将差异值作为图像配准损失。通过计算参考图像与配准图像之间的图像配准损失，在初始图像配准网络更新过程中作为部分损失进行梯度回传，可进一步提高得到的图像配准模型的配准精度。
图3为又一个实施例提供的图像配准模型的训练方法的流程示意图。本实施例涉及的是计算机设备根据n个变形场计算信息约束损失的具体过程。在上述实施例的基础上，可选的，如图3所示，S103可以包括：
S301，将第1个变形场至第n-1个变形场进行变形场叠加，得到叠加后的变形场。
S302，根据叠加后的变形场和第n个变形场计算信息约束损失。
具体的，计算机设备可以先将第2个变形场至第n-1个变形场通过变形场空间变换，都变换到第1个变形场的空间坐标系下，再与第1个变形场进行叠加，得到叠加后的变形场。假设训练样本图像为A、B、C，得到的变形场为fBA、fCB和fCA，可选的，将fCB变换到fBA空间坐标系下的过程可以为：确定fBA中像素点a(x,y)的指向位置b(x+dx,y+dy)，利用插值算法计算fCB在指向位置b处的值c；将c的值作为空间变换后的变形场fCB’在像素点a处的值，最终可得到空间变换后的变形场fCB’。然后将fCB’与fBA进行叠加，得到叠加后的变形场fCA’。
然后，计算机设备根据叠加后的变形场和第n个变形场计算信息约束损失，如计算fCA’与fCA之间的损失。可选的，可以将叠加后的变形场与第n个变形场中对应位置的像素点值进行求差，得到求差结果，将求差结果作为信息约束损失，即将fCA’与fCA进行求差，根据fCA'-fCA的关系式计算该信息约束损失。关于变形场叠加过程的示意图可以参见图3a所示。关于n＝3时图像配准模型的训练过程可以参见图3b所示的示意图。
本实施例提供的图像配准模型的训练方法，计算机设备将第1个变形场至第n-1个变形场进行变形场叠加，得到叠加后的变形场，根据叠加后的变形场和第n个变形场计算信息约束损失。通过将第1个变形场至第n-1个变形场进行叠加来计算与第n个变形场之间的信息约束损失，可提高计算的信息约束损失的准确性，在初始图像配准网络更新过程中作为部分损失进行梯度回传，可加快模型的收敛，大大提高图像配准模型的变形场准确性及形变一致性。
在上述图像配准模型训练完成之后，便可以利用该模型进行图像配准，图4为一个实施例提供的图像配准方法的流程示意图，该方法包括：
S401，获取待配准图像和参考图像。
S402，将待配准图像和参考图像输入图像配准模型进行图像配准，得到待配准图像的配准结果；其中，图像配准模型的训练方式包括：
获取n个训练样本图像；n≥3；
将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场；
根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失，以及根据n个变形场计算信息约束损失；
根据图像配准损失和信息约束损失计算网络损失，利用网络损失对初始图像配准网络进行训练；当网络损失达到满足预设条件时，初始图像配准网络训练完成，得到图像配准模型。
具体的，计算机设备获取待配准图像和参考图像后，将待配准图像作为浮动图像与参考图像输入图像配准模型进行图像配准，可以得到待配准图像的配准结果。可选的，该配准结果可以为待配准图像到参考图像的变形场。而关于该图像配准模型的训练过程，可以参见上述实施例所示的方法，其实现过程与原理类似，在此不再赘述。
本实施例提供的图像配准方法，所采用的图像配准模型具有一定的变形场形变一致性及较高的精度，使得得到的待配准图像的配准结果准确度也较高，同时也具有形变一致性。
应该理解的是，虽然图1-图4的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图1-图4中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
图5为一个实施例提供的图像配准模型的训练装置的结构示意图。如图5所示，该装置包括：第一获取模块11、第一图像配准模块12、计算模块13和训练模块14。
具体的，第一获取模块11，用于获取n个训练样本图像；n≥3。
第一图像配准模块12，用于将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场。
计算模块13，用于根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失，以及根据n个变形场计算信息约束损失。
训练模块14，用于根据图像配准损失和信息约束损失计算网络损失，利用网络损失对初始图像配准网络进行训练；当网络损失达到满足预设条件时，初始图像配准网络训练完成，得到图像配准模型。
本实施例提供的图像配准模型的训练装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。
在其中一个实施例中，计算模块13，还用于根据第n个变形场计算变形场平滑损失；根据图像配准损失、信息约束损失和变形场平滑损失计算网络损失。
在其中一个实施例中，第一图像配准模块12，具体用于将第i个训练样本图像和第i+1个训练样本图像输入初始图像配准网络进行图像配准，得到n-1个变形场；其中，1≤i≤n-1；以及将第1个训练样本图像和第n个训练样本图像输入初始图像配准网络进行图像配准，得到第n个变形场。
在其中一个实施例中，计算模块13包括图像空间变换单元和计算单元；图像空间变换单元，用于利用第n个变形场对第n个变形场对应的第一训练样本图像进行图像空间变换，得到配准图像；计算单元，用于计算配准图像与第n个变形场对应的第二训练样本图像之间的差异值，将差异值作为图像配准损失。
在其中一个实施例中，计算模块13还包括变形场叠加单元；变形场叠加单元，用于将第1个变形场至第n-1个变形场进行变形场叠加，得到叠加后的变形场；计算单元，还用于根据叠加后的变形场和第n个变形场计算信息约束损失。
在其中一个实施例中，计算单元，具体用于将叠加后的变形场与第n个变形场中对应位置的像素点值进行求差，得到求差结果，将求差结果作为信息约束损失。
在其中一个实施例中，当n＝3时，第一图像配准模块12，具体用于将第1个训练样本图像作为参考图像，第2个训练样本图像作为浮动图像输入初始图像配准网络进行图像配准，得到第一个变形场；将第2个训练样本图像作为参考图像，第3个训练样本图像作为浮动图像输入初始图像配准网络进行图像配准，得到第二个变形场；将第1个训练样本图像作为参考图像，第3个训练样本图像作为浮动图像输入初始图像配准网络进行图像配准，得到第三个变形场。
在其中一个实施例中，计算模块13，具体用于计算第n个变形场一阶梯度的所有像素点值的平方和，将平方和作为变形场平滑损失。
图6为一个实施例提供的图像配准装置的结构示意图。如图6所示，该装置包括：第二获取模块21和第二图像配准模块22。
具体的，第二获取模块21，用于获取待配准图像和参考图像。
第二图像配准模块22，用于将待配准图像和参考图像输入图像配准模型进行图像配准，得到待配准图像的配准结果。其中图像配准模型的训练过程可以参见上述图像配准模型的训练装置实施例的实现过程。
本实施例提供的图像配准装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。
关于图像配准模型的训练装置和图像配准装置的具体限定可以参见上文中对于图像配准模型的训练方法和图像配准方法的限定，在此不再赘述。上述图像配准模型的训练装置和图像配准装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图7所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种图像配准模型的训练方法或图像配准方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
本领域技术人员可以理解，图7中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：
获取n个训练样本图像；n≥3；
将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场；
根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失，以及根据n个变形场计算信息约束损失；
根据图像配准损失和信息约束损失计算网络损失，利用网络损失对初始图像配准网络进行训练；当网络损失达到满足预设条件时，初始图像配准网络训练完成，得到图像配准模型。
本实施例提供的计算机设备，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
根据第n个变形场计算变形场平滑损失；
根据图像配准损失、信息约束损失和变形场平滑损失计算网络损失。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
将第i个训练样本图像和第i+1个训练样本图像输入初始图像配准网络进行图像配准，得到n-1个变形场；其中，1≤i≤n-1；
将第1个训练样本图像和第n个训练样本图像输入初始图像配准网络进行图像配准，得到第n个变形场。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
利用第n个变形场对第n个变形场对应的第一训练样本图像进行图像空间变换，得到配准图像；
计算配准图像与第n个变形场对应的第二训练样本图像之间的差异值，将差异值作为图像配准损失。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
将第1个变形场至第n-1个变形场进行变形场叠加，得到叠加后的变形场；
根据叠加后的变形场和第n个变形场计算信息约束损失。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
将叠加后的变形场与第n个变形场中对应位置的像素点值进行求差，得到求差结果，将求差结果作为信息约束损失。
在一个实施例中，当n＝3时，处理器执行计算机程序时还实现以下步骤：
将第1个训练样本图像作为参考图像，第2个训练样本图像作为浮动图像输入初始图像配准网络进行图像配准，得到第一个变形场；
将第2个训练样本图像作为参考图像，第3个训练样本图像作为浮动图像输入初始图像配准网络进行图像配准，得到第二个变形场；
将第1个训练样本图像作为参考图像，第3个训练样本图像作为浮动图像输入初始图像配准网络进行图像配准，得到第三个变形场。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
计算第n个变形场一阶梯度的所有像素点值的平方和，将平方和作为变形场平滑损失。
在一个实施例中，还提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：
获取待配准图像和参考图像；
将待配准图像和参考图像输入图像配准模型进行图像配准，得到待配准图像的配准结果；其中，图像配准模型的训练方式包括：
获取n个训练样本图像；n≥3；
将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场；
根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失，以及根据n个变形场计算信息约束损失；
根据图像配准损失和信息约束损失计算网络损失，利用网络损失对初始图像配准网络进行训练；当网络损失达到满足预设条件时，初始图像配准网络训练完成，得到图像配准模型。
本实施例提供的计算机设备，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
获取n个训练样本图像；n≥3；
将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场；
根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失，以及根据n个变形场计算信息约束损失；
根据图像配准损失和信息约束损失计算网络损失，利用网络损失对初始图像配准网络进行训练；当网络损失达到满足预设条件时，初始图像配准网络训练完成，得到图像配准模型。
本实施例提供的计算机可读存储介质，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
根据第n个变形场计算变形场平滑损失；
根据图像配准损失、信息约束损失和变形场平滑损失计算网络损失。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
将第i个训练样本图像和第i+1个训练样本图像输入初始图像配准网络进行图像配准，得到n-1个变形场；其中，1≤i≤n-1；
将第1个训练样本图像和第n个训练样本图像输入初始图像配准网络进行图像配准，得到第n个变形场。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
利用第n个变形场对第n个变形场对应的第一训练样本图像进行图像空间变换，得到配准图像；
计算配准图像与第n个变形场对应的第二训练样本图像之间的差异值，将差异值作为图像配准损失。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
将第1个变形场至第n-1个变形场进行变形场叠加，得到叠加后的变形场；
根据叠加后的变形场和第n个变形场计算信息约束损失。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
将叠加后的变形场与第n个变形场中对应位置的像素点值进行求差，得到求差结果，将求差结果作为信息约束损失。
在一个实施例中，当n＝3时，计算机程序被处理器执行时还实现以下步骤：
将第1个训练样本图像作为参考图像，第2个训练样本图像作为浮动图像输入初始图像配准网络进行图像配准，得到第一个变形场；
将第2个训练样本图像作为参考图像，第3个训练样本图像作为浮动图像输入初始图像配准网络进行图像配准，得到第二个变形场；
将第1个训练样本图像作为参考图像，第3个训练样本图像作为浮动图像输入初始图像配准网络进行图像配准，得到第三个变形场。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
计算第n个变形场一阶梯度的所有像素点值的平方和，将平方和作为变形场平滑损失。
在一个实施例中，还提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
获取待配准图像和参考图像；
将待配准图像和参考图像输入图像配准模型进行图像配准，得到待配准图像的配准结果；其中，图像配准模型的训练方式包括：
获取n个训练样本图像；n≥3；
将n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场；
根据第n个变形场和第n个变形场对应的训练样本图像计算图像配准损失，以及根据n个变形场计算信息约束损失；
根据图像配准损失和信息约束损失计算网络损失，利用网络损失对初始图像配准网络进行训练；当网络损失达到满足预设条件时，初始图像配准网络训练完成，得到图像配准模型。
本实施例提供的计算机可读存储介质，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。

Claims (10)

1.一种图像配准模型的训练方法，其特征在于，包括：
获取n个训练样本图像；所述n≥3；
将所述n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场；
根据第n个变形场和所述第n个变形场对应的训练样本图像计算图像配准损失，以及将第1个变形场至第n-1个变形场进行变形场叠加，得到叠加后的变形场，将所述叠加后的变形场与所述第n个变形场中对应位置的像素点值进行求差，得到求差结果，将所述求差结果作为信息约束损失；
根据所述图像配准损失和所述信息约束损失计算网络损失，利用所述网络损失对所述初始图像配准网络进行训练；当所述网络损失达到满足预设条件时，所述初始图像配准网络训练完成，得到图像配准模型。
2.根据权利要求1所述的方法，其特征在于，所述方法还包括：
根据所述第n个变形场计算变形场平滑损失；
根据所述图像配准损失、所述信息约束损失和所述变形场平滑损失计算所述网络损失。
3.根据权利要求1或2所述的方法，其特征在于，所述将所述n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场，包括：
将第i个训练样本图像和第i+1个训练样本图像输入所述初始图像配准网络进行图像配准，得到n-1个变形场；其中，1≤i≤n-1；
将第1个训练样本图像和第n个训练样本图像输入所述初始图像配准网络进行图像配准，得到第n个变形场。
4.根据权利要求1所述的方法，其特征在于，所述根据第n个变形场和所述第n个变形场对应的训练样本图像计算图像配准损失，包括：
利用所述第n个变形场对所述第n个变形场对应的第一训练样本图像进行图像空间变换，得到配准图像；
计算所述配准图像与所述第n个变形场对应的第二训练样本图像之间的差异值，将所述差异值作为所述图像配准损失。
5.根据权利要求2所述的方法，其特征在于，所述根据所述图像配准损失、所述信息约束损失和所述变形场平滑损失计算所述网络损失，包括：
将图像配准损失、信息约束损失和变形场平滑损失进行直接求和，或者平均求和，或者加权求和，得到所述网络损失。
6.根据权利要求4所述的方法，其特征在于，所述利用所述第n个变形场对所述第n个变形场对应的第一训练样本图像进行图像空间变换，得到配准图像，包括：
根据所述第n个变形场提供的坐标位置对应关系，利用插值算法对所述第一训练样本图像进行计算，得到所述配准图像。
7.根据权利要求3所述的方法，其特征在于，当n＝3时，所述将所述n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场，包括：
将第1个训练样本图像作为参考图像，第2个训练样本图像作为浮动图像输入所述初始图像配准网络进行图像配准，得到第一个变形场；
将所述第2个训练样本图像作为参考图像，第3个训练样本图像作为浮动图像输入所述初始图像配准网络进行图像配准，得到第二个变形场；
将所述第1个训练样本图像作为参考图像，所述第3个训练样本图像作为浮动图像输入所述初始图像配准网络进行图像配准，得到第三个变形场。
8.根据权利要求2所述的方法，其特征在于，所述根据所述第n个变形场计算变形场平滑损失，包括：
计算所述第n个变形场一阶梯度的所有像素点值的平方和，将所述平方和作为所述变形场平滑损失。
9.一种图像配准方法，其特征在于，包括：
获取待配准图像和参考图像；
将所述待配准图像和所述参考图像输入图像配准模型进行图像配准，得到所述待配准图像的配准结果；其中，所述图像配准模型的训练方式包括：
获取n个训练样本图像；所述n≥3；
将所述n个训练样本图像输入初始图像配准网络进行循环两两图像配准，得到n个变形场；
根据第n个变形场和所述第n个变形场对应的训练样本图像计算图像配准损失，以及将第1个变形场至第n-1个变形场进行变形场叠加，得到叠加后的变形场，将所述叠加后的变形场与所述第n个变形场中对应位置的像素点值进行求差，得到求差结果，将所述求差结果作为信息约束损失；
根据所述图像配准损失和所述信息约束损失计算网络损失，利用所述网络损失对所述初始图像配准网络进行训练；当所述网络损失达到满足预设条件时，所述初始图像配准网络训练完成，得到图像配准模型。
10.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求9所述方法的步骤。
CN201911063804.5A 2019-11-04 2019-11-04 图像配准模型的训练方法、图像配准方法和计算机设备 Active CN110827331B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911063804.5A CN110827331B (zh) 2019-11-04 2019-11-04 图像配准模型的训练方法、图像配准方法和计算机设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911063804.5A CN110827331B (zh) 2019-11-04 2019-11-04 图像配准模型的训练方法、图像配准方法和计算机设备

Publications (2)

Publication Number Publication Date
CN110827331A CN110827331A (zh) 2020-02-21
CN110827331B true CN110827331B (zh) 2022-07-26

Family

ID=69552427

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911063804.5A Active CN110827331B (zh) 2019-11-04 2019-11-04 图像配准模型的训练方法、图像配准方法和计算机设备

Country Status (1)

Country Link
CN (1) CN110827331B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111768436B (zh) * 2020-06-17 2022-10-18 哈尔滨理工大学 一种基于Faster-RCNN改进的图像特征块配准方法
CN114519729A (zh) * 2020-11-20 2022-05-20 腾讯科技（深圳）有限公司 图像配准质量评估模型训练方法、装置和计算机设备
CN114359356A (zh) * 2021-12-28 2022-04-15 上海联影智能医疗科技有限公司 图像配准模型的训练方法、图像配准方法、设备及介质
CN114511599B (zh) * 2022-01-20 2022-09-20 推想医疗科技股份有限公司 模型训练方法及其装置、医学图像配准方法及其装置
CN115393402B (zh) * 2022-08-24 2023-04-18 北京医智影科技有限公司 图像配准网络模型的训练方法、图像配准方法及设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101051386A (zh) * 2007-05-23 2007-10-10 北京航空航天大学 多幅深度图像的精确配准方法
CN104599268A (zh) * 2015-01-06 2015-05-06 广州医科大学附属肿瘤医院 一种联合点配准的局部区域精确变形配准算法
CN109472817A (zh) * 2018-09-27 2019-03-15 浙江工业大学 一种基于循环生成对抗网络的多序列磁共振图像配准方法
CN109801319A (zh) * 2019-01-03 2019-05-24 杭州电子科技大学 基于Hadoop并行加速的分级图分组配准方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170337682A1 (en) * 2016-05-18 2017-11-23 Siemens Healthcare Gmbh Method and System for Image Registration Using an Intelligent Artificial Agent
US11449759B2 (en) * 2018-01-03 2022-09-20 Siemens Heathcare Gmbh Medical imaging diffeomorphic registration based on machine learning

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101051386A (zh) * 2007-05-23 2007-10-10 北京航空航天大学 多幅深度图像的精确配准方法
CN104599268A (zh) * 2015-01-06 2015-05-06 广州医科大学附属肿瘤医院 一种联合点配准的局部区域精确变形配准算法
CN109472817A (zh) * 2018-09-27 2019-03-15 浙江工业大学 一种基于循环生成对抗网络的多序列磁共振图像配准方法
CN109801319A (zh) * 2019-01-03 2019-05-24 杭州电子科技大学 基于Hadoop并行加速的分级图分组配准方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
DGR-Net- Deep Groupwise Registration of Multispectral Images;Tongtong Che等;《International Conference on Information Processing》;20190522;第11492卷;第706-717页 *
Inverse-consistent Deep Networks for Unsupervised Deformable Image Registration;Jun Zhang;《Computer Vision and Pattern Recognition》;20180910;第1-13页 *
基于流形学习的图像配准研究;刘彩凤;《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》;20170515;第I138-902页 *

Also Published As

Publication number Publication date
CN110827331A (zh) 2020-02-21

Similar Documents

Publication Publication Date Title
CN110838139B (zh) 图像配准模型的训练方法、图像配准方法和计算机设备
CN110827331B (zh) 图像配准模型的训练方法、图像配准方法和计算机设备
US12094116B2 (en) Method and system for image registration using an intelligent artificial agent
CN109767461B (zh) 医学影像配准方法、装置、计算机设备和存储介质
Liao et al. An artificial agent for robust image registration
Yang et al. Non-rigid multi-modal medical image registration by combining L-BFGS-B with cat swarm optimization
CN109993726A (zh) 医学图像的检测方法、装置、设备和存储介质
JP2021535482A (ja) 深層学習ベースのコレジストレーション
CN110866909B (zh) 图像生成网络的训练方法、图像预测方法和计算机设备
CN111161270A (zh) 医学影像的血管分割方法、计算机设备和可读存储介质
Pei et al. Non-rigid craniofacial 2D-3D registration using CNN-based regression
CN114298234B (zh) 脑部医学影像分类方法、装置、计算机设备和存储介质
CN110751187B (zh) 异常区域图像生成网络的训练方法和相关产品
Duan et al. Unsupervised learning for deformable registration of thoracic CT and cone‐beam CT based on multiscale features matching with spatially adaptive weighting
CN111223158B (zh) 心脏冠脉图像的伪影校正方法和可读存储介质
Ginn et al. Multislice motion modeling for MRI‐guided radiotherapy gating
Polo et al. Improving hybrid image and structure-based deformable image registration for large internal deformations
CN117078589A (zh) 一种基于多智能体强化学习的自监督神经元分割方法
US20230099906A1 (en) Image registration method, computer device, and storage medium
CN113850710B (zh) 一种跨模态医学影像精准转换方法
Zhang Image registration
Langner Machine Learning Techniques for MRI Data Processing at Expanding Scale
Moreno et al. Combining a breathing model and tumor-specific rigidity constraints for registration of CT-PET thoracic data
Shaikh A survey on medical image analysis using deep reinforcement learning approach
US20240046453A1 (en) Semi-supervised learning leveraging cross-domain data for medical imaging analysis

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