CN108629773B - 建立训练识别心脏血管类型的卷积神经网络数据集的方法 - Google Patents

建立训练识别心脏血管类型的卷积神经网络数据集的方法 Download PDF

Info

Publication number
CN108629773B
CN108629773B CN201810441538.4A CN201810441538A CN108629773B CN 108629773 B CN108629773 B CN 108629773B CN 201810441538 A CN201810441538 A CN 201810441538A CN 108629773 B CN108629773 B CN 108629773B
Authority
CN
China
Prior art keywords
image
data
rough
graph
storing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810441538.4A
Other languages
English (en)
Other versions
CN108629773A (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.)
BEIJING HONGYUN ZHISHENG TECHNOLOGY Co.,Ltd.
Fuwai Hospital of CAMS and PUMC
Original Assignee
Beijing Hongyun Zhisheng Technology Co ltd
Fuwai Hospital of CAMS and PUMC
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 Beijing Hongyun Zhisheng Technology Co ltd, Fuwai Hospital of CAMS and PUMC filed Critical Beijing Hongyun Zhisheng Technology Co ltd
Priority to CN201810441538.4A priority Critical patent/CN108629773B/zh
Publication of CN108629773A publication Critical patent/CN108629773A/zh
Application granted granted Critical
Publication of CN108629773B publication Critical patent/CN108629773B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/241Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • G06F18/2413Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on distances to training or reference patterns
    • 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/10116X-ray image
    • G06T2207/10121Fluoroscopy
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20081Training; Learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20084Artificial neural networks [ANN]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular
    • G06T2207/30104Vascular flow; Blood flow; Perfusion

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Evolutionary Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

本发明公开一种建立训练识别心脏血管类型的卷积神经网络数据集的方法,包括:获取脱敏处理后的心脏冠脉数据,将心脏冠脉数据以单张静态图片形式存储;在单张静态图片中使用对应颜色标注当前体位下的血管,形成粗标图;对粗标图中标注的信息进行像素级标注,形成精标图;识别精标图中血管的颜色,使图片由三通道图像变成单通道图像;将单通道图像存储为二值图;基于脱敏处理后的心脏冠脉数据、粗标图、精标图、单通道图像和二值图建立训练识别心脏血管类型的卷积神经网络数据集。本发明对数据作统一标准的处理,该数据集中的数据真实、多样,数据量大,数据类型丰富,格式标准、错误率较小,可以用于训练不同功能的神经网络,减少人工干预。

Description

建立训练识别心脏血管类型的卷积神经网络数据集的方法
技术领域
本发明涉及的是一种建立训练识别心脏血管类型的卷积神经网络数据集的方法及该数据集,属于通信技术领域。
背景技术
人工神经网络是人工智能领域提出的一种运算模型,它从信息处理角度对人脑神经元网络进行抽象,建立某种简单模型,按不同的连接方式组成不同的网络,以此来解决人工智能领域(如图像识别)中的某些问题。卷积神经网络(Convolutional Neural Network,CNN)是人工神经网络的一种,是一种前馈神经网络,人工神经元可以响应周围单元,可以进行大型图像处理,即输入一张图片,卷积神经网络对此图片进行一定的变换,输出得到的结果。
网络的输出若尽可能的接近预测的值,就需要训练网络。训练网络就是卷积神经网络进行自身权重迭代调整的过程。我们给出一个输入和一个期望得到的输出,网络就会比较目前的输出值和期望值,再根据两者的差异情况来更新每一层的参数权重。训练和测试网络需要用到一个拥有大量数据的数据集。数据集应该主要包括原始图片和每张原始图片对应的标签。
现有技术中,大多数数据集均为用于通用图像分割训练的数据集,这些数据一般来自于日常生活的照片,具有比较高的可识别度。而需要被分割和识别血管的心脏冠状动脉的图像和日常生活的照片具有很大差异,是用数字减影血管造影技术得到的。数字减影血管造影的基本原理是将注入造影剂前后拍摄的两帧X线图像经数字化输入图像计算机,通过减影、增强和再成像过程来获得纯血管影像,同时实时地显现血管影。这种技术得到的图像分辨率较低、噪声很高。同时,心脏冠状动脉血管的类间差异较小、类内差异很大。这两者决定了造影图像和通用数据集中的图像有非常大的区别。用通用图像数据集训练得到的分割网络不能用于心脏冠脉血管的识别和分割。
发明内容
针对上述缺陷,本发明提供了一种建立训练识别心脏血管类型的卷积神经网络数据集的方法,该方法对数据作统一标准的处理,得到一个格式标准、错误率较小的数据集,该数据集中的数据真实、多样,数据量大,数据类型丰富,可以用于训练不同功能的神经网络,做到了尽量减少人工干预。
为达到上述目的,本发明通过以下技术方案来具体实现:
本发明提供了一种建立训练识别心脏血管类型的卷积神经网络数据集的方法,该方法包括:
获取脱敏处理后的心脏冠脉数据,将所述心脏冠脉数据以单张静态图片形式存储,每张所述单张静态图片绑定唯一序号;
按照专业医学的命名规则,在所述单张静态图片中使用对应颜色标注当前体位下的血管,存储形成粗标图;并将不同体位的粗标图分类存储;
对粗标图中标注的信息进行像素级标注,形成精标图;
识别精标图中血管的颜色,将每种颜色存储为对应的唯一标号,使图片由三通道图像变成单通道图像;将单通道图像二值化处理存储为二值图;基于脱敏处理后的心脏冠脉数据、粗标图、精标图、单通道图像和二值图建立训练识别心脏血管类型的卷积神经网络数据集。
进一步的,所述获取脱敏处理后的心脏冠脉数据的步骤包括:
获取包括健康的心脏冠脉数据和病变的心脏冠脉数据;
删除无关数据后形成只含有数字减影血管造影技术(Digital SubtractionAngiography,DSA)影像和拍摄体位信息的dicom文件的心脏冠脉数据。
进一步的,所述将所述心脏冠脉数据以单张静态图片形式存储,每张所述单张静态图片绑定唯一序号的具体步骤包括:
根据心脏冠脉数据中DSA影像的长短,提取并存储影像中预设数量的单张静态图片;
F(N)=f(n)×g(n)×D(n)
Figure BDA0001655962010000021
Figure BDA0001655962010000022
Figure BDA0001655962010000023
其中,F(N)是一组单张静态图片的集合,N表示每张单张静态图片绑定唯一序号;f(n)是脱敏处理后的心脏冠脉数据的原始的DSA影像;g(n)是一个门函数,其中
Figure BDA0001655962010000024
]是向下取整符号,λ是一个可变参数,0≤λ≤0.5,l表示原始的DSA影像的图片帧数,n表示原始的DSA影像的图片帧序号,0≤n≤l;D(n)是冲激函数δ(n)组成的抽样函数,m是可变参数,控制静态图片提取的数量。
进一步的,所述体位包括但不限于头位(CRA),足位(CAU),左前斜(LAO),左前斜足位(LAO_CAU),左前斜头位(LAO_CRA),右前斜足位(RAO_CAU),右前斜头位(RAO_CRA)。
进一步的,所述对粗标图中标注的信息进行像素级标注,形成精标图,包括:
将选取的粗标图中标注的血管,用对应的颜色完全填充,同时单独标出导丝;对于粗标图中没有特别标注的血管,以及导丝附近与背景明显不同的阴影,统一标注为不关注血管或漏液,存储后形成精标图。
进一步的,所述对粗标图中标注的信息进行像素级标注,形成精标图,包括:
按不同体位等比例随机挑选粗标图或顺序选择粗标图。
进一步的,按不同体位等比例随机挑选粗标图的方法,包括:
在脱敏处理后的心脏冠脉数据中按不同体位等比例随机挑选静态图片,根据静态图片绑定唯一序号在粗标图中选取对应图片。
进一步的,所述将单通道图像二值化处理存储为二值图,包括:
将单通道图像中的血管和导丝显示白色,背景显示黑色,包括:
Figure BDA0001655962010000031
其中,G[i,j]为二值化图片中第i行第j列上像素的值,S[i,j]为单通道图片中第i行第j列上像素值;
二值化处理后的图片进行存储后形成二值图。
本发明的有益效果是:
通过本发明提供的技术方案,所处理的数据类型更为复杂,得到的数据也更为多样化。包括dicom(Digital Imaging and Communications in Medicine,即医学数字成像和通信,它定义了质量能满足临床需要的可用于数据交换的医学图像格式)格式的连续图片帧、数字减影血管造影下多体位心脏冠状动脉的图像、血管的黑白二值图标签、血管位置和名称的粗标注、血管位置和名称的像素级标注等,可以对数据作统一标准的处理,得到一个格式标准、错误率较小的数据集,该数据集中的数据真实、多样,数据量大,数据类型丰富,可以用于训练不同功能的神经网络,做到了尽量减少人工干预。
附图说明
图1所示为本发明提供的一种建立训练识别心脏血管类型的卷积神经网络数据集的方法的实施例一流程图。
图2所示为本发明提供的一种建立训练识别心脏血管类型的卷积神经网络数据集的方法的示意图之一。
图3所示为本发明提供的一种建立训练识别心脏血管类型的卷积神经网络数据集的方法的示意图之二。
图4所示为本发明提供的一种建立训练识别心脏血管类型的卷积神经网络数据集的方法的示意图之三。
具体实施方式
下面对本发明的技术方案进行具体阐述,需要指出的是,本发明的技术方案不限于实施例所述的实施方式,本领域的技术人员参考和借鉴本发明技术方案的内容,在本发明的基础上进行的改进和设计,应属于本发明的保护范围。
实施例一
本发明实施例一提供了一种建立训练识别心脏血管类型的卷积神经网络数据集的方法,该方法包括步骤S110-步骤S140:
在步骤S110中,获取脱敏处理后的心脏冠脉数据,将所述心脏冠脉数据以单张静态图片形式存储,每张所述单张静态图片绑定唯一序号。
其中,所述获取脱敏处理后的心脏冠脉数据的步骤包括:
获取包括健康的心脏冠脉数据和病变的心脏冠脉数据。心脏冠脉数据的原始数据均为与医院联系、经过病人同意之后取得的真实数据,包括健康的心脏冠脉数据和病变的心脏冠脉数据。提取的数据均以dicom格式的文件保存。dicom(Digital Imaging andCommunications in Medicine)即医学数字成像和通信,是医学图像和相关信息的国际标准。它定义了质量能满足临床需要的可用于数据交换的医学图像格式。本发明中的dicom文件中除了有DSA影像,还包含病人的姓名、出生年月、性别、拍摄时间、医院信息、拍摄体位(角度)的信息。
删除无关数据后形成只含有DSA影像和拍摄体位信息的dicom文件的心脏冠脉数据。
为了保护病人的隐私,同时减小数据存储量,将与标注无关的信息剔除。使用的脱敏系统具有dicom文件接口,dicom文件读入内存后是一个结构类型的数据,系统自动将脱敏系统自动将姓名、出生年月、性别、拍摄时间和医院信息的词条删除,形成一个新的只含有DSA影像和拍摄体位信息的dicom文件。
进一步的,所述将所述心脏冠脉数据以单张静态图片形式存储,每张所述单张静态图片绑定唯一序号的具体步骤包括:
dicom文件中存储的DSA影像是一组连续图片帧,将其处理成单张静态图片的形式。这一步可以人工提取也可以自动提取,提取后图片保存的统一命名格式为:【序号_体位_src.png】。其中序号是在数据库中唯一的标号,根据图片提取的先后顺序从小到大排列;
自动提取的方法根据心脏冠脉数据中DSA影像的长短,提取并存储影像中预设数量的单张静态图片;如选取视频的中间段的一定数量的图像,原理:
F(N)=f(n)×g(n)×D(n)
Figure BDA0001655962010000051
Figure BDA0001655962010000052
Figure BDA0001655962010000053
其中,F(N)是一组单张静态图片的集合,N表示每张单张静态图片绑定唯一序号;f(n)是脱敏处理后的心脏冠脉数据的原始的DSA影像;g(n)是一个门函数,其中
Figure BDA0001655962010000054
]是向下取整符号,λ是一个可变参数,0≤λ≤0.5,l表示原始的DSA影像的图片帧数,n表示原始的DSA影像的图片帧序号,0≤n≤l;D(n)是冲激函数δ(n)组成的抽样函数,m是可变参数,控制静态图片提取的数量。
获得的图片如图2所示。
在步骤S120中,按照专业医学的命名规则,在所述单张静态图片中使用对应颜色标注当前体位下的血管,存储形成粗标图;并将不同体位的粗标图分类存储。
进一步的,所述体位包括但不限于头位(CRA),足位(CAU),左前斜(LAO),左前斜足位(LAO_CAU),左前斜头位(LAO_CRA),右前斜足位(RAO_CAU),右前斜头位(RAO_CRA)。
按照专业医学的命名规则,在原始静态帧图片上进行简单的描画,标出当前体位下能观察到的主要血管及其走向。其中每个颜色和血管都是一一对应的,即一个颜色在任何体位下都只代表一种血管,每种血管只能由一种颜色表示,如表1所示。
表1
血管 颜色(RGB值)
左主干 中绿(0,159,48)
前降支(近) 土黄(177,135,60)
第一对角支 深桃红(177,63,96)
第一对角支(Ad) 浅桃红(255,39,108)
前降支(中) 橙红(255,111,72)
第二对角支 藏蓝(0,15,120)
第二对角支(Ad) 若草(177,246,132)
前降支(远) 墨绿(0,87,84)
中间支 粉绿(0,198,156)
回旋支(近) 米黄(255,222,144)
回旋支(远) 普鲁士蓝(0,126,192)
第二钝缘支 粉(255,150,180)
第一钝缘支 灰(177,174,168)
左室后侧支 紫(177,102,204)
后降支(回) 品红(255,78,216)
右主干(近) 红(255,0,0)
右主干(中) 浅绿(0,231,12)
右主干(远) 黄绿(177,207,24)
后降支(右) 橙黄(255,183,36)
后侧支 群青蓝(0,54,228)
标注后图片保存的统一命名格式为【序号_体位_mask.png】。其中序号和体位与原图一致。标注结果如图3所示(实际为没有提示框的彩色图)。
根据每张静态图片所携带的体位信息(以文件名的形式),将不同体位的数据分离开来:
使用机器识别每张静态图片的文件名中所包含的体位信息,将数据分为头位(CRA),足位(CAU),左前斜+足位(LAO_CAU),左前斜+头位(LAO_CRA),右前斜+足位(RAO_CAU),右前斜+头位(RAO_CRA),不同体位的图片分别存放在不同的文件夹内。
左前斜和右前斜的图片还分为左冠和右冠。使用机器检测图片内的像素颜色信息,如果图片中包含(255,0,0),(0,231,12),(177,207,24)中任一一个颜色的像素点,就归类为右冠,将文件名改为【序号_体位_R_src.png】和【序号_体位_R_mask.png】,剩下的图片归类为左冠,文件名改为【序号_体位_L_mask.png】和【序号_体位_L_src.png】。
在步骤S130中,对粗标图中标注的信息进行像素级标注,形成精标图。
将选取的粗标图中标注的血管,用对应的颜色完全填充,同时单独标出导丝;对于粗标图中没有特别标注的血管,以及导丝附近与背景明显不同的阴影,统一标注为不关注血管或漏液,存储后形成精标图。
按不同体位等比例随机挑选粗标图或顺序选择粗标图。
其中,按不同体位等比例随机挑选粗标图的方法,包括:
在脱敏处理后的心脏冠脉数据(文件名中带有src的数据)中按不同体位等比例随机挑选静态图片,根据静态图片绑定唯一序号在粗标图中选取对应图片(文件名中带有mask的数据)。
由于相邻编号的图片可能比较相似,顺序选择粗标图工作量大,极大的浪费了时间和精力,可选择按不同体位等比例随机挑选粗标图的方法,保证血管多样性的同时,减少了挑选时间。
根据粗标图,对粗标图按照一定的规则进行像素级标注。即将粗标图中描画过的血管,用对应的颜色完全填充血管。同时单独标出导丝。对于粗标注中没有特别标注的血管,以及导丝附近与背景明显不同的阴影,统一标注为“不关注血管/漏液”。图片保存为【序号_体位_(R/L/无)_label.png】。血管和对应颜色如表2所示:
表2
血管 医生粗标颜色(RGB值) 人工精标颜色(RGB值)
左主干 中绿(0,159,48) (255,0,0)
前降支(近) 土黄(177,135,60) (0,255,0)
第一对角支 深桃红(177,63,96) (0,0,255)
第一对角支(Ad) 浅桃红(255,39,108) (0,128,255)
前降支(中) 橙红(255,111,72) (255,0,255)
第二对角支 藏蓝(0,15,120) (0,255,255)
第二对角支(Ad) 若草(177,246,132) (128,0,255)
前降支(远) 墨绿(0,87,84) (255,255,0)
中间支 粉绿(0,198,156) (128,0,0)
回旋支(近) 米黄(255,222,144) (0,0,128)
回旋支(远) 普鲁士蓝(0,126,192) (0,128,0)
第二钝缘支 粉(255,150,180) (128,128,0)
第一钝缘支 灰(177,174,168) (128,0,128)
左室后侧支 紫(177,102,204) (0,128,128)
后降支(回) 品红(255,78,216) (128,255,0)
右主干(近) 红(255,0,0) (255,128,0)
右主干(中) 浅绿(0,231,12) (255,0,128)
右主干(远) 黄绿(177,207,24) (0,255,128)
后降支(右) 橙黄(255,183,36) (128,128,255)
后侧支 群青蓝(0,54,228) (128,255,128)
不关注血管 (150,150,150)
背景 (0,0,0)
导管 (255,255,255)
像素级标注后的精标图,如图4所示。其中,1为右主干(近),RGB值为(255,0,0)
2为右主干(中),RGB值为(255,0,128)
3为右主干(远),RGB值为(0,255,128)
4为后降支(右),RGB值为(128,128,255)
5为后侧支,RGB值为(128,255,128)
6为导丝,RGB值为(255,255,255)
椭圆圈出的部分为不关注血管/漏液,RGB值为(150,150,150)
所有黑色部分为背景,RGB值为(0,0,0)。
在步骤S140中,识别精标图中血管的颜色,将每种颜色存储为对应的唯一标号,使图片由三通道图像变成单通道图像;将单通道图像二值化处理存储为二值图,基于脱敏处理后的心脏冠脉数据、粗标图、精标图、单通道图像和二值图建立训练识别心脏血管类型的卷积神经网络数据集。
步骤S140之前的标注数据(mask、label)都是有颜色的,因此都是三通道图像,在这一步中使用机器识别每个像素点的RGB颜色信息,将每种颜色对应成一个数字标号存储下来,使图片变成一张单通道图像。文件名为【序号_体位_(R/L/无)_single.png】。如表3所示:
表3
血管 医生粗标颜色(RGB值) 人工精标颜色(RGB值) 单通道标号
左主干 中绿(0,159,48) (255,0,0) 1
前降支(近) 土黄(177,135,60) (0,255,0) 2
第一对角支 深桃红(177,63,96) (0,0,255) 3
第一对角支(Ad) 浅桃红(255,39,108) (0,128,255) 4
前降支(中) 橙红(255,111,72) (255,0,255) 5
第二对角支 藏蓝(0,15,120) (0,255,255) 6
第二对角支(Ad) 若草(177,246,132) (128,0,255) 7
前降支(远) 墨绿(0,87,84) (255,255,0) 8
中间支 粉绿(0,198,156) (128,0,0) 9
回旋支(近) 米黄(255,222,144) (0,0,128) 10
回旋支(远) 普鲁士蓝(0,126,192) (0,128,0) 11
第二钝缘支 粉(255,150,180) (128,128,0) 12
第一钝缘支 灰(177,174,168) (128,0,128) 13
左室后侧支 紫(177,102,204) (0,128,128) 14
后降支(回) 品红(255,78,216) (128,255,0) 15
右主干(近) 红(255,0,0) (255,128,0) 16
右主干(中) 浅绿(0,231,12) (255,0,128) 17
右主干(远) 黄绿(177,207,24) (0,255,128) 18
后降支(右) 橙黄(255,183,36) (128,128,255) 19
后侧支 群青蓝(0,54,228) (128,255,128) 20
不关注血管 (150,150,150) 21
背景 (0,0,0) 0
导管 (255,255,255) 255
将单通道图像二值化处理存储为二值图包括:
将单通道图像中的血管和导丝显示白色,背景显示黑色,包括:
Figure BDA0001655962010000091
其中,G[i,j]为二值化图片中第i行第j列上像素的值,S[i,j]为单通道图片中第i行第j列上像素值;
二值化处理后的图片进行存储后形成二值图,存储的文件名为【序号_体位_(R/L/无)_grey.png】,包括使用数字减影血管造影技术得到的真实人体在不同体位下心脏冠脉的原始视频、原始图片、粗标注图片、精标注图片以及二值化图片等等,可以用来作为后续病变检测的训练数据。
本发明的有益效果是:
本发明提供的技术方案,所处理的数据类型更为复杂,得到的数据也更为多样化。包括dicom(Digital Imaging and Communications in Medicine,即医学数字成像和通信,它定义了质量能满足临床需要的可用于数据交换的医学图像格式)格式的连续图片帧、数字减影血管造影下多体位心脏冠状动脉的图像、血管的黑白二值图标签、血管位置和名称的粗标注、血管位置和名称的像素级标注等,可以对数据作统一标准的处理,得到一个格式标准、错误率较小的数据集,该数据集中的数据真实、多样,数据量大,数据类型丰富,可以用于训练不同功能的神经网络,做到了尽量减少人工干预。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到本发明可以通过硬件实现,也可以借助软件加必要的通用硬件平台的方式来实现。基于这样的理解,本发明的技术方案可以以软件产品的形式体现出来,该软件产品可以存储在一个非易失性存储介质(可以是CD-ROM,U盘,移动硬盘等)中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述的方法。
以上公开的仅为本发明的几个具体实施例,但是,本发明并非局限于上述实施例,任何本领域的技术人员能思之的变化都应落入本发明的保护范围。

Claims (8)

1.一种建立训练识别心脏血管类型的卷积神经网络数据集的方法,其特征在于,该方法包括:
获取脱敏处理后的心脏冠脉数据,将所述心脏冠脉数据以单张静态图片形式存储,每张所述单张静态图片绑定唯一序号;
按照专业医学的命名规则,在所述单张静态图片中使用对应颜色标注当前体位下的血管,存储形成粗标图;并将不同体位的粗标图分类存储;
对粗标图中标注的信息进行像素级标注,形成精标图;
识别精标图中血管的颜色,将每种颜色存储为对应的唯一标号,使图片由三通道图像变成单通道图像;将单通道图像二值化处理存储为二值图;
基于脱敏处理后的心脏冠脉数据、粗标图、精标图、单通道图像和二值图建立训练识别心脏血管类型的卷积神经网络数据集。
2.如权利要求1所述的方法,其特征在于,所述获取脱敏处理后的心脏冠脉数据的步骤包括:
获取包括健康的心脏冠脉数据和病变的心脏冠脉数据;
删除无关数据后形成只含有DSA影像和拍摄体位信息的dicom文件的心脏冠脉数据。
3.如权利要求1或2所述的方法,其特征在于,所述将所述心脏冠脉数据以单张静态图片形式存储,每张所述单张静态图片绑定唯一序号的具体步骤包括:
根据心脏冠脉数据中DSA影像的长短,提取并存储影像中预设数量的单张静态图片;
F(N)=f(n)×g(n)×D(n)
Figure FDA0002860925180000011
Figure FDA0002860925180000012
Figure FDA0002860925180000013
其中,F(N)是一组单张静态图片的集合,N表示每张单张静态图片绑定唯一序号;f(n)是脱敏处理后的心脏冠脉数据的原始的DSA影像;g(n)是一个门函数,其中
Figure FDA0002860925180000014
是向下取整符号,λ是一个可变参数,0≤λ≤0.5,l表示原始的DSA影像的图片帧数,n表示原始的DSA影像的图片帧序号,0≤n≤l;D(n)是冲激函数δ(n)组成的抽样函数,m是可变参数,控制静态图片提取的数量。
4.如权利要求1所述的方法,其特征在于,所述体位包括但不限于头位(CRA),足位(CAU),左前斜(LAO),左前斜足位(LAO_CAU),左前斜头位(LAO_CRA),右前斜足位(RAO_CAU),右前斜头位(RAO_CRA)。
5.如权利要求1所述的方法,其特征在于,所述对粗标图中标注的信息进行像素级标注,形成精标图,包括:
将选取的粗标图中标注的血管,用对应的颜色完全填充,同时单独标出导丝;对于粗标图中没有特别标注的血管,以及导丝附近与背景明显不同的阴影,统一标注为不关注血管或漏液,存储后形成精标图。
6.如权利要求1或5所述的方法,其特征在于,所述对粗标图中标注的信息进行像素级标注,形成精标图,包括:
按不同体位等比例随机挑选粗标图或顺序选择粗标图。
7.如权利要求6所述的方法,其特征在于,按不同体位等比例随机挑选粗标图的方法,包括:
在脱敏处理后的心脏冠脉数据中按不同体位等比例随机挑选静态图片,根据静态图片绑定唯一序号在粗标图中选取对应图片。
8.如权利要求1所述的方法,其特征在于,将单通道图像二值化处理进行存储,包括:
将单通道图像中的血管和导丝显示白色,背景显示黑色,包括:
Figure FDA0002860925180000021
其中,G[i,j]为二值化图片中第i行第j列上像素的值,S[i,j]为单通道图片中第i行第j列上像素值;
二值化处理后的图片进行存储后形成二值图。
CN201810441538.4A 2018-05-10 2018-05-10 建立训练识别心脏血管类型的卷积神经网络数据集的方法 Active CN108629773B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810441538.4A CN108629773B (zh) 2018-05-10 2018-05-10 建立训练识别心脏血管类型的卷积神经网络数据集的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810441538.4A CN108629773B (zh) 2018-05-10 2018-05-10 建立训练识别心脏血管类型的卷积神经网络数据集的方法

Publications (2)

Publication Number Publication Date
CN108629773A CN108629773A (zh) 2018-10-09
CN108629773B true CN108629773B (zh) 2021-06-18

Family

ID=63692508

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810441538.4A Active CN108629773B (zh) 2018-05-10 2018-05-10 建立训练识别心脏血管类型的卷积神经网络数据集的方法

Country Status (1)

Country Link
CN (1) CN108629773B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112308813A (zh) * 2019-07-26 2021-02-02 宏碁股份有限公司 血管状态评估方法与血管状态评估装置
CN110970117B (zh) * 2019-11-29 2023-06-20 北京邮电大学 一种基于二值骨架图的冠状动脉造影数据集存储方法
CN113822839B (zh) * 2020-06-18 2024-01-23 飞依诺科技股份有限公司 医学图像的处理方法、装置、计算机设备和存储介质
CN115953495B (zh) * 2023-03-14 2023-08-01 北京唯迈医疗设备有限公司 基于二维造影图像的智能路径规划装置、系统和存储介质

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107862300A (zh) * 2017-11-29 2018-03-30 东华大学 一种基于卷积神经网络的监控场景下行人属性识别方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7927275B2 (en) * 2002-08-26 2011-04-19 The Cleveland Clinic Foundation System and method of aquiring blood-vessel data
EP2252973A2 (en) * 2008-03-06 2010-11-24 Philips Intellectual Property & Standards GmbH Method of selectively and interactively processing data sets with fusion
DE102014213408B4 (de) * 2014-07-10 2018-10-25 Siemens Healthcare Gmbh Verfahren zur Ermittlung eines dreidimensionalen Modelldatensatzes eines wenigstens ein Gefäßsegment umfassenden Blutgefäßsystems
US9767557B1 (en) * 2016-06-23 2017-09-19 Siemens Healthcare Gmbh Method and system for vascular disease detection using recurrent neural networks
CN106548213B (zh) * 2016-11-30 2019-04-23 上海联影医疗科技有限公司 血管识别方法和装置
CN107154036A (zh) * 2017-03-24 2017-09-12 中南大学 基于亚像素的血管分割方法及其分割系统
CN107993229B (zh) * 2017-12-15 2021-11-19 中科微光医疗研究中心(西安)有限公司 一种基于心血管ivoct图像的组织分类方法和装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107862300A (zh) * 2017-11-29 2018-03-30 东华大学 一种基于卷积神经网络的监控场景下行人属性识别方法

Also Published As

Publication number Publication date
CN108629773A (zh) 2018-10-09

Similar Documents

Publication Publication Date Title
CN108629773B (zh) 建立训练识别心脏血管类型的卷积神经网络数据集的方法
CN108830155B (zh) 一种基于深度学习的心脏冠状动脉分割及识别的方法
CN108021916B (zh) 基于注意力机制的深度学习糖尿病视网膜病变分类方法
Miki et al. Classification of teeth in cone-beam CT using deep convolutional neural network
CN109146872B (zh) 基于深度学习和光流法的心脏冠状动脉影像分割识别方法
CN112017185B (zh) 病灶分割方法、装置及存储介质
CN108615236A (zh) 一种图像处理方法及电子设备
CN110786824B (zh) 基于边界框修正网络的粗标注眼底照出血病变检测方法及系统
CN108921854B (zh) 一种消化道内镜图像不规则病变区域标注方法及系统
CN111062947B (zh) 一种基于深度学习的x光胸片病灶定位方法及系统
CN109063557B (zh) 快速构建心脏冠脉血管识别数据集的方法
CN112541923B (zh) 一种基于眼底图数据集迁移学习的视杯视盘分割方法
CN110110808A (zh) 一种对图像进行目标标注的方法、装置及计算机记录介质
CN111797901A (zh) 一种基于拓扑结构估计的视网膜动静脉分类方法及装置
Wan et al. Retinal image enhancement using cycle-constraint adversarial network
Reddy et al. Retinal fundus image enhancement using piecewise gamma corrected dominant orientation based histogram equalization
CN111340773B (zh) 一种视网膜图像血管分割方法
CN113889238B (zh) 一种图像识别方法、装置、电子设备及存储介质
CN109509179B (zh) 基于医学图像的眼球和晶状体的自动分割方法及系统
Dong et al. Supervised learning-based retinal vascular segmentation by m-unet full convolutional neural network
Shaharum et al. Automatic detection of diabetic retinopathy retinal images using artificial neural network
CN110910409B (zh) 一种灰度图像处理方法、装置和计算机可读存储介质
CN110428431B (zh) 一种心脏医学图像的分割方法、装置、设备及存储介质
CN116664467A (zh) 一种基于交叉神经网络和eca-s的多模态医学图像融合方法
CN111178271B (zh) 一种人脸图像特征增强方法、人脸识别方法及电子设备

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
TA01 Transfer of patent application right

Effective date of registration: 20210125

Address after: 100086 1704-1705, 17th floor, Qingyun contemporary building, building 9, Manting Fangyuan community, Qingyun Li, Haidian District, Beijing

Applicant after: BEIJING HONGYUN ZHISHENG TECHNOLOGY Co.,Ltd.

Applicant after: FUWAI HOSPITAL, CHINESE ACADEMY OF MEDICAL SCIENCES

Address before: 100086 1704-1705, 17th floor, Qingyun contemporary building, building 9, Manting Fangyuan community, Qingyun Li, Haidian District, Beijing

Applicant before: BEIJING HONGYUN ZHISHENG TECHNOLOGY Co.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant