技术领域
本发明涉及图像处理技术领域,尤其涉及胸腔积液性质分析决策函数的构造方法及基于其的分析方法。
背景技术
胸腔积液是一种致命性疾病,而且胸腔积液是肿瘤、胸膜炎症等多种疾病的并发症,因此快速而又精确地判断胸腔积液的性质是非常有必要的。目前,主流的用于分析胸腔积液性质的方法是在有创引流后进行生化检验,这对患者会造成一定的痛苦,并且会花费较多的时间。超声影像是一种无创的辅助诊断方法,但是目前在业界还没有公认的通过对超声图像的特征进行分析来判断积液性质的方法。
本发明将超声图像特征对应到生化检验的指标,实现无创且快速地定量分析积液性质。该方法的可行性基于以下事实:超声图像的胸腔积液区域普遍为低回声区或无回声区,而且不同性质的胸腔积液在超声图像中的表现也有较大差别。此外,许多超声图像积液区域特征能与医生常用的从超声图像积液区域中获取的用于判断积液性质的信息相对应,比如超声图像积液区域的“均值”和“标准差”对应于医生获取的“亮度”和“对比度”信息。
本发明涉及的既有方法包括两大方面,其一为图像特征提取的相关技术,其二为支持向量分类机的相关技术。其中图像特征提取指的是计算机通过算法提取出图像特征,并且用这些特征来描述图像,实现了对图像的抽象化;而支持向量分类机是借助于最优化方法解决数据挖掘中若干问题的有力工具。
发明内容
本发明的目的在于提供一种胸腔积液性质分析决策函数的构造方法及基于其的分析方法,将超声图像特征对应到生化检验的指标,实现无创且快速的定量分析积液性质。
为了实现上述目的,本发明提供一种胸腔积液性质分析决策函数的构造方法,包括:
基于视频采样,获取待处理的图像帧;
基于直方图匹配的图像处理,对获得的图像帧进行图像预处理;
使用U-Net网络对预处理后的图像进行分割,生成胸腔积液区域分割图像;
基于预处理后的图像和所生成的胸腔积液区域分割图像进行特征提取,获得图像特征;
筛选图像特征,并基于筛选后的图像特征构造决策函数。
一种实现方式中,所述基于视频采样,获取待处理的图像帧的步骤:
从视频中采样胸腔超声影像中不受膈肌或肋骨影响的图像帧,作为待处理的图像帧。
一种实现方式中,所述基于直方图匹配的图像处理,对获得的图像帧进行图像预处理的步骤:
选定基准图像,将基准图像的无效区作为直方图匹配的目标对象对其他图像的无效区分别执行直方图匹配,使其他图像的无效区的直方图与基准图像的无效区的直方图的相似度较高;
获取每一次直方图匹配对应的灰度映射;
将其他图像根据其对应的灰度映射改变,改变后的图像即为预处理之后的图像。
一种实现方式中,所述基于预处理后的图像和所生成的胸腔积液区域分割图像进行特征提取,获得图像特征的步骤:
在预处理后的图像中的胸腔积液区域分割图像确定的区域提取第一数量个图像特征;
在预处理后的图像中的灰度值非零的胸腔积液区域的分割图像确定的区域提取第二数量个图像特征。
一种实现方式中,所述筛选图像特征,并基于筛选后的图像特征构造决策函数的步骤:
将化验指标按照其值进行分类;
使用Relief-F进行特征选择;
基于选择后的特征,使用十次十折交叉验证法为支持向量分类机计算准确率;
更改支持向量机的参数和特征选择步骤保留的特征的比例,重复执行所述使用Relief-F进行特征选择的步骤和所述基于选择后的特征,使用十次十折交叉验证法为支持向量分类机计算准确率的步骤;
使用得分最高的支持向量机参数和保留的特征的比例组合计算决策函数。
以及公开一种基于决策函数的胸腔积液性质分析方法,包括:获取带检测病人的化验指标;
(1)设共L个化验指标,每个化验指标t
(2)设从病人的超声视频中采样共M张图片,把图片分别记为I
(3)令
(4)将I
(5)把
(6)令p=p+1,若p≤M,则从步骤(4)开始执行,否则从步骤(7)开始执行;
(7)可得化验指标t
(8)令正在处理的化验指标编号k=k+1,若k≤L,则从步骤(3)开始执行,否则结束,完成胸腔积液性质分析。
本发明还提供了一种基于决策函数的胸腔积液性质分析方法,其特征在于,包括:获取带检测病人的化验指标;
(1)设共L个化验指标,每个化验指标t
(2)设从病人的超声视频中采样共M张图片,把图片分别记为I
(3)令
(4)将I
(5)把
(6)令p=p+1,若p≤M,则从步骤(4)开始执行,否则从步骤(7)开始执行;
(7)可得化验指标t
(8)令正在处理的化验指标编号k=k+1,若k≤L,则从步骤(3)开始执行,否则结束,完成胸腔积液性质分析。
应用本发明实施例提供的胸腔积液性质分析决策函数的构造方法及基于其的胸腔积液性质分析方法,将超声图像特征对应到生化检验的指标,实现无创且快速的定量分析积液性质。
附图说明
图1是本发明实施例胸腔积液性质分析决策函数的构造方法的一种流程示意图。
图2是本发明实施例基于决策函数的胸腔积液性质分析方法的一种流程示意图。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。
如图1本发明提供一种胸腔积液性质分析决策函数的构造方法,包括:
S110,基于视频采样,获取待处理的图像帧;
需要说明的是,医生会在每个肋间对肺部进行扫描,每扫描一个肋间就生成一个dcm格式的视频,从视频中采样得到不受膈肌或肋骨影响的内容尽可能不同的帧。
具体的,胸腔积液超声影像的某些帧中由于膈肌或肋骨的遮挡,胸腔积液的区域会部分或全部变暗,这会影响积液区域的图像特征,而且积液轮廓也会部分或彻底消失,U-Net难以根据这些帧划分出正确的积液区域,所以这些帧应当被去除。这些帧在下文中将被称为“不可用的帧”,其余的帧将被称为“可用的帧”。当“可用的帧”与“不可用的帧”之间相互转变时,灰度值低于t的像素的数量会迅速地增加或减少,称处于转变状态的帧为过渡帧,且称较暗的帧大幅增加的过渡帧为A类过渡帧,称较暗的帧大幅减少的过渡帧为B类过渡帧。只需要找到所有的A类过渡帧和B类过渡帧即可确定哪些帧是“不可用的帧”。
特殊地,第一组过渡帧之前的帧的类别判断方法是从第一帧开始往后检查,如果最先遇到的过渡帧是A类过渡帧,则认为这组过渡帧之前的帧为“可用的帧”,如果最先遇到的过渡帧是B类过渡帧,则认为这组过渡帧之前的帧为“不可用的帧”,如果没有遇到过渡帧,则认为该视频的所有帧都是可用的帧。
为了进一步提高U-Net划分的准确率,所有的过渡帧也都被算作是“不可用的帧”。选取出“可用的帧”之后,由于“可用的帧”中有许多内容相似的帧,所以应当仅取其中内容尽可能不同的帧。具体步骤如下:
1.提取胸腔超声影像中不受膈肌或肋骨影响的帧;
(1)选择合适的灰度阈值t和像素数阈值n。
(2)设共L帧,每帧记为F
(3)计算相邻两帧之间的灰度值低于t的像素个数的差值。
one
计算相邻三帧之间的灰度值低于t的像素个数的差值。
two
(4)找出过渡帧
(5)令status=1,p=1,i=1
(6)若status==1且Up
若status==1且Down
若status==0且Up
(7)i=i+1,若i≤L-1,则返回(6),否则执行(8)
(8)若status==1且p+9<L+1,将F
2.选取内容尽可能不同的帧
(1)设当前共有M帧,每帧记为F
(2)找出灰度均值最大的帧和灰度均值最小的帧,两帧的下标分别记为Imax和Imin。
(3)MSE为均方误差,求standard=MSE(hist
(4)令p=1,i=2,将F
(5)求e
(6)i=i+1,若i≤M,则转到(5),否则算法结束。
S120,基于直方图匹配的图像处理,对获得的图像帧进行图像预处理;
可以理解的是,由于超声设备参数的不同选取,不同胸腔积液超声图像的差异不仅仅是由于病人的肺部情况不同,比如增益Gain,它指的是回声信号的放大程度,可以使图像变亮或变暗;还有时间增益补偿TGC,它可放大返回的信号,以纠正增加组织深度时所导致的衰减。为了减少超声设备参数的不同选取导致的超声图像不同,本方法提出了基于直方图匹配的图像预处理方案。本方法的有效性基于一个事实:对于一个“可用的帧”来说,如果超声设备的深度Depth参数足够大,则胸腔积液超声图像的左下角和右下角均为“无效区”,其中“无效区”指的是这个区域不能体现任何人体组织信息,是纯粹的噪点。而且“无效区”的内容能被超声设备参数的不同选取影响。
预处理中使用到的一个关键算法是直方图匹配,它的作用是可以使一张图I
预处理算法步骤如下:
(1)设图像集为{G
(2)将直方图匹配记作〈I
(3)将G
(4)计算
S130,使用U-Net网络对预处理后的图像进行分割,生成蒙版图像;
需要说明的是,由于需要确定超声图像中的积液区域,故本方法中使用U-Net网络对图像进行分割,生成积液蒙版。此外,还需要生成无零积液蒙版,原因如下:一般的超声图像的积液区域有大量灰度值为0的像素,这是因为液体内部十分均匀,其声阻抗无多大差别,只有很少的反射界面形成。对于超声图像特征来说,从不同超声图像的反射界面对应的部分提取到的特征值差异较大,而从不同超声图像的灰度值为0的区域提取到的特征值是相同的。因此对于某些图像特征比如标准差来说,考虑灰度值为0的像素会降低从不同图像的积液区域中提取出的特征值的差别。
(1)使用U-Net网络对图像进行分割,生成积液蒙版。
使用U-Net网络对图像进行分割。设超声图像中U-Net认为属于积液区域的像素的标签为Lable1,其余像素的标签为Lable2。积液蒙版的大小和其对应的帧一样,若帧中一个像素的标签为Lable1,则积液蒙版的对应位置为1,若帧中的一个像素的标签为Lable2,则积液蒙版的对应位置为0。
(2)生成无零积液蒙版
若帧中的一个像素的标签为Lable2,则无零积液蒙版的对应位置为0;若帧中的一个像素的标签为Lable1且该像素的灰度值为0,则无零积液蒙版的对应位置为0;若帧中的一个像素的标签为Lable1且该像素的灰度值不为0,则无零积液蒙版的对应位置为1。
S140,基于预处理后的图像和所生成的蒙版图像进行特征提取,获得图像特征;
可以理解的是,本发明在超声图像积液区域提取共80个特征。
具体的,在超声图像和积液蒙版确定的超声图像区域提取40个图像特征;
(1)计算以下6个图像特征:
设积液区域总像素数为N,积液区域像素值为i的像素个数为n
三阶距:
一致性:
熵:
(2)选用四个方向(0,1),(1,0),(1,1),(-1,1),针对所选超声图像区域计算四个灰度共生矩阵。对每个灰度共生矩阵计算以下6个特征:
‘contrast’:
‘dissimilarity’:
‘homogeneity’:
‘ASM’:
‘energy’:
‘correlation’:
其中,G为参与计算的GLCM,n等于G的所有元素之和。
令
其中,i=0,1,…,255,j=0,1,…,255
(3)计算关于超声图像的积液区域的旋转不变的均匀模式LBP。取的邻域包括且仅包括中心像素周围的8个像素。将LBP的计算结果作为10个特征,其中第i个(i=1,2,…,9)特征为邻域中“0”的个数为i个的0-1之间跳变次数小于等于两次的中心像素的个数,第10个特征为0-1之间跳变次数大于两次的中心像素的个数。
需要说明的是,在超声图像和无零积液蒙版确定的超声图像区域同样提取40个图像特征,方法步骤与在超声图像和积液蒙版确定的超声图像区域提取40个图像特征相同,在此不做赘述。
S150,筛选图像特征,并给予筛选后的图像特征构造决策函数。
需要说明的是,由于化验指标是数值型数据,而进行特征选择,支持向量分类机需要的都是类别标签,所以首先需要将化验指标的值映射到对应的类别标签。
由于从图像中提取的特征很多,而其中有一些特征对于要预测的化验指标来说分类能力较弱,这些特征不仅会增加计算时间,而且这些特征会扰乱支持向量机的分类,导致分类正确率降低,因此接下来要进行特征选择来筛选出与预测的化验指标最相关的||V
本发明使用的特征提取方法是Relief-F。Relief-F是著名的过滤式特征选择方法,该方法设计了一个“相关统计量”来度量特征的重要性,该统计量是一个向量,每个分量分别对应于一个初始特征,分量值越大,则对应的初始特征的分类能力就越强。最终,选择相关统计量分量最大的||V
最后则使用支持向量分类机得出决策函数。
由于筛选图像特征和使用支持向量分类机得出决策函数都需要事先决定使用哪个参数值。前者需要决定的参数是保留的特征占总特征的比例p,
由于本方法需要定量地估计多个化验指标,所以每估计一个化验指标都要执行一次步骤S150。
将化验指标按照其值进行分类
每种化验指标都有医学上的正常范围,化验指标的值低于正常范围即为偏低,高于正常范围即为偏高,属于正常范围即为正常。
若某一化验指标值偏低,则其对应的类别标签为“0”;若某一化验指标值正常,则其对应的类别标签为“1”;若某一化验指标值偏高,则其对应的类别标签为“2”。
使用Relief-F进行特征选择
使用当前(p,v)组合中的p值执行一次特征选择。
确定相关统计量的方法如下:
给定数据集D={(x
其中p
最终,选择相关统计量分量最大的||V
使用十次十折交叉验证法为支持向量分类机计算准确率
使用十次十折交叉验证法计算当前(p,v)组合对应的分类准确率。支持向量分类机使用选出的图像特征和当前(p,v)组合中的v值。
v-支持向量分类机算法步骤如下:
(1)给定训练集T={(x
(2)选择线性核函数作为核函数K(x,x′)
(3)选择参数v
(4)求解凸二次规划
得解
(5)计算属
(6)构造决策函数
f(x)=sgn(g(x)),
其中
使用得分最高的(p,v)组合计算决策函数
将得分最高的(p,v)组合记为(p
如图2本发明还提供一种基于决策函数的胸腔积液性质分析方法,其特征在于,包括:获取带检测病人的化验指标;
(1)设共L个化验指标,每个化验指标t
(2)设从病人的超声视频中采样共M张图片,把图片分别记为I
(3)令
(4)将I
(5)把
(6)令p=p+1,若p≤M,则从步骤(4)开始执行,否则从步骤(7)开始执行;
(7)可得化验指标t
(8)令正在处理的化验指标编号k=k+1,若k≤L,则从步骤(3)开始执行,否则结束,完成胸腔积液性质分析。
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。
机译: 结构分析方法,结构分析装置,结构分析程序,用于结构分析的物理性质计算方法,用于结构分析的物理性质计算装置以及用于结构分析的物理性质计算程序
机译: 热物理性质分析方法,使用该方法的热物理性质分析装置以及热物理性质分析程序
机译: 固液界面附近物理性质变化的分析装置和方法固液界面附近物理性质变化的分析方法