WO2022178997A1 - Medical image registration method and apparatus, computer device, and storage medium - Google Patents

Medical image registration method and apparatus, computer device, and storage medium Download PDF

Info

Publication number
WO2022178997A1
WO2022178997A1 PCT/CN2021/096702 CN2021096702W WO2022178997A1 WO 2022178997 A1 WO2022178997 A1 WO 2022178997A1 CN 2021096702 W CN2021096702 W CN 2021096702W WO 2022178997 A1 WO2022178997 A1 WO 2022178997A1
Authority
WO
WIPO (PCT)
Prior art keywords
image
registration
medical
training
transformation model
Prior art date
Application number
PCT/CN2021/096702
Other languages
French (fr)
Chinese (zh)
Inventor
李雷来
王健宗
瞿晓阳
Original Assignee
平安科技(深圳)有限公司
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 平安科技(深圳)有限公司 filed Critical 平安科技(深圳)有限公司
Publication of WO2022178997A1 publication Critical patent/WO2022178997A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • 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
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration by the use of more than one image, e.g. averaging, subtraction
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00ICT specially adapted for the handling or processing of medical images
    • G16H30/20ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
    • 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
    • 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/20212Image combination
    • G06T2207/20221Image fusion; Image merging

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Computing Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Biomedical Technology (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Computational Linguistics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Epidemiology (AREA)
  • Medical Informatics (AREA)
  • Primary Health Care (AREA)
  • Public Health (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Processing (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

A medical image registration method and apparatus, a computer device, and a storage medium. The method comprises: obtaining a first medical image set and a second medical image set (S10); after inputting the first medical image set into a first registration transformation model, obtaining a first registered image outputted by the first registration transformation model for multi-view-angle fusion of the first medical image set (S20); after inputting the second medical image set into a second registration transformation model, obtaining a second registered image outputted by the second registration transformation model for multi-time-sequence fusion of the second medical image set (S30); and fusing the first registered image and the second registered image to obtain a third registered image (S40). The method replaces manual determination by means of an artificial intelligence model, thereby improving the determination rate and the determination accuracy.

Description

医学影像配准方法、装置、计算机设备及存储介质Medical image registration method, device, computer equipment and storage medium
本申请要求于2021年02月25日提交中国专利局、申请号为202110210897.0,发明名称为“医学影像配准方法、装置、计算机设备及存储介质”的中国专利申请的优先权,其全部内容通过引用结合在本申请中。This application claims the priority of the Chinese patent application with the application number 202110210897.0 and the invention titled "Medical Image Registration Method, Device, Computer Equipment and Storage Medium", which was filed with the China Patent Office on February 25, 2021, the entire content of which is approved by Reference is incorporated in this application.
技术领域technical field
本申请涉及数字医疗的AI医疗技术领域,应用于数字医疗中的智慧医疗,尤其揭露了一种医学影像配准方法、装置、计算机设备及存储介质。The present application relates to the field of AI medical technology in digital medical treatment, and is applied to smart medical treatment in digital medical treatment, and particularly discloses a medical image registration method, device, computer equipment and storage medium.
背景技术Background technique
医疗科技作为目前人工智能技术的重要研究领域之一,在医学影像、临床决策支持、药物研发以及病理学等方面发挥出了巨大作用,医疗与科技的融合能够有助于医生快速准确地作出诊断,真正地实现“早诊断,早治疗,早康复”的医学目标。其中,医学影像配准在医疗科技中具有很大的实用价值。As one of the important research fields of artificial intelligence technology, medical technology has played a huge role in medical imaging, clinical decision support, drug research and development, and pathology. The integration of medical technology and technology can help doctors make diagnoses quickly and accurately , and truly achieve the medical goal of "early diagnosis, early treatment, and early recovery". Among them, medical image registration has great practical value in medical technology.
目前,随着医学成像设备的不断更新与发展,对于同一患者,可以通过多种成像技术如CT、MRI等去多角度多方面采集患者的病理细胞结构信息。发明人发现以往的传统诊断手段,往往需要一些经验丰富的医生,通过对不同图像的观察以及多方面综合分析,以及结合主观经验和空间想象去作出相关的预测结果,因此发明人发现现有技术中还是过多依赖于医生的主观经验和空间想象去作出判断,现有技术可能存在判断速率和判断准确率不高的问题。At present, with the continuous updating and development of medical imaging equipment, for the same patient, multiple imaging techniques such as CT and MRI can be used to collect the patient's pathological cell structure information from multiple angles and multiple aspects. The inventor found that the traditional diagnostic methods in the past often required some experienced doctors to make relevant prediction results through the observation of different images and comprehensive analysis of various aspects, as well as combining subjective experience and spatial imagination. Therefore, the inventor discovered the prior art. It still relies too much on the subjective experience and spatial imagination of doctors to make judgments, and the existing technology may have the problem of low judgment speed and judgment accuracy.
发明内容SUMMARY OF THE INVENTION
基于此,有必要针对上述技术问题,提供一种医学影像配准方法、装置、计算机设备及存储介质,用于通过人工智能模型代替人工判断,提高判断速率和判断准确率。Based on this, it is necessary to provide a medical image registration method, device, computer equipment and storage medium for the above-mentioned technical problems, which can be used to replace manual judgment with an artificial intelligence model and improve the judgment rate and judgment accuracy.
一种医学影像配准方法,包括:A medical image registration method, comprising:
获取第一医学影像集和第二医学影像集;所述第一医学影像集中包含同一时刻不同视角的多张第一医学影像;所述第二医学影像集中包含不同时刻同一视角的多张第二医学影像;Acquiring a first medical image set and a second medical image set; the first medical image set includes multiple first medical images from different viewing angles at the same time; the second medical image set includes multiple second medical images from the same viewing angle at different times Medical Imaging;
将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像;After the first medical image set is input into the first registration transformation model, a first registration image output by the first registration transformation model for performing multi-view fusion on the first medical image set is obtained;
将所述第二医学影像集输入至第二配准变换模型后,获取所述第二配准变换模型输出的对所述第二医学影像集进行多时序融合的第二配准图像;After the second medical image set is input into the second registration transformation model, a second registration image output by the second registration transformation model for performing multi-sequence fusion on the second medical image set is obtained;
对所述第一配准图像和所述第二配准图像进行融合,得到第三配准图像;所述第三配准图像反映出所述第一医学影像集和所述第二医学影像集的生长追踪结果。The first registration image and the second registration image are fused to obtain a third registration image; the third registration image reflects the first medical image set and the second medical image set growth tracking results.
一种医学影像配准装置,包括:A medical image registration device, comprising:
第一获取模块,用于获取第一医学影像集和第二医学影像集;所述第一医学影像集中包含同一时刻不同视角的多张第一医学影像;所述第二医学影像集中包含不同时刻同一视角的多张第二医学影像;a first acquisition module, configured to acquire a first medical image set and a second medical image set; the first medical image set includes multiple first medical images from different viewing angles at the same time; the second medical image set includes different times Multiple second medical images from the same viewing angle;
第二获取模块,用于将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像;The second acquisition module is configured to, after the first medical image set is input into the first registration transformation model, obtain the first registration transformation model for performing multi-view fusion on the first medical image set outputted by the first registration transformation model. a registered image;
第三获取模块,用于将所述第二医学影像集输入至第二配准变换模型后,获取所述第二配准变换模型输出的对所述第二医学影像集进行多时序融合的第二配准图像;The third acquiring module is configured to acquire the first multi-sequence fusion of the second medical image set output from the second registration transformation model after the second medical image set is input into the second registration transformation model. Two registered images;
融合模块,用于对所述第一配准图像和所述第二配准图像进行融合,得到第三配准图像;所述第三配准图像反映出所述第一医学影像集和所述第二医学影像集的生长追踪结果。a fusion module, configured to fuse the first registration image and the second registration image to obtain a third registration image; the third registration image reflects the first medical image set and the Growth tracking results for the second set of medical images.
一种计算机设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机可读指令,所述处理器执行所述计算机可读指令时实现如下步骤:A computer device, comprising a memory, a processor, and computer-readable instructions stored in the memory and executable on the processor, the processor implementing the following steps when executing the computer-readable instructions:
获取第一医学影像集和第二医学影像集;所述第一医学影像集中包含同一时刻不同视角的多张第一医学影像;所述第二医学影像集中包含不同时刻同一视角的多张第二医学影像;Acquiring a first medical image set and a second medical image set; the first medical image set includes multiple first medical images from different viewing angles at the same time; the second medical image set includes multiple second medical images from the same viewing angle at different times Medical Imaging;
将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像;After the first medical image set is input into the first registration transformation model, a first registration image output by the first registration transformation model for performing multi-view fusion on the first medical image set is obtained;
将所述第二医学影像集输入至第二配准变换模型后,获取所述第二配准变换模型输出的对所述第二医学影像集进行多时序融合的第二配准图像;After the second medical image set is input into the second registration transformation model, a second registration image output by the second registration transformation model for performing multi-sequence fusion on the second medical image set is obtained;
对所述第一配准图像和所述第二配准图像进行融合,得到第三配准图像;所述第三配准图像反映出所述第一医学影像集和所述第二医学影像集的生长追踪结果。The first registration image and the second registration image are fused to obtain a third registration image; the third registration image reflects the first medical image set and the second medical image set growth tracking results.
一个或多个存储有计算机可读指令的可读存储介质,所述计算机可读指令被一个或多个处理器执行时,使得所述一个或多个处理器执行如下步骤:One or more readable storage media storing computer-readable instructions that, when executed by one or more processors, cause the one or more processors to perform the following steps:
获取第一医学影像集和第二医学影像集;所述第一医学影像集中包含同一时刻不同视角的多张第一医学影像;所述第二医学影像集中包含不同时刻同一视角的多张第二医学影像;Acquiring a first medical image set and a second medical image set; the first medical image set includes multiple first medical images from different viewing angles at the same time; the second medical image set includes multiple second medical images from the same viewing angle at different times Medical Imaging;
将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像;After the first medical image set is input into the first registration transformation model, a first registration image output by the first registration transformation model for performing multi-view fusion on the first medical image set is obtained;
将所述第二医学影像集输入至第二配准变换模型后,获取所述第二配准变换模型输出的对所述第二医学影像集进行多时序融合的第二配准图像;After the second medical image set is input into the second registration transformation model, a second registration image output by the second registration transformation model for performing multi-sequence fusion on the second medical image set is obtained;
对所述第一配准图像和所述第二配准图像进行融合,得到第三配准图像;所述第三配准图像反映出所述第一医学影像集和所述第二医学影像集的生长追踪结果。The first registration image and the second registration image are fused to obtain a third registration image; the third registration image reflects the first medical image set and the second medical image set growth tracking results.
上述医学影像配准方法、装置、计算机设备及存储介质,获取第一医学影像集中包含特定条件(同一时刻不同视角)的第一医学影像以及第二医学影像集中包含特定条件(不同时刻同一视角)的第二医学影像中包含特征的第二医学影像,通过两个配准变换模型分别处理对应的第一医学影像和第二医学影像进行图像配准,可提高模型对特定数据的处理能力,提高数据处理效率;接着通过第一图像插值模块将特定条件下多张第一医学图像融合成一张第一配准图像,通过第二图像插值模块将特定条件下多张第二医学图像融合成一张第二配准图像,对不同视角和不同时刻下的图像进行分别处理,得到每一种图像下的所有图像信息;之后再对第一配准图像和第二配准图像进行融合,得到包含更全面图像信息(生长追踪结果)的第三配准图像;总而言之,本方法借用人工智能模型进行处理,抛弃之前过多依赖于医生的主观经验和空间想象的现象(人工参与工作多,且也过多依赖于人工的经验),辅助医生件进行参考处理,提高判断速率和判断准确率。The above-mentioned medical image registration method, apparatus, computer equipment and storage medium are used to obtain a first medical image in a first medical image set containing specific conditions (different viewing angles at the same time) and a second medical image set containing specific conditions (same viewing angle at different times) The second medical image contains features in the second medical image, and the corresponding first medical image and the second medical image are processed by two registration transformation models respectively to perform image registration, which can improve the model's processing ability for specific data and improve Data processing efficiency; then, a plurality of first medical images under specific conditions are fused into a first registration image through the first image interpolation module, and multiple second medical images under specific conditions are fused into a first registration image through the second image interpolation module. Second, register images, process images from different perspectives and at different times to obtain all image information under each image; then fuse the first and second registration images to obtain a more comprehensive The third registration image of the image information (growth tracking results); all in all, this method uses the artificial intelligence model for processing, and abandons the phenomenon that relies too much on the doctor's subjective experience and spatial imagination before (manual participation in the work is too much, and there is also too much Relying on manual experience), assisting the doctor to perform reference processing, improving the judgment rate and judgment accuracy.
附图说明Description of drawings
为了更清楚地说明本申请实施例的技术方案,下面将对本申请实施例的描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。In order to illustrate the technical solutions of the embodiments of the present application more clearly, the following briefly introduces the drawings that are used in the description of the embodiments of the present application. Obviously, the drawings in the following description are only some embodiments of the present application. , for those of ordinary skill in the art, other drawings can also be obtained from these drawings without creative labor.
图1是本申请一实施例中医学影像配准方法的一应用环境示意图;1 is a schematic diagram of an application environment of a medical image registration method in an embodiment of the present application;
图2是本申请一实施例中医学影像配准方法的一流程图;2 is a flowchart of a medical image registration method in an embodiment of the present application;
图3是本申请一实施例中医学影像配准装置的结构示意图;3 is a schematic structural diagram of a medical image registration device in an embodiment of the present application;
图4是本申请一实施中第一标准训练图像的训练流程图;4 is a training flow chart of the first standard training image in an implementation of the present application;
图5是本申请一实施中医学影像集的处理流程图;Fig. 5 is the processing flow chart of the medical image set in one implementation of the present application;
图6是本申请一实施例中计算机设备的一示意图。FIG. 6 is a schematic diagram of a computer device in an embodiment of the present application.
具体实施方式Detailed ways
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative work fall within the protection scope of the present application.
本申请提供的医学影像配准方法,可应用在如图1的应用环境中,其中,客户端通过网络与服务器进行通信。其中,客户端可以包括但不限于各种个人计算机、笔记本电脑、智能手机、平板电脑和便携式可穿戴设备。服务器可以用独立的服务器或者是多个服务器组成的服务器集群来实现。The medical image registration method provided by the present application can be applied in the application environment as shown in FIG. 1 , wherein the client communicates with the server through the network. Wherein, the client may include but not limited to various personal computers, notebook computers, smart phones, tablet computers and portable wearable devices. The server can be implemented as an independent server or a server cluster composed of multiple servers.
在一实施例中,如图2所示,提供一种医学影像配准方法,以该方法应用在图1中的服务器为例进行说明,包括如下步骤:In one embodiment, as shown in FIG. 2 , a medical image registration method is provided, which is described by taking the method applied to the server in FIG. 1 as an example, including the following steps:
S10,获取第一医学影像集和第二医学影像集;所述第一医学影像集中包含同一时刻不同视角的多张第一医学影像;所述第二医学影像集中包含不同时刻同一视角的多张第二医学影像;S10: Acquire a first medical image set and a second medical image set; the first medical image set includes multiple first medical images from different viewing angles at the same time; the second medical image set includes multiple medical images from the same viewing angle at different times second medical image;
可理解地,第一医学影像集存在多张第一医学影像,第一医学影像可为常见的2D图像,如B超图像,第二医学影像与第一医学影像集一致,其中,第一医学影像和第二医学影像由医学成像设备(如CT设备、MRI设备或B超设备等)采集,第一医学影像反映出同一时刻不同视角的细胞组织结果的成像结果,和第二医学影像反映出不同时刻同一个视角的细胞组织结果的成像结果,且第一医学影像和第二医学影像都是针对同一个病理的生长部位,如肿瘤生长部位,在生长部位为肿瘤生长部位时,第一医学影像和第二医学图像都与该肿瘤生长部位关联。Understandably, there are multiple first medical images in the first medical image set, the first medical image may be a common 2D image, such as a B-ultrasound image, and the second medical image is consistent with the first medical image set, wherein the first medical image is The image and the second medical image are collected by medical imaging equipment (such as CT equipment, MRI equipment or B-ultrasound equipment, etc.), the first medical image reflects the imaging results of cell tissue results from different perspectives at the same time, and the second medical image reflects The imaging results of the cell tissue results from the same viewing angle at different times, and the first medical image and the second medical image are both for the same pathological growth site, such as the tumor growth site, when the growth site is the tumor growth site, the first medical image Both the image and the second medical image are associated with the tumor growth site.
S20,将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像;S20, after the first medical image set is input into the first registration transformation model, obtain a first registration image output by the first registration transformation model for performing multi-view fusion on the first medical image set;
可理解地,第一配准变换模型是由卷积神经网络+图像插值模块组成,卷积神经网络是一种深度前馈人工神经网络,应用于图像识别,第一配准变换模型主要是处理第一医学影像集,其中,卷积神经网络模型中可设置最基本的四个层,分别为卷积层、池化层、线性整流层和全连接层,卷积层用于提取输入的第一医学影像中不同特征,池化层用于对卷积层中的特征进行优化处理,得到新的且维度更小的特征,线性整流层使用活性化函数进行线性整流,全连接层用于把所有的特征转换成全局特征,计算第一医学影像所属类别的分数;具体地,在将多张第一医学影像输入至第一配准变换模型后,通过该网络中的卷积神经网络可将多张第一医学影像生成第一变形场图像,对第一变形场图像进行图像插值后,可得到第一配准图像;本实施例可将第一医学影像集中的第一医学影像进行融合成一张第一配准图像。Understandably, the first registration transformation model is composed of a convolutional neural network + an image interpolation module. The convolutional neural network is a deep feedforward artificial neural network, which is applied to image recognition. The first registration transformation model is mainly used for processing. The first medical image set, in which the most basic four layers can be set in the convolutional neural network model, namely convolutional layer, pooling layer, linear rectification layer and fully connected layer. The convolutional layer is used to extract the first Different features in a medical image, the pooling layer is used to optimize the features in the convolution layer to obtain new features with smaller dimensions, the linear rectification layer uses the activation function to perform linear rectification, and the fully connected layer is used to All the features are converted into global features, and the score of the category to which the first medical image belongs is calculated; A first deformation field image is generated from a plurality of first medical images, and after image interpolation is performed on the first deformation field image, a first registration image can be obtained; in this embodiment, the first medical images in the first medical image set can be fused into one. The first registered image.
S30,将所述第二医学影像集输入至第二配准变换模型后,获取所述第二配准变换模型输出的对所述第二医学影像集进行多时序融合的第二配准图像;S30, after inputting the second medical image set into the second registration transformation model, obtain a second registration image output by the second registration transformation model for performing multi-sequence fusion on the second medical image set;
可理解地,第二配准变换模型的组成与第一配准变换模型是一致,不同之处在于第二配准变换模型是处理第二医学影像集,可见,第二配准变换模型中的卷积神经网络训练的数据与第一配准变换模型中的卷积神经网络不同,而两者之间的网络结构可完全相同;本实施例可将第二医学影像集中的第二医学影像进行融合成一张第二配准图像。Understandably, the composition of the second registration transformation model is the same as that of the first registration transformation model, the difference is that the second registration transformation model processes the second medical image set. It can be seen that in the second registration transformation model, the The data trained by the convolutional neural network is different from the convolutional neural network in the first registration transformation model, and the network structure between the two can be exactly the same; in this embodiment, the second medical image in the second medical image set can be processed. fused into a second registered image.
S40,对所述第一配准图像和所述第二配准图像进行融合,得到第三配准图像;所述第三配准图像反映出所述第一医学影像集和所述第二医学影像集的生长追踪结果。S40, fuse the first registration image and the second registration image to obtain a third registration image; the third registration image reflects the first medical image set and the second medical image Growth tracking results for the image set.
可理解地,根据上述两个配准变换模型求得第一配准图像和第二配准图像后,通过特征级图像融合对第一配准图像和第二配准图像进行图像融合,得到第三配准图像,具体地,提取第一配准图像和第二配准图像中的特征信息(如边缘、形状、轮毂和局部特征信息),将提取到的特征信息进行综合处理,综合处理包括目标状态特征融合或目标特性融合中的任意一种,通过目标状态特征融合将与同一时刻不同视角关联的第一配准图像和不同时刻相同视角的第二配准图像进行目标特征统计,对两者提取到的目标特征进行严格配准,得到包含更多图像特征的第三配准图像,通过目标特性融合对第一配准图像和第二配准图像中的特征向量进行重新组合,增加两者之间的图像特征维度,得到第三配准图像,其中,第三配准图像能对多视角和多时序的图像进行融合,因此第三配准图像可完成显示出病理的生长追踪结果;在本实施例之后,可将生长追踪结果的综合变换情况反馈至预设数据接收方。Understandably, after obtaining the first registration image and the second registration image according to the above two registration transformation models, image fusion is performed on the first registration image and the second registration image through feature-level image fusion to obtain the first registration image and the second registration image. Three registration images, specifically, extract feature information (such as edge, shape, hub and local feature information) in the first registration image and the second registration image, and perform comprehensive processing on the extracted feature information. The comprehensive processing includes Either of target state feature fusion or target feature fusion, through target state feature fusion, target feature statistics are performed on the first registration image associated with different viewing angles at the same time and the second registration image at the same viewing angle at different times. The target features extracted by the user are strictly registered to obtain a third registration image containing more image features, and the feature vectors in the first registration image and the second registration image are recombined through target feature fusion, adding two The image feature dimension between the two is obtained to obtain a third registration image, wherein the third registration image can fuse multi-view and multi-sequence images, so the third registration image can complete the pathological growth tracking result; After this embodiment, the comprehensive transformation of the growth tracking result can be fed back to the preset data receiver.
步骤S10至步骤S40所在的实施例,获取第一医学影像集中包含特定条件(同一时刻不同视角)的第一医学影像以及第二医学影像集中包含特定条件(不同时刻同一视角)的第二医学影像中包含特征的第二医学影像,通过两个配准变换模型分别处理对应的第一医学影像和第二医学影像进行图像配准,可提高模型对特定数据的处理能力,提高数据处理效率;接着通过第一图像插值模块将特定条件下多张第一医学图像融合成一张第一配准图像,通过第二图像插值模块将特定条件下多张第二医学图像融合成一张第二配准图像,对不同视角和不同时刻下的图像进行分别处理,得到每一种图像下的所有图像信息;之后再对第一配准图像和第二配准图像进行融合,得到包含更全面图像信息(生长追踪结果)的第三配准图像;总而言之,本申请借用人工智能模型进行处理,辅助医生件处理,提高判断速率和判断准确率。本申请可应用于智慧医疗领域中,从而推动智慧城市的建设。In the embodiment in which steps S10 to S40 are located, a first medical image containing a specific condition (different views at the same time) in the first medical image set and a second medical image containing a specific condition (same view at different times) in the second medical image set are acquired For the second medical image containing features, two registration transformation models are used to process the corresponding first medical image and the second medical image respectively to perform image registration, which can improve the model's processing capability for specific data and improve data processing efficiency; then The first image interpolation module fuses multiple first medical images under specific conditions into one first registration image, and the second image interpolation module fuses multiple second medical images under specific conditions into one second registration image, The images at different viewing angles and at different times are processed separately to obtain all image information under each image; then the first registration image and the second registration image are fused to obtain more comprehensive image information (growth tracking). Results) of the third registration image; all in all, this application uses an artificial intelligence model for processing, assisting in the processing of doctor's documents, and improving the judgment rate and judgment accuracy. The present application can be applied in the field of smart medical care, thereby promoting the construction of smart cities.
进一步地,如图5所述,在通过医学成像设备获取关于肿瘤患者的医学影像集后,按照视角和时刻的关系划分医学影像集后得到第一医学影像集合第二医学影像集,再将第一医学影像集和第二医学医学集分别输入至第一配准网络和第二配准网络后,得到包含综合情况的第一配准图像和包含变化情况的第二配准图像,最后通过融合两张图像得到包含肿瘤生长追踪情况的第三配准图像。Further, as shown in FIG. 5 , after obtaining a medical image set about a tumor patient through a medical imaging device, the medical image set is divided according to the relationship between the viewing angle and the time to obtain the first medical image set and the second medical image set, and then the first medical image set is obtained. After a medical image set and a second medical medical set are respectively input to the first registration network and the second registration network, the first registration image including the comprehensive situation and the second registration image including the change situation are obtained, and finally through fusion The two images resulted in a third registration image containing tumor growth tracking.
进一步地,所述第一配准变换模型包括第一卷积神经网络和第一图像插值模块,所述将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像,包括:Further, the first registration transformation model includes a first convolutional neural network and a first image interpolation module, and after the first medical image set is input into the first registration transformation model, the first registration transformation model is obtained. The first registration image output by the registration transformation model for performing multi-view fusion on the first medical image set includes:
获取所述第一医学影像集中的多张所述第一医学影像;acquiring a plurality of the first medical images in the first medical image set;
将所有所述第一医学影像输入至第一卷积神经网络后,获取将各所述第一医学影像变形之后生成的多个第一变形场图像;After inputting all the first medical images into the first convolutional neural network, acquiring a plurality of first deformation field images generated after deforming each of the first medical images;
将所有所述第一变形场图像通过第一图像插值模块进行图像插值后,得到一张所述第一配准图像。After performing image interpolation on all the first deformation field images through the first image interpolation module, a first registration image is obtained.
可理解地,本实施例中的第一配准变换模型是包括图像插值模块,其中图像插值的目的是将多个变形场变换到同一个图像空间中,具体公式见下述;本实施例中的一张第一医学影像对应一个变形场,变形场中包括了卷积神经网络对第一医学图像进行处理后所得到的图像特征。Understandably, the first registration transformation model in this embodiment includes an image interpolation module, where the purpose of image interpolation is to transform multiple deformation fields into the same image space, and the specific formula is as follows; A first medical image corresponding to a deformation field, and the deformation field includes image features obtained by processing the first medical image by a convolutional neural network.
进一步地,如图4所示,所述将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像之前,还包括:Further, as shown in FIG. 4 , after the first medical image set is input into the first registration transformation model, the output of the first registration transformation model is obtained to perform multiple operations on the first medical image set. Before the first registration image for perspective fusion, it also includes:
获取第一标准训练图像以及第一参考图像,将一张所述第一标准训练图像以及一张所述第一参考图像作为一对第一配准对,通过第一卷积神经网络对所有所述第一配准对进行变形后,得到多个第一训练变形场;一个所述第一训练变形场对应一张所述第一标准训练图像;Obtain a first standard training image and a first reference image, use a first standard training image and a first reference image as a pair of first registration After the first registration pair is deformed, a plurality of first training deformation fields are obtained; one of the first training deformation fields corresponds to a piece of the first standard training image;
将所述第一训练变形场通过第一图像插值模块进行图像插值后,得到第一插值结果;After the first training deformation field is subjected to image interpolation by the first image interpolation module, a first interpolation result is obtained;
将所述第一插值结果和所述第一参考图像输入至与所述第一卷积神经网络关联的第一预设损失函数后,获取第一训练变形场和第一参考图像的第一损失值;After the first interpolation result and the first reference image are input into the first preset loss function associated with the first convolutional neural network, the first loss of the first training deformation field and the first reference image is obtained value;
在所述第一损失值大于第一预设损失值时,通过反向传播算法不断更新所述第一卷积神经网络的网络参数直至所述第一损失值小于或等于所述第一预设损失值;When the first loss value is greater than a first preset loss value, the network parameters of the first convolutional neural network are continuously updated through a back-propagation algorithm until the first loss value is less than or equal to the first preset loss value loss value;
在所述第一损失值小于或等于第一预设损失值时,确定所述第一配准变换模型训练完成,得到所述第一配准变换模型。When the first loss value is less than or equal to the first preset loss value, it is determined that the training of the first registration transformation model is completed, and the first registration transformation model is obtained.
可理解地,第一标准训练图像可用影像序列进行排列,影像序列为I m(m=1,2...n),每一个序号代表同一时刻不同视角的第一标准训练图像,第一参考图像(可理解为一个视角下的标准图像)用I 0表示,一个第一参考图像对应一个视角的第一标准训练图像;第一训练变形场只与第一标准训练图像存在对应关系;本实施例是模型的训练方式,通过该训练方式先求出损失函数(f m为第一训练变形场,假设第一插值结果为I(f m),I(f m)中包含第一训练变形场与第一参考图像之间的数学关系,为线性插值求出的结果,第一损失函数为L(m)=∑ m|I 0-I(f m)|,通过该第一损失函数计算第一标准训练图像与第一参数图像之间的差别)的值,在确定第一损失值大于第一预设损失值,不断地更新卷积神经网络中隐藏层(全连接层)中的网络参数(权重)直至述第一损失值小于或等于第一预设损失值;本实施例相比于之前传统训练方式,无需先提取图像特征再选取模型去进行不断迭代优化模型的繁琐过程,可利用无监督学习的卷积神经网络模型直接实现训练。 Understandably, the first standard training images can be arranged in an image sequence, and the image sequence is Im ( m =1, 2...n), each serial number represents the first standard training images from different perspectives at the same time, and the first reference The image (which can be understood as a standard image under one viewing angle) is represented by I 0 , and a first reference image corresponds to a first standard training image of one viewing angle; the first training deformation field only has a corresponding relationship with the first standard training image; this implementation The example is the training method of the model, through which the loss function (f m is the first training deformation field, assuming that the first interpolation result is I(f m ), I(f m ) contains the first training deformation field The mathematical relationship with the first reference image is the result obtained by linear interpolation, the first loss function is L(m)=∑ m |I 0 -I(f m )|, and the first loss function is used to calculate the first loss function. The value of the difference between a standard training image and the first parameter image), when it is determined that the first loss value is greater than the first preset loss value, the network parameters in the hidden layer (full connection layer) in the convolutional neural network are continuously updated (weight) until the first loss value is less than or equal to the first preset loss value; compared with the previous traditional training method, this embodiment does not need to extract image features first and then select the model to perform the cumbersome process of iteratively optimizing the model. Unsupervised learning of convolutional neural network models directly implements training.
进一步地,所述第二配准变换模型包括第二卷积神经网络和第二图像插值模块,所述将所述第二医学影像集输入至第二配准变换模型后,获取所述第二配准变换模型输出的对所述第二医学影像集进行多时序融合的第二配准图像,包括:Further, the second registration transformation model includes a second convolutional neural network and a second image interpolation module, and after the second medical image set is input into the second registration transformation model, the second registration transformation model is obtained. The second registration image output by the registration transformation model for performing multi-sequence fusion on the second medical image set includes:
获取所述第二医学影像集中的多张第二配准对;一张所述第二配准对包括一张所述第二医学影像和一张第二参考图像;acquiring a plurality of second registration pairs in the second medical image set; a second registration pair includes a second medical image and a second reference image;
将所有所述第二配准对输入至第二卷积神经网络后,获取将各所述第二配准对变形之后生成的多个第二变形场;After inputting all the second registration pairs to the second convolutional neural network, acquiring a plurality of second deformation fields generated after deforming each of the second registration pairs;
将所有所述第二变形场通过第二图像插值模块进行图像插值后,得到所述第二配准图像。After all the second deformation fields are subjected to image interpolation by the second image interpolation module, the second registration image is obtained.
可理解地,该实施例的具体内容可参照上述第一配准图像所在实施例,在此不在赘述。Understandably, for the specific content of this embodiment, reference may be made to the above-mentioned embodiment where the first registered image is located, and details are not described here.
进一步地,所述将所有所述第二配准对输入至由卷积神经网络构建成的第一配准变换模型之前,还包括:Further, before the inputting all the second registration pairs to the first registration transformation model constructed by the convolutional neural network, the method further includes:
获取第二标准训练图像以及第二参考图像,将一张所述第二标准训练图像以及一张所述第二参考图像作为一对第二配准对,通过第二卷积神经网络对所有所述第二配准对进行变形后,得到多个第二训练变形场;一个所述第二训练变形场对应一张所述第二标准训练图像;Acquire a second standard training image and a second reference image, use one of the second standard training image and one of the second reference images as a pair of second registration pairs, and use the second convolutional neural network for all the After the second registration pair is deformed, a plurality of second training deformation fields are obtained; one of the second training deformation fields corresponds to a piece of the second standard training image;
将所述第二训练变形场通过第二图像插值模块进行图像插值后,得到第二插值结果;After the second training deformation field is subjected to image interpolation by the second image interpolation module, a second interpolation result is obtained;
将所述第二插值结果和所述第二参考图像输入至与所述第二卷积神经网络关联的第二预设损失函数后,获取第二训练变形场和第二参考图像的第二损失值;After the second interpolation result and the second reference image are input to a second preset loss function associated with the second convolutional neural network, a second loss of the second training deformation field and the second reference image is obtained value;
在所述第二损失值大于第二预设损失值时,通过反向传播算法不断更新所述第二卷积神经网络的网络参数直至所述第二损失值小于或等于所述第二预设损失值;When the second loss value is greater than a second preset loss value, the network parameters of the second convolutional neural network are continuously updated through a back-propagation algorithm until the second loss value is less than or equal to the second preset loss value loss value;
在所述第二损失值小于或等于第二预设损失值时,确定所述第二配准变换模型训练完成,得到所述第二配准变换模型。When the second loss value is less than or equal to the second preset loss value, it is determined that the training of the second registration transformation model is completed, and the second registration transformation model is obtained.
可理解地,第二标准训练图像可用影像序列进行排列,影像序列为I m'(m'=1,2...n),每一个序号代表不同时刻同一视角的第二标准训练图像,第二参考图像(可理解为一个时刻下的标准图像)用I 0'表示,一个第二参考图像对应一个时刻的第二标准训练图像;本实施例的的具体内容可参照上述第一配准变换模型一致,在此不在赘述。 Understandably, the second standard training images can be arranged in an image sequence, where the image sequence is Im ' ( m '=1, 2...n), each serial number represents the second standard training images of the same viewing angle at different times, the first Two reference images (which can be understood as a standard image at a moment) are represented by I 0 ′, and a second reference image corresponds to a second standard training image at a moment; for the specific content of this embodiment, refer to the above-mentioned first registration transformation The model is consistent and will not be repeated here.
进一步地,所述将所述第一训练变形场通过第一图像插值模块进行图像插值后,得到第一插值结果,包括:Further, after performing image interpolation on the first training deformation field through the first image interpolation module, a first interpolation result is obtained, including:
通过第一变换公式将与所述第一训练变形场中的第一标准训练图像变换到所述第一训练变形场中的第一参考图像的图像空间中,得到包含在同一个图像空间的第一目标配准图像的第一插值结果;其中,所述第一变换公式为:Transform the first standard training image in the first training deformation field into the image space of the first reference image in the first training deformation field by using the first transformation formula, so as to obtain the first standard training image contained in the same image space. A first interpolation result of a target registration image; wherein, the first transformation formula is:
X'=ax'+by',Y'=cx'+dy'X'=ax'+by', Y'=cx'+dy'
其中:in:
x'为所述第一标准训练图像中像素的横坐标;x' is the abscissa of the pixel in the first standard training image;
y'为所述第一标准训练图像中像素的纵坐标;y' is the ordinate of the pixel in the first standard training image;
X'为所述第一参考图像或所述第一目标配准图像中像素的横坐标;X' is the abscissa of the pixel in the first reference image or the first target registration image;
Y'为所述第一参考图像或所述第一目标配准图像中像素的纵坐标;Y' is the ordinate of the pixel in the first reference image or the first target registration image;
a、b、c和d均为变换参数值,a、b、c和d构成的变换矩阵
Figure PCTCN2021096702-appb-000001
可理解地,变换矩阵可预先通过多组测到的变换之前和变换之后的横坐标和纵坐标数据求出。
a, b, c and d are all transformation parameter values, and the transformation matrix formed by a, b, c and d
Figure PCTCN2021096702-appb-000001
Understandably, the transformation matrix can be obtained in advance through multiple sets of measured abscissa and ordinate data before and after transformation.
进一步地,所述通过第一卷积神经网络对所有所述第一配准对进行变形后,得到多个第一训练变形场,包括:Further, after all the first registration pairs are deformed by the first convolutional neural network, a plurality of first training deformation fields are obtained, including:
将所有所述第一配准对输入至所述卷积神经网络后,通过与所述卷积神经网络对应的第一预设表达式对所有所述第一配准对进行变形,得到与所述第一配准对对应的多个第一训练变形场;所述第一预设表达式包括:After all the first registration pairs are input into the convolutional neural network, all the first registration pairs are deformed by the first preset expression corresponding to the convolutional neural network to obtain the same Multiple first training deformation fields corresponding to the first registration pair; the first preset expression includes:
Figure PCTCN2021096702-appb-000002
Figure PCTCN2021096702-appb-000002
Figure PCTCN2021096702-appb-000003
Figure PCTCN2021096702-appb-000003
其中:in:
x为所述第一标准训练图像的强度;x is the intensity of the first standard training image;
y为所述第一参考图像的强度;y is the intensity of the first reference image;
X为x的取值范围;X如x1、x2...xn等;X is the value range of x; X is x1, x2...xn, etc.;
Y为y的取值范围;y如y1、y2...yn等;Y is the value range of y; y such as y1, y2...yn, etc.;
D m为所述第一标准训练图像和所述第一参考图像像素之间的相似度; D m is the similarity between the first standard training image and the first reference image pixel;
x i为所述第一标准训练图像第i个像素的强度; x i is the intensity of the ith pixel of the first standard training image;
y i为所述第一参考图像第i个像素的强度; y i is the intensity of the ith pixel of the first reference image;
x m为所述第一标准训练图像的平均强度; x m is the average intensity of the first standard training image;
y m为所述第一参考图像的平均强度; y m is the average intensity of the first reference image;
p(x)为所述第一标准训练图像概率密度分布函数;p(x) is the probability density distribution function of the first standard training image;
p(y)为所述第一参考图像概率密度分布函数;p(y) is the probability density distribution function of the first reference image;
p(x,y)为所述第一标准训练图像和所述第一参考图像的联合概率密度函数。p(x,y) is the joint probability density function of the first standard training image and the first reference image.
进一步地,所述将所述第二训练变形场通过第二图像插值模块进行图像插值后,得到第二插值结果,包括:Further, after performing image interpolation on the second training deformation field through the second image interpolation module, a second interpolation result is obtained, including:
通过第二变换公式将与所述第二训练变形场中的第二标准训练图像变换到所述第二训练变形场中的第二参考图像的图像空间中,得到包含在同一个图像空间的第二目标配准图像的第二插值结果;其中,所述第二变换公式为:Transform the second standard training image in the second training deformation field into the image space of the second reference image in the second training deformation field by using the second transformation formula, so as to obtain the first standard training image contained in the same image space. The second interpolation result of the two-target registration image; wherein, the second transformation formula is:
X”=a'x”+b'y”,Y”=c'x”+d'y”X"=a'x"+b'y", Y"=c'x"+d'y"
其中:in:
x”为所述第二标准训练图像中像素的横坐标;x" is the abscissa of the pixel in the second standard training image;
y”为所述第二标准训练图像中像素的纵坐标;y" is the ordinate of the pixel in the second standard training image;
X”为所述第二参考图像或所述第二目标配准图像中像素的横坐标;X" is the abscissa of the pixel in the second reference image or the second target registration image;
Y"为所述第二参考图像或第二目标配准图像中像素的纵坐标;Y" is the ordinate of the pixel in the second reference image or the second target registration image;
a’、b’、c’和d’均为变换参数值,a’、b’、c’和d’构成的变换矩阵
Figure PCTCN2021096702-appb-000004
a', b', c' and d' are all transformation parameter values, and the transformation matrix composed of a', b', c' and d'
Figure PCTCN2021096702-appb-000004
进一步地,所述通过第二卷积神经网络对所有所述第二配准对进行变形后,得到多个第二训练变形场,包括:Further, after all the second registration pairs are deformed by the second convolutional neural network, a plurality of second training deformation fields are obtained, including:
将所有所述第二配准对输入至所述卷积神经网络后,通过与所述卷积神经网络对应的第二预设表达式对所有所述第二配准对进行变形,得到与所述第二配准对对应的多个第二训练变形场;所述第二预设表达式包括:After all the second registration pairs are input into the convolutional neural network, all the second registration pairs are deformed by a second preset expression corresponding to the convolutional neural network to obtain the same Multiple second training deformation fields corresponding to the second registration pair; the second preset expression includes:
Figure PCTCN2021096702-appb-000005
Figure PCTCN2021096702-appb-000005
Figure PCTCN2021096702-appb-000006
Figure PCTCN2021096702-appb-000006
其中:in:
x”'为所述第二标准训练图像的强度;x"' is the intensity of the second standard training image;
y”'为所述第二参考图像的强度;y"' is the intensity of the second reference image;
X”'为x”'的取值范围;X"' is the value range of x"';
Y”'为y”'的取值范围;Y"' is the value range of y"';
D m'为所述第二标准训练图像和所述第二参考图像像素之间的相似度; D m ' is the similarity between the second standard training image and the second reference image pixel;
x i'为所述第二标准训练图像第i个像素的强度; x i ' is the intensity of the i-th pixel of the second standard training image;
y i'为所述第二参考图像第i个像素的强度; y i ' is the intensity of the ith pixel of the second reference image;
x m'为所述第二标准训练图像的平均强度; x m ' is the average intensity of the second standard training image;
y m'为所述第二参考图像的平均强度; y m ' is the average intensity of the second reference image;
p(x”')为所述第二标准训练图像概率密度分布函数;p(x"') is the probability density distribution function of the second standard training image;
p(y”')为所述第二参考图像概率密度分布函数;p(y"') is the probability density distribution function of the second reference image;
p(x”',y”')为所述第二标准训练图像和所述第二参考图像的联合概率密度函数。p(x"', y"') is the joint probability density function of the second standard training image and the second reference image.
综上所述,上述提供了一种医学影像配准方法,获取第一医学影像集中包含特定条件(同一时刻不同视角)的第一医学影像以及第二医学影像集中包含特定条件(不同时刻同一视角)的第二医学影像中包含特征的第二医学影像,通过两个配准变换模型分别处理对应的第一医学影像和第二医学影像进行图像配准,可提高模型对特定数据的处理能力,提高数据处理效率;接着通过第一图像插值模块将特定条件下多张第一医学图像融合成一张第一配准图像,通过第二图像插值模块将特定条件下多张第二医学图像融合成一张第二配准图像,对不同视角和不同时刻下的图像进行分别处理,得到每一种图像下的所有图像信息;之后再对第一配准图像和第二配准图像进行融合,得到包含更全面图像信息(生长追踪结果)的第三配准图像;总而言之,本方法借用人工智能模型进行处理,抛弃之前过多依赖于医生的主观经验和空间想象的现象,辅助医生件进行参考处理,提高判断速率和判断准确率。本方法可应用于智慧医疗领域中,从而推动智慧城市的建设。To sum up, a medical image registration method is provided above, and a first medical image containing a specific condition (different viewing angles at the same time) in the first medical image set and a second medical image set containing specific conditions (the same viewing angle at different times) are obtained. ) in the second medical image of the second medical image containing features, and processing the corresponding first medical image and the second medical image through two registration transformation models respectively to perform image registration, which can improve the model's processing capability for specific data, Improve data processing efficiency; then fuse multiple first medical images under specific conditions into a first registration image through a first image interpolation module, and fuse multiple second medical images under specific conditions into a single image through a second image interpolation module For the second registration image, the images from different perspectives and at different times are processed separately to obtain all the image information under each image; then the first registration image and the second registration image are fused to obtain a more The third registration image of comprehensive image information (growth tracking results); in a word, this method uses artificial intelligence model for processing, abandons the phenomenon that relies too much on the doctor's subjective experience and spatial imagination, and assists the doctor's software in reference processing, improving Judgment rate and judgment accuracy. The method can be applied in the field of smart medical care, thereby promoting the construction of smart cities.
应理解,上述实施例中各步骤的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本申请实施例的实施过程构成任何限定。It should be understood that the size of the sequence numbers of the steps in the above embodiments does not mean the sequence of execution, and the execution sequence of each process should be determined by its function and internal logic, and should not constitute any limitation to the implementation process of the embodiments of the present application.
在一实施例中,提供一种医学影像配准装置,该医学影像配准装置与上述实施例中医学影像配准方法一一对应。如图3所示,该医学影像配准装置包括第一获取模块11、第二获取模块12、第三获取模块13和融合模块14。各功能模块详细说明如下:In one embodiment, a medical image registration apparatus is provided, and the medical image registration apparatus is in one-to-one correspondence with the medical image registration method in the above embodiment. As shown in FIG. 3 , the medical image registration device includes a first acquisition module 11 , a second acquisition module 12 , a third acquisition module 13 and a fusion module 14 . The detailed description of each functional module is as follows:
第一获取模块11,用于获取第一医学影像集和第二医学影像集;所述第一医学影像集中包含同一时刻不同视角的多张第一医学影像;所述第二医学影像集中包含不同时刻同一视角的多张第二医学影像;The first acquisition module 11 is configured to acquire a first medical image set and a second medical image set; the first medical image set includes a plurality of first medical images from different viewing angles at the same time; the second medical image set includes different Multiple second medical images at the same viewing angle at the same time;
第二获取模块12,用于将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像;The second acquisition module 12 is configured to acquire, after the first medical image set is input into the first registration transformation model, the output of the first registration transformation model for performing multi-view fusion on the first medical image set the first registered image;
第三获取模块13,用于将所述第二医学影像集输入至第二配准变换模型后,获取所述第二配准变换模型输出的对所述第二医学影像集进行多时序融合的第二配准图像;The third acquisition module 13 is configured to acquire, after the second medical image set is input into the second registration transformation model, the output of the second registration transformation model for performing multi-sequence fusion on the second medical image set the second registration image;
融合模块14,用于对所述第一配准图像和所述第二配准图像进行融合,得到第三配 准图像;所述第三配准图像反映出所述第一医学影像集和所述第二医学影像集的生长追踪结果。The fusion module 14 is configured to fuse the first registration image and the second registration image to obtain a third registration image; the third registration image reflects the first medical image set and the Growth tracking results for the second set of medical images described.
进一步地,所述第二获取模块包括:Further, the second acquisition module includes:
第一获取子模块,用于获取所述第一医学影像集中的多张所述第一医学影像;a first acquisition submodule, configured to acquire a plurality of the first medical images in the first medical image set;
第二获取子模块,用于将所有所述第一医学影像输入至第一卷积神经网络后,获取将各所述第一医学影像变形之后生成的多个第一变形场图像;a second acquisition sub-module, configured to acquire a plurality of first deformation field images generated by deforming each of the first medical images after inputting all the first medical images into the first convolutional neural network;
第一图像插值模块子模块,用于将所有所述第一变形场图像通过第一图像插值模块进行图像插值后,得到一张所述第一配准图像。The first image interpolation module sub-module is configured to obtain a first registration image after performing image interpolation on all the first deformation field images through the first image interpolation module.
进一步地,所述医学影像配准装置还包括:Further, the medical image registration device further includes:
第一变形模块,用于获取第一标准训练图像以及第一参考图像,将一张所述第一标准训练图像以及一张所述第一参考图像作为一对第一配准对,通过第一卷积神经网络对所有所述第一配准对进行变形后,得到多个第一训练变形场;一个所述第一训练变形场对应一张所述第一标准训练图像;The first deformation module is used to obtain a first standard training image and a first reference image, and use one of the first standard training image and one of the first reference images as a pair of first registration pairs, through the first After the convolutional neural network deforms all the first registration pairs, a plurality of first training deformation fields are obtained; one of the first training deformation fields corresponds to one of the first standard training images;
第一图像插值模块模块,用于将所述第一训练变形场通过第一图像插值模块进行图像插值后,得到第一插值结果;a first image interpolation module, configured to obtain a first interpolation result after performing image interpolation on the first training deformation field through the first image interpolation module;
第四获取模块,用于将所述第一插值结果和所述第一参考图像输入至与所述第一卷积神经网络关联的第一预设损失函数后,获取第一训练变形场和第一参考图像的第一损失值;The fourth obtaining module is configured to obtain the first training deformation field and the first training deformation field after inputting the first interpolation result and the first reference image into the first preset loss function associated with the first convolutional neural network. a first loss value of the reference image;
第一更新模块,用于在所述第一损失值大于第一预设损失值时,通过反向传播算法不断更新所述第一卷积神经网络的网络参数直至所述第一损失值小于或等于所述第一预设损失值;A first update module, configured to continuously update the network parameters of the first convolutional neural network through a back-propagation algorithm when the first loss value is greater than a first preset loss value until the first loss value is less than or is equal to the first preset loss value;
第一确定子模块,用于在所述第一损失值小于或等于第一预设损失值时,确定所述第一配准变换模型训练完成,得到所述第一配准变换模型。A first determination submodule, configured to determine that the training of the first registration transformation model is completed when the first loss value is less than or equal to a first preset loss value, and obtain the first registration transformation model.
进一步地,所述第三获取模块包括:Further, the third acquisition module includes:
第三获取子模块,用于获取所述第二医学影像集中的多张第二配准对;一张所述第二配准对包括一张所述第二医学影像和一张第二参考图像;The third acquisition sub-module is configured to acquire a plurality of second registration pairs in the second medical image set; a second registration pair includes a second medical image and a second reference image ;
第四获取子模块,用于将所有所述第二配准对输入至第二卷积神经网络后,获取将各所述第二配准对变形之后生成的多个第二变形场;a fourth acquisition sub-module, configured to acquire a plurality of second deformation fields generated after deforming each of the second registration pairs after inputting all the second registration pairs into the second convolutional neural network;
第二图像插值模块子模块,用于将所有所述第二变形场通过第二图像插值模块进行图像插值后,得到所述第二配准图像。The second image interpolation module sub-module is configured to obtain the second registration image after performing image interpolation on all the second deformation fields through the second image interpolation module.
进一步地,所述医学影像配准装置还包括:Further, the medical image registration device further includes:
第二变形模块,用于获取第二标准训练图像以及第二参考图像,将一张所述第二标准训练图像以及一张所述第二参考图像作为一对第二配准对,通过第二卷积神经网络对所有所述第二配准对进行变形后,得到多个第二训练变形场;一个所述第二训练变形场对应一张所述第二标准训练图像;The second deformation module is configured to obtain a second standard training image and a second reference image, and use a second standard training image and a second reference image as a pair of second registration pairs. After the convolutional neural network deforms all the second registration pairs, a plurality of second training deformation fields are obtained; one of the second training deformation fields corresponds to a second standard training image;
第二图像插值模块模块,用于将所述第二训练变形场通过第二图像插值模块进行图像插值后,得到第二插值结果;a second image interpolation module, configured to obtain a second interpolation result after performing image interpolation on the second training deformation field through the second image interpolation module;
第五获取模块,用于将所述第二插值结果和所述第二参考图像输入至与所述第二卷积神经网络关联的第二预设损失函数后,获取第二训练变形场和第二参考图像的第二损失值;A fifth acquisition module, configured to acquire the second training deformation field and the first the second loss value of the second reference image;
第二更新模块,用于在所述第二损失值大于第二预设损失值时,通过反向传播算法不断更新所述第二卷积神经网络的网络参数直至所述第二损失值小于或等于所述第二预设损失值;a second update module, configured to continuously update the network parameters of the second convolutional neural network through a back-propagation algorithm when the second loss value is greater than a second preset loss value until the second loss value is less than or is equal to the second preset loss value;
第二确定模块,用于在所述第二损失值小于或等于第二预设损失值时,确定所述第二配准变换模型训练完成,得到所述第二配准变换模型。A second determination module, configured to determine that the training of the second registration transformation model is completed when the second loss value is less than or equal to a second preset loss value, and obtain the second registration transformation model.
进一步地,所述第一图像插值模块模块包括:Further, the first image interpolation module module includes:
变换子模块,用于通过第一变换公式将与所述第一训练变形场中的第一标准训练图像变换到所述第一训练变形场中的第一参考图像的图像空间中,得到包含在同一个图像空间的第一目标配准图像的第一插值结果;其中,所述第一变换公式为:The transformation sub-module is used to transform the first standard training image in the first training deformation field into the image space of the first reference image in the first training deformation field through the first transformation formula, and obtain the image space contained in the first training deformation field. The first interpolation result of the first target registration image in the same image space; wherein, the first transformation formula is:
X'=ax'+by',Y'=cx'+dy'X'=ax'+by', Y'=cx'+dy'
其中:in:
x'为所述第一标准训练图像中像素的横坐标;x' is the abscissa of the pixel in the first standard training image;
y'为所述第一标准训练图像中像素的纵坐标;y' is the ordinate of the pixel in the first standard training image;
X'为所述第一参考图像中像素的横坐标;X' is the abscissa of the pixel in the first reference image;
Y'为所述第一参考图像中像素的纵坐标;Y' is the ordinate of the pixel in the first reference image;
a、b、c和d均为变换参数值,a、b、c和d构成的变换矩阵
Figure PCTCN2021096702-appb-000007
a, b, c and d are all transformation parameter values, and the transformation matrix formed by a, b, c and d
Figure PCTCN2021096702-appb-000007
进一步地,所述第一变形模块包括:Further, the first deformation module includes:
第五获取子模块,包括将所有所述第一配准对输入至所述卷积神经网络后,通过与所述卷积神经网络对应的第一预设表达式对所有所述第一配准对进行变形,得到与所述第一配准对对应的多个第一训练变形场;所述第一预设表达式包括:The fifth acquisition sub-module includes, after inputting all the first registration pairs into the convolutional neural network, performing registration on all the first registrations through a first preset expression corresponding to the convolutional neural network Deformation is performed on the pair to obtain a plurality of first training deformation fields corresponding to the first registration pair; the first preset expression includes:
Figure PCTCN2021096702-appb-000008
Figure PCTCN2021096702-appb-000008
Figure PCTCN2021096702-appb-000009
Figure PCTCN2021096702-appb-000009
其中:in:
x为所述第一标准训练图像的强度;x is the intensity of the first standard training image;
y为所述第一标准训练图像的强度;y is the intensity of the first standard training image;
X为x的取值范围;X is the value range of x;
Y为y的取值范围;Y is the value range of y;
D m为所述第一标准训练图像和所述第一参考图像像素之间的相似度; D m is the similarity between the first standard training image and the first reference image pixel;
x i为所述第一标准训练图像第i个像素的强度; x i is the intensity of the ith pixel of the first standard training image;
y i为所述第一标准训练图像第i个像素的强度; y i is the intensity of the ith pixel of the first standard training image;
x m为所述第一参考图像的平均强度; x m is the average intensity of the first reference image;
y m为所述第一参考图像的平均强度; y m is the average intensity of the first reference image;
p(x)为所述第一标准训练图像概率密度分布函数;p(x) is the probability density distribution function of the first standard training image;
p(y)为所述第一参考图像概率密度分布函数;p(y) is the probability density distribution function of the first reference image;
p(x,y)为所述第一标准训练图像和所述第一参考图像的联合概率密度函数。p(x,y) is the joint probability density function of the first standard training image and the first reference image.
关于医学影像配准装置的具体限定可以参见上文中对于医学影像配准方法的限定,在此不再赘述。上述医学影像配准装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。For the specific limitations of the medical image registration device, reference may be made to the above limitations on the medical image registration method, which will not be repeated here. Each module in the above-mentioned medical image registration apparatus may be implemented in whole or in part by software, hardware and combinations thereof. The above modules can be embedded in or independent of the processor in the computer device in the form of hardware, or stored in the memory in the computer device in the form of software, so that the processor can call and execute the operations corresponding to the above modules.
在一个实施例中,提供了一种计算机设备,该计算机设备可以是服务器,其内部结构图可以如图4所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口和数据库。其中,该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括计算机可读指令、内存储器。该计算机可读指令存储有操作系统、计算机可读指令和数据库。该内存储器为计算机可读指令中的操作系统和计算机可读指令的运行提供环境。该计算机设备的数据库用于存储多条历史测试数据,每条历史测试数据对应有测试问题记录。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机可读指令被处理器执行时以实现一种医学影像配准方法。In one embodiment, a computer device is provided, and the computer device may be a server, and its internal structure diagram may be as shown in FIG. 4 . The computer device includes a processor, memory, a network interface, and a database connected by a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes computer readable instructions, internal memory. The computer readable instructions store an operating system, computer readable instructions and a database. The internal memory provides an environment for the operating system in the computer readable instructions and the execution of the computer readable instructions. The database of the computer equipment is used to store multiple pieces of historical test data, and each piece of historical test data corresponds to a test problem record. The network interface of the computer device is used to communicate with an external terminal through a network connection. The computer readable instructions, when executed by a processor, implement a medical image registration method.
在一个实施例中,提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机可读指令,所述处理器执行所述计算机可读指令时实现如下步骤:In one embodiment, a computer device is provided, comprising a memory, a processor, and computer-readable instructions stored on the memory and executable on the processor, the processor performing the following when executing the computer-readable instructions step:
获取第一医学影像集和第二医学影像集;所述第一医学影像集中包含同一时刻不同视角的多张第一医学影像;所述第二医学影像集中包含不同时刻同一视角的多张第二医学影像;Acquiring a first medical image set and a second medical image set; the first medical image set includes multiple first medical images from different viewing angles at the same time; the second medical image set includes multiple second medical images from the same viewing angle at different times Medical Imaging;
将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像;After the first medical image set is input into the first registration transformation model, a first registration image output by the first registration transformation model for performing multi-view fusion on the first medical image set is obtained;
将所述第二医学影像集输入至第二配准变换模型后,获取所述第二配准变换模型输出的对所述第二医学影像集进行多时序融合的第二配准图像;After the second medical image set is input into the second registration transformation model, a second registration image output by the second registration transformation model for performing multi-sequence fusion on the second medical image set is obtained;
对所述第一配准图像和所述第二配准图像进行融合,得到第三配准图像;所述第三配准图像反映出所述第一医学影像集和所述第二医学影像集的生长追踪结果。The first registration image and the second registration image are fused to obtain a third registration image; the third registration image reflects the first medical image set and the second medical image set growth tracking results.
在一个实施例中,提供了一个或多个存储有计算机可读指令的可读存储介质,本实施例所提供的可读存储介质包括非易失性可读存储介质和易失性可读存储介质;该可读存储介质上存储有计算机可读指令,该计算机可读指令被一个或多个处理器执行时,使得所述一个或多个处理器执行如下步骤:In one embodiment, one or more readable storage media storing computer-readable instructions are provided, and the readable storage media provided in this embodiment include non-volatile readable storage media and volatile readable storage media medium; computer-readable instructions are stored on the readable storage medium, and when the computer-readable instructions are executed by one or more processors, the one or more processors perform the following steps:
获取第一医学影像集和第二医学影像集;所述第一医学影像集中包含同一时刻不同视角的多张第一医学影像;所述第二医学影像集中包含不同时刻同一视角的多张第二医学影像;Acquiring a first medical image set and a second medical image set; the first medical image set includes multiple first medical images from different viewing angles at the same time; the second medical image set includes multiple second medical images from the same viewing angle at different times Medical Imaging;
将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像;After the first medical image set is input into the first registration transformation model, a first registration image output by the first registration transformation model for performing multi-view fusion on the first medical image set is obtained;
将所述第二医学影像集输入至第二配准变换模型后,获取所述第二配准变换模型输出的对所述第二医学影像集进行多时序融合的第二配准图像;After the second medical image set is input into the second registration transformation model, a second registration image output by the second registration transformation model for performing multi-sequence fusion on the second medical image set is obtained;
对所述第一配准图像和所述第二配准图像进行融合,得到第三配准图像;所述第三配准图像反映出所述第一医学影像集和所述第二医学影像集的生长追踪结果。The first registration image and the second registration image are fused to obtain a third registration image; the third registration image reflects the first medical image set and the second medical image set growth tracking results.
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机可读指令来指令相关的硬件来完成,所述的计算机可读指令可存储于一非易失性计算机可读取存储介质或易失性可读存储介质中,该计算机可读指令在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可 包括只读存储器(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)等。Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be implemented by instructing relevant hardware through computer-readable instructions, and the computer-readable instructions can be stored in a non-volatile computer. In a readable storage medium or a volatile readable storage medium, the computer-readable instructions, when executed, may include the processes of the foregoing method embodiments. Wherein, any reference to memory, storage, database or other medium used in the various embodiments provided in this application may include non-volatile and/or volatile memory. Nonvolatile memory may include read only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory may include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous chain Road (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.
所属领域的技术人员可以清楚地了解到,为了描述的方便和简洁,仅以上述各功能单元、模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能单元、模块完成,即将所述装置的内部结构划分成不同的功能单元或模块,以完成以上描述的全部或者部分功能。Those skilled in the art can clearly understand that, for the convenience and simplicity of description, only the division of the above-mentioned functional units and modules is used as an example. Module completion, that is, dividing the internal structure of the device into different functional units or modules to complete all or part of the functions described above.
以上所述实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围,均应包含在本申请的保护范围之内。The above-mentioned embodiments are only used to illustrate the technical solutions of the present application, but not to limit them; although the present application has been described in detail with reference to the above-mentioned embodiments, those of ordinary skill in the art should understand that: it can still be used for the above-mentioned implementations. The technical solutions described in the examples are modified, or some technical features thereof are equivalently replaced; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions in the embodiments of the application, and should be included in the within the scope of protection of this application.

Claims (20)

  1. 一种医学影像配准方法,其中,包括:A medical image registration method, comprising:
    获取第一医学影像集和第二医学影像集;所述第一医学影像集中包含同一时刻不同视角的多张第一医学影像;所述第二医学影像集中包含不同时刻同一视角的多张第二医学影像;Acquiring a first medical image set and a second medical image set; the first medical image set includes multiple first medical images from different viewing angles at the same time; the second medical image set includes multiple second medical images from the same viewing angle at different times Medical Imaging;
    将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像;After the first medical image set is input into the first registration transformation model, a first registration image output by the first registration transformation model for performing multi-view fusion on the first medical image set is obtained;
    将所述第二医学影像集输入至第二配准变换模型后,获取所述第二配准变换模型输出的对所述第二医学影像集进行多时序融合的第二配准图像;After the second medical image set is input into the second registration transformation model, a second registration image output by the second registration transformation model for performing multi-sequence fusion on the second medical image set is obtained;
    对所述第一配准图像和所述第二配准图像进行融合,得到第三配准图像;所述第三配准图像反映出所述第一医学影像集和所述第二医学影像集的生长追踪结果。The first registration image and the second registration image are fused to obtain a third registration image; the third registration image reflects the first medical image set and the second medical image set growth tracking results.
  2. 如权利要求1所述的医学影像配准方法,其中,所述第一配准变换模型包括第一卷积神经网络和第一图像插值模块,所述将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像,包括:The medical image registration method according to claim 1, wherein the first registration transformation model comprises a first convolutional neural network and a first image interpolation module, and the first medical image set is input to the first After a registration transformation model, acquiring a first registration image output by the first registration transformation model for performing multi-view fusion on the first medical image set, including:
    获取所述第一医学影像集中的多张所述第一医学影像;acquiring a plurality of the first medical images in the first medical image set;
    将所有所述第一医学影像输入至第一卷积神经网络后,获取将各所述第一医学影像变形之后生成的多个第一变形场图像;After inputting all the first medical images into the first convolutional neural network, acquiring a plurality of first deformation field images generated after deforming each of the first medical images;
    将所有所述第一变形场图像通过第一图像插值模块进行图像插值后,得到一张所述第一配准图像。After performing image interpolation on all the first deformation field images through the first image interpolation module, a first registration image is obtained.
  3. 如权利要求1所述的医学影像配准方法,其中,所述将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像之前,还包括:The medical image registration method according to claim 1, wherein after the first medical image set is input into the first registration transformation model, a pair of the first registration transformation model outputted by the first registration transformation model is obtained. Before performing the first registration image for multi-view fusion, a medical image set further includes:
    获取第一标准训练图像以及第一参考图像,将一张所述第一标准训练图像以及一张所述第一参考图像作为一对第一配准对,通过第一卷积神经网络对所有所述第一配准对进行变形后,得到多个第一训练变形场;一个所述第一训练变形场对应一张所述第一标准训练图像;Obtain a first standard training image and a first reference image, use a first standard training image and a first reference image as a pair of first registration After the first registration pair is deformed, a plurality of first training deformation fields are obtained; one of the first training deformation fields corresponds to a piece of the first standard training image;
    将所述第一训练变形场通过第一图像插值模块进行图像插值后,得到第一插值结果;After the first training deformation field is subjected to image interpolation by the first image interpolation module, a first interpolation result is obtained;
    将所述第一插值结果和所述第一参考图像输入至与所述第一卷积神经网络关联的第一预设损失函数后,获取第一训练变形场和第一参考图像的第一损失值;After the first interpolation result and the first reference image are input into the first preset loss function associated with the first convolutional neural network, the first loss of the first training deformation field and the first reference image is obtained value;
    在所述第一损失值大于第一预设损失值时,通过反向传播算法不断更新所述第一卷积神经网络的网络参数直至所述第一损失值小于或等于所述第一预设损失值;When the first loss value is greater than a first preset loss value, the network parameters of the first convolutional neural network are continuously updated through a back-propagation algorithm until the first loss value is less than or equal to the first preset loss value loss value;
    在所述第一损失值小于或等于第一预设损失值时,确定所述第一配准变换模型训练完成,得到所述第一配准变换模型。When the first loss value is less than or equal to the first preset loss value, it is determined that the training of the first registration transformation model is completed, and the first registration transformation model is obtained.
  4. 如权利要求1所述的医学影像配准方法,其中,所述第二配准变换模型包括第二卷积神经网络和第二图像插值模块,所述将所述第二医学影像集输入至第二配准变换模型后,获取所述第二配准变换模型输出的对所述第二医学影像集进行多时序融合的第二配准图像,包括:The medical image registration method according to claim 1, wherein the second registration transformation model comprises a second convolutional neural network and a second image interpolation module, and the second medical image set is input to the first After the second registration transformation model, acquiring the second registration image output by the second registration transformation model for performing multi-sequence fusion on the second medical image set, including:
    获取所述第二医学影像集中的多张第二配准对;一张所述第二配准对包括一张所述第二医学影像和一张第二参考图像;acquiring a plurality of second registration pairs in the second medical image set; a second registration pair includes a second medical image and a second reference image;
    将所有所述第二配准对输入至第二卷积神经网络后,获取将各所述第二配准对变形之后生成的多个第二变形场;After inputting all the second registration pairs to the second convolutional neural network, acquiring a plurality of second deformation fields generated after deforming each of the second registration pairs;
    将所有所述第二变形场通过第二图像插值模块进行图像插值后,得到所述第二配准图像。After all the second deformation fields are subjected to image interpolation by the second image interpolation module, the second registration image is obtained.
  5. 如权利要求4所述的医学影像配准方法,其中,所述将所有所述第二配准对输入至由卷积神经网络构建成的第一配准变换模型之前,还包括:The medical image registration method according to claim 4, wherein before inputting all the second registration pairs to the first registration transformation model constructed by a convolutional neural network, the method further comprises:
    获取第二标准训练图像以及第二参考图像,将一张所述第二标准训练图像以及一张所述第二参考图像作为一对第二配准对,通过第二卷积神经网络对所有所述第二配准对进行变形后,得到多个第二训练变形场;一个所述第二训练变形场对应一张所述第二标准训练图像;Acquire a second standard training image and a second reference image, use one of the second standard training image and one of the second reference images as a pair of second registration pairs, and use the second convolutional neural network for all the After the second registration pair is deformed, a plurality of second training deformation fields are obtained; one of the second training deformation fields corresponds to a piece of the second standard training image;
    将所述第二训练变形场通过第二图像插值模块进行图像插值后,得到第二插值结果;After the second training deformation field is subjected to image interpolation by the second image interpolation module, a second interpolation result is obtained;
    将所述第二插值结果和所述第二参考图像输入至与所述第二卷积神经网络关联的第二预设损失函数后,获取第二训练变形场和第二参考图像的第二损失值;After the second interpolation result and the second reference image are input to a second preset loss function associated with the second convolutional neural network, a second loss of the second training deformation field and the second reference image is obtained value;
    在所述第二损失值大于第二预设损失值时,通过反向传播算法不断更新所述第二卷积神经网络的网络参数直至所述第二损失值小于或等于所述第二预设损失值;When the second loss value is greater than a second preset loss value, the network parameters of the second convolutional neural network are continuously updated through a back-propagation algorithm until the second loss value is less than or equal to the second preset loss value loss value;
    在所述第二损失值小于或等于第二预设损失值时,确定所述第二配准变换模型训练完成,得到所述第二配准变换模型。When the second loss value is less than or equal to the second preset loss value, it is determined that the training of the second registration transformation model is completed, and the second registration transformation model is obtained.
  6. 如权利要求3所述的医学影像配准方法,其中,所述将所述第一训练变形场通过第一图像插值模块进行图像插值后,得到第一插值结果,包括:The medical image registration method according to claim 3, wherein the obtaining a first interpolation result after performing image interpolation on the first training deformation field by a first image interpolation module, comprising:
    通过第一变换公式将与所述第一训练变形场中的第一标准训练图像变换到所述第一训练变形场中的第一参考图像的图像空间中,得到包含在同一个图像空间的第一目标配准图像的第一插值结果;其中,所述第一变换公式为:Transform the first standard training image in the first training deformation field into the image space of the first reference image in the first training deformation field by using the first transformation formula, so as to obtain the first standard training image contained in the same image space. A first interpolation result of a target registration image; wherein, the first transformation formula is:
    X'=ax'+by',Y'=cx'+dy'X'=ax'+by', Y'=cx'+dy'
    其中:in:
    x'为所述第一标准训练图像中像素的横坐标;x' is the abscissa of the pixel in the first standard training image;
    y'为所述第一标准训练图像中像素的纵坐标;y' is the ordinate of the pixel in the first standard training image;
    X'为所述第一参考图像中像素的横坐标;X' is the abscissa of the pixel in the first reference image;
    Y'为所述第一参考图像中像素的纵坐标;Y' is the ordinate of the pixel in the first reference image;
    a、b、c和d均为变换参数值,a、b、c和d构成的变换矩阵
    Figure PCTCN2021096702-appb-100001
    a, b, c and d are all transformation parameter values, and the transformation matrix formed by a, b, c and d
    Figure PCTCN2021096702-appb-100001
  7. 如权利要求3所述的医学影像配准方法,其中,所述通过第一卷积神经网络对所有所述第一配准对进行变形后,得到多个第一训练变形场,包括:The medical image registration method according to claim 3, wherein after deforming all the first registration pairs through the first convolutional neural network, a plurality of first training deformation fields are obtained, comprising:
    将所有所述第一配准对输入至所述卷积神经网络后,通过与所述卷积神经网络对应的第一预设表达式对所有所述第一配准对进行变形,得到与所述第一配准对对应的多个第一训练变形场;所述第一预设表达式包括:After all the first registration pairs are input into the convolutional neural network, all the first registration pairs are deformed by the first preset expression corresponding to the convolutional neural network to obtain the same Multiple first training deformation fields corresponding to the first registration pair; the first preset expression includes:
    Figure PCTCN2021096702-appb-100002
    Figure PCTCN2021096702-appb-100002
    Figure PCTCN2021096702-appb-100003
    Figure PCTCN2021096702-appb-100003
    其中:in:
    x为所述第一标准训练图像的强度;x is the intensity of the first standard training image;
    y为所述第一参考图像的强度;y is the intensity of the first reference image;
    X为x的取值范围;X is the value range of x;
    Y为y的取值范围;Y is the value range of y;
    D m为所述第一标准训练图像和所述第一参考图像像素之间的相似度; D m is the similarity between the first standard training image and the first reference image pixel;
    x i为所述第一标准训练图像第i个像素的强度; x i is the intensity of the ith pixel of the first standard training image;
    y i为所述第一参考图像第i个像素的强度; y i is the intensity of the ith pixel of the first reference image;
    x m为所述第一标准训练图像的平均强度; x m is the average intensity of the first standard training image;
    y m为所述第一参考图像的平均强度; y m is the average intensity of the first reference image;
    p(x)为所述第一标准训练图像概率密度分布函数;p(x) is the probability density distribution function of the first standard training image;
    p(y)为所述第一参考图像概率密度分布函数;p(y) is the probability density distribution function of the first reference image;
    p(x,y)为所述第一标准训练图像和所述第一参考图像的联合概率密度函数。p(x,y) is the joint probability density function of the first standard training image and the first reference image.
  8. 一种医学影像配准装置,其中,包括:A medical image registration device, comprising:
    第一获取模块,用于获取第一医学影像集和第二医学影像集;所述第一医学影像集中包含同一时刻不同视角的多张第一医学影像;所述第二医学影像集中包含不同时刻同一视角的多张第二医学影像;a first acquisition module, configured to acquire a first medical image set and a second medical image set; the first medical image set includes multiple first medical images from different viewing angles at the same time; the second medical image set includes different times Multiple second medical images from the same viewing angle;
    第二获取模块,用于将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像;The second acquisition module is configured to, after the first medical image set is input into the first registration transformation model, obtain the first registration transformation model for performing multi-view fusion on the first medical image set outputted by the first registration transformation model. a registered image;
    第三获取模块,用于将所述第二医学影像集输入至第二配准变换模型后,获取所述第二配准变换模型输出的对所述第二医学影像集进行多时序融合的第二配准图像;The third acquiring module is configured to acquire the first multi-sequence fusion of the second medical image set output from the second registration transformation model after the second medical image set is input into the second registration transformation model. Two registered images;
    融合模块,用于对所述第一配准图像和所述第二配准图像进行融合,得到第三配准图像;所述第三配准图像反映出所述第一医学影像集和所述第二医学影像集的生长追踪结果。a fusion module, configured to fuse the first registration image and the second registration image to obtain a third registration image; the third registration image reflects the first medical image set and the Growth tracking results for the second set of medical images.
  9. 一种计算机设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机可读指令,其中,所述处理器执行所述计算机可读指令时实现如下步骤:A computer device comprising a memory, a processor, and computer-readable instructions stored in the memory and executable on the processor, wherein the processor implements the following steps when executing the computer-readable instructions:
    获取第一医学影像集和第二医学影像集;所述第一医学影像集中包含同一时刻不同视角的多张第一医学影像;所述第二医学影像集中包含不同时刻同一视角的多张第二医学影像;Acquiring a first medical image set and a second medical image set; the first medical image set includes multiple first medical images from different viewing angles at the same time; the second medical image set includes multiple second medical images from the same viewing angle at different times Medical Imaging;
    将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像;After the first medical image set is input into the first registration transformation model, a first registration image output by the first registration transformation model for performing multi-view fusion on the first medical image set is obtained;
    将所述第二医学影像集输入至第二配准变换模型后,获取所述第二配准变换模型输出的对所述第二医学影像集进行多时序融合的第二配准图像;After the second medical image set is input into the second registration transformation model, a second registration image output by the second registration transformation model for performing multi-sequence fusion on the second medical image set is obtained;
    对所述第一配准图像和所述第二配准图像进行融合,得到第三配准图像;所述第三配准图像反映出所述第一医学影像集和所述第二医学影像集的生长追踪结果。The first registration image and the second registration image are fused to obtain a third registration image; the third registration image reflects the first medical image set and the second medical image set growth tracking results.
  10. 如权利要求9所述的计算机设备,其中,所述第一配准变换模型包括第一卷积神经网络和第一图像插值模块,所述将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像,包括:10. The computer device of claim 9, wherein the first registration transformation model includes a first convolutional neural network and a first image interpolation module, the first set of medical images being input to the first registration After transforming the model, acquiring the first registration image output by the first registration transformation model for performing multi-view fusion on the first medical image set, including:
    获取所述第一医学影像集中的多张所述第一医学影像;acquiring a plurality of the first medical images in the first medical image set;
    将所有所述第一医学影像输入至第一卷积神经网络后,获取将各所述第一医学影像变形之后生成的多个第一变形场图像;After inputting all the first medical images into the first convolutional neural network, acquiring a plurality of first deformation field images generated after deforming each of the first medical images;
    将所有所述第一变形场图像通过第一图像插值模块进行图像插值后,得到一张所述第一配准图像。After performing image interpolation on all the first deformation field images through the first image interpolation module, a first registration image is obtained.
  11. 如权利要求9所述的计算机设备,其中,所述将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像之前,所述处理器执行所述计算机可读指令时实现如下步骤:The computer device according to claim 9, wherein, after the first medical image set is input into the first registration transformation model, a pair of the first medical images output by the first registration transformation model is acquired. Before collecting the first registration image for multi-view fusion, the processor implements the following steps when executing the computer-readable instructions:
    获取第一标准训练图像以及第一参考图像,将一张所述第一标准训练图像以及一张所述第一参考图像作为一对第一配准对,通过第一卷积神经网络对所有所述第一配准对进行变形后,得到多个第一训练变形场;一个所述第一训练变形场对应一张所述第一标准训练图像;Obtain a first standard training image and a first reference image, use a first standard training image and a first reference image as a pair of first registration After the first registration pair is deformed, a plurality of first training deformation fields are obtained; one of the first training deformation fields corresponds to a piece of the first standard training image;
    将所述第一训练变形场通过第一图像插值模块进行图像插值后,得到第一插值结果;After the first training deformation field is subjected to image interpolation by the first image interpolation module, a first interpolation result is obtained;
    将所述第一插值结果和所述第一参考图像输入至与所述第一卷积神经网络关联的第一预设损失函数后,获取第一训练变形场和第一参考图像的第一损失值;After the first interpolation result and the first reference image are input into the first preset loss function associated with the first convolutional neural network, the first loss of the first training deformation field and the first reference image is obtained value;
    在所述第一损失值大于第一预设损失值时,通过反向传播算法不断更新所述第一卷积神经网络的网络参数直至所述第一损失值小于或等于所述第一预设损失值;When the first loss value is greater than a first preset loss value, the network parameters of the first convolutional neural network are continuously updated through a back-propagation algorithm until the first loss value is less than or equal to the first preset loss value loss value;
    在所述第一损失值小于或等于第一预设损失值时,确定所述第一配准变换模型训练完成,得到所述第一配准变换模型。When the first loss value is less than or equal to the first preset loss value, it is determined that the training of the first registration transformation model is completed, and the first registration transformation model is obtained.
  12. 如权利要求9所述的计算机设备,其中,所述第二配准变换模型包括第二卷积神经网络和第二图像插值模块,所述将所述第二医学影像集输入至第二配准变换模型后,获取所述第二配准变换模型输出的对所述第二医学影像集进行多时序融合的第二配准图像,包括:9. The computer device of claim 9, wherein the second registration transformation model includes a second convolutional neural network and a second image interpolation module, the second set of medical images being input to the second registration After transforming the model, acquiring a second registration image output by the second registration transformation model for performing multi-sequence fusion on the second medical image set, including:
    获取所述第二医学影像集中的多张第二配准对;一张所述第二配准对包括一张所述第二医学影像和一张第二参考图像;acquiring a plurality of second registration pairs in the second medical image set; a second registration pair includes a second medical image and a second reference image;
    将所有所述第二配准对输入至第二卷积神经网络后,获取将各所述第二配准对变形之后生成的多个第二变形场;After inputting all the second registration pairs to the second convolutional neural network, acquiring a plurality of second deformation fields generated after deforming each of the second registration pairs;
    将所有所述第二变形场通过第二图像插值模块进行图像插值后,得到所述第二配准图像。After all the second deformation fields are subjected to image interpolation by the second image interpolation module, the second registration image is obtained.
  13. 如权利要求12所述的计算机设备,其中,所述将所有所述第二配准对输入至由卷积神经网络构建成的第一配准变换模型之前,所述处理器执行所述计算机可读指令时实现如下步骤:13. The computer device of claim 12, wherein before said inputting all said second registration pairs to a first registration transformation model constructed by a convolutional neural network, said processor executes said computer-readable The following steps are implemented when reading the command:
    获取第二标准训练图像以及第二参考图像,将一张所述第二标准训练图像以及一张所述第二参考图像作为一对第二配准对,通过第二卷积神经网络对所有所述第二配准对进行变形后,得到多个第二训练变形场;一个所述第二训练变形场对应一张所述第二标准训练图像;Acquire a second standard training image and a second reference image, use one of the second standard training image and one of the second reference images as a pair of second registration pairs, and use the second convolutional neural network for all the After the second registration pair is deformed, a plurality of second training deformation fields are obtained; one of the second training deformation fields corresponds to a piece of the second standard training image;
    将所述第二训练变形场通过第二图像插值模块进行图像插值后,得到第二插值结果;After the second training deformation field is subjected to image interpolation by the second image interpolation module, a second interpolation result is obtained;
    将所述第二插值结果和所述第二参考图像输入至与所述第二卷积神经网络关联的第二预设损失函数后,获取第二训练变形场和第二参考图像的第二损失值;After the second interpolation result and the second reference image are input to a second preset loss function associated with the second convolutional neural network, a second loss of the second training deformation field and the second reference image is obtained value;
    在所述第二损失值大于第二预设损失值时,通过反向传播算法不断更新所述第二卷积神经网络的网络参数直至所述第二损失值小于或等于所述第二预设损失值;When the second loss value is greater than a second preset loss value, the network parameters of the second convolutional neural network are continuously updated through a back-propagation algorithm until the second loss value is less than or equal to the second preset loss value loss value;
    在所述第二损失值小于或等于第二预设损失值时,确定所述第二配准变换模型训练完成,得到所述第二配准变换模型。When the second loss value is less than or equal to the second preset loss value, it is determined that the training of the second registration transformation model is completed, and the second registration transformation model is obtained.
  14. 如权利要求11所述的计算机设备,其中,所述将所述第一训练变形场通过第一图像插值模块进行图像插值后,得到第一插值结果,包括:The computer device according to claim 11, wherein, after performing image interpolation on the first training deformation field through a first image interpolation module, obtaining a first interpolation result, comprising:
    通过第一变换公式将与所述第一训练变形场中的第一标准训练图像变换到所述第一 训练变形场中的第一参考图像的图像空间中,得到包含在同一个图像空间的第一目标配准图像的第一插值结果;其中,所述第一变换公式为:Transform the first standard training image in the first training deformation field into the image space of the first reference image in the first training deformation field by using the first transformation formula, so as to obtain the first standard training image contained in the same image space. A first interpolation result of a target registration image; wherein, the first transformation formula is:
    X'=ax'+by',Y'=cx'+dy'X'=ax'+by', Y'=cx'+dy'
    其中:in:
    x'为所述第一标准训练图像中像素的横坐标;x' is the abscissa of the pixel in the first standard training image;
    y'为所述第一标准训练图像中像素的纵坐标;y' is the ordinate of the pixel in the first standard training image;
    X'为所述第一参考图像中像素的横坐标;X' is the abscissa of the pixel in the first reference image;
    Y'为所述第一参考图像中像素的纵坐标;Y' is the ordinate of the pixel in the first reference image;
    a、b、c和d均为变换参数值,a、b、c和d构成的变换矩阵
    Figure PCTCN2021096702-appb-100004
    a, b, c and d are all transformation parameter values, and the transformation matrix formed by a, b, c and d
    Figure PCTCN2021096702-appb-100004
  15. 一个或多个存储有计算机可读指令的可读存储介质,其中,所述计算机可读指令被一个或多个处理器执行时,使得所述一个或多个处理器执行如下步骤:One or more readable storage media storing computer-readable instructions, wherein the computer-readable instructions, when executed by one or more processors, cause the one or more processors to perform the following steps:
    获取第一医学影像集和第二医学影像集;所述第一医学影像集中包含同一时刻不同视角的多张第一医学影像;所述第二医学影像集中包含不同时刻同一视角的多张第二医学影像;Acquiring a first medical image set and a second medical image set; the first medical image set includes multiple first medical images from different viewing angles at the same time; the second medical image set includes multiple second medical images from the same viewing angle at different times Medical Imaging;
    将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像;After the first medical image set is input into the first registration transformation model, a first registration image output by the first registration transformation model for performing multi-view fusion on the first medical image set is obtained;
    将所述第二医学影像集输入至第二配准变换模型后,获取所述第二配准变换模型输出的对所述第二医学影像集进行多时序融合的第二配准图像;After the second medical image set is input into the second registration transformation model, a second registration image output by the second registration transformation model for performing multi-sequence fusion on the second medical image set is obtained;
    对所述第一配准图像和所述第二配准图像进行融合,得到第三配准图像;所述第三配准图像反映出所述第一医学影像集和所述第二医学影像集的生长追踪结果。The first registration image and the second registration image are fused to obtain a third registration image; the third registration image reflects the first medical image set and the second medical image set growth tracking results.
  16. 如权利要求15所述的可读存储介质,其中,所述第一配准变换模型包括第一卷积神经网络和第一图像插值模块,所述将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像,包括:16. The readable storage medium of claim 15, wherein the first registration transformation model comprises a first convolutional neural network and a first image interpolation module, and wherein the first set of medical images is input to the first After registering the transformation model, acquiring the first registration image output by the first registration transformation model for performing multi-view fusion on the first medical image set, including:
    获取所述第一医学影像集中的多张所述第一医学影像;acquiring a plurality of the first medical images in the first medical image set;
    将所有所述第一医学影像输入至第一卷积神经网络后,获取将各所述第一医学影像变形之后生成的多个第一变形场图像;After inputting all the first medical images into the first convolutional neural network, acquiring a plurality of first deformation field images generated by deforming each of the first medical images;
    将所有所述第一变形场图像通过第一图像插值模块进行图像插值后,得到一张所述第一配准图像。After performing image interpolation on all the first deformation field images through the first image interpolation module, a first registration image is obtained.
  17. 如权利要求15所述的可读存储介质,其中,所述将所述第一医学影像集输入至第一配准变换模型后,获取所述第一配准变换模型输出的对所述第一医学影像集进行多视角融合的第一配准图像之前,所述计算机可读指令被一个或多个处理器执行时,使得所述一个或多个处理器执行如下步骤:The readable storage medium according to claim 15, wherein after the first medical image set is input into the first registration transformation model, a pair of the first registration transformation model outputted by the first registration transformation model is obtained. Before the medical image set performs the first registration image for multi-view fusion, when the computer-readable instructions are executed by one or more processors, the one or more processors perform the following steps:
    获取第一标准训练图像以及第一参考图像,将一张所述第一标准训练图像以及一张所述第一参考图像作为一对第一配准对,通过第一卷积神经网络对所有所述第一配准对进行变形后,得到多个第一训练变形场;一个所述第一训练变形场对应一张所述第一标准训练图像;Obtain a first standard training image and a first reference image, use a first standard training image and a first reference image as a pair of first registration After the first registration pair is deformed, a plurality of first training deformation fields are obtained; one of the first training deformation fields corresponds to a piece of the first standard training image;
    将所述第一训练变形场通过第一图像插值模块进行图像插值后,得到第一插值结果;After the first training deformation field is subjected to image interpolation by the first image interpolation module, a first interpolation result is obtained;
    将所述第一插值结果和所述第一参考图像输入至与所述第一卷积神经网络关联的第一预设损失函数后,获取第一训练变形场和第一参考图像的第一损失值;After the first interpolation result and the first reference image are input into the first preset loss function associated with the first convolutional neural network, the first loss of the first training deformation field and the first reference image is obtained value;
    在所述第一损失值大于第一预设损失值时,通过反向传播算法不断更新所述第一卷积神经网络的网络参数直至所述第一损失值小于或等于所述第一预设损失值;When the first loss value is greater than a first preset loss value, the network parameters of the first convolutional neural network are continuously updated through a back-propagation algorithm until the first loss value is less than or equal to the first preset loss value loss value;
    在所述第一损失值小于或等于第一预设损失值时,确定所述第一配准变换模型训练完成,得到所述第一配准变换模型。When the first loss value is less than or equal to the first preset loss value, it is determined that the training of the first registration transformation model is completed, and the first registration transformation model is obtained.
  18. 如权利要求15所述的可读存储介质,其中,所述第二配准变换模型包括第二卷积神经网络和第二图像插值模块,所述将所述第二医学影像集输入至第二配准变换模型后,获取所述第二配准变换模型输出的对所述第二医学影像集进行多时序融合的第二配准图像,包括:16. The readable storage medium of claim 15, wherein the second registration transformation model comprises a second convolutional neural network and a second image interpolation module, the second set of medical images being input to the second After registering the transformation model, acquiring the second registration image output by the second registration transformation model for performing multi-sequence fusion on the second medical image set, including:
    获取所述第二医学影像集中的多张第二配准对;一张所述第二配准对包括一张所述第二医学影像和一张第二参考图像;acquiring a plurality of second registration pairs in the second medical image set; a second registration pair includes a second medical image and a second reference image;
    将所有所述第二配准对输入至第二卷积神经网络后,获取将各所述第二配准对变形之后生成的多个第二变形场;After inputting all the second registration pairs to the second convolutional neural network, acquiring a plurality of second deformation fields generated after deforming each of the second registration pairs;
    将所有所述第二变形场通过第二图像插值模块进行图像插值后,得到所述第二配准图像。After all the second deformation fields are subjected to image interpolation by the second image interpolation module, the second registration image is obtained.
  19. 如权利要求18所述的可读存储介质,其中,所述将所有所述第二配准对输入至由卷积神经网络构建成的第一配准变换模型之前,所述计算机可读指令被一个或多个处理器执行时,使得所述一个或多个处理器执行如下步骤:19. The readable storage medium of claim 18, wherein the computer readable instructions are executed by a prior to inputting all of the second registration pairs to a first registration transformation model constructed by a convolutional neural network. When executed by one or more processors, the one or more processors are caused to perform the following steps:
    获取第二标准训练图像以及第二参考图像,将一张所述第二标准训练图像以及一张所述第二参考图像作为一对第二配准对,通过第二卷积神经网络对所有所述第二配准对进行变形后,得到多个第二训练变形场;一个所述第二训练变形场对应一张所述第二标准训练图像;Obtain a second standard training image and a second reference image, use one of the second standard training image and one of the second reference images as a pair of second registration pairs, and use the second convolutional neural network to register all the After the second registration pair is deformed, a plurality of second training deformation fields are obtained; one of the second training deformation fields corresponds to a piece of the second standard training image;
    将所述第二训练变形场通过第二图像插值模块进行图像插值后,得到第二插值结果;After the second training deformation field is subjected to image interpolation by the second image interpolation module, a second interpolation result is obtained;
    将所述第二插值结果和所述第二参考图像输入至与所述第二卷积神经网络关联的第二预设损失函数后,获取第二训练变形场和第二参考图像的第二损失值;After the second interpolation result and the second reference image are input to a second preset loss function associated with the second convolutional neural network, a second loss of the second training deformation field and the second reference image is obtained value;
    在所述第二损失值大于第二预设损失值时,通过反向传播算法不断更新所述第二卷积神经网络的网络参数直至所述第二损失值小于或等于所述第二预设损失值;When the second loss value is greater than a second preset loss value, the network parameters of the second convolutional neural network are continuously updated through a back-propagation algorithm until the second loss value is less than or equal to the second preset loss value loss value;
    在所述第二损失值小于或等于第二预设损失值时,确定所述第二配准变换模型训练完成,得到所述第二配准变换模型。When the second loss value is less than or equal to the second preset loss value, it is determined that the training of the second registration transformation model is completed, and the second registration transformation model is obtained.
  20. 如权利要求17所述的可读存储介质,其中,所述将所述第一训练变形场通过第一图像插值模块进行图像插值后,得到第一插值结果,包括:The readable storage medium according to claim 17, wherein, after performing image interpolation on the first training deformation field by a first image interpolation module, obtaining a first interpolation result, comprising:
    通过第一变换公式将与所述第一训练变形场中的第一标准训练图像变换到所述第一训练变形场中的第一参考图像的图像空间中,得到包含在同一个图像空间的第一目标配准图像的第一插值结果;其中,所述第一变换公式为:Transform the first standard training image in the first training deformation field into the image space of the first reference image in the first training deformation field by using the first transformation formula, so as to obtain the first standard training image contained in the same image space. A first interpolation result of a target registration image; wherein, the first transformation formula is:
    X'=ax'+by',Y'=cx'+dy'X'=ax'+by', Y'=cx'+dy'
    其中:in:
    x'为所述第一标准训练图像中像素的横坐标;x' is the abscissa of the pixel in the first standard training image;
    y'为所述第一标准训练图像中像素的纵坐标;y' is the ordinate of the pixel in the first standard training image;
    X'为所述第一参考图像中像素的横坐标;X' is the abscissa of the pixel in the first reference image;
    Y'为所述第一参考图像中像素的纵坐标;Y' is the ordinate of the pixel in the first reference image;
    a、b、c和d均为变换参数值,a、b、c和d构成的变换矩阵
    Figure PCTCN2021096702-appb-100005
    a, b, c and d are all transformation parameter values, and the transformation matrix formed by a, b, c and d
    Figure PCTCN2021096702-appb-100005
PCT/CN2021/096702 2021-02-25 2021-05-28 Medical image registration method and apparatus, computer device, and storage medium WO2022178997A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202110210897.0 2021-02-25
CN202110210897.0A CN112837357B (en) 2021-02-25 2021-02-25 Medical image registration method, device, computer equipment and storage medium

Publications (1)

Publication Number Publication Date
WO2022178997A1 true WO2022178997A1 (en) 2022-09-01

Family

ID=75933377

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2021/096702 WO2022178997A1 (en) 2021-02-25 2021-05-28 Medical image registration method and apparatus, computer device, and storage medium

Country Status (2)

Country Link
CN (1) CN112837357B (en)
WO (1) WO2022178997A1 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112837357B (en) * 2021-02-25 2024-03-05 平安科技(深圳)有限公司 Medical image registration method, device, computer equipment and storage medium
CN114373004B (en) * 2022-01-13 2024-04-02 强联智创(北京)科技有限公司 Dynamic image registration method
CN116402865B (en) * 2023-06-06 2023-09-15 之江实验室 Multi-mode image registration method, device and medium using diffusion model

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108053428A (en) * 2017-12-28 2018-05-18 上海联影医疗科技有限公司 A kind of method for registering of medical image, device and equipment
US20180360313A1 (en) * 2017-06-20 2018-12-20 Siemens Healthcare Gmbh Deep-learnt tissue deformation for medical imaging
CN109658377A (en) * 2018-10-31 2019-04-19 泰格麦迪(北京)医疗科技有限公司 A kind of breast MRI lesion region detection method based on multi-dimension information fusion
CN109801251A (en) * 2019-01-28 2019-05-24 李�杰 The fusion method of medical image and the image detecting method learnt based on fusion medical image
CN110599528A (en) * 2019-09-03 2019-12-20 济南大学 Unsupervised three-dimensional medical image registration method and system based on neural network
CN112001889A (en) * 2020-07-22 2020-11-27 杭州依图医疗技术有限公司 Medical image processing method and device and medical image display method
CN112837357A (en) * 2021-02-25 2021-05-25 平安科技(深圳)有限公司 Medical image registration method and device, computer equipment and storage medium

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10043272B2 (en) * 2014-09-16 2018-08-07 Esaote S.P.A. Method and apparatus for acquiring and fusing ultrasound images with pre-acquired images
CN109767461B (en) * 2018-12-28 2021-10-22 上海联影智能医疗科技有限公司 Medical image registration method and device, computer equipment and storage medium

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180360313A1 (en) * 2017-06-20 2018-12-20 Siemens Healthcare Gmbh Deep-learnt tissue deformation for medical imaging
CN108053428A (en) * 2017-12-28 2018-05-18 上海联影医疗科技有限公司 A kind of method for registering of medical image, device and equipment
CN109658377A (en) * 2018-10-31 2019-04-19 泰格麦迪(北京)医疗科技有限公司 A kind of breast MRI lesion region detection method based on multi-dimension information fusion
CN109801251A (en) * 2019-01-28 2019-05-24 李�杰 The fusion method of medical image and the image detecting method learnt based on fusion medical image
CN110599528A (en) * 2019-09-03 2019-12-20 济南大学 Unsupervised three-dimensional medical image registration method and system based on neural network
CN112001889A (en) * 2020-07-22 2020-11-27 杭州依图医疗技术有限公司 Medical image processing method and device and medical image display method
CN112837357A (en) * 2021-02-25 2021-05-25 平安科技(深圳)有限公司 Medical image registration method and device, computer equipment and storage medium

Also Published As

Publication number Publication date
CN112837357A (en) 2021-05-25
CN112837357B (en) 2024-03-05

Similar Documents

Publication Publication Date Title
WO2020238734A1 (en) Image segmentation model training method and apparatus, computer device, and storage medium
WO2022178997A1 (en) Medical image registration method and apparatus, computer device, and storage medium
US11270169B2 (en) Image recognition method, storage medium and computer device
US11842487B2 (en) Detection model training method and apparatus, computer device and storage medium
US10706333B2 (en) Medical image analysis method, medical image analysis system and storage medium
Jin et al. 3D PBV-Net: an automated prostate MRI data segmentation method
CN110599528B (en) Unsupervised three-dimensional medical image registration method and system based on neural network
WO2021164306A1 (en) Image classification model training method, apparatus, computer device, and storage medium
CN110276745B (en) Pathological image detection algorithm based on generation countermeasure network
WO2021233017A1 (en) Image processing method and apparatus, and device and computer-readable storage medium
US20220076052A1 (en) Similarity determining method and device, network training method and device, search method and device, and electronic device and storage medium
CN110751187B (en) Training method of abnormal area image generation network and related product
CN112509119A (en) Spatial data processing and positioning method and device for temporal bone and electronic equipment
CN115393361B (en) Skin disease image segmentation method, device, equipment and medium with low annotation cost
CN110866909A (en) Training method of image generation network, image prediction method and computer equipment
CN112530550A (en) Image report generation method and device, computer equipment and storage medium
CN116469545A (en) Method, device and equipment for carrying out auxiliary disease prediction by using medical image
Tang et al. Lesion segmentation and RECIST diameter prediction via click-driven attention and dual-path connection
Li et al. Retinal vessel segmentation by using AFNet
CN111582449B (en) Training method, device, equipment and storage medium of target domain detection network
CN112164447A (en) Image processing method, device, equipment and storage medium
CN111192252B (en) Image segmentation result optimization method and device, intelligent terminal and storage medium
CN111598144B (en) Training method and device for image recognition model
CN114169467A (en) Image annotation method, electronic device and storage medium
Ezzat et al. Automatic liver segmentation in computed tomography scans using deep semantic segmentation

Legal Events

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

Ref document number: 21927430

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 21927430

Country of ref document: EP

Kind code of ref document: A1