首页 > 专利 > 安庆师范大学 > 粗精调焦的医学显微图像清晰度评价方法专利详情

粗精调焦的医学显微图像清晰度评价方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2016-10-18
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2017-04-26
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2019-04-05
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2036-10-18
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN201610922544.2 申请日 2016-10-18
公开/公告号 CN106548471B 公开/公告日 2019-04-05
授权日 2019-04-05 预估到期日 2036-10-18
申请年 2016年 公开/公告年 2019年
缴费截止日
分类号 G06T7/00 主分类号 G06T7/00
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 1
权利要求数量 2 非专利引证数量 1
引用专利数量 3 被引证专利数量 0
非专利引证 1、Gergely Windisch and MiklósKozlovszky.Evaluating Sharpness Metricsfor HE Stained《.CINTI 2014 •15th IEEEInternational Symposium on ComputationalIntelligence and Informatics》.2014,;
引用专利 CN101902366A、CN102496157A、CN103985108A 被引证专利
专利权维持 2 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 安庆师范大学 当前专利权人 安庆师范大学
发明人 郑馨、程玉胜、江伟 第一发明人 郑馨
地址 安徽省安庆市集贤北路1318号安庆师范大学 邮编 246133
申请人数量 1 发明人数量 3
申请人所在省 安徽省 申请人所在市 安徽省安庆市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
摘要
本发明公开了一种粗精调焦的医学显微图像清晰度评价方法,先将输入的RGB彩色图像转换为Lab颜色空间;接着通过统计亮度图IL的灰度方差获取灰度变化值,利用颜色通道图Ia和颜色通道图Ib获取颜色显著图,并得到颜色显著图的局部梯度值;最后,将灰度变化值和局部梯度变化值相乘,获取最终的图像清晰度。本发明提出的方法综合考虑图像的全局灰度变化和局部梯度变化,既具有较宽的陡峭区又具有较高的灵敏度,仅使用一种清晰度评价函数就可以同时用于大步长粗调焦和小步长精细调焦。
  • 摘要附图
    粗精调焦的医学显微图像清晰度评价方法
  • 说明书附图:其中,
    粗精调焦的医学显微图像清晰度评价方法
  • 说明书附图:其中,-1
    粗精调焦的医学显微图像清晰度评价方法
  • 说明书附图:[0012]
    粗精调焦的医学显微图像清晰度评价方法
  • 说明书附图:[0025]
    粗精调焦的医学显微图像清晰度评价方法
  • 说明书附图:[0026]
    粗精调焦的医学显微图像清晰度评价方法
  • 说明书附图:[0032]
    粗精调焦的医学显微图像清晰度评价方法
  • 说明书附图:[0035]
    粗精调焦的医学显微图像清晰度评价方法
  • 说明书附图:[0042]
    粗精调焦的医学显微图像清晰度评价方法
  • 说明书附图:[0047]
    粗精调焦的医学显微图像清晰度评价方法
  • 说明书附图:[0048]
    粗精调焦的医学显微图像清晰度评价方法
  • 说明书附图:[0054]
    粗精调焦的医学显微图像清晰度评价方法
  • 说明书附图:[0057]
    粗精调焦的医学显微图像清晰度评价方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2019-04-05 授权
2 2017-04-26 实质审查的生效 IPC(主分类): G06T 7/00 专利申请号: 201610922544.2 申请日: 2016.10.18
3 2017-03-29 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
1.一种粗精调焦的医学显微图像清晰度评价方法,其特征在于,所述方法包括以下几个步骤:
步骤一、输入RGB格式的原始图像,并将其转换为Lab颜色空间,分别得到亮度图IL和两个颜色通道图:颜色通道图Ia和颜色通道图Ib;
步骤二、通过以下函数获取原始图像的灰度变化值FVar:
其中, 为IL的灰度均值,(x,y)为像素的位置;
步骤三、获取原始图像的局部梯度值,具体过程如下:
(3.1)分别对颜色通道图Ia和颜色通道图Ib中每个像素的灰度进行去均值和归一化操作,再将操作结果线性叠加得到颜色显著图S;
(3.2)对以上步骤所得的颜色显著图S中每个像素计算局部梯度值,具体地,统计当前像素横向、纵向和两个斜方向上的像素灰度差的绝对值之积;
G(x,y)=[|S(x+2,y)-S(x,y)|×|S(x,y+2)-S(x,y)|]×[|S(x+1,y+1)-S(x,y)|×|S(x-1,y+1)-S(x,y)|]
(3.3)仅统计大于门限T的函数值之和作为输入图像的局部梯度值:
FGrad=∑[G(x,y)>T]
本步骤中,T设定为10;
步骤四、将图像灰度变化值FVar和局部梯度值FGrad相乘,获取最终的图像清晰度F。

2.如权利要求1所述的粗精调焦的医学显微图像清晰度评价方法,其特征在于,所述的步骤(3.1)具体包括:
(3.1.1)使用以下公式进行去均值操作:
其中, 和 分别为颜色通道图Ia和颜色通道图Ib的灰度均值,||为绝对值符号,I′a和I′b分别为颜色通道图Ia和颜色通道图Ib去均值操作后的图像;
(3.1.2)分别对I′a和I′b进行归一化操作,得到归一化图像Na和归一化图像Nb;
Na(x,y)=N(I′a(x,y))
Nb(x,y)=N(I′b(x,y))
其中,N()为归一化函数,其公式如下:
其中,max()和min()分别为最大值和最小值函数;
(3.1.3)通过以下公式将归一化图像Na和归一化图像Nb线性叠加得到颜色显著图S:
(x,y)为图像中像素的位置。
说明书

技术领域

[0001] 本发明属于图像处理技术领域,更具体地,涉及一种医学显微图像清晰度评价方法。

背景技术

[0002] 医学显微图像广泛应用于血细胞分析、染色体分析、尿沉渣分析等医学研究和临床医疗诊断中。通过自动聚焦获取清晰的显微图像是后续医学图像处理、目标识别与病理分析的前提。显微镜通过光学系统的作用将肉眼不可见的观察对象放大,从微观形态研究和认识对象的特性,因此,显微镜自动聚焦是非常精密复杂的过程。
[0003] 现有的显微镜自动聚焦方法通常采用基于图像反馈的被动式聚焦方式,即利用清晰度评价函数对不同离焦位置采集到的图像进行清晰度评价,以此找到具有最佳清晰度的聚焦位置。因此,自动聚焦的关键在于图像清晰度评价函数的设计。一个理想的医学显微图像清晰度评价函数除了应当具备无偏性、单峰性和抗噪能力之外,还应当能够同时实现快速大步长粗调焦和小步长精细调焦。
[0004] 显微镜自动聚焦过程一般分为粗调焦和精细调焦两个阶段。在粗调焦阶段,显微镜载物台的位置往往远离聚焦位置,为了快速靠近聚焦位置,通常需要采用大步长粗调焦。为了保证能够到达聚焦位置,此时的清晰度评价函数必须具有足够宽的陡峭区来避免因步长过大越过陡峭区而导致的聚焦失败。在精细调焦阶段,一般需要采用小步长获取精确的聚焦位置,此时的清晰度评价函数应当具备足够高的灵敏度才能获取最精确的聚焦位置。
[0005] 目前,应用最广泛的清晰度评价函数主要有两类:梯度函数和统计函数。基于梯度函数的清晰度评价函数,如灰度差分绝对值之和(SMD)、Brenner函数、图像灰度梯度能量函数(Tenengrad)、拉普拉斯能量(EOL)等函数,它们的灵敏度好,但陡峭区过窄,不适用于大范围粗调焦;灰度方差等统计函数的计算量小,拥有较宽的陡峭区,且在平缓区域的单调性好,然而,其灵敏度差,不适用于小范围精调焦。为了解决该问题,一些粗精结合的图像自动聚焦方法在粗、细两个调焦阶段分别使用不同的清晰度评价函数。然而,这些聚焦方法的性能会受到聚焦策略的影响,如粗调焦阶段的步长的选取、搜索策略等。因此,目前尚没有一种清晰度评价函数能够同时实现快速大步长粗调焦和小步长精细调焦,即,既具备足够宽的陡峭区又拥有足够高的灵敏度。

发明内容

[0006] 针对现有技术的以上缺陷或改进需求,本发明提供了一种粗精调焦的医学显微图像清晰度评价方法,其目的在于设计一种全面评价图像清晰度的函数,使其既具备足够宽的陡峭区又拥有足够高的灵敏度,无论在粗调焦阶段还是在精细调焦阶段中都能表现出较好的性能。
[0007] 本发明提出的技术方案如下:
[0008] 一种粗精调焦的医学显微图像清晰度评价方法,其特征在于,所述方法包括以下几个步骤:
[0009] 步骤一、输入RGB格式的原始图像,并将其转换为Lab颜色空间,分别得到亮度图IL和两个颜色通道图:颜色通道图Ia和颜色通道图Ib;
[0010] 步骤二、通过以下函数获取原始图像的灰度变化值FVar:
[0011]
[0012] 其中, 为IL的灰度均值,(x,y)为像素的位置;
[0013] 步骤三、获取原始图像的局部梯度值,具体过程如下:
[0014] (3.1)分别对颜色通道图Ia和颜色通道图Ib中每个像素的灰度进行去均值和归一化操作,再将操作结果线性叠加得到颜色显著图S;
[0015] (3.2)对以上步骤所得的颜色显著图S中每个像素计算局部梯度值,具体地,统计当前像素横向、纵向和两个斜方向上的像素灰度差的绝对值之积;
[0016] G(x,y)=[|S(x+2,y)-S(x,y)|×|S(x,y+2)-S(x,y)|]×
[0017] [|S(x+1,y+1)-S(x,y)|×|S(x-1,y+1)-S(x,y)|]
[0018] (3.3)仅统计大于门限T的函数值之和作为输入图像的局部梯度值:
[0019] FGrad=∑[G(x,y)>T]
[0020] 本步骤中,优选地,T设定为10;
[0021] 步骤四、将图像灰度变化值FVar和局部梯度值FGrad相乘,获取最终的图像清晰度F。
[0022] 优选的,所述的步骤(3.1)具体包括:
[0023] (3.1.1)使用以下公式进行去均值操作:
[0024]
[0025]
[0026] 其中, 和 分别为颜色通道图Ia和颜色通道图Ib的灰度均值,||为绝对值符号,I′a和I′b分别为颜色通道图Ia和颜色通道图Ib去均值操作后的图像;
[0027] (3.1.2)分别对I′a和I′b进行归一化操作,得到归一化图像Na和归一化图像Nb。
[0028] Na(x,y)=N(I′a(x,y))
[0029] Nb(x,y)=N(I′b(x,y))
[0030] 其中,N()为归一化函数,其公式如下:
[0031]
[0032] 其中,max()和min()分别为最大值和最小值函数;
[0033] (3.1.3)通过以下公式将归一化图像Na和归一化图像Nb线性叠加得到颜色显著图S:
[0034]
[0035] (x,y)为图像中像素的位置。
[0036] 本发明能够达到的有益效果如下:
[0037] 医学显微图像在粗、细两个不同的调焦阶段对清晰度评价函数又不同的要求,即在大步长粗调焦阶段需要具有较宽的陡峭区,在小步长精细调焦阶段需要具有较高的灵敏度。本发明提出的方法综合考虑了图像的全局灰度变化和局部梯度变化,仅使用一种清晰度评价函数就可以同时用于大步长粗调焦和小步长精细调焦。

实施方案

[0038] 为了使本发明的目的、技术方案及优点更加清楚明白,以下对本发明的具体实施方式进行进一步详细说明。
[0039] 步骤一、输入RGB格式的原始图像,并将其转换为Lab颜色空间,分别得到亮度图IL和两个颜色通道图:颜色通道图Ia和颜色通道图Ib;
[0040] 步骤二、获取原始图像的灰度变化值,具体地,通过以下函数计算亮度图IL的灰度方差FVar:
[0041]
[0042] 其中, 为IL的灰度均值,(x,y)为像素的位置。
[0043] 步骤三、获取原始图像的局部梯度值,具体过程如下:
[0044] (3.1)分别对颜色通道图Ia和颜色通道图Ib中每个像素的灰度进行去均值和归一化操作,再将操作结果线性叠加得到颜色显著图S。
[0045] (3.1.1)去均值操作是为了消除医学显微图像中的大面积的背景干扰,其公式如下:
[0046]
[0047]
[0048] 其中, 和 分别为颜色通道图Ia和颜色通道图Ib的灰度均值,||为绝对值符号,I′a和I′b分别为颜色通道图Ia和颜色通道图Ib去均值操作后的图像。当医学显微图像中背景区域面积远大于目标区域面积时,去均值操作后的颜色通道图像中,背景区域将接近于零,而目标区域因与背景区域具有显著的颜色差异而远大于零。因此,去均值操作有利于消除背景区域对聚焦的干扰。
[0049] (3.1.2)分别对I′a和I′b进行归一化操作,得到归一化图像Na和归一化图像Nb。
[0050] Na(x,y)=N(I′a(x,y))
[0051] Nb(x,y)=N(I′b(x,y))
[0052] 其中,N()为归一化函数,其公式如下:
[0053]
[0054] 其中,max()和min()分别为最大值和最小值函数。归一化后的图像Na和Nb的像素灰度区间将扩展至[0,255],使目标区域与背景区域的颜色对比度增大,从而进一步消除背景干扰。
[0055] (3.1.3)通过以下公式将归一化图像Na和归一化图像Nb线性叠加得到颜色显著图S:
[0056]
[0057] (3.2)对以上步骤所得的颜色显著图中每个像素计算局部梯度值,具体地,统计当前像素横向、纵向和两个斜方向上的像素灰度差的绝对值之积;
[0058] G(x,y)=[|S(x+2,y)-S(x,y)|×|S(x,y+2)-S(x,y)|]×
[0059] [|S(x+1,y+1)-S(x,y)|×|S(x-1,y+1)-S(x,y)|]
[0060] 其中,|S(x+2,y)-S(x,y)|和|S(x,y+2)-S(x,y)|分别用于提取水平方向和垂直方向的梯度值,|S(x+1,y+1)-S(x,y)|和|S(x-1,y+1)-S(x,y)|分别用于提取斜向和反斜向的梯度值。因此,G(x,y)综合了四个方向的梯度值,能更全面反映每个像素的局部梯度变化。
[0061] (3.3)为了减少背景噪声的干扰,仅统计大于门限T的函数值之和作为输入图像的局部梯度值:
[0062] FGrad=∑[G(x,y)>T]
[0063] 本步骤中,优选地,T设定为10。
[0064] 步骤四、将图像灰度变化值FVar和局部梯度值FGrad相乘,得到最终的图像清晰度F:
[0065] F=FVar×FGrad
[0066] 当离聚焦位置较远时,图像的局部梯度变化不明显,呈现轻微波动状态;而全局灰度变化较为明显,呈明显的单调性。此时F的变化率由FVar主导,因此,F的陡峭区宽度接近于FVar的陡峭区宽度。当距离聚焦位置越来越近时,图像的边缘、纹理等细节信息逐渐丰富,局部梯度变化远大于全局灰度变化,此时,F的变化率主要由FGrad决定,且F的变化率大于FGrad,因此F的灵敏度更好。
[0067] 本发明的粗精调焦的医学显微图像清晰度评价方法,先将输入的RGB彩色图像转换为Lab颜色空间;接着通过统计亮度图L的灰度方差获取灰度变化值,利用颜色通道图a和颜色通道图b获取显著图,并得到显著图的局部梯度值;最后,将灰度变化值和局部梯度值相乘,获取最终的图像清晰度。本发明提出的方法综合考虑图像的全局灰度变化和局部梯度变化,既具有较宽的陡峭区又具有较高的灵敏度,仅使用一种清晰度评价函数就可以同时用于大步长粗调焦和小步长精细调焦。
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号