CN109903269A - 确定脊柱横断面图像的异常类型的方法及计算设备 - Google Patents

确定脊柱横断面图像的异常类型的方法及计算设备 Download PDF

Info

Publication number
CN109903269A
CN109903269A CN201910069196.2A CN201910069196A CN109903269A CN 109903269 A CN109903269 A CN 109903269A CN 201910069196 A CN201910069196 A CN 201910069196A CN 109903269 A CN109903269 A CN 109903269A
Authority
CN
China
Prior art keywords
sectional image
backbone cross
image
convolution
backbone
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910069196.2A
Other languages
English (en)
Other versions
CN109903269B (zh
Inventor
刘星宇
张逸凌
冯世庆
张蒂
刘慧烨
张云东
安奕成
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Liu Xingyu
Longwood Valley Medtech Co Ltd
Original Assignee
Longwood Valley Medtech Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Longwood Valley Medtech Co Ltd filed Critical Longwood Valley Medtech Co Ltd
Priority to CN201910069196.2A priority Critical patent/CN109903269B/zh
Publication of CN109903269A publication Critical patent/CN109903269A/zh
Application granted granted Critical
Publication of CN109903269B publication Critical patent/CN109903269B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Analysis (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

本发明公开了确定脊柱横断面图像的异常类型的方法及计算设备。其中方法包括步骤:从脊柱横断面图像中获取异常区域;确定该脊柱横断面图像中的至少一个关键点;以及结合异常区域和关键点以确定该脊柱横断面图像的异常类型。

Description

确定脊柱横断面图像的异常类型的方法及计算设备
技术领域
本发明涉及图像处理技术领域,尤其涉及一种确定脊柱横断面图像的异常类型的方法及计算设备。
背景技术
随着计算机视觉的迅速发展,利用计算机技术来辅助骨科精准手术的技术也逐渐增多。最常见的应用如通过磁共振(Magnetic Resonance Imaging,MRI)成像,从磁共振图像中可以得到物质的多种物理特性参数,如质子密度、自旋-晶格驰豫时间T1、自旋-自旋驰豫时间T2、扩散系数、磁化系数、化学位移等等,以将其应用于全身各系统的成像诊断。其中主要的应用之一是脊柱成像,以获得人体的脊柱矢状面图像和脊柱横断面图像,进而应用在各种脊柱疾病诊断中,最具代表性的应用如腰椎间盘突出的辅助诊断。
然而,目前在利用MRI图像进行辅助诊断时多采用人工判读的方式,其不可避免地引入以下问题:1)存在误差,由于医务人员的业务水平及不同区域诊疗水平的参差不齐,同一MRI图像会得出不同结果;2)腰椎间盘突出的分型较复杂,不易在短时间内进行判读,致使医务工作者负担较重,且效率不高。
鉴于此,需要新的对脊柱的相关图像进行处理的方案,以更好地辅助专业医生实现对各种脊柱疾病的诊断。
发明内容
为此,本发明提供了一种确定脊柱横断面图像的异常类型的方法及计算设备,以力图解决或者至少缓解上面存在的至少一个问题。
根据本发明的一个方面,提供了一种确定脊柱横断面图像的异常类型的方法,该方法适于在计算设备中执行,包括步骤:从脊柱横断面图像中获取异常区域;确定脊柱横断面图像中的至少一个关键点;以及结合异常区域和关键点以确定脊柱横断面图像的异常类型。
可选地,在根据本发明的方法中,结合异常区域和关键点以确定脊柱横断面图像的异常类型的步骤还包括:根据所确定的至少一个关键点生成网格;以及根据异常区域在网格中的位置,确定脊柱横断面图像的异常类型。
可选地,在根据本发明的方法中,从脊柱横断面图像中获取异常区域的步骤还包括:利用第一神经网络来处理脊柱横断面图像,以判断脊柱横断面图像是否异常;以及若确定该脊柱横断面图像异常,则从该脊柱横断面图像中获取异常区域。
可选地,根据本发明的方法还包括步骤:从脊柱矢状面图像中确定出存在异常的感兴趣区域;以及获取与存在异常的感兴趣区域相对应的脊柱横断面图像作为脊柱横断面图像。
可选地,在根据本发明的方法中,若确定脊柱横断面图像异常,则从脊柱横断面图像中获取异常区域的步骤包括:利用第二神经网络来处理脊柱横断面图像,以输出第一图像,在第一图像中具有特定像素值的区域为异常区域;以及对第一图像进行插值处理,以获得第二图像,第二图像与该脊柱横断面图像具有相同尺寸。
可选地,在根据本发明的方法中,确定脊柱横断面图像中的至少一个关键点的步骤包括:利用第三神经网络来处理脊柱横断面图像,以输出包含至少一个关键点位置的至少一个热图;以及将这至少一个热图合并生成预测图,其中预测图包含至少一个关键点。
可选地,在根据本发明的方法中,第三神经网络通过卷积处理层和空间融合层耦接而成。
可选地,在根据本发明的方法中,卷积处理层包含8层卷积,且各层卷积中采用的卷积核大小依次为:5*5、5*5、5*5、5*5、9*9、1*1、1*1、1*1;以及空间融合层包含5层卷积,且各层卷积中采用的卷积核大小依次为:7*7、13*13、13*13、1*1、1*1;并且,将卷积处理层中第三层卷积和第七层卷积的输出结合后,作为空间融合层的输入。
可选地,根据本发明的方法还包括预先训练生成第三神经网络的步骤:构造预训练的第三神经网络,并设置初始的网络参数;将训练图像输入预训练的第三神经网络,通过损失函数计算网络的输出所对应的损失值;根据损失值调整网络参数;以及重复迭代计算损失值和调整网络参数的步骤,直到满足预定条件时训练结束。
可选地,在根据本发明的方法中,将训练图像输入预训练的第三神经网络,利用损失函数计算网络的输出所对应的损失值的步骤包括:将训练图像输入预训练的第三神经网络,通过损失函数分别计算网络中预定层的输出所对应的损失;以及结合各层的损失得到该网络所对应的损失值,其中,预定层包括卷积处理层中第4、6、8层卷积、以及空间融合层中第5层卷积。
可选地,在根据本发明的方法中,第一神经网络通过至少一个深度残差网络耦合生成,如,第一神经网络采用34层的ResNet网络结构。
可选地,在根据本发明的方法中,第二神经网络通过收缩路径和扩张路径耦接生成;其中,收缩路径包含多个第一卷积处理块,且每个第一卷积处理块包含两个3*3的卷积核、激活函数和一个最大池化单元,扩张路径包含多个第二卷积处理块和输出块,且每个第二卷积处理块包含两个3*3的卷积核、激活函数和一个上采样单元,以及输出块包含1*1的卷积核,同时,收缩路径还适于分别对经各第一卷积处理块处理后所生成的特征图进行裁剪,再将裁剪后的特征图输入到对应的扩张路径的第二卷积处理块中。
根据本发明的又一方面,提供了一种计算设备,包括:一个或多个处理器;和存储器;一个或多个程序,其中一个或多个程序存储在存储器中并被配置为由一个或多个处理器执行,一个或多个程序包括用于执行如上所述方法中的任一方法的指令。
根据本发明的再一方面,提供了一种存储一个或多个程序的计算机可读存储介质,一个或多个程序包括指令,指令当计算设备执行时,使得计算设备执行如上所述的方法中的任一方法。
综上,根据本发明的方案,首先从脊柱横断面图像中获取异常区域;同时,对脊柱横断面图像进行关键点检测,确定出其中的多个关键点;之后,结合异常区域和关键点,就可以确定出该脊柱横断面图像的异常类型。根据本发明的方案基于卷积神经网络,大大提高了处理脊柱横断面图像的准确率,且节省了确定异常类型的时间,能够很好地辅助专业医生实现对各种脊柱疾病的精确诊断。
附图说明
为了实现上述以及相关目的,本文结合下面的描述和附图来描述某些说明性方面,这些方面指示了可以实践本文所公开的原理的各种方式,并且所有方面及其等效方面旨在落入所要求保护的主题的范围内。通过结合附图阅读下面的详细描述,本公开的上述以及其它目的、特征和优势将变得更加明显。遍及本公开,相同的附图标记通常指代相同的部件或元素。
图1示出了根据本发明一个实施例的计算设备100的构造示意图;
图2示出了根据本发明一个实施例的确定脊柱横断面图像的异常类型的方法200的流程示意图;
图3A和图3B分别示出了脊柱矢状面图像和脊柱横断面图像的示意图;
图4A和图4B分别示出了根据本发明实施方式的正常的脊柱横断面图像和异常的脊柱横断面图像的示意图;
图5示出了根据本发明一个实施例的脊柱横断面图像中的关键点的示意图;
图6示出了根据本发明一个实施例的第三神经网络600的结构示意图;
图7示出了根据本发明一个实施例的脊柱横断面图像的异常类型的示意图;
图8示出了根据本发明实施例的应用方法200来确定脊柱横断面图像的异常类型的过程示意图;以及
图9示出了根据本发明实施例的确定脊柱横断面图像的异常类型的示意性示意图。
具体实施方式
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例,然而应当理解,可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本公开,并且能够将本公开的范围完整的传达给本领域的技术人员。
图1是示例计算设备100的框图。在基本的配置102中,计算设备100典型地包括系统存储器106和一个或者多个处理器104。存储器总线108可以用于在处理器104和系统存储器106之间的通信。
取决于期望的配置,处理器104可以是任何类型的处理器,包括但不限于:微处理器(μP)、微控制器(μC)、数字信息处理器(DSP)或者它们的任何组合。处理器104可以包括诸如一级高速缓存110和二级高速缓存112之类的一个或者多个级别的高速缓存、处理器核心114和寄存器116。示例的处理器核心114可以包括运算逻辑单元(ALU)、浮点数单元(FPU)、数字信号处理核心(DSP核心)或者它们的任何组合。示例的存储器控制器118可以与处理器104一起使用,或者在一些实现中,存储器控制器118可以是处理器104的一个内部部分。
取决于期望的配置,系统存储器106可以是任意类型的存储器,包括但不限于:易失性存储器(诸如RAM)、非易失性存储器(诸如ROM、闪存等)或者它们的任何组合。系统存储器106可以包括操作系统120、一个或者多个应用122以及程序数据124。在一些实施方式中,应用122可以布置为在操作系统上利用程序数据124进行操作。
计算设备100还可以包括有助于从各种接口设备(例如,输出设备142、外设接口144和通信设备146)到基本配置102经由总线/接口控制器130的通信的接口总线140。示例的输出设备142包括图形处理单元148和音频处理单元150。它们可以被配置为有助于经由一个或者多个A/V端口152与诸如显示器或者扬声器之类的各种外部设备进行通信。示例外设接口144可以包括串行接口控制器154和并行接口控制器156,它们可以被配置为有助于经由一个或者多个I/O端口158和诸如输入设备(例如,键盘、鼠标、笔、语音输入设备、图像输入设备)或者其他外设(例如打印机、扫描仪等)之类的外部设备进行通信。示例的通信设备146可以包括网络控制器160,其可以被布置为便于经由一个或者多个通信端口164与一个或者多个其他计算设备162通过网络通信链路的通信。
网络通信链路可以是通信介质的一个示例。通信介质通常可以体现为在诸如载波或者其他传输机制之类的调制数据信号中的计算机可读指令、数据结构、程序模块,并且可以包括任何信息递送介质。“调制数据信号”可以是这样的信号,它的数据集中的一个或者多个或者它的改变可以在信号中以编码信息的方式进行。作为非限制性的示例,通信介质可以包括诸如有线网络或者专线网络之类的有线介质,以及诸如声音、射频(RF)、微波、红外(IR)或者其它无线介质在内的各种无线介质。这里使用的术语计算机可读介质可以包括存储介质和通信介质二者。在一些实施例中,计算机可读介质中存储一个或多个程序,这一个或多个程序中包括执行某些方法的指令。
计算设备100可以实现为小尺寸便携(或者移动)电子设备的一部分,这些电子设备可以是诸如蜂窝电话、数码照相机、个人数字助理(PDA)、个人媒体播放器设备、无线网络浏览设备、个人头戴设备、应用专用设备、或者可以包括上面任何功能的混合设备。当然,计算设备100也可以实现为包括桌面计算机和笔记本计算机配置的个人计算机,或者是具有上述配置的服务器。本发明的实施方式对此均不作限制。
在根据本发明的一些实施例中,计算设备100被配置为执行确定脊柱横断面图像的异常类型的方法200,程序数据124中包含了用于执行该方法200的相关指令。
图2示出了根据本发明一个实施例的确定脊柱横断面图像的异常类型的方法200的流程图。
在结合图2对方法200进行详细阐述之前,首先需要说明的是,根据本发明的实施例,关于脊柱的MRI图像有两种,一种是脊柱矢状面图像,另一种是脊柱横断面图像。如图3A示出的是脊柱矢状面图像,图3B示出的是横断面图像,可以将图3B看成是图3A的一个水平横切视角。
如图2所示,方法200始于步骤S210。在步骤S210中,从脊柱横断面图像中获取异常区域。
根据本发明的实施方式,在从脊柱横断面图像中获取异常区域之前,先判断该脊柱横断面图像是否异常。在确定该脊柱横断面图像异常时,再从中获取异常区域;否则,不对其进行分割。根据本发明的实施方式,脊柱图像的异常主要针对椎间盘突出,若脊柱横断面图像中包含的间盘区域不存在突出,则认为该图像正常;反之,若脊柱横断面图像中包含的间盘区域存在突出,则认为该图像异常(医学上通常认为这种情况可能存在病变)。如图4A和图4B,分别示出了根据本发明实施方式的正常的脊柱横断面图像和异常的脊柱横断面图像。
例如,可以采用传统的图像处理算法或是深度学习方法来对脊柱横断面图像进行初步筛选,以判断脊柱横断面图像是否异常。在一种实施例中,通过构建深度学习网络来判断脊柱横断面图像是否异常。具体地,利用第一神经网络来处理该脊柱横断面图像,以判断脊柱横断面图像是否异常;若经判断后确定该脊柱横断面图像异常,则从该脊柱横断面图像中获取异常区域。
根据本发明的一种实施例,基于卷积神经网络,第一神经网络通过至少一个深度残差网络耦合生成。通过采用残差网络的结构,在不显著增加网络深度的前提下,提升网络的准确率。在本发明的实施例中,第一神经网络采用34层的ResNet网络结构,输入脊柱横断面图像,输出两个类别,例如,若输出为“1”,就表示输入的脊柱横断面图像异常,若输出为“0”,就表示输入的脊柱横断面图像正常。当然不限于此。第一神经网络的结构可参考ResNet网络,鉴于ResNet网络属于本领域技术人员已知内容,故此处不再赘述。
此外,方法200还包含了训练生成第一神经网络的步骤。以下示出根据本发明一种实施例的训练生成第一神经网络的过程。
首先,采集一些脊柱横断面图像作为训练图像(例如,采集1500张脊柱横断面图像),并在专业医师的指导下对这些训练图像进行标注,标注方式例如是:若脊柱横断面图像正常,则标注为0;若脊柱横断面图像异常,则标注为1。这样,每个训练图像都具有相应的标注数据。同时,构造预训练的第一神经网络,并设置初始的网络参数。之后,将训练图像输入到预训练的第一神经网络中,并根据输出的结果和训练图像对应的标注数据对网络参数进行微调,以生成新的第一神经网络,重复上述步骤,直到第一神经网络的输出与标注数据的损失满足预定条件(本发明的实施例对损失函数的选取不做限制),即,第一神经网络的输出与标注数据接近一致,训练结束。此时生成的第一神经网络就作为第一神经网络。
接下来,在确定脊柱横断面图像异常后,利用第二神经网络来处理该脊柱横断面图像,以输出第一图像。在一些实施例中,在第一图像中具有特定像素值的区域就是异常区域,例如,以二值化的方式来表示第一图像,当像素值为1时,表示该像素属于异常区域,当像素值为0时,表示该像素不属于异常区域,但不限于此。根据本发明的实施例,第二神经网络采用卷积神经网络,如U-net网络结构,输入脊柱横断面图像,输出分割出异常区域的图像,即第一图像。应当了解,经U-net网络结构输出的第一图像的尺寸要小于输入图像的尺寸,故在得到第一图像后,还要对第一图像进行插值处理,以获得第二图像,以保证第二图像与输入的脊柱横断面图像具有相同尺寸。插值处理例如可以是双线性插值处理,但不限于此。
以下对第二神经网络的网络结构进行进一步阐述。参考U-net网络,第二神经网络通过一个收缩路径和一个扩张路径耦接生成。其中,收缩路径遵循典型的卷积网络结构,其包含多个第一卷积处理块,各第一卷积处理块包含两个重复的3*3的卷积核(无填充卷积,unpadded convolution)、激活函数和一个最大池化单元,其中激活函数均使用修正线性单元(rectified linear unit,ReLU)激活函数,最大池化单元包含一个用于下采样(downsample)的步长为2的2*2最大池化操作。另外,在每一个下采样的步骤中,特征通道数量都加倍。在扩张路径中,包含多个第二卷积处理块和输出块。其中,每个第二卷积处理块包含两个3*3的卷积核、激活函数和一个上采样单元。收缩路径还适于分别对经各第一卷积处理块处理后所生成的特征图进行裁剪(由于在每次卷积操作中,边界像素存在缺失问题,因此有必要对特征图进行裁剪),再将裁剪后的特征图输入到对应的扩张路径的第二卷积处理块中。假设收缩路径中包含两个第一卷积处理块,扩张路径中包含两个第二卷积处理块,那么,将第二个第一卷积处理块输出的特征图裁剪后,输入到扩张路径的第一个第二卷积处理块,经第一个第一卷积处理块输出的特征图裁剪后,输入到扩张路径的第二个第二卷积处理块。这样,每个第二卷积处理块都包含对由上一第二卷积处理块输出的特征图进行上采样(upsample)——即用2*2的卷积核进行卷积运算(上卷积,up-convolution)的操作,用于减少一半的特征通道数量;接着级联收缩路径中相应的裁剪后的特征图;再用两个3*3的卷积核进行卷积运算,同样均使用ReLU激活函数。在最后一个输出块,利用1*1的卷积核进行卷积运算,将每个特征向量映射到网络的输出块。
同样地,方法200还包含了训练生成第二神经网络的步骤。以下示出根据本发明一种实施例的训练生成第二神经网络的过程。
在上述用于训练第一神经网络的脊柱横断面图像中,将被标注为异常图像的脊柱横断面图像,作为训练生成第二神经网络的训练图像。在专业医师的指导下、利用标注软件标注出各训练图像中的异常区域。这样,每个训练图像均具有相应的标注数据(标注数据例如是异常区域的位置)。或者,对训练图像进行二值化处理,根据专业医师的指导,将训练图像中属于异常区域的像素的像素值设为1,不属于异常区域的像素的像素值设为0。这样,各像素的像素值就可以作为各训练图像的标注数据,不限于此。同时,构造预训练的第二神经网络,并设置初始的网络参数。之后,将训练图像输入到预训练的第二神经网络中,并根据输出的结果和训练图像对应的标注数据对网络参数进行微调,以生成新的第二神经网络,重复上述步骤,直到第二神经网络的输出与标注数据的损失满足预定条件(本发明的一个实施例中,损失函数采用pixel wise softmax cross entropy,像素级别的交叉熵,但不限于此),即,第二神经网络的输出与标注数据接近一致,训练结束。此时生成的第二神经网络就作为第二神经网络。
至此,经步骤S210就获取了脊柱横断面图像中的异常区域。在另一些实施方式中,会先对脊柱矢状面图像进行判断识别,当确认脊柱矢状面图像中存在异常时,再执行步骤S210,来对与该脊柱矢状面图像中异常部位相对应的脊柱横断面图像做进一步判断,以得到更加精确的预测结果。
具体地,在利用第一神经网络对脊柱横断面图像进行处理(来判断该脊柱横断面图像是否异常)的步骤之前,还包括步骤:先从脊柱矢状面图像中确定存在异常的区域;再获取与该存在异常的区域相对应的脊柱横断面图像,作为要处理的脊柱横断面图像。从脊柱矢状面图像中识别出存在异常的区域的过程,可以通过人工手段来实现,也可以采用计算机视觉方法来实现,以下示出根据本发明一些实施例的实现过程。
第一步,从脊柱矢状面图像中截取出每块椎骨的感兴趣区域,并生成感兴趣区域图像。
在本发明的实施例中,感兴趣区域定义为包含间盘和骨髓交汇处的区域。在一个实施例中,先定位出脊柱矢状面图像中的每一块椎骨,将定位到的椎骨通过矩形框标注出来,矩形框应当包含整块椎骨区域,那么,根据矩形框的四个顶点就可以确定出椎骨的中心点。接着,以相邻椎骨的中心点的连线为边生成对应的正方形,该正方形所包含的区域就是感兴趣区域。
第二步,分别提取每个感兴趣区域图像的至少一个特征,并根据所提取的特征拼接生成每个感兴趣区域图像对应的特征向量。在一种实施例中,所提取的特征应当尽可能包含图像的形状特征、纹理特征和亮度特征。故,提取的特征可以是以下特征中的一种或多种:局部纹理特征、方向梯度的金字塔直方图特征、亮度直方图特征和胡矩特征,且不限于此。
第三步,根据特征向量确定对应的感兴趣区域图像是否存在异常。根据一种实施例,将感兴趣区域图像输入到预先训练生成的识别模型中,经识别模型处理后,输出表征该感兴趣区域图像是否存在异常的值。例如,当输出值为0时,表示对应的感兴趣区域图像正常;当输出值为1时,表示对应的感兴趣区域图像异常。识别模型可以是机器学习模型(如支持向量机)、也可以是其他深度卷积模型,本发明的实施例对此不做限制。
如前文所述,当确定脊柱矢状面图像中包含存在异常的感兴趣区域时,再获取该感兴趣区域所对应的水平横切视图,即脊柱横断面图像,并执行步骤S210,从中进一步分割出异常区域。
随后在步骤S220中,确定该脊柱横断面图像中的至少一个关键点。为便于后期确定异常类型,在本申请的实施方式中,利用关键点检测技术,从脊柱横断面图像中提取出6个关键点,这6个关键点的位置如图5所示,分别位于上下关节突、椎体正中及两椎板汇聚处内缘,由这6个关键点形成一个菱形,且为突出显示,在图5中用红色显示这6个关键点。
以下示出根据本发明一个实施例的检测关键点的过程。首先,利用第三神经网络来处理脊柱横断面图像,以输出包含至少一个关键点位置的至少一个热图;接着,将所输出的这至少一个热图合并,生成预测图,其中预测图包含了所述关键点。在根据本发明的实施例中,要检测的关键点数目为6,故经第三神经网络输出6个热图。
图6示出了根据本发明一个实施例的第三神经网络600的结构示意图。如图6所示,第三神经网络600通过卷积处理层610和空间融合层620耦接而成。其中,卷积处理层610采用典型的卷积网络结构,其由8层卷积构成,分别记作“卷积1”“卷积2”“卷积3”“卷积4”“卷积52”“卷积6”“卷积7”“卷积8”,且每一层卷积中采用的卷积核大小依次为:5*5、5*5、5*5、5*5、9*9、1*1、1*1、1*1。另外,在前两层卷积中,除了卷积处理外,还包括池化处理(池化处理例如是最大值池化,不限于此)。在本实施例中,通过空间融合层620来进一步学习关键点之间的内在联系。如图6所示,将第三层卷积和第七层卷积(即,卷积3和卷积7)的输出结合后,输入到空间融合层620,再经过5层卷积处理,最终输出包含了关键点的热图(heatmap)。空间融合层620中包含5层卷积,分别记作“卷积1_f”“卷积2_f”“卷积3_f”“卷积4_f”“卷积5_f”,且各层卷积中所采用的卷积核大小依次为:7*7、13*13、13*13、1*1、1*1。应当指出,在卷积处理层610和空间融合层620中,各层卷积中还可以包含激活函数(如ReLU),激活函数可以通过设置单独的激活层实现,也可以在构造层对象时通过传递activation参数实现,本发明的实施例对此不做限制。关于卷积神经网络的基本构成属于本领域技术人员已知内容,此处不予赘述。
方法200还包含了预先训练生成第三神经网络600的步骤。以下示出根据本发明一个实施例的训练生成第三神经网络600的过程。
同样地,将上述用于训练生成第二神经网络的脊柱横断面图像,再作为训练生成第三神经网络的训练图像。在专业医师的指导下、利用标注软件标注出各训练图像中的6个关键点,并根据标注出的6个关键点的位置生成高斯热力图(作为目标图像)。同时,构造预训练的第三神经网络(预训练的第三神经网络的网络结构可参考上文关于第三神经网络600的相关描述),并设置初始的网络参数。之后,将训练图像输入到预训练的第三神经网络中,利用损失函数计算第三神经网络的输出相对于目标图像的损失值。在根据本发明的一些实施例中,将训练图像输入到预训练的第三神经网络时,通过损失函数分别计算第三神经网络中预定层的输出所对应的损失,再结合各层的损失得到该检测网络所对应的损失值。如前文所述,在第三神经网络中,卷积处理层由8层卷积构成,空间融合层由5层卷积构成,预定层设为卷积处理层中第4、6、8层、以及空间融合层中第5层(即,图6所示的“卷积4”“卷积6”“卷积8”“卷积5_f”)。通过损失函数分别计算这4层所输出的热图中关键点的位置,相对于目标图像中关键点的位置的损失,得到4个损失,如图6,示意性地示出了根据这4个预定层的输出计算出4个损失(分别记作第一损失、第二损失、第三损失、第四损失),再将这4个损失相结合,就得到了本次训练中第三神经网络的损失值。需要说明的是,可以通过加权求和的方式计算出这4个损失对应的损失值,本发明的实施例对此不做过多限制。之后,根据损失值对网络参数进行微调,以生成新的第三神经网络。重复迭代上述计算损失值的步骤和调整网络参数的步骤,直到满足预定条件(预定条件可以是损失值小于某个阈值,也可以是训练的次数达到某个值,不限于此),总之,当第三神经网络最终输出的预测图与目标图像中关键点的位置接近一致时,训练结束。此时生成的第三神经网络就作为第三神经网络。
随后在步骤S230中,结合异常区域和关键点以确定脊柱横断面图像的异常类型。
根据MSU医学定义,结合异常区域所处位置、及异常区域的大小,将脊柱横断面图像的异常分成10种异常类型。根据一种实施例,这10种异常类型主要是针对椎间盘突出而进行的分类,通过T2加权像上对突出最严重节段与关节突关节间线关系的测量,用1、2、3表示异常区域的大小(即,间盘突出程度),用A、B、C表示异常区域所处位置(即,间盘突出位置)。图7示出了根据本发明实施例的脊柱横断面图像的异常类型的示意图。如图7所示,这10种异常类型分别记作:1-A、1-B、1-C、2-A、2-B、2-C、3-A、3-B、2-AB、3-AB。图7中黑色填充部分表征了脊柱横断面图像中的异常区域。
在根据本发明的实施例中,根据所确定的至少一个关键点生成网格,例如对确定的6个关键点进行编号,按编号顺序连接这6个关键点,最终生成一个3*3的网格。而后,根据异常区域在该网格中的位置,就可以确定出脊柱横断面图像的异常类型。结合图7,当脊柱横断面图像中的异常区域大部分处于网格的第1行第2列中时,就确定该脊柱横断面图像的异常类型是1-A;同样,当脊柱横断面图像的异常区域大部分处于网格第1行第3列中时,就确定该脊柱横断面图像的异常类型是1-B;依次类推,此处不再一一赘述。
通过确定出的异常类型,更加直观地将临床症状体征和脊柱横断面图像中异常区域(即,椎间盘突出的部位)联系起来,这对于临床手术入路和手术方式的选择具有很大的指导意义。
为进一步说明根据方法200确定脊柱横断面图像的异常类型的过程,图8示出了根据本发明实施例的应用方法200来确定脊柱横断面图像的异常类型的过程示意图。
如图8所示,首先,从脊柱矢状面图像801中确定出存在异常的感兴趣区域,而后获取与该存在异常的感兴趣区域相对应的脊柱横断面图像802,以利用该脊柱横断面图像802作更加精确地判断。接着,将脊柱横断面图像802输入第一神经网络,以判断该脊柱横断面图像802是否异常。当确认脊柱横断面图像802异常时,再将该脊柱横断面图像802输入到第二神经网络中,以输出指示了异常区域的脊柱横断面图像803。另外,将脊柱横断面图像802输入到第三神经网络中,以输出确定了关键点的脊柱横断面图像804,并根据所确定的关键点生成网格。最终,将脊柱横断面图像803中的异常区域映射到脊柱横断面图像804中,生成脊柱横断面图像805。
进一步地,图9示出了根据本发明实施例的确定脊柱横断面图像的异常类型的示意性示意图。将指示了异常区域的脊柱横断面图像803(为便于突出显示,用绿色标记了图像中的异常区域)映射到确定出关键点且生成网格的脊柱横断面图像804(为便于突出显示,用红色标记了图像中的关键点,网格如图9中虚线所示)中,生成脊柱横断面图像805。基于脊柱横断面图像805,就可以得出异常区域处于网格中的具体位置,进而确定出异常类型。如图9中所示出的脊柱横断面图像805,可以预测出异常类型为2-AB和1-B。根据所确定的异常类型,可以为专业医生提供参考,以辅助医生进行诊断。
综上,根据本发明的方案,首先从脊柱横断面图像中获取异常区域;同时,对脊柱横断面图像进行关键点检测,确定出其中的多个关键点;之后,结合异常区域和关键点,就可以确定出该脊柱横断面图像的异常类型。本方案基于卷积神经网络,通过训练生成了第一神经网络、第二神经网络和第三神经网络,以完成对脊柱横断面图像的各种处理,大大提高了处理脊柱横断面图像的准确率,减小了由于不同医务人员水平或不同地区诊疗水平的差异而带来的诊疗差异。最终确定的结果可以作为参考,医务人员可以依据该结果进行后续的脊髓受压节段定位及精确地疾病诊断,并由此选择合适的临床手术入路和手术方式。
进一步地,根据本发明的方案,能够节省确定异常类型的时间,提高诊疗效率。辅助专业医生实现对各种脊柱疾病的精确诊断,大大缓解了医务人员负担。
应当理解,为了精简本公开并帮助理解各个发明方面中的一个或多个,在上面对本发明的示例性实施例的描述中,本发明的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而,并不应将该公开的方法解释成反映如下意图:即所要求保护的本发明要求比在每个权利要求中所明确记载的特征更多特征。更确切地说,如下面的权利要求书所反映的那样,发明方面在于少于前面公开的单个实施例的所有特征。因此,遵循具体实施方式的权利要求书由此明确地并入该具体实施方式,其中每个权利要求本身都作为本发明的单独实施例。
本领域那些技术人员应当理解在本文所公开的示例中的设备的模块或单元或组件可以布置在如该实施例中所描述的设备中,或者可替换地可以定位在与该示例中的设备不同的一个或多个设备中。前述示例中的模块可以组合为一个模块或者此外可以分成多个子模块。
本领域那些技术人员可以理解,可以对实施例中的设备中的模块进行自适应性地改变并且把它们设置在与该实施例不同的一个或多个设备中。可以把实施例中的模块或单元或组件组合成一个模块或单元或组件,以及此外可以把它们分成多个子模块或子单元或子组件。除了这样的特征和/或过程或者单元中的至少一些是相互排斥之外,可以采用任何组合对本说明书(包括伴随的权利要求、摘要和附图)中公开的所有特征以及如此公开的任何方法或者设备的所有过程或单元进行组合。除非另外明确陈述,本说明书(包括伴随的权利要求、摘要和附图)中公开的每个特征可以由提供相同、等同或相似目的的替代特征来代替。
本发明一并公开了:
A9、如A8所述的方法,其中,卷积处理层还适于将其中第三层卷积和第七层卷积的输出结合后,作为空间融合层的输入。
A10、如A6-9中任一项所述的方法,还包括预先训练生成第三神经网络的步骤:构造预训练的第三神经网络,并设置初始的网络参数;将训练图像输入预训练的第三神经网络,通过损失函数计算网络的输出所对应的损失值;根据损失值调整网络参数;以及重复迭代计算损失值和调整网络参数的步骤,直到满足预定条件时训练结束。
A11、如A10所述的方法,其中,将训练图像输入预训练的第三神经网络,利用损失函数计算网络的输出所对应的损失值的步骤包括:将训练图像输入预训练的第三神经网络,通过损失函数分别计算网络中预定层的输出所对应的损失;以及结合各层的损失得到该网络所对应的损失值。
A12、如A11所述的方法,其中,预定层包括卷积处理层中第4、6、8层卷积、以及空间融合层中第5层卷积。
A13、如A2-12中任一项所述的方法,其中,第一神经网络通过至少一个深度残差网络耦合生成。
A14、如A13所述的方法,其中,第一神经网络采用34层的ResNet网络结构。
A15、如A5-14中任一项所述的方法,其中,第二神经网络通过收缩路径和扩张路径耦接生成。
A16、如A15所述的方法,其中,在第二神经网络中,收缩路径包含多个第一卷积处理块,且每个第一卷积处理块包含两个3*3的卷积核、激活函数和一个最大池化单元;扩张路径包含多个第二卷积处理块和输出块,且每个第二卷积处理块包含两个3*3的卷积核、激活函数和一个上采样单元,以及所述输出块包含1*1的卷积核。
A17、如A16所述的方法,其中,收缩路径还适于分别对经各第一卷积处理块处理后所生成的特征图进行裁剪,再将裁剪后的特征图输入到对应的扩张路径的第二卷积处理块中。
此外,本领域的技术人员能够理解,尽管在此所述的一些实施例包括其它实施例中所包括的某些特征而不是其它特征,但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。例如,在下面的权利要求书中,所要求保护的实施例的任意之一都可以以任意的组合方式来使用。
这里描述的各种技术可结合硬件或软件,或者它们的组合一起实现。从而,本发明的方法和设备,或者本发明的方法和设备的某些方面或部分可采取嵌入有形媒介,例如软盘、CD-ROM、硬盘驱动器或者其它任意机器可读的存储介质中的程序代码(即指令)的形式,其中当程序被载入诸如计算机之类的机器,并被所述机器执行时,所述机器变成实践本发明的设备。
在程序代码在可编程计算机上执行的情况下,计算设备一般包括处理器、处理器可读的存储介质(包括易失性和非易失性存储器和/或存储元件),至少一个输入装置,和至少一个输出装置。其中,存储器被配置用于存储程序代码;处理器被配置用于根据该存储器中存储的所述程序代码中的指令,执行本发明所述的方法。
以示例而非限制的方式,计算机可读介质包括计算机存储介质和通信介质。计算机可读介质包括计算机存储介质和通信介质。计算机存储介质存储诸如计算机可读指令、数据结构、程序模块或其它数据等信息。通信介质一般以诸如载波或其它传输机制等已调制数据信号来体现计算机可读指令、数据结构、程序模块或其它数据,并且包括任何信息传递介质。以上的任一种的组合也包括在计算机可读介质的范围之内。
此外,所述实施例中的一些在此被描述成可以由计算机系统的处理器或者由执行所述功能的其它装置实施的方法或方法元素的组合。因此,具有用于实施所述方法或方法元素的必要指令的处理器形成用于实施该方法或方法元素的装置。此外,装置实施例的在此所述的元素是如下装置的例子:该装置用于实施由为了实施该发明的目的的元素所执行的功能。
如在此所使用的那样,除非另行规定,使用序数词“第一”、“第二”、“第三”等等来描述普通对象仅仅表示涉及类似对象的不同实例,并且并不意图暗示这样被描述的对象必须具有时间上、空间上、排序方面或者以任意其它方式的给定顺序。
尽管根据有限数量的实施例描述了本发明,但是受益于上面的描述,本技术领域内的技术人员明白,在由此描述的本发明的范围内,可以设想其它实施例。此外,应当注意,本说明书中使用的语言主要是为了可读性和教导的目的而选择的,而不是为了解释或者限定本发明的主题而选择的。因此,在不偏离所附权利要求书的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。对于本发明的范围,对本发明所做的公开是说明性的,而非限制性的,本发明的范围由所附权利要求书限定。

Claims (10)

1.一种确定脊柱横断面图像的异常类型的方法,所述方法适于在计算设备中执行,所述方法包括步骤:
从脊柱横断面图像中获取异常区域;
确定所述脊柱横断面图像中的至少一个关键点;以及
结合所述异常区域和关键点以确定所述脊柱横断面图像的异常类型。
2.如权利要求1所述的方法,其中,所述结合异常区域和关键点以确定脊柱横断面图像的异常类型的步骤还包括:
根据所确定的至少一个关键点生成网格;以及
根据所述异常区域在所述网格中的位置,确定所述脊柱横断面图像的异常类型。
3.如权利要求1或2所述的方法,其中,所述从脊柱横断面图像中获取异常区域的步骤还包括:
利用第一神经网络来处理所述脊柱横断面图像,以判断所述脊柱横断面图像是否异常;以及
若确定所述脊柱横断面图像异常,则从所述脊柱横断面图像中获取异常区域。
4.如权利要求3所述的方法,其中,在所述利用第一神经网络来处理所述脊柱横断面图像的步骤之前,还包括步骤:
从脊柱矢状面图像中确定存在异常的感兴趣区域;以及
获取与所述存在异常的感兴趣区域相对应的脊柱横断面图像作为所述脊柱横断面图像。
5.如权利要求3或4所述的方法,其中,所述若确定脊柱横断面图像异常,则从所述脊柱横断面图像中获取异常区域的步骤包括:
利用第二神经网络来处理所述脊柱横断面图像,以输出第一图像,其中在所述第一图像中具有特定像素值的区域为异常区域;以及
对所述第一图像进行插值处理,以获得第二图像,所述第二图像与所述脊柱横断面图像具有相同尺寸。
6.如权利要求1-5中任一项所述的方法,其中,所述确定脊柱横断面图像中的至少一个关键点的步骤包括:
利用第三神经网络来处理所述脊柱横断面图像,以输出包含至少一个关键点位置的至少一个热图;以及
将所述至少一个热图合并生成预测图,其中所述预测图包含所述至少一个关键点。
7.如权利要求6所述的方法,其中,所述第三神经网络通过卷积处理层和空间融合层耦接而成。
8.如权利要求7所述的方法,其中,
所述卷积处理层包含8层卷积,且各层卷积中采用的卷积核大小依次为:5*5、5*5、5*5、5*5、9*9、1*1、1*1、1*1;以及
所述空间融合层包含5层卷积,且各层卷积中采用的卷积核大小依次为:7*7、13*13、13*13、1*1、1*1。
9.一种计算设备,包括:
一个或多个处理器;和
存储器;
一个或多个程序,其中所述一个或多个程序存储在所述存储器中并被配置为由所述一个或多个处理器执行,所述一个或多个程序包括用于执行根据权利要求1-8所述方法中的任一方法的指令。
10.一种存储一个或多个程序的计算机可读存储介质,所述一个或多个程序包括指令,所述指令当计算设备执行时,使得所述计算设备执行根据权利要求1-8所述的方法中的任一方法。
CN201910069196.2A 2019-01-24 2019-01-24 确定脊柱横断面图像的异常类型的方法及计算设备 Active CN109903269B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910069196.2A CN109903269B (zh) 2019-01-24 2019-01-24 确定脊柱横断面图像的异常类型的方法及计算设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910069196.2A CN109903269B (zh) 2019-01-24 2019-01-24 确定脊柱横断面图像的异常类型的方法及计算设备

Publications (2)

Publication Number Publication Date
CN109903269A true CN109903269A (zh) 2019-06-18
CN109903269B CN109903269B (zh) 2021-05-04

Family

ID=66944167

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910069196.2A Active CN109903269B (zh) 2019-01-24 2019-01-24 确定脊柱横断面图像的异常类型的方法及计算设备

Country Status (1)

Country Link
CN (1) CN109903269B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112071422A (zh) * 2020-09-04 2020-12-11 平安科技(深圳)有限公司 基于神经网络的腰椎病变诊断系统
CN112754458A (zh) * 2019-11-01 2021-05-07 上海联影医疗科技股份有限公司 磁共振成像方法、系统及存储介质
CN113516614A (zh) * 2020-07-06 2021-10-19 阿里巴巴集团控股有限公司 脊柱影像的处理方法、模型训练方法、装置及存储介质
CN113962927A (zh) * 2021-09-01 2022-01-21 北京长木谷医疗科技有限公司 基于强化学习的髋臼杯位置调整方法、装置及存储介质
CN115861346A (zh) * 2023-02-16 2023-03-28 邦世科技(南京)有限公司 一种基于场景感知融合网络的脊柱核磁共振图像分割方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1868162A2 (en) * 2006-06-16 2007-12-19 National Institute of Advanced Industrial Science and Technology Unusual action detection based on image analysis involving CHLAC feature extraction
CN101933045A (zh) * 2008-01-31 2010-12-29 皇家飞利浦电子股份有限公司 短轴后期增强心脏mri的自动3d分割
US20110144482A1 (en) * 2009-12-11 2011-06-16 Fujifilm Corporation Image display device and method, as well as program
CN102129699A (zh) * 2010-01-12 2011-07-20 富士胶片株式会社 断层图像生成设备、断层图像生成方法和生成断层的程序
CN102737250A (zh) * 2011-01-14 2012-10-17 西门子公司 3d医学图像数据中对脊椎骨损伤自动检测的方法和系统
US8447074B2 (en) * 2009-08-17 2013-05-21 Sony Corporation Image processing apparatus, image processing method, and program
EP2879094A2 (en) * 2013-11-28 2015-06-03 Konica Minolta, Inc. Medical image processing system and program
CN106923779A (zh) * 2010-04-28 2017-07-07 基文影像公司 用于显示体内图像部分的系统和方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1868162A2 (en) * 2006-06-16 2007-12-19 National Institute of Advanced Industrial Science and Technology Unusual action detection based on image analysis involving CHLAC feature extraction
CN101933045A (zh) * 2008-01-31 2010-12-29 皇家飞利浦电子股份有限公司 短轴后期增强心脏mri的自动3d分割
US8447074B2 (en) * 2009-08-17 2013-05-21 Sony Corporation Image processing apparatus, image processing method, and program
US20110144482A1 (en) * 2009-12-11 2011-06-16 Fujifilm Corporation Image display device and method, as well as program
CN102129699A (zh) * 2010-01-12 2011-07-20 富士胶片株式会社 断层图像生成设备、断层图像生成方法和生成断层的程序
CN106923779A (zh) * 2010-04-28 2017-07-07 基文影像公司 用于显示体内图像部分的系统和方法
CN102737250A (zh) * 2011-01-14 2012-10-17 西门子公司 3d医学图像数据中对脊椎骨损伤自动检测的方法和系统
EP2879094A2 (en) * 2013-11-28 2015-06-03 Konica Minolta, Inc. Medical image processing system and program

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
MATTOX, T.,ET.AL: "Abnormal spinal curvature and its relationship to pelvic organ prolapse", 《AMERICAN JOURNAL OF OBSTETRICS AND GYNECOLOGY 》 *
王博亮,等: "数字医学在精准骨科手术中的应用", 《 厦门大学学报(自然科学版) 》 *
王立军: "基于关联规则的医学图像分类研究", 《中国优秀硕士学位论文全文数据库》 *
胡如雪: "三维医学图像快速交互分割方法研究", 《中国优秀硕士学位论文全文数据库》 *
郭艳蓉: "医学图像分割与配准若干关键问题研究", 《中国优秀硕士学位论文全文数据库》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112754458A (zh) * 2019-11-01 2021-05-07 上海联影医疗科技股份有限公司 磁共振成像方法、系统及存储介质
CN113516614A (zh) * 2020-07-06 2021-10-19 阿里巴巴集团控股有限公司 脊柱影像的处理方法、模型训练方法、装置及存储介质
CN112071422A (zh) * 2020-09-04 2020-12-11 平安科技(深圳)有限公司 基于神经网络的腰椎病变诊断系统
WO2021121159A1 (zh) * 2020-09-04 2021-06-24 平安科技(深圳)有限公司 基于神经网络的腰椎病变诊断结果输出系统以及方法
CN112071422B (zh) * 2020-09-04 2024-05-03 平安科技(深圳)有限公司 基于神经网络的腰椎病变诊断系统
CN113962927A (zh) * 2021-09-01 2022-01-21 北京长木谷医疗科技有限公司 基于强化学习的髋臼杯位置调整方法、装置及存储介质
CN115861346A (zh) * 2023-02-16 2023-03-28 邦世科技(南京)有限公司 一种基于场景感知融合网络的脊柱核磁共振图像分割方法

Also Published As

Publication number Publication date
CN109903269B (zh) 2021-05-04

Similar Documents

Publication Publication Date Title
CN109903269A (zh) 确定脊柱横断面图像的异常类型的方法及计算设备
Gibson et al. NiftyNet: a deep-learning platform for medical imaging
US20210225027A1 (en) Image region localization method, image region localization apparatus, and medical image processing device
CN110110617B (zh) 医学影像分割方法、装置、电子设备和存储介质
CN109872306A (zh) 医学图像分割方法、装置和存储介质
CN109903268A (zh) 确定脊柱图像集的异常类型的方法及计算设备
CN109993726A (zh) 医学图像的检测方法、装置、设备和存储介质
CN109493347A (zh) 在图像中对稀疏分布的对象进行分割的方法和系统
CN108603922A (zh) 自动心脏体积分割
CN109559300A (zh) 图像处理方法、电子设备及计算机可读存储介质
EP3921799A1 (en) Method and apparatus for automated target and tissue segmentation using multi-modal imaging and ensemble machine learning models
CN108596904A (zh) 生成定位模型的方法及对脊柱矢状位图像进行处理的方法
US10275909B2 (en) Systems and methods for an integrated system for visualizing, simulating, modifying and 3D printing 3D objects
CN111145147B (zh) 多模态医学图像的分割方法及终端设备
JP7239701B2 (ja) 人工神経網を利用する臓器の体積測定方法及びその装置
CN109003270A (zh) 一种图像处理方法及电子设备
CN114298234B (zh) 脑部医学影像分类方法、装置、计算机设备和存储介质
CN108846829A (zh) 病变部位识别方法及装置、计算机装置及可读存储介质
CN106062782A (zh) 针对基于图集的配准的无监督的训练
Ruppert et al. Medical image registration based on watershed transform from greyscale marker and multi-scale parameter search
CN114998247A (zh) 异常预测方法、装置、计算机设备、存储介质
CN114998374A (zh) 基于位置先验的图像分割方法、装置、设备和存储介质
Li et al. S 3 egANet: 3D spinal structures segmentation via adversarial nets
CN108634934A (zh) 对脊柱矢状位图像进行处理的方法和设备
Aldemir et al. Chain code strategy for lossless storage and transfer of segmented binary medical data

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Liu Xingyu

Inventor after: Zhang Yiling

Inventor after: An Yicheng

Inventor before: Liu Xingyu

Inventor before: Zhang Yiling

Inventor before: Feng Shiqing

Inventor before: Zhang Di

Inventor before: Liu Huiye

Inventor before: Zhang Yundong

Inventor before: An Yicheng

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant
CP03 Change of name, title or address

Address after: Room 1904, unit 2, building 9, yard 2, Simiao Road, Daxing District, Beijing 100176

Patentee after: Liu Xingyu

Patentee after: Beijing Changmugu Medical Technology Co.,Ltd.

Address before: Room 602-5, 6th Floor, No. 52 North Fourth Ring West Road, Haidian District, Beijing, 100190

Patentee before: Liu Xingyu

Patentee before: BEIJING CHANGMUGU MEDICAL TECHNOLOGY Co.,Ltd.

CP03 Change of name, title or address