一种眼底图像微动脉瘤自动检测方法
【技术领域】
[0001] 本发明涉及眼底图像微动脉瘤自动检测方法。
【背景技术】
[0002] 糖尿病患者的血糖长期高于正常值，导致视网膜毛细血管的形状和结构发生改 变，毛细血管壁膨大，便会形成眼底视网膜上的微动脉瘤（Microaneurysm，MA)。微动脉瘤是 糖尿病性视网膜病变最早出现且能被观察到的症状之一，在眼底图像上显示为极小的暗红 色点，在绿色通道上则显示为暗色的点。微动脉瘤的直径大小约为10-100 ym, 10-20 μπι的 最为常见，在正常尺寸的眼底图像中直径约为1像素，因此肉眼很难发现。眼底微动脉瘤极 小，难以通过医生肉眼发现，只有受过多年专业训练的眼科专家医生才能凭经验做出判断。 这就给眼底微动脉瘤的筛查带来不便，检测效率低下，同时使得误诊的概率大大增加，而且 不利于医疗条件欠佳的地区展开大规模眼底微动脉瘤筛查。
[0003] 解决方法便是通过眼底图像的处理自动检测微动脉瘤，目前提出的方法多是基于 微动脉瘤的形状、颜色等特点，通过形态学、匹配滤波等技术提取出疑似微动脉瘤的区域。 再进一步通过微动脉瘤的大小、对比度、形状等相关特征对疑似区域进行排除，去除非微动 脉瘤的出血等病变点，选取出真实的微动脉瘤。但实际操作中。眼底图像往往受到患者个 体差异以及环境光源的影响而导致亮度不均匀，不便于之后的检测；微动脉瘤的尺寸也非 常小，尤其是在发病早期，直径很小难以判断；图像中的血管与微动脉瘤亮度相似，也是重 要的干扰因素。基于这些原因，类似算法会检测出大量的疑似区域，在排除疑似时大幅增加 了系统复杂度。因此，目前的微动脉瘤自动检测还会存在漏判和误判的问题。
【发明内容】
[0004] 本发明所要解决的技术问题就是提供一种眼底图像微动脉瘤自动检测方法，能够 最大限度地获得完整的血管和微动脉瘤结构，避免漏判和误判。
[0005] 为解决上述技术问题，本发明采用如下技术方案：一种眼底图像微动脉瘤自动检 测方法，首先，使用高斯匹配滤波技术增强眼底图像中的血管和微动脉瘤像素；然后，利用 局部熵确定阈值，从而将血管和微动脉瘤对象从眼底图像中分割出来；最后，利用长度筛选 的方法，将血管排除在疑似区域之外，保留微动脉瘤像素，从而实现了对眼底图像中微动脉 瘤的分割和检测。
[0006] 优选的，在检测前对眼底图像进行预处理，通过图像平滑处理减少图像的噪声。
[0007] 优选的，使用二维高斯匹配滤波增强血管和微动脉瘤像素，其采用的函数定义如
[0008] 其中L为片段长度，在这个片段内血管方向固定不变，视网膜血管的走向不固定， 所以核函数需要不断地旋转，只有当滤波器方向与血管走向相同才会产生峰值响应，然后 利用一组12个16X15像素的内核对眼底图像进行卷积，在每个像素点上只保留最大的响 应，得到响应图像，σ 2是高斯函数拟合二项式系数的最小方差。
[0009] 优选的，设阈值为Τ，0彡T彡L-1，将滤波处理后的图像LXL灰度共生矩阵分为 A、B、C、D四个部分，A部分即灰度值较低的取值，为图像血管和微动脉瘤的目标元素的灰度 分布；而C部分即灰度值较高的取值，为图像中背景的灰度分布；B、D两部分则是目标元素 与背景之间的像素点或元素边界点的灰度分布；由局部熵确定阈值时只需要考虑这四个部 分中的A、C两个部分;A部分目标元素的整体灰度分布概率可定义为：
[0011] 而C部分背景元素的概率则定义为：
[0013] 对这两个概率进行归一化，得到以下两个关于向量（T，T)的函数：
[0016] 则A部分目标元素的局部熵是一个关于T的函数，可由下式得出：
[0020] 则目标元素和背景元素两部分的总局部熵为：
[0021] Ht(T) =Ha (T)+Hc (T)
[0022] 则由此得出阈值T的取值应使Ht取到最大值，此时有：
[0024] 此时得到理想的阈值T，以进行目标区域和背景区域的分类；其中L是整个图像的 灰度级数，Plj是i和j两个灰度值共生的概率，即LXL灰度共生矩阵中的（i，j)元素。
[0025] 整个图像的灰度级数为L，则共有LXL种数对组合。相邻两个像素的灰度值数对 (i，j)在整个图像的所有相邻像素对的灰度数对中的出现概率记为P,，即是矩阵的（i，j) 元素。
[0026] 优选的，长度筛选方法以如下方式操作：对于每个目标区域像素点，比较其周围相 邻的八个像素点，相邻的八个像素中，同样处在目标区域之内的像素少于两个，便可视为不 连通的单像素小点，即为所需的真实微动脉瘤点。
[0027] 本发明采用了基于高斯匹配滤波和局部熵确定阈值实现图像分割的方法，能够最 大限度地获得完整的血管和微动脉瘤结构，避免漏判和误判。
【附图说明】
[0028] 下面结合附图和【具体实施方式】对本发明作进一步描述：
[0029] 图1为灰度共生矩阵的结构示意图；
[0030] 图2为对图像进行分割之后产生的二值图像；
[0031] 图3为长度筛选后确定的微动脉瘤位置。
【具体实施方式】
[0032] 本发明眼底图像微动脉瘤自动检测的算法分为三个步骤：首先对图像进行匹配滤 波以增强图像中的血管和微动脉瘤像素；之后利用局部熵确定阈值，分割出图像的血管、微 动脉瘤等目标区域，这种方法能反映图像的灰度空间分布，保留二值图像的空间结构；最后 加入长度筛选的步骤，在筛选出的暗色区域中去除掉长度较长且相互连通的血管结构，即 可检测出真实的微动脉瘤结构。
[0033] 首先，对基于匹配滤波的血管增强方法做出说明。
[0034] 眼底图像中，如果背景较暗，则血管、微动脉瘤和出血等区域在图中与暗色背 景区分不明显。图像中的背景纹理、噪声，以及黄斑、视盘等区域也会对微动脉瘤检测 造成干扰，降低检测的准确率。为了增强血管和微动脉瘤等目标区域，需要对图像进 行匹配滤波处理。血管和微动脉瘤截面的灰度分布都可以用一条高斯分布曲线来近 似。在这里使用二维高斯匹配滤波可以增强血管和微动脉瘤像素。函数定义如下：
[0035] 其中L为片段长度，在这个片段内血管方向固定不变，视网