placido盘图像中特征点的提取方法及系统

1.本发明涉及图像处理技术领域，特别涉及一种placido盘图像中特征点的提取方法及系统。


背景技术：

2.角膜地形是借用了地形勘察中的概念用来反映角膜表面的形态。角膜地形图就是对整个角膜进行分析，使整个角膜的曲率在角膜地形图上客观的展现出来。在临床上，角膜地形图用于诊断角膜散光和圆锥角膜，使得早期圆锥角膜的诊断成为可能，角膜地形图问世以来，圆锥角膜早期诊断准确率高达96％。此外，角膜地形图还可以用于角膜屈光手术的术前检查和术后疗效评价，用于指导临床手术。
3.而角膜地形图是根据计算机处理过的伪彩色图像，在地形图中，使用等高线的形式来表示角膜形状的变化，可以在二维图像中表示三维的信息。进行角膜地形图绘制的第一步就是进行角膜placido图像的采集，并进行特征点的提取。根据placido环计算还原的眼角膜的曲率半径，进而计算屈光度达成绘制角膜地形图的目的。
4.如申请号为cn108734701a的发明专利中，提出了一种placido环图像形态变化识别算法，用于进行placido环图像的处理，但是其方法复杂，步骤繁琐，运算速度有待提高。再如专利cn108734701a公开了一种圆形拟合的方法去计算圆心，但其运算量复杂并且不具有足够的鲁棒性。
5.所以现在有必要提供一种可靠的提取placido盘图像中特征点的方案。


技术实现要素：

6.本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种placido盘图像中特征点的提取方法及系统。
7.为解决上述技术问题，本发明采用的技术方案是：一种placido盘图像中特征点的提取方法，包括以下步骤：
8.1)获取placido图像；
9.2)对获取的placido图像进行预处理；
10.3)计算placido图像的中心；
11.4)利用极坐标展开placido图像，并进行形态学处理和高斯滤波处理；
12.5)对步骤4)的结果分别进行顶帽变换和黑帽变换，然后将得到的两个结果相减；
13.6)对步骤5)的结果进行笛卡尔坐标转换，然后利用方差法得到角膜区域；
14.7)对步骤6)的结果进极坐标化，然后进行形态学处理、直方图均衡处理和自适应二值化处理；
15.8)对步骤7)的结果进行边缘提取，然后再进行笛卡尔坐标系化，得到含有placido圆环的图像；
16.9)在步骤8)得到的图像中的每个placido圆环上进行特征点的提取。
17.优选的是，所述步骤1)的具体方法为：在placido盘的中心两侧设置四个圆形的小灯，placido盘的中心处于该四个圆形的小灯所围绕形成的矩形框内，然后利用该placido盘进行角膜placido图像的采集。
18.优选的是，所述步骤2)中进行预处理的步骤具体包括：对placido图像进行归一化处理，然后使用进行高斯滤波处理。
19.优选的是，所述步骤2)中进行归一化处理的方法为：先寻找到placido图中的像素最大值和像素最小值，将二者做差，得到两者的差值res，然后将每个像素点像素值的和像素最小值的差去除以res，取整数，这样将整个图像的像素值都映射到了0和1之间。
20.优选的是，所述所述步骤2)中使用10*10的高斯核进行高斯滤波处理。
21.优选的是，所述步骤3)具体包括：
22.1)对于整个placido图像进行灰度化；
23.2)对灰度化之后的placido图像进行roi区域提取：灰度化后的placido图像中，四个圆形的小灯会生成四个过曝的亮斑，将四个过曝的亮斑连接形成的矩形框作为roi区域，placido盘的中心处于该roi区域内；
24.3)对于roi区域进行全局的质心提取，提取到的质心即为placido图像的中心。
25.优选的是，所述步骤7)中形态学处理的具体方法为：使用10*10的椭圆核进行开操作。
26.本发明还提供一种placido盘图像中特征点的提取系统，其采用如上所述的方法进行placido盘图像中特征点的提取。
27.本发明还提供一种存储介质，其上存储有计算机程序，其特征在于，该程序被执行时用于实现如上所述的方法。
28.本发明还提供一种计算机设备，包括存储器、处理器以及存储在所述存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权上所述的方法。
29.本发明的有益效果是：
30.本发明中通过利用placido盘上的小灯辅助进行placido盘图像中心的定位，能减小计算量、提高运算速度和准确度，图像处理速度可达到每秒钟30fps；提高了placido环上特征点的提取速度和效率，定位到的特征点可以用于角膜地形图的绘制。
31.本发明中对placido图像上的睫毛区域进行了智能抹去，有效避免了因为睫毛的干扰而造成的误识别。
附图说明
32.图1为本发明的placido盘图像中特征点的提取方法的流程图；
33.图2为本发明的实施例1中的原始placido图像；
34.图3
‑
8为本发明的实施例1中placido图像的处理过程图；
35.图9为本发明的实施例1中的最终得到的placido图像。
具体实施方式
36.下面结合实施例对本发明做进一步的详细说明，以令本领域技术人员参照说明书
文字能够据以实施。
37.应当理解，本文所使用的诸如“具有”、“包含”以及“包括”术语并不排除一个或多个其它元件或其组合的存在或添加。
38.实施例1
39.参照图1，本实施例提供的一种placido盘图像中特征点的提取方法，包括以下步骤：
40.s1、获取placido图像。
41.具体为：在placido盘的中心两侧设置四个圆形的小灯，placido盘的中心处于该四个圆形的小灯所围绕形成的矩形框内，然后利用该placido盘进行角膜placido图像的采集。获取的原始placido图像如图2所示，四个圆形的小灯在成像后会在placido图像上生成四个过曝的亮斑，通过这四个亮斑可以辅助进行后续的placido图像中心的定位。
42.s2、对获取的placido图像进行预处理。
43.具体包括：
44.对placido图像进行归一化处理：先寻找到placido图中的像素最大值和像素最小值，将二者做差，得到两者的差值res，然后将每个像素点像素值的和像素最小值的差去除以res，取整数，这样将整个图像的像素值都映射到了0和1之间，可以有效避免图像仿射变换带来的噪声，便于后续的处理。
45.然后使用进行高斯滤波处理：使用10*10的高斯核进行高斯滤波处理，从而进行全图的平滑。
46.s3、计算placido图像的中心。
47.具体包括：
48.1)对于整个placido图像进行灰度化；
49.2)对灰度化之后的placido图像进行roi区域提取：灰度化后的placido图像中，四个圆形的小灯会生成四个过曝的亮斑，将四个过曝的亮斑连接形成的矩形框作为roi区域，placido盘的中心处于该roi区域内；
50.3)对于roi区域进行全局的质心提取，提取到的质心即为placido图像的中心。
51.通过四个过曝的亮斑作为辅助，再计算placido盘的中心，能够显著减小计算量，提高运算速度和定位准确度。
52.s4、利用极坐标展开placido图像，如图3；再进行形态学处理和高斯滤波处理，如图4。
53.s5、对步骤4)的结果分别进行顶帽变换和黑帽变换，然后将得到的两个结果相减，从而能够很好的去除睫毛的影响，有效避免了因为睫毛的干扰而造成的误识别，如图5。
54.s6、对步骤5)的结果进行笛卡尔坐标转换，然后利用方差法得到角膜区域，如图6。使用方差法可以得到含有placido环映射的roi_placido区域，因为在存在placido环的区域和不存在placido环的区域，图像方差有着明显的差距，所以通过方差的差别能精准定位出roi_placido区域。
55.s7、对步骤6)的结果进极坐标化，然后进行形态学处理(使用10*10的椭圆核进行开操作，用于平滑表面)、直方图均衡处理和自适应二值化处理，如图7。
56.s8、对步骤7)的结果进行边缘提取，然后再进行笛卡尔坐标系化，得到含有
placido圆环的图像，如图8。
57.s9、在步骤8)得到的图像中的每个placido圆环上进行特征点的提取，如图9。定位到的特征点可以用于角膜地形图的绘制，例如通过采取插值的方式进行角膜地形图的绘制可以得到完整的角膜地形图，在临床和圆锥角膜筛查方面有重要的意义。
58.实施例2
59.本实施例提供一种placido盘图像中特征点的提取系统，其采用实施例1的方法进行placido盘图像中特征点的提取。
60.本实施例还提供一种存储介质，其上存储有计算机程序，该程序被执行时用于实现实施例1的方法。
61.本实施例还提供一种计算机设备，包括存储器、处理器以及存储在所述存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现实施例1的方法。
62.尽管本发明的实施方案已公开如上，但其并不仅仅限于说明书和实施方式中所列运用，它完全可以被适用于各种适合本发明的领域，对于熟悉本领域的人员而言，可容易地实现另外的修改，因此在不背离权利要求及等同范围所限定的一般概念下，本发明并不限于特定的细节。