一种个体人眼大视场入射波前像差的获取方法与流程

1.本发明涉及眼视光学领域，特别是涉及一种个体人眼大视场入射波前像差的获取方法。


背景技术：

2.目前，近视已成为一个世界性的问题，近视人群基数的增大以及近视的早发，会导致白内障、青光眼、视网膜脱落等致盲性疾病的高发。因此，近视的有效防控引起了科学界以及政府部门的高度关注。
3.光学方法是近视防控的有效手段之一，通过视网膜周边的近视性离焦能够有效地控制近视的发展。因此，针对视网膜周边离焦的研究是光学方法防控近视的重要研究内容。其中，波前像差仪是研究周边视场成像性质的有效手段。
4.目前，临床上测量人眼波前像差的仪器可以分为入射型和出射型两类。(1)对于出射型波前像差仪，测量的是光线从人眼出射时的波前像差。该像差仪的主要代表是基于hartmann
‑
shack原理的波前像差仪，其原理为视网膜上的点光源发出的光线通过人眼出射，并经过一系列光学器件最终成像在ccd上面。由于人眼的像差的存在，ccd所成的像，会偏离理想成像位置。根据实际像点相对理想像点的偏离程度，得到人眼波前像差。(2)入射型波前像差仪的典型代表为主观型的波前像差仪，其测量原理为：每次测量，有两条光线入射至人眼，第一条光线与人眼的光轴重合，与视网膜相较交于a点；另外一条光线与光轴存在一定的距离，当其与第一条光线平行时，由于人眼存在像差，其与视网膜的交点b与a点并不重合。通过调节第二条光线的角度，使b点与a点重合。重复以上过程，根据第二条光线的倾斜角度，通过一系列的计算可以得出人眼的波前像差。该类型像差仪，虽命名为入射型像差仪，但根据光路可逆的原理可知，该像差仪相当于将点光源置于视网膜上，然后逐一测量每条光线对应的波前的斜率，进而重构出波前像差。综上，目前临床上采用的出射型和入射型的波前像差仪，实际上测量的是光线从人眼出射时人眼的波前像差，即人眼的出射波前像差。
5.然而，对于近视防控具有意义的波前像差应为光线进入人眼时对应的波前像差，即人眼的入射波前像差。由于临床上测量人眼的入射波前像差存在诸多的实际困难，虽有一些方法被提出，但主要测量的是0
°
视场的波前像差；此外，这些方法也存在重复性、准确性不高的问题。
6.因此，找到一种能够得到大视场下人眼的入射波前像差的方法成为现今需要解决的问题。


技术实现要素：

7.本发明的目的是提供一种个体人眼大视场入射波前像差的获取方法，利用出射波前像差数据和其他临床测量数据，并采用构建人眼模型的方式，来得到大视场人眼入射波前像差，该方法具有可供分析的视场角度大、不增加设备购置成本、易于推广等特点。
8.为实现上述目的，本发明提供一种个体人眼大视场入射波前像差的获取方法，具体包括以下步骤：
9.s1、根据角膜地形图数据、人眼出射波前像差数据和navarro人眼模型参数，构建半个性化出射眼模型；
10.s2、对所述半个性化出射眼模型的晶状体进行优化，得到大视场下出射个性化眼模型；
11.s3、对所述大视场下出射个性化眼模型进行反转操作，得到大视场下入射个性化眼模型；基于所述大视场下入射个性化眼模型，得到人眼大视场入射波前像差。
12.优选的，所述s1具体为：
13.s1.1、建立角膜面型；
14.s1.2、采集离散的角膜地形图数据，并利用所述角膜地形图数据进行拟合，得到所述角膜面型的曲率半径、圆锥系数和各项zernike系数；
15.s1.3、采用波前像差仪测量大视场下的人眼出射波前像差数据；
16.s1.4、将navarro人眼模型参数代入zemax中，得到初始人眼模型；基于所述s1.2和所述初始人眼模型，得到半个性化出射眼模型。
17.优选的，所述角膜面型，具体表达式为：
[0018][0019]
式中，x、y表示角膜的坐标，z表示角膜高度，r为曲面的曲率半径，k为曲面的圆锥系数，z
i
(x,y)为第i项zernike多项式，a
i
为第z
i
(x,y)项的系数，i＝1,2,3...30。
[0020]
优选的，所述角膜地形图数据通过pentacam三维眼前节分析仪进行采集。
[0021]
优选的，所述s1.4具体为：
[0022]
将所述角膜面型的曲率半径、圆锥系数和各项zernike系数替代所述初始人眼模型中的角膜前表面和后表面；将临床测量的眼轴数据，替换所述初始人眼模型中的眼轴数据；然后，将所述初始人眼模型的入射光瞳大小、波长、视场参数设置为与测量人眼出射波前像差数据时的设置相一致，进而得到半个性化出射眼模型。
[0023]
优选的，所述s2具体为：
[0024]
s2.1、将半个性化出射眼模型的角膜前表面的设置为所述角膜面型，并将所述角膜面型的曲率半径、圆锥系数和各项zernike系数代入所述角膜面型中；
[0025]
s2.2、将半个性化出射眼模型的晶状体的前后表面设置为所述角膜面型，并将其各参数设为变量；
[0026]
s2.3、将0
°
视场的波前像差数据作为优化目标，依次优化所述晶状体的前后表面的曲率半径、圆锥系数和各阶zernike系数，得到0
°
视场的出射型眼模型；
[0027]
s2.4、基于所述0
°
视场出射型眼模型，增加水平视场和竖直视场，并根据不同视场设置相应的瞳孔直径；然后将所述0
°
视场、水平视场和所述竖直视场下的波前像差数据作为优化目标，并且对优化函数的权重进行相应的调整；
[0028]
s2.5、继续优化晶状体的参数，直到所述navarro出射眼模型在各个视场下的波前像差与实测值基本一致为止，优化完成，得到大视场下出射个性化眼模型。
[0029]
优选的，所述s2.3具体为：
[0030]
将0
°
视场的波前像差数据作为优化目标，依次优化晶状体前后表面的曲率半径、圆锥系数和各阶zernike系数，直至所述半个性化出射眼模型的0
°
视场波前像差与实际测量的波前像差值基本一致，优化完成，从而构建出0
°
视场的出射型眼模型。
[0031]
与现有技术相比，本发明具有以下技术效果：
[0032]
本发明采用了个体人眼的临床数据，最终得到人眼模型为个性化的人眼模型，在临床应用中具有巨大的价值；本发明解决了实际测量大视场人眼入射波前像差存在的困难，借助临床上容易获得的出射波前像差，通过一系列的处理，得出入射波前像差。
附图说明
[0033]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0034]
图1为本发明实施例的方法流程图；
[0035]
图2为本发明实施例的人眼波前像差0
°
测量结果图；
[0036]
图3为本发明实施例的人眼数据构建的大视场出射个性化眼模型示意图；
[0037]
图4为本发明实施例的人眼数据的大视场入射个性化眼模型示意图。
具体实施方式
[0038]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0039]
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
[0040]
实施例1
[0041]
参照图1所示，本发明提出一种个体人眼大视场入射波前像差的获取方法(本实施例以一只人眼的数据为例，对大视场入射人眼模型的构建过程做详细说明)，具体包括以下步骤：
[0042]
s1、根据角膜地形图数据、实际测量得到的大视场(包括水平视场、竖直视场以及其他方向的视场)的人眼出射波前像差数据和navarro人眼模型参数，构建半个性化出射眼模型，具体为：
[0043]
s1.1、建立角膜面型，其表达式为：(该角膜面型能够更精准的描述角膜形态。)
[0044][0045]
式中，x、y表示角膜的坐标，z表示角膜高度，r为曲面的曲率半径，k为曲面的圆锥系数，z
i
(x,y)为第i项zernike多项式，a
i
为第z
i
(x,y)项的系数，在角膜前表面的拟合中，i
＝1,2,3...30。
[0046]
s1.2、通过pentacam三维眼前节分析仪(oculus，德国)采集离散的角膜地形图数据，并将采集的角膜地形图数据存储于excel表格中。以上数据在matlab中拟合出(1)式中的各个参数r,k,a
i
，其中i＝1,2,3...30。
[0047]
其中，通过角膜地形图数据拟合出的曲率半径r、圆锥系数k和各项zernike系数a
i
的结果如表1所示：
[0048]
表1
[0049]
[0050][0051]
s1.3、利用波前像差仪wavescan测量0
°
视场、(
‑
15
°
，
‑
10
°
，
‑5°
，+5
°
，+10
°
，+15
°
)水平视场和(
‑
15
°
，
‑
10
°
，
‑5°
，+5
°
，+10
°
，+15
°
)竖直视场的波前像差。如图2所示，显示了该眼在0
°
视场的波前像差数据。将波前像差数据输入zemax中的相应参数中。
[0052]
s1.4、将表2的navarro人眼模型参数输入到zemax中，得到出射人眼模型的初始结
构(即，初始人眼模型)。
[0053]
表2
[0054][0055]
然后，将根据角膜地形图数据拟合得出的曲率半径r、圆锥系数k、第1项到第30项zernike系数a
i
替代所述初始人眼模型中的角膜前表面和后表面。如果测量了人眼的眼轴数据，则可用临床测量得到的眼轴数据，替换所述初始人眼模型中的眼轴数据。同时将初始人眼模型的入射光瞳大小、波长、视场等参数设置为与测量人眼波前像差时的设置相一致，进而得到半个性化出射眼模型。在本实施例中，将视场角设置为0
°
，出射光瞳直径与实测相同为7.99mm，波长为780nm。
[0056]
s2、对半个性化出射眼模型的晶状体进行优化，得到大视场下出射个性化眼模型；
[0057]
将半个性化出射眼模型的角膜前表面的设置为式(1)表示的面型，并将前面拟合得出的曲率半径r、圆锥系数k和各项zernike系数代入。将半个性化出射眼模型的晶状体的前后表面设置为式(1)表示的面型，并将各参数设为变量。将0
°
视场的波前像差数据作为优化目标，依次优化晶状体前后表面的曲率半径、圆锥系数和各阶zernike系数，直至人眼模型的0
°
视场波前像差与实际测量的波前像差值基本一致。此时优化完成，从而构建出0
°
视场的出射型眼模型。
[0058]
在上述优化完成的0
°
视场出射型眼模型的基础上，增加
‑
15
°
、
‑
10
°
、
‑5°
、+5
°
、+10
°
、+15
°
水平视场和
‑
15
°
、
‑
10
°
、
‑5°
、+5
°
、+10
°
、+15
°
竖直视场，并根据不同视场设置相应的瞳孔直径。然后将这些视场下的波前像差数据作为优化目标，并且对优化函数的权重进行相应的调整。
[0059]
继续优化晶状体的参数，直到人眼模型在各个视场下的波前像差与实测值基本一致为止，优化完成。此时，就确定了眼模型晶状体的面型参数，进而得到大视场下出射个性化眼模型，如图3所示。
[0060]
其中，最终优化完的晶状体面型参数见表3：
[0061]
表3
[0062]
[0063][0064]
s3、将上述大视场下出射个性化眼模型进行反转操作，使其变为大视场下入射个性化眼模型，如图4所示、
[0065]
基于该人眼模型，可以分析其在6mm瞳孔、555nm下的某特定视场下的像差。表4给出了大视场下入射个性化眼模型在6mm瞳孔、555nm波长、0
°
和水平10
°
视场下的波前像差(单位波数)。
[0066]
表4
[0067][0068]
以上所述的实施例仅是对本发明的优选方式进行描述，并非对本发明的范围进行限定，在不脱离本发明设计精神的前提下，本领域普通技术人员对本发明的技术方案做出的各种变形和改进，均应落入本发明权利要求书确定的保护范围内。