留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

高精度角膜曲率半径测量系统

李华建 肖作江 刘颖 赵媛媛 王睿智 贺肖影

李华建, 肖作江, 刘颖, 赵媛媛, 王睿智, 贺肖影. 高精度角膜曲率半径测量系统[J]. 中国光学, 2020, 13(3): 501-509. doi: 10.3788/CO.2019-0174
引用本文: 李华建, 肖作江, 刘颖, 赵媛媛, 王睿智, 贺肖影. 高精度角膜曲率半径测量系统[J]. 中国光学, 2020, 13(3): 501-509. doi: 10.3788/CO.2019-0174
LI Hua-jian, XIAO Zuo-jiang, LIU Ying, ZHAO Yuan-yuan, WANG Rui-zhi, HE Xiao-ying. High precision corneal curvature radius measurement system[J]. Chinese Optics, 2020, 13(3): 501-509. doi: 10.3788/CO.2019-0174
Citation: LI Hua-jian, XIAO Zuo-jiang, LIU Ying, ZHAO Yuan-yuan, WANG Rui-zhi, HE Xiao-ying. High precision corneal curvature radius measurement system[J]. Chinese Optics, 2020, 13(3): 501-509. doi: 10.3788/CO.2019-0174

高精度角膜曲率半径测量系统

doi: 10.3788/CO.2019-0174
基金项目: 吉林省重点研发项目(No. 20180201025GX)
详细信息
    作者简介:

    李华建(1995—),男,江苏泰州人,硕士研究生,2018年于长春理工大学光电信息学院光电信息科学与工程专业获得学士学位,现为长春理工大学硕士研究生,主要研究方向为光学仪器设计、光电检测。E-mail:lhj1663126998@163.com

    肖作江(1973—),男,吉林长春人,博士,副研究员,硕士生导师,2003年于长春理工大学获得机械设计理论硕士学位,2009年于长春理工大学获得光学工程博士学位,主要研究方向为精密测量、光电检测。E-mail:xiao2706@163.com

    通讯作者: 肖作江(通讯作者)(1973—),吉林长春人,2003年获得长春理工大学机械设计理论硕士学位,2009年获得长春理工大学光学工程博士学位,现为长春理工大学副研究员、硕士生导师,主要研究方向为精密测量、光电检测。E-mail:xiao2706@163.com
  • 中图分类号: TP394.1;TH691.9

High precision corneal curvature radius measurement system

Funds: Supported by Jilin Province Key Research and Development Project (No. 20180201025GX)
More Information
  • 摘要: 为了实现成像角膜曲率计沿光轴方向的精确对准,提高角膜曲率的测量精度,设计了一种高精度成像角膜曲率测量系统。对该系统的成像光源、成像光学系统、干涉测量系统等进行研究。通过LED阵列均匀照射靶环形成光源;成像物镜采用双远心镜头,使景深变大,有利于对准测量,同时保证成像物镜的放大率不会受景深的影响;引用低相干干涉测量技术,利用低相干干涉信号定位角膜顶点和测量光源的位置,利用光栅尺监测扫描反射镜的位置,实现了角膜顶点到测量光源距离的精确测量。最后,分析了该系统成像物镜放大率的稳定性和角膜曲率测量误差,在理论基础上,做出了实验样机,用所设计的样机对标准角膜模拟眼进行测试,系统精度达到±0.02 mm,基本满足角膜曲率的测量要求。
  • 图  1  角膜曲率半径测量原理图

    Figure  1.  Measurement principle of corneal curvature radius

    图  2  误差分析图

    Figure  2.  Error analysis chart

    图  3  高精度角膜曲率半径测量系统原理图

    Figure  3.  Schematic diagram of measurement system for high precision corneal curvature radius

    图  4  靶环装置

    Figure  4.  Target ring device

    图  5  靶环装置成像原理图

    Figure  5.  Schematic diagram of imaging principle of target ring device

    图  6  靶环装置组成图

    Figure  6.  Composition diagram of target ring device

    图  7  靶环装置光源仿真图

    Figure  7.  Simulation diagram of light source for target ring device

    图  8  模拟面的辐照度

    Figure  8.  Irradiances of the simulated surface

    图  9  双远心镜头成像原理图

    Figure  9.  Imaging principle of double telecentric lens

    图  10  干涉测量系统

    Figure  10.  Interferometry measurement system

    图  11  干涉图

    Figure  11.  Interferogram

    图  12  $b$为100 μm时角膜曲率半径的误差

    Figure  12.  The error of corneal curvature radius when $b$ is 100 μm

    图  13  $b$为−100 μm时角膜曲率半径的误差

    Figure  13.  The error of corneal curvature radius when $b$ is −100 μm

    图  14  实验样机

    Figure  14.  Experimental prototype

    图  15  图像处理结果

    Figure  15.  Image processing results

    表  1  景深对应的$b$$\beta $

    Table  1.   The $b$ and $\beta $ values corresponding to different depthes of field

    景深/mm$b$/mm$\beta $/mm
    −367−0.464
    −268−0.459
    −169−0.456
    070−0.450
    171−0.446
    272−0.441
    373−0.437
    下载: 导出CSV

    表  2  双远心镜头参数

    Table  2.   Double telecentric lens parameters

    参数数值
    物距/mm100
    前景深/mm4.68
    后景深/mm5.47
    放大率–0.45
    下载: 导出CSV

    表  3  不同物距的放大率

    Table  3.   Magnifications corresponding to different object distances

    物距/mm放大率
    97−0.450 19
    98−0.450 15
    99−0.450 12
    100−0.450 00
    101−0.449 92
    102−0.449 86
    103−0.449 82
    104−0.449 78
    下载: 导出CSV

    表  4  实验结果

    Table  4.   Experimental results

    标准值测量值/mm误差/mm标准值测量值/mm误差/mm
    6.656${r_{\max }}$6.675﹢0.0197.988${r_{\max}}$7.9970.009
    ${r_{\min }}$6.641−0.015${r_{\min}}$7.975−0.013
    6.656${r_{\max}}$6.673+0.0179.458${r_{\max}}$9.466+0.008
    ${r_{\min }}$6.637−0.019${r_{\min }}$9.443−0.015
    6.656${r_{\max}}$6.669+0.0139.458${r_{\max}}$9.473+0.015
    ${r_{\min }}$6.647+0.011${r_{\min }}$9.455−0.003
    7.988${r_{\max}}$7.999+0.0119.458${r_{\max}}$9.477+0.019
    ${r_{\min }}$7.974−0.014${r_{\min }}$9.459+0.001
    7.988${r_{\max}}$7.998+0.010
    ${r_{\min }}$7.980−0.008
    下载: 导出CSV
  • [1] HOWLAND H C, HOWLAND B. Photorefraction: a technique for study of refractive state at a distance[J]. Journal of the Optical Society of America, 1974, 64(2): 240-249. doi:  10.1364/JOSA.64.000240
    [2] CHOI M, WEISS S, SCHAEFFEL F, et al. Laboratory, clinical, and kindergarten test of a new eccentric infrared photorefractor (Power Refractor)[J]. Optometry and Vision Science, 2000, 77(10): 537-548. doi:  10.1097/00006324-200010000-00008
    [3] MOUROULIS P. Visual Instrumentation: Optical Design and Engineering Principles[M]. New York: McGraw-Hill Professional, 1999.
    [4] 闫洁, 孟鹏花, 赵俊奇. 人眼角膜曲率测量系统的研究[J]. 应用基础与工程科学学报,2011,19(S1):254-261.

    YAN J, MENG P H, ZHAO J Q. Research of curvature measuring system of eyes cornea[J]. Journal of Basic Science and Engineering, 2011, 19(S1): 254-261. (in Chinese)
    [5] 潘兵, 俞立平, 吴大方. 使用双远心镜头的高精度二维数字图像相关测量系统[J]. 光学学报,2013,33(4):0412004. doi:  10.3788/AOS201333.0412004

    PAN B, YU L P, WU D F. High-accuracy two-dimensional digital image correlation measurement system using a bilateral telecentric lens[J]. Acta Optica Sinica, 2013, 33(4): 0412004. (in Chinese) doi:  10.3788/AOS201333.0412004
    [6] 李明东. 基于机器视觉的高分辨率双远心物镜的设计[D]. 桂林: 桂林电子科技大学, 2016: 226-232.

    LI M D. The design of high resolution double telecentric lens based on machine vision[D]. Guilin: Guilin University of Electronic Technology, 2016: 226-232. (in Chinese)
    [7] 马森, 谢芳, 刘义秦, 等. 光纤双干涉在线绝对测量技术研究[J]. 仪器仪表学报,2013,34(2):268-274. doi:  10.3969/j.issn.0254-3087.2013.02.005

    MA S, XIE F, LIU Y Q, et al. Research on optical fiber dual-interferometry for on-line and absolute measurement[J]. Chinese Journal of Scientific Instrument, 2013, 34(2): 268-274. (in Chinese) doi:  10.3969/j.issn.0254-3087.2013.02.005
    [8] 刘富国, 查学军, 杨波, 等. 基于光纤低相干干涉技术的透镜中心厚度测量方法研究[J]. 应用激光,2016,36(5):605-610.

    LIU F G, ZHA X J, YANG B, et al. Study on the method of measuring the center thickness of the lenses based on low coherence interferometry of optical fiber[J]. Applied Laser, 2016, 36(5): 605-610. (in Chinese)
    [9] SEIDEMANN A, SCHAEFFEL F, GUIRAO A, et al. Peripheral refractive errors in myopic, emmetropic, and hyperopic young subjects[J]. Journal of the Optical Society of America A, 2002, 19(12): 2363-2373. doi:  10.1364/JOSAA.19.002363
    [10] SCHAEFFEL F, FARKAS L, HOWLAND H C. Infrared photoretinoscope[J]. Applied Optics, 1987, 26(8): 1505-1509. doi:  10.1364/AO.26.001505
    [11] 闫蓓, 王斌, 李媛. 基于最小二乘法的椭圆拟合改进算法[J]. 北京航空航天大学学报,2008,34(3):295-298.

    YAN B, WANG B, LI Y. Optimal ellipse fitting method based on least-square principle[J]. Journal of Beijing University of Aeronautics and Astronautics, 2008, 34(3): 295-298. (in Chinese)
  • [1] 于海, 万秋华, 孙莹, 卢新然, 贾兴丹.  一种自适应安装的高精度图像式角位移测量装置 . 中国光学, 2020, 13(3): 510-516. doi: 10.3788/CO.2019-0107
    [2] 刘军号, 李瑞辰.  高精度干涉式光纤陀螺热漂移分析 . 中国光学, 2020, 13(2): 333-343. doi: 10.3788/CO.20201302.0333
    [3] 张洁, 倪小龙, 刘智, 齐冀, 姚海峰.  高精度连续变倍率激光扩束系统设计 . 中国光学, 2019, 12(3): 693-700. doi: 10.3788/CO.20191203.0693
    [4] 高旭, 李舒航, 马庆林, 陈伟.  光栅精密位移测量技术发展综述 . 中国光学, 2019, 12(4): 741-752. doi: 10.3788/CO.20191204.0741
    [5] 吕强, 李文昊, 巴音贺希格, 柏杨, 刘兆武, 王玮.  基于衍射光栅的干涉式精密位移测量系统 . 中国光学, 2017, 10(1): 39-50. doi: 10.3788/CO.20171001.0039
    [6] 徐乐, 张春雷, 代雷, 张健.  高精度非回转对称非球面加工方法研究 . 中国光学, 2016, 9(3): 364-370. doi: 10.3788/CO.20160903.0364
    [7] 马占龙, 彭利荣, 王高文, 谷勇强.  应用离子束修正高精度CGH基底 . 中国光学, 2016, 9(2): 270-276. doi: 10.3788/CO.20160902.0270
    [8] 陈凯, 雷枫, 伊藤雅英.  白光干涉测量超薄透明电极ITO薄膜的厚度 . 中国光学, 2015, 8(4): 567-573. doi: 10.3788/CO.20150804.0567
    [9] 董磊.  光场空间相干性的测量方法及比较 . 中国光学, 2015, 8(6): 1020-1026. doi: 10.3788/CO.20150806.1020
    [10] 王永伟, 艾华, 卓仁善, 曹艳波.  干涉测量波面重构DCT算法研究 . 中国光学, 2014, 7(6): 1012-1018. doi: 10.3788/CO.20140706.1012
    [11] 陈宝刚, 明名, 吕天宇.  大口径球面反射镜曲率半径的精确测量 . 中国光学, 2014, 7(1): 163-168. doi: 10.3788/CO.20140701.0163
    [12] 何远清, 刘永基, 翟奕.  成像角膜曲率计的光学设计 . 中国光学, 2014, 7(6): 956-961. doi: 10.3788/CO.20140706.0956
    [13] 吴威, 许廷发, 王亚伟, 闫辉, 徐磊.  高精度全景补偿电子稳像 . 中国光学, 2013, 6(3): 378-385. doi: 10.3788/CO.20130603.0378
    [14] 张春雷, 向阳, 于长淞.  高精度干涉仪针孔空间滤波器研制 . 中国光学, 2013, 6(6): 952-957. doi: 10.3788/CO.20130606.952
    [15] 张晓龙, 刘英, 孙强.  高精度非致冷长波红外热像仪的辐射标定 . 中国光学, 2012, 5(3): 235-241. doi: 10.3788/CO.20120503.0235
    [16] 张景和, 张顺国, 董跃华.  乌氏干涉测量的信息读取 . 中国光学, 2011, 4(2): 147-153.
    [17] 齐向东, 滕丽华, 于海利.  丝杠精度双频激光干涉测量中的 阿贝误差实时补偿 . 中国光学, 2010, 3(3): 279-284.
    [18] 刘艳, 苏东奇, 杨怀江, 隋永新.  高精度干涉检验移相算法对振动误差的免疫能力 . 中国光学, 2010, 3(5): 500-508.
    [19] 耿天文, 刘建红, 刘绍锦, 刘畅.  高精度光学对准测量装置的设计 . 中国光学, 2010, 3(5): 467-473.
    [20] 李东景, 于平, 李先锋, 赵昶宇.  星上定标系统探测器高精度信号处理电路的实现 . 中国光学, 2009, 2(3): 211-217.
  • 加载中
图(15) / 表ll (4)
计量
  • 文章访问数:  698
  • HTML全文浏览量:  548
  • PDF下载量:  16
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-08-27
  • 修回日期:  2019-11-08
  • 网络出版日期:  2020-06-30
  • 刊出日期:  2020-06-01

高精度角膜曲率半径测量系统

doi: 10.3788/CO.2019-0174
    基金项目:  吉林省重点研发项目(No. 20180201025GX)
    作者简介:

    李华建(1995—),男,江苏泰州人,硕士研究生,2018年于长春理工大学光电信息学院光电信息科学与工程专业获得学士学位,现为长春理工大学硕士研究生,主要研究方向为光学仪器设计、光电检测。E-mail:lhj1663126998@163.com

    肖作江(1973—),男,吉林长春人,博士,副研究员,硕士生导师,2003年于长春理工大学获得机械设计理论硕士学位,2009年于长春理工大学获得光学工程博士学位,主要研究方向为精密测量、光电检测。E-mail:xiao2706@163.com

    通讯作者: 肖作江(通讯作者)(1973—),吉林长春人,2003年获得长春理工大学机械设计理论硕士学位,2009年获得长春理工大学光学工程博士学位,现为长春理工大学副研究员、硕士生导师,主要研究方向为精密测量、光电检测。E-mail:xiao2706@163.com
  • 中图分类号: TP394.1;TH691.9

摘要: 为了实现成像角膜曲率计沿光轴方向的精确对准,提高角膜曲率的测量精度,设计了一种高精度成像角膜曲率测量系统。对该系统的成像光源、成像光学系统、干涉测量系统等进行研究。通过LED阵列均匀照射靶环形成光源;成像物镜采用双远心镜头,使景深变大,有利于对准测量,同时保证成像物镜的放大率不会受景深的影响;引用低相干干涉测量技术,利用低相干干涉信号定位角膜顶点和测量光源的位置,利用光栅尺监测扫描反射镜的位置,实现了角膜顶点到测量光源距离的精确测量。最后,分析了该系统成像物镜放大率的稳定性和角膜曲率测量误差,在理论基础上,做出了实验样机,用所设计的样机对标准角膜模拟眼进行测试,系统精度达到±0.02 mm,基本满足角膜曲率的测量要求。

English Abstract

李华建, 肖作江, 刘颖, 赵媛媛, 王睿智, 贺肖影. 高精度角膜曲率半径测量系统[J]. 中国光学, 2020, 13(3): 501-509. doi: 10.3788/CO.2019-0174
引用本文: 李华建, 肖作江, 刘颖, 赵媛媛, 王睿智, 贺肖影. 高精度角膜曲率半径测量系统[J]. 中国光学, 2020, 13(3): 501-509. doi: 10.3788/CO.2019-0174
LI Hua-jian, XIAO Zuo-jiang, LIU Ying, ZHAO Yuan-yuan, WANG Rui-zhi, HE Xiao-ying. High precision corneal curvature radius measurement system[J]. Chinese Optics, 2020, 13(3): 501-509. doi: 10.3788/CO.2019-0174
Citation: LI Hua-jian, XIAO Zuo-jiang, LIU Ying, ZHAO Yuan-yuan, WANG Rui-zhi, HE Xiao-ying. High precision corneal curvature radius measurement system[J]. Chinese Optics, 2020, 13(3): 501-509. doi: 10.3788/CO.2019-0174
    • 角膜是人眼最前面的透明部分,人眼大部分的屈光力由其提供。通过角膜曲率检测可以判断人眼有无杂散光,也可以指导角膜屈光度矫正手术以及为佩戴隐形眼镜提供科学验配依据等[1-2]。因此,精确测量人眼的角膜曲率是非常必要的。最早的角膜曲率测量装置是1856年由Helmhotlz设计的,它是一种基于目视手动调焦原理的传统刻度结果显示的仪器[3]。这种仪器结构简单,由人眼进行光学对准,对准精度依赖操作经验,测量速度慢。20世纪90年代,日本Topcon公司和Nidek公司相继研制出了光机电式角膜曲率测量装置,能够快速测出人眼角膜曲率。这种光机电式角膜曲率测量装置根据成像CCD上同心圆环的大小、光学成像系统的放大率、同心圆环的实际大小以及同心圆环到人眼角膜顶点的距离就能通过反射原理计算出角膜曲率的大小[4]。上述测量仪器测量时很容易实现人眼角膜的水平$x$方向和垂直$y$方向的对准,但在沿光轴的z方向,由于受光学成像物镜景深的限制,因而很难实现高精度的z向对焦,会引起角膜顶点到测量光源距离的微小变化和CCD上测量圆环像大小的微小变化,从而使得计算出的角膜曲率半径值不够精确。

      针对上述情况,在光机电式角膜曲率测量装置的基础上研究了一种高精度成像角膜曲率测量系统,光学成像物镜采用双远心镜头,放大率在一定景深范围内恒定不变,从而保证CCD上测量光源像的大小不会受镜头景深的影响[5-6]。同时,在对焦完成后,利用低相干干涉测量技术测量角膜顶点到测量光源的距离[7-8],以提高角膜曲率半径的测量精度。

    • 角膜是非球面的,但中间表面可以看作是一个凸面镜,角膜曲率半径测量就是利用这一特点。角膜曲率半径的测量原理图如图1所示[9-10]

      图  1  角膜曲率半径测量原理图

      Figure 1.  Measurement principle of corneal curvature radius

      图1中,$AB$的高度为$h$,从$A$发出两束光,一束射向角膜曲率中心 $C$,经角膜反射后原路返回,一束射向角膜主焦点$F'$经角膜反射后平行于光轴出射,两反射光线的交点$A'$$A$的像。根据几何关系可知:

      $$\frac{{h'}}{h} = \frac{{r - a}}{{r + b}}.$$ (1)

      $A'$经光学成像物镜DE投影后成像在CCD上,在CCD上像点A''的高度为${h_0}$,则光学系统的放大率$\beta $为:

      $$\beta = \frac{{{h_0}}}{{h'}}.$$ (2)

      由物像位置关系可知:

      $$\frac{{\rm{1}}}{a} - \frac{1}{b} = \frac{2}{r}.$$ (3)

      $b$>>$r$时,可得:

      $$r = \frac{{2b{h_0}}}{{\beta h - {h_0}}},$$ (4)

      式中$r$为角膜曲率半径,$\;\beta $为光学系统的放大率,$b$为角膜顶点到AB的距离,${h_0}$为CCD上像的大小。

    • 由式(4)可知,影响角膜曲率测量精度的有角膜顶点到测量光源的轴向距离$b$、CCD上像的大小${h_0}$、成像物镜的垂轴放大率$\;\beta $和测量光源到光轴的距离$h$

      ${h_0}$产生误差的因素主要有角膜中心的对准和CCD的图像处理,在如今的角膜曲率计中这些误差已经控制到测量精度要求内。$h$是一个固定值,测量时一般会用一个标准眼来校准这个值。而$b$$\beta $的大小受光学系统景深的影响。$b$的大小是物距$L$减去测量光源到成像物镜的距离,在几何光学中,光学系统的横向放大率$\beta $可以由下式表示。

      $$\beta = \frac{{L'}}{L},$$ (5)

      式中$L'$为像距,$L$为物距。在一种现有的角膜曲率计中,$\beta $=−0.45、$h$=20 mm,物距是100 mm,测量光源到成像物镜的距离是30 mm,则不同景深对应的$b$$\beta $值如表1所示(前景深为负,后景深为正)。

      表 1  景深对应的$b$$\beta $

      Table 1.  The $b$ and $\beta $ values corresponding to different depthes of field

      景深/mm$b$/mm$\beta $/mm
      −367−0.464
      −268−0.459
      −169−0.456
      070−0.450
      171−0.446
      272−0.441
      373−0.437

      标定的$b$值和$\,\beta $值分别是70 mm和−0.45,当景深为零时测量的${r_0}$值是理想值,设不同景深对应的测量值为${r_i}$,则角膜曲率半径的测量误差$\Delta r = $${r_i} - {r_0} $。角膜曲率半径的测量误差标准是±0.02 mm,由公式(3)可知,假设${h_0}$=0.5 mm,不同景深下,角膜曲率的测量误差如图2所示。

      图  2  误差分析图

      Figure 2.  Error analysis chart

      图2可知,当成像物镜的景深在−0.137~0.111 mm之间时,测量结果才能满足测量要求。因此,采用传统的图像清晰度方法很难实现角膜曲率的高精度测量。

    • 高精度角膜曲率半径测量系统的测量原理如图3所示,系统包括成像系统和干涉测量系统。成像系统中测量光源发出的光投射到人眼角膜前表面,反射后经成像物镜成像在CCD上。

      图  3  高精度角膜曲率半径测量系统原理图

      Figure 3.  Schematic diagram of measurement system for high precision corneal curvature radius

      干涉测量系统中,激光光源发出的光经过光纤耦合器后分成两束:一束光通过分光镜进入成像系统中,经过测量光源中心的半反半透镜和人眼角膜中心反射后原路返回;另一束光经过光学扫描延迟线原路返回,并与角膜顶点和测量光源中心反射回来的光发生干涉,利用光栅尺记录两次干涉时扫描反射镜的位置,就能精确得到角膜顶点到测量光源的距离。通过CCD上像的大小、测量光源的实际高度、成像物镜的放大率和角膜顶点到测量光源的距离就能计算出角膜曲率半径值。

    • 成像系统包括测量光源、人眼角膜、成像物镜和CCD。利用红外LED灯均匀照射如图4所示的靶环装置,形成测量光源,接着投射到人眼,若人眼是正常眼或单纯的近视眼及远视眼,则在CCD上会形成清晰的同心圆环,若是散光眼,则形成清晰的同心椭圆。

      图  4  靶环装置

      Figure 4.  Target ring device

    • 需要计算的靶环装置参数主要是靶环上光阑的高度以及靶环到角膜顶点的最优距离。图5为靶环装置成像原理图。

      图  5  靶环装置成像原理图

      Figure 5.  Schematic diagram of imaging principle of target ring device

      设靶环上光阑$AB$的高度为h,光阑像高$A'B'$$h'$,靶环到角膜顶点的距离为b,角膜曲率半径为r,根据角膜曲率测量原理可得:

      $$h' = \frac{r}{{2r + 2b}} \times h.$$ (6)

      角膜曲率半径的测量范围是 5.5~11 mm,角膜曲率的计算是在角膜中心大约4 mm区域内进行的,靶环投影在角膜表面的像高不能大于 2 mm。因此,可以得出下列关系式:

      $$\left\{ \begin{aligned} & 11h {\text{≤}} 44 + 4b \\ & 5.5h {\text{≤}} 22 + 4b \\ \end{aligned} \right..$$ (7)

      人眼的平均角膜曲率为 8 mm,根据公式(6)将靶环到角膜顶点的最优距离设计成32 mm,相应的光阑高度为20 mm。当角膜曲率半径为11 mm时,h≤15.6 mm,当角膜曲率半径为5.5 mm时,h≤27.27 mm。在满足条件的情况下,角膜曲率半径越大,对应的光阑高度就越小。因此,为保证角膜曲率的测量范围,设计高度分别为15 mm和20 mm的两组光阑。

    • 整个靶环装置如图6所示,用反光罩和光阑将光源封装在内部,同时,根据积分球漫反射原理,在反光罩内部和光阑的内表面镀上漫反射材料,有利于确保光强均匀分布,反光罩和光阑的制作材料都是不透光的。用LED光源照射光阑,从而形成靶环。

      图  6  靶环装置组成图

      Figure 6.  Composition diagram of target ring device

      利用光学仿真软件TracePro对靶环装置进行仿真,将LED发光面模拟成光源面,每个 LED的光照度模拟为使用的正常范围,且发散角为120°,并将反光罩的内部和光阑的背面均模拟成漫反射面。在靶环装置的前面32 mm处模拟一个直径为30 mm的圆面,作为光源接收面,充当人眼。靶环装置的仿真结果如图7所示,光线由光源发出,在反光罩内部和光阑的背面漫反射后从光阑的两组环形圆孔出射后照射接收面。

      图  7  靶环装置光源仿真图

      Figure 7.  Simulation diagram of light source for target ring device

      图8(彩图见期刊电子版)给出了接收面的辐照度分布图。左边为辐照度的等高线图,颜色越浅,表明辐照度越大,由图8可看出,整个接收面上辐照度的分布比较均匀。右边为接收面中间剖面的辐照度变化图,蓝色线表示水平方向的辐照度,绿色线表示竖直方向的辐照度,辐照度为旋转对称分布,因此中间剖面每一点处水平方向的辐照度和竖直方向的辐照度基本一样,并且最大辐照度小于230 W/m,根据IEC60825-1可知,基本不会对人眼造成伤害。综上所述,由靶环装置投射出来到人眼的光不会损伤人眼,且光源是均匀分布的,使CCD接收面得到的图像更易于处理与分析。

      图  8  模拟面的辐照度

      Figure 8.  Irradiances of the simulated surface

    • 成像物镜采用双远心镜头,图9为双远心物镜的成像原理图。双远心镜头的景深很大,并且双远心镜头在景深范围内放大率保持恒定,从而保证了CCD上像的大小不会受景深的影响。

      图  9  双远心镜头成像原理图

      Figure 9.  Imaging principle of double telecentric lens

      所用双远心镜头的参数如表2所示。为了使测量数据具有可比性,所选用的双远心镜头的物距和放大率同上述现有角膜曲率计中的成像物镜的物距和放大率一样。

      表 2  双远心镜头参数

      Table 2.  Double telecentric lens parameters

      参数数值
      物距/mm100
      前景深/mm4.68
      后景深/mm5.47
      放大率–0.45

      在双远心镜头的景深范围内,放大率和物距的关系如表3所示。由表3可以看出,双远心镜头的放大率的稳定性非常好,在景深范围内几乎保持不变,这有助于减小测量误差。

      表 3  不同物距的放大率

      Table 3.  Magnifications corresponding to different object distances

      物距/mm放大率
      97−0.450 19
      98−0.450 15
      99−0.450 12
      100−0.450 00
      101−0.449 92
      102−0.449 86
      103−0.449 82
      104−0.449 78
    • 图10为干涉测量系统的原理图,为了实现测量光源到角膜顶点距离的测量,在测量光源处放一个半透半反镜。干涉系统采用迈克尔逊干涉仪,光线经反射镜后原路返回,所以扫描延迟线的扫描距离比角膜顶点到测量靶环距离的两倍还大。

      图  10  干涉测量系统

      Figure 10.  Interferometry measurement system

      为提高测量精度、缩小干涉长度,测量光源采用中心波长为850 nm、光谱宽度为30 nm的宽带光源。宽带光源发出的低相干光经过光纤耦合器后分成两束:一束为测量光,经光纤准直器后进入测量臂,测量光分别在半反半透镜和人眼角膜顶点产生反射,反射光再经准直器回到光纤耦合器中;另一束为参考光,经光纤准直器后进入参考臂,经扫描反射镜反射后原路返回。测量臂中由两次反射表面反射回来的光具有不同的光程,参考臂中反射光的光程随着扫描反射镜位置的改变而不断变化,当参考反射光与测量臂中某一反射面反射光的光程差满足干涉条件时,会产生干涉信号。当光程差为零时,干涉信号达到极大值,由探测器接收并处理。为了精确获得扫描反射镜的移动距离,利用高精度可读数光栅尺监测扫描反射镜的位置,当出现干涉信号的极大值时,通过光栅尺和信号处理技术记录扫描反射镜的位置,根据产生两次干涉信号极大值时反射镜的位置,就能确定角膜顶点到测量靶环的距离,如图11所示。

      图  11  干涉图

      Figure 11.  Interferogram

    • 由上述分析可知,在高精度角膜曲率测量系统中,影响角膜曲率半径测量的有干涉仪的定位精度和成像物镜的放大率稳定性。干涉的定位精度直接影响到角膜顶点到测量光源的距离$b$,所用光源的干涉长度为100 μm,则干涉的定位精度是±50 μm,角膜顶点和测量光源距离的两次干涉定位精度为±100 μm,即$b$的测量误差为±100 μm。放大率的大小与景深有关,因此分析了当$b$的误差分别为100 μm和−100 μm时,不同景深对测量造成的误差(前景深为负,后景深为正),如图12图13所示(为保证结结果的可比性,假设$h$值和${h_0}$值分别为20 mm和0.5 mm)。

      图  12  $b$为100 μm时角膜曲率半径的误差

      Figure 12.  The error of corneal curvature radius when $b$ is 100 μm

      图  13  $b$为−100 μm时角膜曲率半径的误差

      Figure 13.  The error of corneal curvature radius when $b$ is −100 μm

      图12图13可以看出,当$b$的误差取最大值,景深从−4 mm变化到5 mm时,高精度角膜曲率测量的理论总误差符合测量标准。由此可知,成像物镜采用双远心镜头时放大率的稳定性和低相干干涉测量出的角膜顶点到测量光源的距离的精度均满足角膜曲率的测量要求。

    • 在上述的理论基础上做出了实验样机,如图14所示。由于人眼中心区域的角膜曲率半径不同,故测量的同心圆灯组经人眼反射后在CCD上成的像不再是一个标准的同心圆,而是近似为一个椭圆。图15是对CCD上获取的图像的处理结果,15(a)为直接获取的图像,15(b)为二值化处理图像,15(c)为Sobel边缘检测图像。这样便可通过提取图像中关键区域的白点,再用最小二乘曲线拟合计算椭圆长轴${h_{0\max }}$和短轴${h_{0\min }}$大小[11],即可求出相对应的角膜曲率的最小值${r_{\min }}$和最大值${r_{\max }}$

      图  14  实验样机

      Figure 14.  Experimental prototype

      图  15  图像处理结果

      Figure 15.  Image processing results

      在二维平面上,椭圆方程可以表示为:

      $$A{x^2} + Bxy + C{y^2} + Dx + Ey + F = 0.$$ (8)

      即6个点就能确定一个椭圆,但实际上由于人眼表面的曲率不是规则性的,不能直接用6个点求出椭圆的大小,本文采用了基于最小二乘法的椭圆拟合方法。最小二乘法的曲线拟合是使测量误差的平方和最小,即对于目标函数:

      $$ \begin{split} f(A,B,C,D,E,F) = \sum\nolimits_{i = 1}^n ( &Ax_i^2 + B{x_i}{y_i} + Cy_i^2 +\\ & D{x_i} + E{y_i} + F), \end{split} \!\!\!\!$$

      使:

      $$\frac{{\partial f}}{{\partial A}} = \frac{{\partial f}}{{\partial B}} = \frac{{\partial f}}{{\partial C}} = \frac{{\partial f}}{{\partial D}} = \frac{{\partial f}}{{\partial E}} = \frac{{\partial f}}{{\partial F}} = 0,$$ (10)

      成立,由此可得到一个线性方程组。根据图像中得到的点,应用求解方程组的算法即可求出系数的大小,从而可根据椭圆定义求出椭圆长轴、短轴和轴位角的大小,进而就能算出角膜中心处曲率半径的最大值${r_{\max }}$和最小值${r_{\min }}$。根据公式(11)就可算出曲率半径所对应的屈光度的值(其中n为人眼的折射率,一般取1.337 5)。

      $$F = \frac{{(n - 1)1\;000}}{r}.$$ (11)

      样机能对角膜曲率半径、角膜屈光度进行测量,通过编制程序即可完成转换。用样机对标准模拟眼进行测量,对3种不同的标准眼重复测量3次,测量结果如表4所示。测量结果使用图像处理拟合椭圆方法计算得到,所以角膜曲率测量时计算的是最大值和最小值。从表4可以看出:系统的测量精度在0.02 mm以内,说明成像物镜放大率的稳定性非常好,同时低相干干涉测量精度也非常高;重复性在0.03 mm以内,说明测量结果不受景深的影响。

      表 4  实验结果

      Table 4.  Experimental results

      标准值测量值/mm误差/mm标准值测量值/mm误差/mm
      6.656${r_{\max }}$6.675﹢0.0197.988${r_{\max}}$7.9970.009
      ${r_{\min }}$6.641−0.015${r_{\min}}$7.975−0.013
      6.656${r_{\max}}$6.673+0.0179.458${r_{\max}}$9.466+0.008
      ${r_{\min }}$6.637−0.019${r_{\min }}$9.443−0.015
      6.656${r_{\max}}$6.669+0.0139.458${r_{\max}}$9.473+0.015
      ${r_{\min }}$6.647+0.011${r_{\min }}$9.455−0.003
      7.988${r_{\max}}$7.999+0.0119.458${r_{\max}}$9.477+0.019
      ${r_{\min }}$7.974−0.014${r_{\min }}$9.459+0.001
      7.988${r_{\max}}$7.998+0.010
      ${r_{\min }}$7.980−0.008
    • 本文在传统角膜曲率计的基础上,提出一种高精度角膜曲率测量系统,其成像物镜采用双远心光路镜头,以使放大率在景深范围内保持恒定,从而降低了由于CCD上所成像大小的微小变化对测量结果的影响。利用一种低相干干涉测量技术测量出角膜顶点到测量光源的距离,距离测量精度达到±100 μm,相比于通过成像清晰度法判断角膜顶点到测量光源的距离,这种方法大大提高了角膜曲率半径系统的测量精度,同时由于成像物镜采用的是双远心镜头,也避免了无法利用成像清晰度法判断角膜顶点到测量光源距离的缺陷。用实验样机对标准眼进行测试,测量结果表明,高精度角膜曲率测量系统对角膜曲率测量的精度满足测量要求,系统测量精度达到±0.02 mm。

参考文献 (11)

目录

    /

    返回文章
    返回