CN106778829A - 一种主动学习的肝脏损伤类别的图像检测方法 - Google Patents
一种主动学习的肝脏损伤类别的图像检测方法 Download PDFInfo
- Publication number
- CN106778829A CN106778829A CN201611062637.9A CN201611062637A CN106778829A CN 106778829 A CN106778829 A CN 106778829A CN 201611062637 A CN201611062637 A CN 201611062637A CN 106778829 A CN106778829 A CN 106778829A
- Authority
- CN
- China
- Prior art keywords
- pixel
- gray scale
- scale difference
- liver
- density
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30056—Liver; Hepatic
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V2201/00—Indexing scheme relating to image or video recognition or understanding
- G06V2201/03—Recognition of patterns in medical or anatomical images
- G06V2201/031—Recognition of patterns in medical or anatomical images of internal organs
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种主动学习的肝脏损伤类别的图像检测方法,具体步骤包括:步骤1,使用肝脏主动形状模型从腹部医学图像集中定位并提取肝脏目标区域;步骤2,定义和提取肝脏目标区域中像素点的方向模式特征向量和灰度幅度变化模式特征向量;步骤3,基于肝脏目标区域中像素点局部密度和到更大密度像素点最小距离二维空间的肝脏像素密度聚类中心定义和聚类划分;步骤4,基于肝脏像素密度聚类边缘、聚类中心和聚类划分的主动学习方法进行肝脏损伤类别的图像检测。
Description
技术领域
本发明属于计算机图像处理领域,特别涉及一种主动学习的肝脏损伤类别的图像检测方法。
背景技术
随着交通事业的快速发展,交通事故等引起的腹部损伤呈不断上升的趋势,特别是肝脏因弹性较差,容易导致暴力性损伤。此外,临床中有许多病理性损伤,如甲肝、乙肝等病毒或一些药物和化学物质引起的肝脏损伤,或肝脏肿瘤、心功能不全导致肝脏淤血,或某些先天性肝脏疾病等引起的肝脏损伤。暴力性或病理性损伤会让肝脏受累,而且可能会累及更多器官,对人体造成更大的伤害,严重威胁到患者的生命安全和身体健康。肝脏损伤检测为后续治疗方案的制定和选择提供重要依据,在临床中具有重要的意义。
临床中医生借助医学图像肉眼观察和临床经验检测肝脏损伤。Kakihara利用增强型CT(Computerized Tomography)图像人工分析损伤的大小、形状和位置,检测由蛔虫感染引起的内脏幼虫移动症肝脏损伤。为了克服灌注CT图像的高辐射缺陷,Wang在双能量CT上分析了虚拟光谱曲线的肝脏损伤定量特征。Purysko利用增强型MRI(Magnetic ResonanceImaging)图像人工分析了良性、肿瘤和肝细胞癌的肝脏损伤特征。Matthias Eiber比较了扩散加权的MRI和多层螺旋CT两种模态图像对人工检测肝脏损伤准确率的影响。目前,医学图像人工肝脏损伤检测会面临几个方面的问题:
(1)暴力性肝脏损伤,常伴有腹腔积血、积液、胃肠内容物及气体等干扰,加之患者***不能配合,这些干扰信号增加了肝脏损伤检测的难度。而病理性肝脏损伤,临床症状及患者体征表现通常不够明显,导致损伤检测面临挑战。
(2)人工检测方式存在信息利用率不高(如CT有4096级灰度,而人眼的灰度分辨率仅约40级),且人工长时间大量读图容易产生疲倦而导致准确率降低等问题。
(3)肝脏损伤呈现多样性,人工检测受个人经验限制,检测存在错误或遗漏。此外,研究表明当肝脏损伤小于20毫米时,人工检测的准确率急剧降低。
临床中肝脏损伤病人越来越多,建立一个健壮、客观和可重复的计算机辅助工具以提高肝脏损伤检测的效率和准确率,具有重要的临床意义和价值。
近年来,国内外开始探索基于计算机辅助的肝脏损伤自动检测研究,取得了一些阶段性成果。Bilello使用阈值方法分割肝脏实质,研究支持向量机(Support VectorMachine,SVM)方法分类肝脏实质的空洞和缺口检测损伤。Kumar使用模糊C均值方法分割肝脏肿瘤区域,研究神经网络分类器识别恶性和良性肿瘤。Masuda用EM(Expectationmaximization)算法和高斯混合模型检测肝脏损伤。Milizer使用三级连接概率Boosting树分类肝脏实质和损伤。Pescia使用AdaBoost学习技术,线性组合3个弱分类器构成肝脏损伤检测的强分类器。Wu使用深度学习方法分类肝脏良性和恶性损伤。尽管这些方法在肝脏损伤检测中取得到了阶段性成果,但是还存在几个问题:
(1)大多数方法需要手工标注大量肝脏损伤样本,在临床中很难满足。且肝脏损伤检测需要在肝脏目标区域上提取高维特征,高维空间的数据非常稀少,需要的样本量呈指数级增长。
(2)大多数方法需要损伤类别数作为先验知识,只能把肝脏损伤分为有或无损伤,恶性或良性损伤,低密度或高密度损伤等2类。事实上,暴力性或病理性肝脏损伤存在多种可能,固定损伤的类别数限制了肝脏损伤自动检测的使用范围。
发明内容
为解决上述技术问题,本发明提出一种主动学习的肝脏损伤类别的图像检测方法。其特征在于该方法利用利用肝脏主动形状模式,最优匹配,特征定义和提取、密度聚类和主动学习等技术实现了肝脏损伤的图像检测方法。该方法只需要专家手工标注少量代表性的肝脏损伤,且可以检测多种类型的损伤。具体步骤包括:
步骤1,使用肝脏主动形状模型从医院采集或网上公开的腹部医学图像测试集中定位并提取肝脏目标区域,腹部医学图像测试集中包含腹部医学图像;
步骤2,定义和提取肝脏目标区域中像素点的方向模式特征向量和灰度幅度变化模式特征向量;
步骤3,基于肝脏目标区域中像素点局部密度和到更大密度像素点最小距离二维空间的肝脏像素密度聚类中心定义和聚类划分;
步骤4,基于肝脏像素密度聚类边缘、聚类中心和聚类划分的主动学习方法进行肝脏损伤类别的图像检测。
本发明步骤1包括:
肝脏主动形状模型的几何形状参数包括平移量t、旋转量r和缩放量s,以及纹理缩放参数,所述纹理缩放参数包括全局增缩参数u和偏差参数w;将上述几何形状参数记为模型参数p,p={t,r,s,u,w},模型参数对应的增量为Δp,Δp={Δt,Δr,Δs,Δu,Δw},Δt表示平移量t的增量,Δr表示旋转量r的增量,Δs表示缩放量s的增量,Δu表示全局增缩参数u的增量,Δw表示偏差参数w的增量,通过如下公式对模型参数p进行k次线性模型迭代更新:
其中,i取值范围为1~k,当经过k次线性模型迭代更新后,肝脏模板(由TimCootes提出的主动形状模型(Active Shape Model,ASM)从医院采集或网上公开的训练集中训练出肝脏模型,软件可以从网站http://personalpages.manchester.ac.uk/staff/timothy.f.cootes/下载)与腹部医学图像的部分区域均方差最小时,定位并提取此腹部医学图像区域为肝脏目标区域
本发明步骤2包括以下步骤:
步骤2-1,假设提取的肝脏目标区域中像素点Ic及其8邻域像素点的空间关系如表1所示:
表1
I6 | I7 | I8 |
I5 | Ic | I1 |
I4 | I3 | I2 |
其中,I1、I2、I3、I4、I5、I6、I7和I8为像素点Ic的邻域像素点,分别计算以像素点Ic的4邻域像素点I1、I3、I5和I7为中心的五个方向的灰度差;
步骤2-2,分别计算以像素点Ic,的4邻域像素点I1、I3、I5和I7为中心的方向模式数特征,构造像素点Ic的方向模式特征向量{fc,1,fc,2,fc,3,fc,4},即统计4邻域像素点在5个方向中灰度差小于0的个数:
fc,1=Num(D1<0,D2<0,D3<0,D4<0,D5<0),
其中D1为像素点I8和像素点I1的灰度差,D2为像素点I7和像素点I1的灰度差,D3为像素点Ic和像素点I1的灰度差,D4为像素点I3和像素点I1的灰度差,D5为像素点I2和像素点I1的灰度差;
fc,2=Num(D6<0,D7<0,D8<0,D9<0,D10<0),
其中D6为像素点I2和像素点I3的灰度差,D7为像素点I1和像素点I3的灰度差,D8为像素点Ic和像素点I3的灰度差,D9为像素点I5和像素点I3的灰度差,D10为像素点I4和像素点I3的灰度差;
fc,3=Num(D11<0,D12<0,D13<0,D14<0,D15<0),
其中,D11为像素点I4和像素点I5的灰度差,D12为像素点I3和像素点I5的灰度差,D13为像素点Ic和像素点I5的灰度差,D14为像素点I7和像素点I5的灰度差,D15为像素点I6和像素点I5的灰度差;
fc,4=Num(D16<0,D17<0,D18<0,D19<0,D20<0),
其中,D16为像素点I6和像素点I7的灰度差,D17为像素点I5和像素点I7的灰度差,D18为像素点Ic和像素点I5的灰度差,D19为像素点I1和像素点I7的灰度差,D20为像素点I8和像素点I7的灰度差。
步骤2-3,计算像素点Ic的4邻域像素点I1、I3、I5和I7为中心的五个方向幅度变化模式数特征构造像素点Ic的灰度幅度变化模式特征向量{gc,1,gc,2,gc,3,gc,4},即统计4邻域数据点I1、I3、I5和I7五个方向的幅度变化值大于像素点Ic五个方向的幅度变化值的个数:
其中,M1,1表示像素点I1与相邻的像素点I8、I7、Ic、I3和I2的灰度幅度变化值,计算公式为:
M1,2表示像素点Ic与相邻的像素点I8、I7、I1、I3和I2的灰度幅度变化值,计算公式为:
其中,E1为像素点I8和像素点Ic的灰度差,E2为像素点I7和像素点Ic的灰度差,E3为像素点Ic和像素点I1的灰度差,E4为像素点I3和像素点Ic的灰度差,E5为像素点I2和像素点Ic的灰度差。
其中M3,1表示像素点I3与相邻的像素点I2、I1、Ic、I5和I4的灰度幅度变化值,计算公式为:
其中M3,2表示像素点Ic与相邻的像素点I2、I1、I3、I5和I4的灰度幅度变化值,计算公式为
其中E6为像素点I5和像素点Ic的灰度差,E7为像素点I4和像素点Ic的灰度差。
其中M5,1表示像素点I5与相邻的像素点I4、I3、Ic、I7和I6的灰度幅度变化值,计算公式为:
其中M5,2表示像素点Ic与相邻的像素点I4、I3、I5、I7和I6的灰度幅度变化值,计算公式为:
其中E8为像素点I6和像素点Ic的灰度差。
其中,M7,1表示像素点I7与相邻的像素点I6、I5、Ic、I1和I8的幅度变化值,计算公式为:
其中M7,2表示像素点Ic与相邻的像素点I6、I5、I7、I1和I8的幅度变化值,计算公式为:
本发明步骤3包括以下步骤:
步骤3-1,计算肝脏目标区域中所有像素点高维特征之间的欧氏距离:
设步骤1中定位和提取的的肝脏目标区域共有N个像素点v1,v2,…,vN,根据步骤2提取每个像素点的方向模式特征向量和灰度幅度变化模式特征向量,计算任意两个像素点vm和vn之间的欧氏距离dm,n为:
其中fm,l表示像素点vm的第l个方向模式特征,fn,l表示像素点vn的第l个方向模式特征,gm,l表示像素点vm的第l个灰度幅度变化模式特征,gn,l表示像素点vn的第l个灰度幅度变化模式特征;
步骤3-2,计算每个像素点的局部密度和每个像素点到比其密度更大的像素点的最小距离:
局部密度ρi定义为像素点vi到周围相邻像素点的距离小于阈值dc(本发明,阈值dc取所有像素点的相互距离中由小到大排列的第3%的位置)的像素点个数,通过如下公式计算像素点vi的局部密度ρi:
其中函数di,j表示像素点vi到vj的欧氏距离。
通过如下公式计算像素点vi到任意比它密度大的像素点vj的最小距离δi:
其中ρj表示像素点vj的局部密度。
步骤3-3,基于两维参数空间的肝脏像素密度聚类中心定义和聚类划分:
肝脏像素密度聚类中心Clusters定义为局部密度大于密度阈值ρc,且与比其密度更大的像素点的距离大于最小距离阈值δc的像素点:
Clusters={i|(ρi>ρc)&(δi>δc)},
本发明中,密度阈值ρc取最大密度的25%,即最小距离阈值δc取最大距离的25%,即
确定聚类中心后,聚类的类别数为密度聚类中心的个数,除聚类中心像素点以外的像素点划分到距离聚类中心最近的聚类簇中。
本发明步骤4包括:聚类边缘Edges定义为局部密度小于密度阈值ρc,且与比其密度更大的像素点的距离大于最小距离阈值δc的像素点:
Edges={i|(ρi<ρc)&(δi>δc)},
选择聚类边缘的像素点作为最不确定性样本集,让影像专家手工标注这些像素点,并添加到标注过的数据集;选择步骤3-3定义的聚类中心的像素点作为最有代表性的样本集,让影像专家手工标注这些像素点,标注聚类中心的像素点的损伤类别,并添加到标注过的数据集;与聚类中心同簇的没有标注过的像素点直接标注为聚类中心的损伤类别。
肝脏图像的形状变化较大,一般图像分割方法很难定位和提取。本发明采用基于形状先验知识的肝脏主动形状模型可以准确地定位肝脏在腹部图像中的位置。通过肝脏形状的平移、旋转和缩放,以及纹理的增缩和偏差分析,可以精确定的分割出肝脏。肝脏的实质和损伤区别可以通过局部纹理特征表现出来,传统的局部特征主要利用相邻像素的灰度值,忽略了方向模式和幅度变化模式等局部特征。为此,本发明采用了肝脏像素的局部五个方向模式和幅度变化模式的纹理特征。现有的肝脏损伤检测方法需要手工标注大量的训练样本,且损伤类别有限。本发明通过把像素点的高维空间转换到二维参数空间,实现了基于肝脏像素密度聚类主动学习方法的肝脏损伤类别的图像检测。只需要专家手工标注最有代表性的聚类中心样本,以及最不确定性的聚类边缘样本,且能检测多种类别的损伤。
综上,本发明利用肝脏图像的客观成像规律,在处理步骤中,采用了一些符合这些规律的计算机图像处理方法进行组合。相比以往手工检测的方法,得到了更加客观和准确的损伤检测结果,相比以外计算机辅助检测方法,大大地降低了训练样本量,提高了适用性和效率,属于专利法的保护范围。
有益效果:(1)本发明提出了肝脏像素的方向模式特征向量和灰度幅度变化模式特征向量。一方面,避免了现有局部特征一般都考虑所有相邻像素且仅考虑像素的灰度值,没有考虑方向的变化模式。另一方面,幅度变化模式体现了更多的纹理细节变化,更加适合肝脏微小损伤的特征。
(2)本发明提出了把肝脏高维特征映射到二维参数空间的肝脏像素密度聚类方法。统计一定距离范围内像素点的个数作为局部密度,简单有效。最小距离是局部密度的函数,适合任意形状的聚类。基于二维参数的密度聚类定义和划分,可以快速地把所有的肝脏像素聚类成不同的模式。
(3)基于肝脏像素密度聚类的主动学习方法,提出了新的主动学习代表性元素抽取方法和最不确定性的度量策略。基于主动学习的肝脏损伤只需要专家手工标注少量代表性样本的损伤类别,节约了大量的人力和财力。本发明结合了主动可变模型,医学图像分割,特征定义与提取,肝脏像素密度聚类和主动学习等技术,实现了只需要少量标注样本的肝脏损伤图像检测方法,对医学图像的临床应用有着很高的应用价值。
附图说明
下面结合附图和具体实施方式对本发明做更进一步的具体说明,本发明的上述或其他方面的优点将会变得更加清楚。
图1为本发明流程图。
图2为数据点Ic及其8邻域数据点的空间关系。
具体实施方式
下面结合附图对发明的技术内容进行详细说明:
本方法分为肝脏目标区域提取、特征定义和提取、肝脏像素密度聚类和主动学习检测四个部分,具体的工作流程如图1所示。
(一)使用肝脏主动形状模型从腹部医学图像集中定位并提取肝脏目标区域。
假设肝脏主动形状模型的几何形状参数:平移量t、旋转量r和缩放量s,以及纹理缩放参数:全局增缩u和偏差w,这些模型参数记为p={t,r,s,u,w},模型参数的增量Δp={Δt,Δr,Δs,Δu,Δw},经过k次线性模型迭代更新:
调整参数后的肝脏模板与腹部医学图像的区域均方差最小时,定位并提取腹部图像的对应区域为肝脏。
(二)定义和提取肝脏像素点的方向模式特征向量和灰度幅度变化模式特征向量。
假设提取的像素点Ic及其8邻域数据点的空间关系如图2所示,方向模式特征向量的详细步骤如下:
(21)对于每个以像素点Ic为中心的3×3方块,计算Ic的4邻域像素点(I1,I3,I5,I7)为中心的五个方向的灰度差。
像素点I1与相邻的I8、I7、Ic、I3和I2五个像素点的灰度差:即D1为像素点I8和像素点I1的灰度差,D2为像素点I7和像素点I1的灰度差,D3为像素点Ic和像素点I1的灰度差,D4为像素点I3和像素点I1的灰度差,D5为像素点I2和像素点I1的灰度差;
像素点I3与相邻的I2、I1、Ic、I5和I4五个像素点的灰度差,即D6为像素点I2和像素点I3的灰度差,D7为像素点I1和像素点I3的灰度差,D8为像素点Ic和像素点I3的灰度差,D9为像素点I5和像素点I3的灰度差,D10为像素点I4和像素点I3的灰度差;
像素点I5与相邻的I4、I3、Ic、I7和I6五个像素点的灰度差,即D11为像素点I4和像素点I5的灰度差,D12为像素点I3和像素点I5的灰度差,D13为像素点Ic和像素点I5的灰度差,D14为像素点I7和像素点I5的灰度差,D15为像素点I6和像素点I5的灰度差;
像素点I7与相邻的I6、I5、Ic、I1和I8五个像素点的灰度差,即D16为像素点I6和像素点I7的灰度差,D17为像素点I5和像素点I7的灰度差,D18为像素点Ic和像素点I5的灰度差,D19为像素点I1和像素点I7的灰度差,D20为像素点I8和像素点I7的灰度差。
(22)定义和计算Ic的4邻域像素点(I1,I3,I5,I7)为中心的五个方向模式数特征,构造像素点Ic的方向模式特征向量{fc,1,fc,2,fc,3,fc,4},即统计4邻域像素点在5个方向中灰度差小于0的个数:
fc,1=Num(D1<0,D2<0,D3<0,D4<0,D5<0)
fc,2=Num(D6<0,D7<0,D8<0,D9<0,D10<0)
fc,3=Num(D11<0,D12<0,D13<0,D14<0,D15<0)
fc,4=Num(D16<0,D17<0,D18<0,D19<0,D20<0)
(23)计算Ic的4邻域像素点(I1,I3,I5,I7)为中心的五个方向的幅度变化模式数特征:
像素点I1与相邻的像素点I8、I7、Ic、I3和I2的幅度变化值M1,1:
像素点Ic与相邻的像素点I8、I7、I1、I3和I2的幅度变化值M1,2:
其中,E1为像素点I8和像素点Ic的灰度差,E2为像素点I7和像素点Ic的灰度差,E3为像素点Ic和像素点I1的灰度差,E4为像素点I3和像素点Ic的灰度差,E5为像素点I2和像素点Ic的灰度差。
像素点I3与相邻的像素点I2、I1、Ic、I5和I4的幅度变化值M3,1:
计算像素点Ic与相邻的像素点I2、I1、I3、I5和I4的幅度变化值M3,2:
计算像素点I5与相邻的像素点I4、I3、Ic、I7和I6的幅度变化值M5,1:
计算像素点Ic与相邻的像素点I4、I3、I5、I7和I6的幅度变化值M5,2
其中E8为像素点I6和像素点Ic的灰度差。
计算像素点I7与相邻的像素点I6、I5、Ic、I1和I8的幅度变化值M7,1:
计算像素点Ic与相邻的像素点I6、I5、I7、I1和I8的幅度变化值M7,2:
(24)定义和计算Ic的4邻域像素点(I1,I3,I5,I7)为中心的五个方向幅度变化模式数特征,构造像素点Ic的灰度幅度变化模式特征向量{gc,1,gc,2,gc,3,gc,4},4邻域数据点(I1,I3,I5,I7)五个方向的幅度变化值大于Ic五个方向的幅度变化值的个数。
(三)基于肝脏像素局部密度和到更大密度数据点的最小距离的二维空间的肝脏像素密度聚类中心定义和聚类划分。
步骤31,计算所有像素点高维特征之间的欧氏距离。
设步骤1从医院采集或网上公开的腹部医学图像中定位和提取的肝脏图像集共有N个像素点v1,v2,…,vN,根据步骤2提取每个像素点的方向模式特征向量和灰度幅度变化模式特征向量。任意两个像素点vm和vn之间的欧氏距离为:
其中fm,l表示像素点vm的第l个方向模式特征,fn,l表示像素点vn的第l个方向模式特征,gm,l表示像素点vm的第l个灰度幅度变化模式特征,gn,l表示像素点vn的第l个灰度幅度变化模式特征;
步骤32,计算每个数据点的参数:局部密度和到更大密度数据点的最小距离。
局部密度ρi定义为像素点vi到周围相邻像素点的距离小于阈值dc的像素点个数:
其中函数di,j表示像素点vi到vj的欧氏距离。本发明中,阈值dc为所有像素点的相互距离中由小到大排列的第3%的位置。
参数δi表示像素点vi到任何比它密度大的像素点vj的最小距离:
步骤33,基于两维参数空间的肝脏像素密度聚类中心定义和聚类划分。
密度聚类中心定义为,局部密度大于密度阈值ρc,且距离比其密度更大的像素点大于距离阈值δc的像素点。
Clusters={i|(ρi>ρc)&(δi>δc)} (8)
本发明中密度阈值ρc取最大密度的25%,即距离阈值δc取最大距离的25%,即
确定聚类中心后,聚类的类别数为密度聚类中心的个数,除聚类中心像素点以外的像素点划分到距离聚类中心最近的聚类簇中。
(四)基于肝脏像素密度聚类边缘、中心和聚类划分的主动学习方法进行肝脏损伤类别的图像检测。
聚类边缘定义为:局部密度小于密度阈值ρc,且距离比其密度更大的像素点大于距离阈值δc的像素点。
Edges={i|(ρi<ρc)&(δi>δc)} (9)
选择聚类边缘的像素点作为最不确定性样本集,标注这些像素点的损伤,并添加到标注过的数据集。选择聚类中心的像素点作为最有代表性的样本集,标注这些像素点的损伤,并添加到标注过的数据集。与聚类中心同簇的没有标注过的像素点直接标注为聚类中心的肝脏损伤。
本发明提供了一种主动学习的肝脏损伤图像检测方法,具体实现该技术方案的方法和途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
Claims (5)
1.一种主动学习的肝脏损伤类别的图像检测方法,其特征在于,包括以下步骤:
步骤1,使用肝脏主动形状模型从腹部医学图像测试集中定位并提取肝脏目标区域,腹部医学图像测试集中包含腹部医学图像;
步骤2,定义和提取肝脏目标区域中像素点的方向模式特征向量和灰度幅度变化模式特征向量;
步骤3,基于肝脏目标区域中像素点局部密度和到更大密度像素点最小距离二维空间的肝脏像素密度聚类中心定义和聚类划分;
步骤4,基于肝脏像素密度聚类边缘、聚类中心和聚类划分的主动学习方法进行肝脏损伤类别的图像检测。
2.根据权利要求1所述的方法,其特征在于,步骤1包括:
肝脏主动形状模型的几何形状参数包括平移量t、旋转量r和缩放量s,以及纹理缩放参数,所述纹理缩放参数包括全局增缩参数u和偏差参数w;将上述几何形状参数记为模型参数p,p={t,r,s,u,w},模型参数对应的增量为Δp,Δp={Δt,Δr,Δs,Δu,Δw},通过如下公式对模型参数p进行k次线性模型迭代更新:
其中,i取值范围为1~k,当经过k次线性模型迭代更新后,肝脏模板与腹部医学图像的部分区域均方差最小时,定位并提取此腹部医学图像区域为肝脏目标区域。
3.根据权利要求2所述的方法,其特征在于,步骤2包括以下步骤:
步骤2-1,假设提取的肝脏目标区域中像素点Ic及其8邻域像素点的空间关系如下表所示:
其中,I1、I2、I3、I4、I5、I6、I7和I8为像素点Ic的邻域像素点,分别计算以像素点Ic的4邻域像素点I1、I3、I5和I7为中心的五个方向的灰度差;
步骤2-2,分别计算以像素点Ic的4邻域像素点I1、I3、I5和I7为中心的方向模式数特征,构造像素点Ic的方向模式特征向量{fc,1,fc,2,fc,3,fc,4},即统计4邻域像素点在5个方向中灰度差小于0的个数:
fc,1=Num(D1<0,D2<0,D3<0,D4<0,D5<0),
其中D1为像素点I8和像素点I1的灰度差,D2为像素点I7和像素点I1的灰度差,D3为像素点Ic和像素点I1的灰度差,D4为像素点I3和像素点I1的灰度差,D5为像素点I2和像素点I1的灰度差;
fc,2=Num(D6<0,D7<0,D8<0,D9<0,D10<0),
其中D6为像素点I2和像素点I3的灰度差,D7为像素点I1和像素点I3的灰度差,D8为像素点Ic和像素点I3的灰度差,D9为像素点I5和像素点I3的灰度差,D10为像素点I4和像素点I3的灰度差;
fc,3=Num(D11<0,D12<0,D13<0,D14<0,D15<0),
其中,D11为像素点I4和像素点I5的灰度差,D12为像素点I3和像素点I5的灰度差,D13为像素点Ic和像素点I5的灰度差,D14为像素点I7和像素点I5的灰度差,D15为像素点I6和像素点I5的灰度差;
fc,4=Num(D16<0,D17<0,D18<0,D19<0,D20<0),
其中,D16为像素点I6和像素点I7的灰度差,D17为像素点I5和像素点I7的灰度差,D18为像素点Ic和像素点I5的灰度差,D19为像素点I1和像素点I7的灰度差,D20为像素点I8和像素点I7的灰度差;
步骤2-3,计算像素点Ic的4邻域像素点I1、I3、I5和I7为中心的五个方向幅度变化模式数特征,构造像素点Ic的灰度幅度变化模式特征向量{gc,1,gc,2,gc,3,gc,4},即统计4邻域数据点I1、I3、I5和I7五个方向的幅度变化值大于像素点Ic五个方向的幅度变化值的个数:
其中,M1,1表示像素点I1与相邻的像素点I8、I7、Ic、I3和I2的灰度幅度变化值,计算公式为:
M1,2表示像素点Ic与相邻的像素点I8、I7、I1、I3和I2的灰度幅度变化值,计算公式为:
其中,E1为像素点I8和像素点Ic的灰度差,E2为像素点I7和像素点Ic的灰度差,E3为像素点Ic和像素点I1的灰度差,E4为像素点I3和像素点Ic的灰度差,E5为像素点I2和像素点Ic的灰度差;
其中,M3,1表示像素点I3与相邻的像素点I2、I1、Ic、I5和I4的灰度幅度变化值,计算公式为:
其中,M3,2表示像素点Ic与相邻的像素点I2、I1、I3、I5和I4的灰度幅度变化值,计算公式为:
其中,E6为像素点I5和像素点Ic的灰度差,E7为像素点I4和像素点Ic的灰度差;
其中M5,1表示像素点I5与相邻的像素点I4、I3、Ic、I7和I6的灰度幅度变化值,计算公式为:
其中M5,2表示像素点Ic与相邻的像素点I4、I3、I5、I7和I6的灰度幅度变化值,计算公式为:
其中,E8为像素点I6和像素点Ic的灰度差;
其中,M7,1表示像素点I7与相邻的像素点I6、I5、Ic、I1和I8的幅度变化值,计算公式为:
其中M7,2表示像素点Ic与相邻的像素点I6、I5、I7、I1和I8的幅度变化值,计算公式为:
4.根据权利要求3所述的方法,其特征在于,步骤3包括以下步骤:
步骤3-1,计算肝脏目标区域中所有像素点高维特征之间的欧氏距离:
设步骤1中定位和提取的肝脏目标区域共有N个像素点v1,v2,…,vN,根据步骤2提取每个像素点的方向模式特征向量和灰度幅度变化模式特征向量,计算任意两个像素点vm和vn之间的欧氏距离dm,n为:
其中fm,l表示像素点vm的第l个方向模式特征,fn,l表示像素点vn的第l个方向模式特征,gm,l表示像素点vm的第l个灰度幅度变化模式特征,gn,l表示像素点vn的第l个灰度幅度变化模式特征;
步骤3-2,计算每个像素点的局部密度和每个像素点到比其密度更大的像素点的最小距离:
局部密度ρi定义为像素点vi到周围相邻像素点的距离小于阈值dc的像素点个数,通过如下公式计算像素点vi的局部密度ρi:
其中函数di,j表示像素点vi到vj的欧氏距离;
通过如下公式计算像素点vi到任意比它密度大的像素点vj的最小距离δi:
其中ρj表示像素点vj的局部密度;
步骤3-3,基于两维参数空间的肝脏像素密度聚类中心定义和聚类划分:
肝脏像素密度聚类中心Clusters定义为局部密度大于密度阈值ρc,且与比其密度更大的像素点的距离大于最小距离阈值δc的像素点:
Clusters={i|(ρi>ρc)&(δi>δc)},
确定聚类中心后,聚类的类别数为密度聚类中心的个数,除聚类中心像素点以外的像素点划分到距离聚类中心最近的聚类簇中。
5.根据权利要求4所述的方法,其特征在于,步骤4包括:聚类边缘Edges定义为局部密度小于密度阈值ρc,且与比其密度更大的像素点的距离大于最小距离阈值δc的像素点:
Edges={i|(ρi<ρc)&(δi>δc)},
选择聚类边缘的像素点作为最不确定性样本集,标注聚类边缘的像素点的损失类别,并添加到标注过的数据集;
选择步骤3-3定义的聚类中心的像素点作为最有代表性的样本集,标注聚类中心的像素点的损伤类别,并添加到标注过的数据集;
与聚类中心同簇的没有标注过的像素点直接标注为聚类中心的损伤类别。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611062637.9A CN106778829B (zh) | 2016-11-28 | 2016-11-28 | 一种主动学习的肝脏损伤类别的图像检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611062637.9A CN106778829B (zh) | 2016-11-28 | 2016-11-28 | 一种主动学习的肝脏损伤类别的图像检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106778829A true CN106778829A (zh) | 2017-05-31 |
CN106778829B CN106778829B (zh) | 2019-04-30 |
Family
ID=58901703
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611062637.9A Active CN106778829B (zh) | 2016-11-28 | 2016-11-28 | 一种主动学习的肝脏损伤类别的图像检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106778829B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108875734A (zh) * | 2018-05-23 | 2018-11-23 | 平安科技(深圳)有限公司 | 肝脏癌变定位方法、装置及存储介质 |
CN109034190A (zh) * | 2018-06-15 | 2018-12-18 | 广州深域信息科技有限公司 | 一种动态选择策略的主动样本挖掘的物体检测***及方法 |
CN111274425A (zh) * | 2020-01-20 | 2020-06-12 | 平安科技(深圳)有限公司 | 医疗影像分类方法、装置、介质及电子设备 |
CN111539354A (zh) * | 2020-04-27 | 2020-08-14 | 易普森智慧健康科技(深圳)有限公司 | 一种液基细胞玻片扫描区域识别方法 |
US20210035296A1 (en) * | 2018-04-27 | 2021-02-04 | Delphinus Medical Technologies, Inc. | System and method for feature extraction and classification on ultrasound tomography images |
CN113408670A (zh) * | 2021-08-17 | 2021-09-17 | 深圳电通信息技术有限公司 | 图像分类方法和装置 |
US20230147029A1 (en) * | 2019-03-07 | 2023-05-11 | Beautai Ltd. | Systems and methods for automated makeup application |
CN117854623A (zh) * | 2024-01-04 | 2024-04-09 | 山东瑞谱检测技术有限公司 | 一种自适应网格的医药辅料浓度分布快速检测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101477628A (zh) * | 2009-01-06 | 2009-07-08 | 青岛海信电子产业控股股份有限公司 | 车辆阴影去除方法和装置 |
CN104299185A (zh) * | 2014-09-26 | 2015-01-21 | 京东方科技集团股份有限公司 | 一种图像放大方法、图像放大装置及显示设备 |
-
2016
- 2016-11-28 CN CN201611062637.9A patent/CN106778829B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101477628A (zh) * | 2009-01-06 | 2009-07-08 | 青岛海信电子产业控股股份有限公司 | 车辆阴影去除方法和装置 |
CN104299185A (zh) * | 2014-09-26 | 2015-01-21 | 京东方科技集团股份有限公司 | 一种图像放大方法、图像放大装置及显示设备 |
Non-Patent Citations (2)
Title |
---|
凌华强: "主动形体模型法在肝脏CT图像分割中的应用", 《浙江工业大学学报》 * |
陶辉: "基于密度峰值选取聚类中心的优化", 《内江科技》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20210035296A1 (en) * | 2018-04-27 | 2021-02-04 | Delphinus Medical Technologies, Inc. | System and method for feature extraction and classification on ultrasound tomography images |
CN108875734A (zh) * | 2018-05-23 | 2018-11-23 | 平安科技(深圳)有限公司 | 肝脏癌变定位方法、装置及存储介质 |
CN109034190A (zh) * | 2018-06-15 | 2018-12-18 | 广州深域信息科技有限公司 | 一种动态选择策略的主动样本挖掘的物体检测***及方法 |
CN109034190B (zh) * | 2018-06-15 | 2022-04-12 | 拓元(广州)智慧科技有限公司 | 一种动态选择策略的主动样本挖掘的物体检测***及方法 |
US20230147029A1 (en) * | 2019-03-07 | 2023-05-11 | Beautai Ltd. | Systems and methods for automated makeup application |
US11983958B2 (en) * | 2019-03-07 | 2024-05-14 | Gemma Robotics, Inc. | Systems and methods for automated makeup application |
CN111274425A (zh) * | 2020-01-20 | 2020-06-12 | 平安科技(深圳)有限公司 | 医疗影像分类方法、装置、介质及电子设备 |
WO2021147217A1 (zh) * | 2020-01-20 | 2021-07-29 | 平安科技(深圳)有限公司 | 医疗影像分类方法、装置、介质及电子设备 |
CN111274425B (zh) * | 2020-01-20 | 2023-08-22 | 平安科技(深圳)有限公司 | 医疗影像分类方法、装置、介质及电子设备 |
CN111539354A (zh) * | 2020-04-27 | 2020-08-14 | 易普森智慧健康科技(深圳)有限公司 | 一种液基细胞玻片扫描区域识别方法 |
CN113408670A (zh) * | 2021-08-17 | 2021-09-17 | 深圳电通信息技术有限公司 | 图像分类方法和装置 |
CN117854623A (zh) * | 2024-01-04 | 2024-04-09 | 山东瑞谱检测技术有限公司 | 一种自适应网格的医药辅料浓度分布快速检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106778829B (zh) | 2019-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106778829B (zh) | 一种主动学习的肝脏损伤类别的图像检测方法 | |
Soltaninejad et al. | Automated brain tumour detection and segmentation using superpixel-based extremely randomized trees in FLAIR MRI | |
Zhao et al. | Automatic nasopharyngeal carcinoma segmentation using fully convolutional networks with auxiliary paths on dual-modality PET-CT images | |
US10667778B2 (en) | Accurate detection and assessment of radiation induced lung injury based on a computational model and computed tomography imaging | |
US9092691B1 (en) | System for computing quantitative biomarkers of texture features in tomographic images | |
Wajid et al. | Local energy-based shape histogram feature extraction technique for breast cancer diagnosis | |
Yu et al. | An automatic computer-aided detection scheme for pneumoconiosis on digital chest radiographs | |
Rehman et al. | Lung cancer detection and classification from chest CT scans using machine learning techniques | |
Kakar et al. | Automatic segmentation and recognition of lungs and lesion from CT scans of thorax | |
CN109035160A (zh) | 医学影像的融合方法及基于融合医学影像学习的图像检测方法 | |
Peng et al. | Segmentation of lung in chest radiographs using hull and closed polygonal line method | |
Li et al. | Learning image context for segmentation of the prostate in CT-guided radiotherapy | |
CN108776962A (zh) | 一种构建肺肿瘤良恶性预测模型的方法 | |
CN109753997B (zh) | 一种ct图像中的肝脏肿瘤自动精确鲁棒分割方法 | |
Li et al. | Learning image context for segmentation of prostate in CT-guided radiotherapy | |
CN104268873A (zh) | 基于核磁共振图像的乳腺肿瘤分割方法 | |
Li et al. | Brain tumor segmentation based on region of interest-aided localization and segmentation U-Net | |
Jony et al. | Detection of lung cancer from CT scan images using GLCM and SVM | |
Yao et al. | Computer-assisted diagnosis of lung cancer using quantitative topology features | |
Hoogi et al. | A fully-automated pipeline for detection and segmentation of liver lesions and pathological lymph nodes | |
Bahadure et al. | Feature extraction and selection with optimization technique for brain tumor detection from MR images | |
Chandra et al. | Localization of the suspected abnormal region in chest radiograph images | |
CN104835155A (zh) | 基于分形的早期乳腺癌钙化点计算机辅助检测方法 | |
CN109191452A (zh) | 一种基于主动学习的腹腔ct图像腹膜转移自动标记方法 | |
Pawar et al. | Tuberculosis screening using digital image processing techniques |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |