首页 / 专利库 / 动物学 / 哺乳动物 / 有蹄类动物 / 马科动物 / / 睑板腺图像处理方法

睑板腺图像处理方法

阅读:279发布:2020-10-28

专利汇可以提供睑板腺图像处理方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种睑板腺 图像处理 方法,包含以下步骤:获取睑板腺图像;将睑板腺图像进行灰度化;对灰度化后的睑板腺图像进行降噪处理;计算分割 阈值 ;根据分割阈值对降噪处理后的睑板腺图像进行二值化处理得到二值图;对二值图进行闭运算后保留最大区域并对最大区域进行填充;提取二值图中的最大区域的边界;将边界映射到降噪处理后的睑板腺图像。本发明提供的睑板腺图像处理方法首先计算出一个科学的分割阈值,通过该分割阈值对睑板腺图像进行二值化,再对二值图进行闭运算,区域选取以及区域填充等操作得到睑板腺区域,找出睑板腺区域的边界后将边界映射至睑板腺图像,快速将睑板腺图像中的睑板腺区域分割出来。,下面是睑板腺图像处理方法专利的具体信息内容。

1.一种睑板腺图像处理方法,其特征在于,包含以下步骤:
获取睑板腺图像;
将所述睑板腺图像进行灰度化;
对灰度化后的所述睑板腺图像进行降噪处理;
计算分割阈值
根据所述分割阈值对降噪处理后的所述睑板腺图像进行二值化处理得到二值图;
对所述二值图进行闭运算后保留最大区域并对所述最大区域进行填充;
提取所述二值图中的所述最大区域的边界;
将所述边界映射到降噪处理后的所述睑板腺图像。
2.根据权利要求1所述的睑板腺图像处理方法,其特征在于,
所述获取睑板腺图像的具体方法为:
通过照明装置照明;
图像采集装置的成像显示器的显示界面设置粗选框;
调整图像采集装置的采集头的位置以将睑板腺区域完全包含在所述粗选框内后采集图像得到所述睑板腺图像。
3.根据权利要求2所述的睑板腺图像处理方法,其特征在于,
所述照明装置照射红外光。
4.根据权利要求1所述的睑板腺图像处理方法,其特征在于,
所述对灰度化的所述睑板腺图像进行降噪处理的具体方法为;
对灰度化的所述睑板腺图像进行均值滤波降噪。
5.根据权利要求1所述的睑板腺图像处理方法,其特征在于,
所述计算分割阈值的具体方法为:
使用最大类间方法计算降噪处理后的所述睑板腺图像的灰度阈值;
统计降噪处理后的所述睑板腺图像的灰度值得到直方分布图;
根据下述公式计算分割阈值,

其中,Th为所述分割阈值,a为所述灰度阈值,b为255,i为灰度值,f(i)为所述直方分布图中灰度值为i的像素的个数。
6.根据权利要求5所述的睑板腺图像处理方法,其特征在于,
所述根据所述分割阈值对降噪处理后的所述睑板腺图像进行二值化处理得到二值图的具体方法为:
将降噪处理后的所述睑板腺图像中灰度值在0.9Th-1.1Th之间的区域的灰度值置为
255且其他区域的灰度值置为0得到所述二值图。
7.根据权利要求1所述的睑板腺图像处理方法,其特征在于,
所述提取所述二值图中的所述最大区域的边界的具体方法是:
通过canny算子提取所述二值图中的所述最大区域的边界。
8.根据权利要求1所述的睑板腺图像处理方法,其特征在于,
在所述将所述边界映射到降噪处理后的所述睑板腺图像之后,所述睑板腺图像处理方法还包含步骤:
通过所述边界内的睑板腺区域的灰度值范围对全图进行灰度拉伸;
对灰度拉伸后的所述睑板腺图像进行图像增强。
9.根据权利要求8所述的睑板腺图像处理方法,其特征在于,
所述通过所述边界内的睑板腺区域的灰度值范围对全图进行灰度拉伸的具体方法为:
获取所述边界内的睑板腺区域的最小灰度值和最大灰度值;
将降噪处理后的所述睑板腺图像中的灰度值小于所述最小灰度值的区域的灰度值置为0;
将降噪处理后的所述睑板腺图像中的灰度值大于所述最大灰度值的区域的灰度值置为255;
将降噪处理后的所述睑板腺图像中的灰度值介于所述最小灰度值和所述最灰度值范围内的区域通过下述伽变换公式进行灰度转换,

其中,s为变换后的灰度值,r为变换前的灰度值,c为变换常数,γ为伽马系数,γ的取值范围为大于等于1.2且小于等于2。
10.根据权利要求9所述的睑板腺图像处理方法,其特征在于,
所述对灰度拉伸后的所述睑板腺图像进行图像增强的具体方法为:
对灰度拉伸后的所述睑板腺图像经过高斯高通滤波器滤波,设定第一截止频率,得到第一高频图后将所述第一高频图叠加至灰度拉伸后的所述睑板腺图像得到第一增强图;
对所述第一增强图经过高斯高通滤波器滤波,设定第二截止频率,得到第二高频图后将所述第二高频图叠加至灰度拉伸后的所述睑板腺图像得到第二增强图;
对所述第二增强图经过高斯高通滤波器滤波,设定第三截止频率,得到第三高频图后将所述第三高频图叠加至灰度拉伸后的所述睑板腺图像得到最终增强图;
其中,所述第一截止频率小于第二截止频率,所述第二截止频率小于第三截止频率。

说明书全文

睑板腺图像处理方法

技术领域

[0001] 本发明涉及一种睑板腺图像处理方法。

背景技术

[0002] 睑板腺埋藏于上下睑板之中,其开口位于睑缘,分泌多种脂质成分,构成泪膜的脂质层,影响泪液蒸发及泪膜的稳定性,因此,睑板腺对维持眼表健康至关重要。睑板腺形态与其功能相关,炎症、堵塞、激素平变化、环境因素等多种原因均可以影响睑板腺的数量和形态,进而改变睑板腺油脂分泌功能,导致泪膜稳定性下降,发生睑板腺功能障碍(MGD)。MGD是蒸发过强型干眼症的主要病因,蒸发过强型干眼症约占据干眼构成的86%,是一组困扰人们视觉质量乃至生活质量的常见病。睑板腺功能障碍是蒸发过强型干眼症的主要原因,是判断是否干眼的一个重要指标。
[0003] 现有方案拍摄睑板腺图像观察睑板腺的形态结构,主要缺陷在于:人工判断睑板腺区域一方面不够精确,且花费的时间较长。另一方面,从拍摄的睑板腺图像中很难观察到清晰的睑板腺结构。

发明内容

[0004] 本发明提供了一种睑板腺图像处理方法,采用如下的技术方案:一种睑板腺图像处理方法,包含以下步骤:
获取睑板腺图像;
将睑板腺图像进行灰度化;
对灰度化后的睑板腺图像进行降噪处理;
计算分割阈值
根据分割阈值对降噪处理后的睑板腺图像进行二值化处理得到二值图;
对二值图进行闭运算后保留最大区域并对最大区域进行填充;
提取二值图中的最大区域的边界;
将边界映射到降噪处理后的睑板腺图像。
[0005] 进一步地,获取睑板腺图像的具体方法为:通过照明装置照明;
图像采集装置的成像显示器的显示界面设置粗选框;
调整图像采集装置的采集头的位置以将睑板腺区域完全包含在粗选框内后采集图像得到睑板腺图像。
[0006] 进一步地,照明装置照射红外光。
[0007] 进一步地,对灰度化的睑板腺图像进行降噪处理的具体方法为;对灰度化的睑板腺图像进行均值滤波降噪。
[0008] 进一步地,计算分割阈值的具体方法为:使用最大类间方法计算降噪处理后的睑板腺图像的灰度阈值;
统计降噪处理后的睑板腺图像的灰度值得到直方分布图;
根据下述公式计算分割阈值,
 ,
其中,Th为分割阈值,a为灰度阈值,b为255,i为灰度值,f(i)为直方分布图中灰度值为i的像素的个数。
[0009] 进一步地,根据分割阈值对降噪处理后的睑板腺图像进行二值化处理得到二值图的具体方法为:将降噪处理后的睑板腺图像中灰度值在0.9Th-1.1Th之间的区域的灰度值置为255且其他区域的灰度值置为0得到二值图。
[0010] 进一步地,提取二值图中的最大区域的边界的具体方法是:通过canny算子提取二值图中的最大区域的边界。
[0011] 进一步地,在将边界映射到降噪处理后的睑板腺图像之后,睑板腺图像处理方法还包含步骤:通过边界内的睑板腺区域的灰度值范围对全图进行灰度拉伸;
对灰度拉伸后的睑板腺图像进行图像增强。
[0012] 进一步地,通过边界内的睑板腺区域的灰度值范围对全图进行灰度拉伸的具体方法为:获取边界内的睑板腺区域的最小灰度值和最大灰度值;
将降噪处理后的睑板腺图像中的灰度值小于最小灰度值的区域的灰度值置为0;
将降噪处理后的睑板腺图像中的灰度值大于最大灰度值的区域的灰度值置为255;
将降噪处理后的睑板腺图像中的灰度值介于最小灰度值和最灰度值范围内的区域通过下述伽变换公式进行灰度转换,

其中,s为变换后的灰度值,r为变换前的灰度值,c为变换常数,γ为伽马系数,γ的取值范围为大于等于1.2且小于等于2。
[0013] 进一步地,对灰度拉伸后的睑板腺图像进行图像增强的具体方法为:对灰度拉伸后的睑板腺图像经过高斯高通滤波器滤波,设定第一截止频率,得到第一高频图后将第一高频图叠加至灰度拉伸后的睑板腺图像得到第一增强图;
对第一增强图经过高斯高通滤波器滤波,设定第二截止频率,得到第二高频图后将第二高频图叠加至灰度拉伸后的睑板腺图像得到第二增强图;
对第二增强图经过高斯高通滤波器滤波,设定第三截止频率,得到第三高频图后将第三高频图叠加至灰度拉伸后的睑板腺图像得到最终增强图;
其中,第一截止频率小于第二截止频率,第二截止频率小于第三截止频率。
[0014] 本发明的有益之处在于所提供的睑板腺图像处理方法首先计算出一个科学的分割阈值,通过该分割阈值对睑板腺图像进行二值化,再对二值图进行闭运算,区域选取以及区域填充等操作得到睑板腺区域,找出睑板腺区域的边界后将边界映射至睑板腺图像,快速将睑板腺图像中的睑板腺区域分割出来。
[0015] 本发明的有益之处还在于所提供的睑板腺图像处理方法使用最大类间方法计算降噪处理后的睑板腺图像的灰度阈值,再统计出降噪处理后的睑板腺图像的灰度值得到直方分布图,最后根据计算公式计算分割阈值,经过上述步骤计算得到的分割阈值更加科学。同时,以该分割阈值为中心,将降噪处理后的睑板腺图像中灰度值在0.9Th-1.1Th之间的区域的灰度值置为255且其他区域的灰度值置为0得到二值图,对二值图进行调整最后得到较为准确的睑板腺区域。
[0016] 本发明的有益之处还在于所提供的睑板腺图像处理方法根据分割出的睑板腺区域的灰度值范围对图像进行拉伸,能够避免全局灰度拉伸造成的拉伸效果不明显等问题。附图说明
[0017] 图1是本发明的睑板腺图像处理方法的示意图;图2是本发明的图像采集装置的显示界面的示意图;
图3是本发明的灰度转换后的睑板腺图像的示意图;
图4是本发明的睑板腺图像的灰度值统计直方分布图示意图;
图5是本发明的二值图的示意图;
图6是对图5中的二值图进行闭运算区域选取以及填充后得到的二值图的示意图;
图7是对图6进行边界提取得到的边界的示意图;
图8是将边界映射至睑板腺图像的示意图;
图9是本发明的伽马变换的示意图;
图10是最终增强后的睑板腺图像的示意图。

具体实施方式

[0018] 以下结合附图和具体实施例对本发明作具体的介绍。如图1所示,本发明揭示一种睑板腺图像处理方法,用于提取睑板腺图像中的睑板腺区域,包含以下步骤:S1:获取睑板腺图像。S2:将睑板腺图像进行灰度化。S3:对灰度化后的睑板腺图像进行降噪处理。S4:计算分割阈值。S5:根据分割阈值对降噪处理后的睑板腺图像进行二值化处理得到二值图。S6:对二值图进行闭运算后保留最大区域并对最大区域进行填充。S7:提取二值图中的最大区域的边界。S8:将边界映射到降噪处理后的睑板腺图像。通过上述方法,根据计算出的分割阈值对图像进行二值化,再经过闭运算和图像填充,能够快速得到睑板腺图像中的睑板腺区域。以下具体介绍上述步骤。
[0019] 对于步骤S1:获取睑板腺图像。
[0020] 具体而言,首先通过照明装置照明,将图像采集装置的采集头对准被检测对象的眼部,从图像采集装置的成像显示器中观察采集头的图像采集范围,在本发明中,通过计算机软件在显示界面中设置粗选框,调整图像采集装置的采集头的位置以将睑板腺区域完全包含在粗选框内后采集图像得到睑板腺图像,在采集到图像的同时也完成了睑板腺区域的粗分割,如图2所示,显示界面中的方框为粗选框。
[0021] 其中,照明装置为红外照明装置,照射红外光,获取到的睑板腺图像为红外图像。
[0022] 对于步骤S2:将睑板腺图像进行灰度化。
[0023] 将获取到的睑板腺图像转换为灰度图,如图3所示。
[0024] 对于步骤S3:对灰度化后的睑板腺图像进行降噪处理。
[0025] 通过降噪手段对灰度化后的睑板腺图像进行降噪处理。在本发明中,对灰度化后的睑板腺图像进行均值滤波,具体的,使用3*3的模板对图像进行均值滤波去除噪声。
[0026] 对于步骤S4:计算分割阈值。
[0027] 具体而言,首先使用最大类间方法计算降噪处理后的睑板腺图像的灰度阈值,统计降噪处理后的睑板腺图像的灰度值得到直方分布图,如图4所示,根据下述公式(1)计算分割阈值,,                                            (1)
其中,Th为分割阈值,a为灰度阈值,b为255,i为灰度值,f(i)为直方分布图中灰度值为i的像素的个数。
[0028] 对于步骤S5:根据分割阈值对降噪处理后的睑板腺图像进行二值化处理得到二值图。
[0029] 具体而言,将降噪处理后的睑板腺图像中灰度值在0.9Th-1.1Th之间的区域的灰度值置为255且其他区域的灰度值置为0得到二值图。如图5所示。
[0030] 对于步骤S6:对二值图进行闭运算后保留最大区域并对最大区域进行填充。
[0031] 首先对二值图进行闭运算,使各区域边界更为光滑。此时,二值图中包含多个区域,计算每个区域的面积,保留面积最大的区域,再对该最大区域进行填充,如图6所示。
[0032] 对于步骤S7:提取二值图中的最大区域的边界。
[0033] 具体的,在本发明中,通过canny算子提取二值图中的最大区域的边界。如图7所示。
[0034] 对于步骤S8:将边界映射到降噪处理后的睑板腺图像。
[0035] 在本发明中,将提取出的边界映射到降噪处理后的睑板腺图像中如图8所示,边界内的区域为提取出的睑板腺区域。
[0036] 作为一种优选的实施方式,在将边界映射到降噪处理后的睑板腺图像之后,睑板腺图像处理方法还包含步骤:通过边界内的睑板腺区域的灰度值范围对全图进行灰度拉伸。对灰度拉伸后的睑板腺图像进行图像增强。经过上述步骤,能够使睑板腺区域中的腺体更为明显且便于观察。
[0037] 进一步地,通过边界内的睑板腺区域的灰度值范围对全图进行灰度拉伸的具体方法为:获取边界内的睑板腺区域的最小灰度值和最大灰度值。将降噪处理后的睑板腺图像中的灰度值小于最小灰度值的区域的灰度值置为0。将降噪处理后的睑板腺图像中的灰度值大于最大灰度值的区域的灰度值置为255。将降噪处理后的睑板腺图像中的灰度值介于最小灰度值和最灰度值范围内的区域通过下述伽马变换公式(2)进行灰度转换,,                                                     (2)其中,s为变换后的灰度值,r为变换前的灰度值,c为变换常数,γ为伽马系数。
[0038] 如图9所示,对于不同的γ值,伽马变换有不同的拉伸曲线,因为需要增强亮的区域且降低暗的区域,所以选择一个大于1的γ值进行变换。γ的取值范围优选为大于等于1.2且小于等于2。在本发明中,γ取值为1.5。
[0039] 而对灰度拉伸后的睑板腺图像进行图像增强的具体方法为:对灰度拉伸后的睑板腺图像经过高斯高通滤波器滤波,设定第一截止频率,得到第一高频图后将第一高频图叠加至灰度拉伸后的睑板腺图像得到第一增强图。对第一增强图经过高斯高通滤波器滤波,设定第二截止频率,得到第二高频图后将第二高频图叠加至灰度拉伸后的睑板腺图像得到第二增强图。对第二增强图经过高斯高通滤波器滤波,设定第三截止频率,得到第三高频图后将第三高频图叠加至灰度拉伸后的睑板腺图像得到最终增强图,如图10所示。其中,第一截止频率小于第二截止频率,第二截止频率小于第三截止频率。在本发明中,第一截止频率为8Hz,第二截止频率为16Hz,第三截止频率为32Hz。
[0040] 以上显示和描述了本发明的基本原理、主要特征和优点。本行业的技术人员应该了解,上述实施例不以任何形式限制本发明,凡采用等同替换或等效变换的方式所获得的技术方案,均落在本发明的保护范围内。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈