技术领域
[0001] 本
发明涉及人眼色差测量的技术领域，具体涉及一种哈特曼人眼色差测量系统。
背景技术
[0002] 人类的眼睛结构复杂，它的成像系统主要由
角膜、前房、虹膜、晶状体和玻璃体组成。每一部分的屈光程度不一样，而大多数情况我们希望研究的
眼底视网膜却在人眼的成像系统之后。多
波长光经过这一透光系统进入眼底视网膜的过程，必然伴随色差的产生，有研究表明，角膜的不规则是横向色差产生的主要原因。换言之，无论是希望进行多波长的眼底视网膜成像，还是
白内障复明手术所需要的人工晶体(IOL)眼的评价，还有其他术前、术后眼的色差评估都可能要对人眼色差进行测量、校正和补偿。因此，研究人眼色差的测量有极其重大的意义，是视网膜成像
质量提高和白内障患者在可见光下清晰视物的重要前提。
[0003] 由几何光学的知识，可以将人眼色差可以分为轴向色差和横向色差。轴向色差(Longitudinal or Axial Chromatic Aberration)简称LCA，描述两种色光对轴上物点成像
位置的差异，造成视网膜成像的离焦；横向色差(Transverse or Lateral Chromatic Aberration)简称TCA，主要是由人眼光学系统本身色散性质造成，在视网膜成像中表现为图像放大率的不同和空间位移的改变(见图1)。
[0004] 若将人眼由角膜、前房、虹膜、晶状体和玻璃体组成的成像系统，看成一个透镜组合系统，又由于物方为实物，物方色差为0，可以定义轴上色差为(几何光学、像差、光学设计，李晓彤，岑兆丰，范世福，2012年12月第二版)：
[0005]
[0006] 其中ΣC1为初级轴上色差系数，也定义为第一色差和数；n′k为不同组成部分的折射率，uk则为每一部分的视场角。
[0007] 类似的，横向色差可以定义为：
[0008]
[0009] 其中ΣCΠ为初级横向色差系数，也定义为第二色差和数。从横向色差的公式可以看出，初级横向色差仅与视场的一次方成比例，表明光学系统在视场不大时，就会有横向色差的有害影响。
[0010] 单色球面波(或平面波)经过光学系统后，将由于像差而发生
变形。如果物方球面波时复色的，那么各色波面经系统后，将因各自像差的不同而又不同程度的变形。而不同波长的两光波面间的偏离量，可用来表征色差，称之为波色差。按照习惯，轴向色差一般用光焦度作为单位，而横向色差则采用角度作为单位。
[0011] 从二十世纪40年代开始，已经有人开展了人眼色差测量的研究，George Wald等人采用
光谱细孔屈光镜测量得到14只人眼在365nm至750nm波段的LCA为3.2D(George Wald,Donald R,Griffin.May 1947,Vol.37,No.5：321～336)。一直到今天，色差测量的研究从未停止过，色差测量的方法可分为直接测量法和间接测量法，其中直接测量法由主观测量发展为更为准确、快捷的客观测量，而间接测量法则通过
软件建立模拟人眼进行仿真研究。
[0012] 若按照色差种类划分，也就是分为轴向色差和横向色差。轴向色差随着哈特曼人眼像差技术的发展实现了测量，通过将不同波长光的泽尼克多项式第四离焦项系数相减可得轴向色差。
[0013] 2008年，Silvestre Manzanera等人采用哈特曼波前
传感器对人眼LCA进行了客观测量。(A wavelength tunable wavefront sensor for the human eye.Silvestre Manzanera,Carmen Canovas,Pedro M.Prieto.et al.OPTICS EXPRESS.26May 2008.Vol.16,No.11：7748～7755)该方法具体采用了一个Xe白光灯和一个用于波长选择的干涉
滤波器作为
光源，一台哈特曼波前传感器用于人眼波前信息的采集，再采用一个可动变焦模
块适应波长改变焦距的调整。该研究方法没有实现人眼横向色差的测量。
[0014] 对于横向色差的测量，1987年，Youmay U.Ogboso和Harold E.Bedell通过双色视标法从主观上测量横向色差(Magnitude of lateral chromatic aberration across the retina of the human eye.Youmay U.Ogboso,Harold E.Bedell.Optical Society of America.1987.August.Vol.4,No.8.1666～1672)。2012年，Wolf M.Harmening通过不同波长光的眼底成像图之间的偏移量，客观测量横向色差，并与主观方法做了比较，验证了准确性。(Measurement and correction of transverse chromatic offsets for multi-wavelength retinal microscopy in the living eye.Wolf M.Harmening,Pavan Tiruveedhula,Austin Roorda.BIOMEDICAL OPTICS EXPRESS.2012september.Vol.3,No.9.2066～2077)2013年，饶丰等人在其
专利(基于哈特曼传感器的人眼色差测量装置，授权公告号：CN103230254A)中提出了另一种采用哈特曼测量人眼色差的方法，可分时测量人眼的横向色差与轴向色差。该方法采用交替发光的单色光光源，通过一台哈特曼波前传感器分时测量每一波长的波前信息，再通过计算得到人眼的色差。该方法有不足之处：
[0015] 1、在实际的哈特曼单色相差测量中，固视状态下，活体人眼自身不可避免的有震颤tremor、微快速扫视microsaccade、和偏移drift的存在(Ralf Engbert,Reinhold Kliegl.Vision Research.432003：1035–1045)，这些视网膜抖动使得瞳孔对应的微透(或棱)镜上光线发生变化，有可能影响横向色差的准确性，必须予以考虑。
[0016] 在这些固视视网膜抖动中，不规则的、高
频率(50～100Hz)的震颤幅度极小，约为视锥细胞直径大小，仅约20秒视角；微快速扫视速度约每秒几分视角；由大量震颤运动构成的偏移可达6分视角，是缓慢的曲线运动。可知微快速扫视是造成视网膜抖动的主要原因。与此同时，人眼内在的像差的动态变化也给测量带来影响。
[0017] 2、由于哈特曼的量程限制，对于屈光不正的人眼可能无法实现完美的视网膜聚焦，造成测量结果不准确。
发明内容
[0018] 本发明克服了目前哈特曼色差测量系统的不足，提供一种基于哈特曼波前传感器的人眼色差测量系统，是一种人眼横向色差与轴向色差可同时测量的光学仪器，旨在实现准确的人眼两种色差(横向色差、轴向色差)的测量，从而为眼底成像诊断眼科
疾病和病眼
治疗提供了更为可靠的数据。通过三个哈特曼波前传感器对三波长光进行同时
数据采集，避免了视网膜抖动和人眼像差随时间动态变化造成的误差，简化了实验；可通过在人眼前加入试镜片trial，对屈光不正的被测人眼进行预校正，并在系统中加入离焦补偿系统进行补偿。
[0019] 本发明采用的技术方案是：哈特曼人眼色差测量系统，包括可预调焦的可见光或
近红外光信标，光阑，第一分光镜，活体人眼，调焦系统，第二分光镜，口径匹配系统，第一滤波器，第二滤波器，第二反射镜，三个哈特曼波前传感器，观察目标系统和计算机组成；可预调焦的可见光或近红外光信标包括第一波长信标，第二波长信标，第三波长信标，第一反射镜，第一合束镜和第二合束镜；三个哈特曼波前传感器为第一哈特曼波前传感器，第二哈特曼波前传感器和第三哈特曼波前传感器；
[0020] 可见光或近红外光三波长信标的第一波长信标发出信标光，信标光预先
准直后，经第一反射镜反射后再经第一合束镜透射到达第二合束镜后被第二合束镜反射，可见光或近红外光三波长信标的第二波长信标发出信标光，信标光预先准直后，经第一合束镜反射后到达第二合束镜后被第二合束镜反射，可见光或近红外光三波长信标的第三波长信标发出信标光，信标光预先准直后，经第二合束镜透射，该第二合束镜将三束光波汇聚成一束不相干的离散光，经光阑、第一分光镜反射进入人眼瞳孔；人眼眼底视网膜漫反射的后向散射光，通过第一分光镜、调焦系统、第二分光镜，经口径匹配系统，第一滤波器选择出第一波长光透射进入第一哈特曼波前传感器，同理，由第一滤波器反射的剩余的两束光，经第二滤波器选择出第二波长光反射到第二哈特曼波前传感器，最后一束光透过第二滤波器由第二反射镜反射，进入第三哈特曼波前传感器，三个哈特曼波前传感器将采集到的光斑图像送至计算机，计算机根据测得的三波长人眼波相差经控制软件转化为泽尼克多项式，通过各波长泽尼克多项式系数第四项(即离焦)之差计算人眼轴向色差；通过各波长在相应哈特曼波前传感器的每一微透镜或棱镜或某几个微透镜或棱镜区域的位置偏差，计算对应瞳孔位置的局部横向色差；
[0021] 通过改变观察目标系统的位置，让人眼主动进行眼球的转动，从而改变视轴与测量轴的角度，用于测量不同视轴角度下的人眼色差；通过前后移动调焦系统，配合观察目标系统进行测量前的调焦。
[0022] 其中，所述的可见光或近红外光三波长信标是考虑多个离散波长组合的
叠加效应，同时入射人眼，满足人眼安全剂量，可以是可见光、近红外
激光器，或可见光、近红外
半导体激光器，或可见光、近红外超
辐射半导体器件。
[0023] 其中，所述的三个哈特曼波前传感器实现同时采集不同波长的波前数据。
[0024] 其中，所述的局部横向色差测量，可以通过增减微透镜或棱镜阵列的单元数或者选择某几个临近的微透镜或棱镜来缩小或扩大每次测量的瞳孔局部范围。
[0025] 其中，所述的三个哈特曼波前传感器均是基于微棱镜阵列的哈特曼波前传感器，或基于微透镜阵列的哈特曼波前传感器。
[0026] 其中，所述调焦系统可以是一个双透镜4F系统、Badal调焦系统，用于补偿系统无法测量的较大离焦。
[0027] 本发明与
现有技术相比所具有的优点是：
[0028] 1、本发明实现了人眼轴向色差与横向色差的同时测量，且为客观测量法，不需要被测量者接受训练；
[0029] 2、本发明采用三台哈特曼同时测量不同波长的人眼波像差，克服了微快速扫视、人眼像差
波动带来的实验误差，即避免了时间因素带来的误差。
[0030] 3、本发明有调焦系统，可以补偿哈特曼无法测量的较大离焦。
附图说明
[0031] 图1为几何光学概念上的轴向色差与横向色差示意图；
[0032] 图中虚线和实线分别表示不同波长的两色光，左上图为像空间的轴向色差，右上图为物空间的轴向色差，左下图为某一瞳孔位置像空间的横向色差，右下图为某一瞳孔位置物空间的横向色差；
[0033] 图2为三台哈特曼人眼色差测量系统结构示意图；
[0034] 图中1、3、5为可预调焦的可见或近红外光三波长信标的第一波长信标、第二波长信标、第三波长信标，2为第一反射镜，4为第一合束镜，6为第二合束镜，7为光阑，8为第一分光镜，9为活体人眼，10为调焦系统，11为第二分光镜，12为口径匹配系统，13为第一滤波器，15为第二滤波器，17为第二反射镜，14、16、18为三个哈特曼波前传感器，即第一哈特曼波前传感器14，第二哈特曼波前传感器16，第三哈特曼波前传感器18，19为观察目标系统和20为计算机；
[0035] 图3为哈特曼波前传感器微透镜阵列，某一微小透镜的光束汇聚到CCD面示意图；
[0036] 图中实线和虚线分别表示波长不同的两束光，实心点与空心点分别对应实线和虚线的光，左图为两束不同波长光在Y方向上的偏差，右图为左视图，可见两束不同波长光在X，Y方向上的偏差。
具体实施方式
[0037] 下面结合附图及具体实施方式详细介绍本发明。
[0038] 如图2所示，本发明一种哈特曼人眼色差测量系统包括可预调焦的可见光或近红外光信标，光阑7，第一分光镜8，活体人眼9，调焦系统10，第二分光镜11，口径匹配系统12，第一滤波器13，第二滤波器15，第二反射镜17，三个哈特曼波前传感器14、16、18，观察目标系统19和计算机20。可预调焦的可见光或近红外光信标包括第一波长信标1，第二波长信标3，第三波长信标5，第一反射镜2，第一合束镜4和第二合束镜6，其中可见或近红外光信标可以是激光器laser、半导体激光器laser diode和超辐射半导体器件super luminescent diode—SLD；第一、第二、第三哈特曼波前传感器14、16、18可以是基于微棱镜阵列的哈特曼波前传感器，或是基于微透镜阵列的哈特曼波前传感器；合束镜、滤波器可以是热镜、冷镜或者二向色滤波器；分光镜可以是玻璃分光镜、
薄膜分光镜；调焦系统可以是一个双透镜4F系统、Badal调焦系统。
[0039] 本实例的哈特曼人眼色差测量系统工作过程如下：可见或近红外三波长信标的第一波长信标1、第二波长信标3、第三波长信标5发出信标光，信标光预先准直后，经第一反射镜2，第一合束镜4，第二合束镜6，将三束光波汇聚成一束不相干的离散光，经光阑7、第一分光镜8反射进入人眼9瞳孔；人眼9眼底视网膜漫反射的后向散射光，通过第一分光镜8、调焦系统10、第二分光镜11，经口径匹配系统12，由第一滤波器13选择出第一波长光透射进入第一哈特曼14，同理，由第一滤波器13反射的剩余的两束光，经第二滤波器15选择出第二波长光反射到第二哈特曼16，最后一束光透过第二滤波器15由第二反射镜17反射，进入第三哈特曼18，三台哈特曼将采集到的光斑图像送至计算机20，计算机20根据测得的三波长人眼波相差经控制软件转化为泽尼克多项式，通过各波长泽尼克多项式系数第四项(离焦)之差计算人眼轴向色差；通过各波长在相应哈特曼波前传感器的每一微透镜(或棱镜)或某几个微透镜(或棱镜)区域的位置偏差，计算对应瞳孔位置的横向色差。通过改变观察目标系统19的位置，让人眼9主动进行眼球的转动，从而改变视轴与测量轴的角度，用于测量不同视轴角度下人眼9的色差；通过前后移动调焦系统10，配合观察目标系统19进行测量前的调焦。
[0040] 被测波前为圆域波前时，通常可采用一组泽尼克多项式来描述：
[0041]
[0042] 式中， 为哈特曼波前传感器的入射光波前(即包含人眼像差信息的反射光波前)，a0为平均位相波面；ak为第k项泽尼克多项式系数；Zk为第k项泽尼克多项式。
[0043] 泽尼克多项式是一组在圆域上
正交的多项式，它在单位圆上定义为：
[0044]
[0045] 其中，
[0046]
[0047] θ表示单位圆不同角度，m和n分别为角向频率和径向频率，它们恒为整数并且满足：
[0048] m≤n，n-|m|＝even(4)
[0049] 第一、第二、第三哈特曼波前传感器14、16、18可以是基于微棱镜阵列的哈特曼波前传感器，也可以是基于微透镜阵列的哈特曼波前传感器，两者的区别主要在采集波前信息的方式，但对采集到的波前信息的处理过程是一致的(如图2)。在微棱镜或微透镜阵列焦平面采集到的光斑偏移信息，通过计算机的处理，采用质心
算法，假设光斑的位置为(xi，yi)，则探测到的全孔径的波面误差信息可表示为：
[0050]
[0051]
[0052] 式中，m＝1～M，n＝1～N为子孔径映射到CCD光敏靶面上对应的
像素区域，M和N分别为子孔径映射到光敏靶面上对应区域的横向和纵向像素数，Inm是CCD光敏靶面上第(n,m)个像素接收到的
信号，xnm，ynm分别为第(n,m)个像素的x坐标和y坐标。
[0053] 再根据以下公式计算入射波前的波前斜率Gx，Gy：
[0054]
[0055]
[0056] 式中，s为子孔径面积；Φ(x,y)为入射光束波前
相位；f为微透镜焦距。
[0057] 像差波前转换为可用于计算的泽尼克多项式，一个完整的波前Φ(x,y)可以用泽尼克多项式(1)来描述。子孔径内的斜率数据与泽尼克多项式系数的关系为：
[0058]
[0059] 其中εx、εy为波前相位测量误差，n为模式阶数，Zxk(i)和Zyk(i)是第k项泽尼克多项式在第i子孔径内的平均斜率，Si为子孔径的归一化面积。m个子孔径斜率n项泽尼克系数的关系用矩阵表示为：
[0060]
[0061] 记为：
[0062] G＝DA+ε (11)
[0063] 对于任意的2m和n，上述方程的最小二乘和最小范数解可用广义逆D+表示：
[0064] A＝D+G (12)
[0065] 此时得到泽尼克系数A。
[0066] 泽尼克系数A用Cn表示，其中n表示泽尼克系数第n项。
[0067] 通过波前拟合出的不同波长的波平面之差，可以得到不同波长间的波色差。
[0068] 通过泽尼克系数A求解轴向色差的过程如下：
[0069] Y方向和X方向散光项C3、C5，可计算散光的角度：
[0070]
[0071] 通过计算两个中间变量B和E，其中C4为离焦项：
[0072] or
[0073]
[0074] 通过(11)、(12)计算以屈光度为单位的两波长λ1、λ2的轴向色差，其中R为瞳孔半径, 为λ1的E中间变量， 为λ2的E中间变量：
[0075]
[0076] 通过微透镜或棱镜阵列的每一微透镜或棱镜或某几个微透镜或棱镜区域中各波长的偏离量求解横向色差过程如下(以下以某一微透镜为例)：
[0077] 已知某一微透镜阵列的每一个微透镜的半径为r，焦距为f，波长λ1汇聚于该微透镜区域(x1，y1)，波长λ2汇聚于(x2，y2)，两波长汇聚点在X方向距离为Δx，在Y方向距离为Δy，则可以得到该微透镜对应于人眼瞳孔位置的在X、Y两个方向上的局部横向色差βy、βx为：
[0078]
[0079]
