基于矢量像差理论的屈光手术后视觉质量评价方法和系统

本发明属于视光学领域，特别是涉及基于矢量像差理论的屈光手术后视觉质量评价方法和系统。

背景技术：

随着物理光学和医学技术的快速发展，屈光手术逐渐成为矫正屈光不正的重要手段。目前屈光手术技术成熟度较高、治疗高效、安全性较高，与传统佩戴框架眼镜相比，屈光手术可以使患者获得较好的裸眼视力，摆脱了眼镜带来的不便，也满足了部分职业对于裸眼视力的要求，屈光手术已成为目前屈光不正的主流矫正方式之一。随着近年来近视等屈光不正患者的快速增加以及屈光手术技术的发展，屈光手术在全世界得到了广泛的应用，目前屈光手术在全球的手术量已超过2亿。

屈光手术的目的是矫正患者的屈光不正，提升患者的裸眼视力，一般采用视力来评价手术的有效性。随着社会发展以及人们生活水平的提高，患者对于视觉质量的要求也相应提高。而视力仅可以反映在高空间频率下患者的视觉质量，无法全面反映日常生活中患者的视觉效果。

对比敏感度曲线(contrastsensitivityfunction,csf)可以反映多种空间频率下的视觉质量，但是这种方式依赖患者主观感知，受心理状态等主观因素影响，无法准确评估人眼这一光学系统实际的成像质量。由于屈光手术后像差处于较低水平，此时像差的均方根(rootmeansquare，rms)值与视觉质量的对应关系较弱，因此像差均方根也无法准确评价屈光手术后的视觉质量。本发明能够提供一种基于矢量像差理论的新型的能够全面、准确、客观地评价屈光手术后视觉质量的评价方法。

技术实现要素：

为了有效克服对于屈光手术后像差处于较低水平的人眼无法全面、准确、客观地评估其视觉质量的问题，本发明提出基于矢量像差理论的屈光手术后视觉质量评价方法，对屈光手术后患者视觉质量进行综合的客观评价。

为实现上述目的，本发明提供了如下方案：

基于矢量像差理论的屈光手术后视觉质量评价方法，包括如下步骤：

s1.建立标准眼模型，采用zernike自由曲面表征所述标准眼模型的角膜前表面，对标准眼模型全眼像差进行定量表述，建立标准眼模型角膜改变量与标准眼模型全眼像差的定量关系；

s2.获取屈光手术后的临床患者全眼像差系数，所述临床患者全眼像差系数采用临床标准zernike多项式表示；

利用所述s1建立的定量关系，反向求解患者眼模型的角膜改变量及对应的zernike系数；

s3.利用s2得到的zernike系数表征所述标准眼模型的角膜前表面，获得术后患者眼模型；

所述术后患者眼模型具有眼波前像差分布；

s4.使用光学成像系统像质评价方法得到视觉质量客观评价指标，完成对所述术后患者眼模型进行患者手术后视觉质量进行综合的客观评价；所述光学成像系统像质评价方法包括但不限于调制传递函数、点扩散函数、三维图像仿真。

优选的，所述zernike自由曲面的表达式为：

其中，表示二次曲面，c和k分别为二次曲面的顶点曲率和圆锥系数，为zernike多项式各项子项的叠加式，其中ki代表zernike多项式系数，zi表示第i项zernike多项式，i对应不同波前像差，ρ、分别为zernike多项式的径向坐标和方位角，z(x,y)为自由曲面在坐标(x,y)处的矢高。

优选的，对所述标准眼模型全眼像差进行定量表述的过程为：利用光学软件对所述标准眼模型进行真实光线追迹，获取边缘光线在所述标准眼模型的角膜前表面的入射高度h，以及边缘视场主光线在所述标准眼模型的角膜前表面的入射高度将h、和共轴zernike自由曲面像差结合，完成对所述标准眼模型全眼像差的定量表述；

所述共轴zernike自由曲面像差，为在共轴zernike自由曲面光学系统的非光阑处，其zernike项对应系数k产生的光学系统像差。

优选的，所述标准眼模型全眼像差通过所述标准眼模型的瞳孔放缩比例因子a、光束孔径离心因子b、视场矢量h、zernike项对应系数k进行定量表述，其中：

r为共轴zernike自由曲面的归一化半径。

优选的，所述步骤s2还包括：

在获取由上述临床标准zernike多项式表示的所述临床患者全眼像差系数后,将所述临床标准zernike多项式转换为矢量像差理论使用的fringezernike多项式，获取屈光手术后的术后患者全眼像差系数。

优选的，所述步骤s4包括如下步骤：

s4.1获取所述术后患者眼模型的调制传递函数，求取调制传递函数曲线下的归一化面积；将所述归一化面积作为屈光手术后视觉质量评价的量化指标；

s4.2获取所述术后患者眼模型的点扩散函数分析图；

s4.3获取所述术后患者眼模型的二维图像仿真分析图；

s4.4基于所述调制传递函数及所述归一化面积、所述点扩散函数分析图、所述二维图像仿真分析图，完成患者屈光手术后视觉质量综合的客观评价。

优选的，所述步骤s4.1还包括：

s4.1.1通过光学软件，得到所述术后患者眼模型调制传递函数的曲线参数；

s4.1.2利用梯形积分法求取调制传递函数曲线下归一化面积。

优选的，利用梯形积分法求取调制传递函数曲线下归一化面积的过程为：将调制传递函数曲线f(x)沿着离散等间隔空间频率xm(m＝1,2,...,n)进行积分，调制传递函数曲线下归一化面积s的公式为：

其中，δh为采样空间频率间隔，smax＝nδh，为调制传递函数曲线最大归一化面积。

本发明还提供了基于矢量像差理论的屈光手术后视觉质量评价系统，所述系统包括：标准眼模块、患者临床眼模块、患者眼模块、指标模块；

所述标准眼模块用于使用zernike自由曲面表征所述标准眼模型的角膜前表面，对所述标准眼模型全眼像差进行定量表述，建立标准眼模型角膜改变量与标准眼模型全眼像差的定量关系；

所述患者临床眼模块用于获取采用临床标准zernike多项式表示的临床患者全眼像差系数，根据所述标准眼模块建立的定量关系，得到患者眼模型的角膜改变量及对应的zernike系数；

所述患者眼模块用于根据所述患者临床眼模块得到的zernike系数，建立术后患者眼模型；

所述指标模块用于根据所述患者眼模块建立的所述术后患者眼模型，得到调制传递函数曲线以及曲线下的归一化面积、点扩散函数分析图和二维图像仿真分析图。

本发明的有益效果为：

本发明基于矢量像差理论，从光学经典像质评价的角度出发，通过标准眼模型光线追迹参数和共轴zernike自由曲面像差，建立标准眼模型角膜改变量与标准眼模型全眼像差的定量关系，再通过患者眼模型角膜改变量，建立术后患者眼模型，最后基于调制传递函数曲线下的归一化面积、点扩散函数分析图和二维图像仿真分析图，对屈光手术后患者视觉质量进行综合的客观评价。本发明有效克服了因过多的主观因素影响导致的无法准确评价患者术后视觉质量改善程度，以及无法明确手术是否达到预期效果等问题，对分析患者术后视觉质量具有重要的意义和广阔的应用前景。

附图说明

为了更清楚地说明本发明的技术方案，下面对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例的方法流程图；

图2为本发明实施例中标准眼模型示意图；

图3为本发明实施例中标准眼模型角膜改变量与标准眼模型全眼像差值分布；

图4为本发明实施例中临床全眼像差分布；

图5为本发明实施例中临床标准zernike多项式与fringezernike多项式的对应关系；

图6为本发明实施例中反向求解出的眼模型角膜改变量对应的fringezernike系数；

图7为本发明实施例中患者眼模型的调制传递函数mtf分析图；

图8为本发明实施例中梯形积分法原理图；

图9为本发明实施例中术后患者眼模型的调制传递函数mtf参数以及求解出的归一化mtf曲线下面积；

图10为本发明实施例中术后患者眼模型的点扩散函数psf分析图；

图11为本发明实施例中术后患者眼模型的二维图像仿真分析图；

图12为本发明实施例系统组成示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。

本实施例以一位接收飞秒激光辅助的准分子激光原位角膜磨镶术(femtosecondlaser-assistedlaserinsitukeratomileusis,fs-lasik)的患者为例，说明术后视觉质量评价的过程。

在本实施例中，采用的光学软件为codev。

如图1所示，本发明提供基于矢量像差理论的屈光手术后视觉质量评价方法，包括如下步骤：

s1.建立标准眼模型，采用zernike自由曲面表征所述标准眼模型的角膜前表面，对所述标准眼模型全眼像差进行定量表述，建立标准眼模型角膜改变量与标准眼模型全眼像差的定量关系；

在光学软件codev中，按照表1所示的结构参数建立标准眼模型。

表1

对于如图2所示的标准人眼模型，对其复杂的角膜前表面，采用由zernike多项式表示的自由曲面(zernike自由曲面)进行表征，zernike自由曲面的面型表达式如下，

其中，表示二次曲面，c和k分别为二次曲面的顶点曲率和圆锥系数，为zernike多项式各项子项的叠加式，其中ki代表zernike多项式系数，zi表示第i项zernike多项式，i对应不同波前像差，ρ、分别为zernike多项式的径向坐标和方位角，z(x,y)为自由曲面在坐标(x,y)处的矢高。

在本实施例中，利用codev光学软件对标准眼模型进行真实光线追迹，获取边缘光线在角膜前表面的入射高度h为3mm，边缘视场主光线在角膜前表面的入射高度为0，结合h与和共轴zernike自由曲面像差，对标准眼模型角膜改变量产生的眼模型全眼像差定量表述，得到标准眼模型角膜改变量与产生的标准眼模型全眼像差的定量关系，如图3所示。

其中，标准眼模型全眼像差为标准眼模型的角膜发生改变后产生的对应像差；共轴zernike自由曲面像差，为在共轴zernike自由曲面光学系统的非光阑处，其zernike项对应系数k产生的光学系统像差。

图3中，为眼模型视场矢量，a和b分别为标准眼模型的瞳孔放缩比例因子和光束孔径离心因子，标准眼模型全眼像差通过所述标准眼模型的瞳孔放缩比例因子a、光束孔径离心因子b、视场矢量h、zernike项对应系数k进行定量表述，其中：

r为自由曲面的归一化半径。

定量关系建立过程具体为：设置zernike自由曲面的归一化半径为3mm，则根据式(2)和式(3)可得图3中的a和b的大小分别为1和0，带入共轴zernike自由曲面光学系统非光阑处自由曲面的zernike项产生的光学系统像差表达式中，得到眼模型角膜改变量与产生的眼模型全眼像差的定量关系。

s2.获取屈光手术后的临床患者全眼像差系数，所述临床患者全眼像差系数采用临床标准zernike多项式表示；

在本实施例中，通过波前像差测量获取由临床标准zernike多项式表示的临床患者全眼像差系数keye，全眼波前像差数据的测量与收集采用美国强生idesign波前像差测量系统完成，瞳孔直径均标化为6mm。临床认为高阶像差中对于视觉质量影响较大的垂直彗差、水平彗差和初级球差的全眼像差分布如图4所示。

根据图4中的全眼像差，将临床标准zernike多项式转换为矢量像差理论利用的fringezernike多项式，从而获取屈光手术术后患者的全眼fringezernike像差系数kfringe。两种zernike多项式的对应关系如图5所示。在本实施例中，这里的像差类型为临床认为对视觉质量最具影响的离焦、像散、初级彗差、初级球差和三叶像差。

本实施例以初级球差为例，说明临床标准zernike多项式与fringezernike多项式之间的转换方法。临床标准zernike项和fringezernike多项式表示的初级球差表达式分别为

zfringesphere＝6ρ⁴-6ρ²+1………………(5)

其中ρ为光瞳径向坐标，因此为保证总测量值不变，即，

weyesphere＝wfringesphere…………………(6)

其中weyesphere和wfringesphere分别表示初级球差在临床标准zernike多项式和fringezernike多项式表征下的总测量值，他们可表示为以下形式，

weyesphere＝keyespherezeyesphere……………(7)

wfringesphere＝kfringespherezfringesphere………(8)

其中keyesphere和kfringesphere分别为临床测量的临床标准zernike初级球差系数和对应的fringezernike初级球差系数，因此根据式(6)，在转换过程中，keyesphere与kfringesphere之间需满足以下关系，

keyespherezeyesphere＝kfringespherezfringesphere

因此可根据临床测量得到的初级球差系数keyesphere计算得出其对应的fringezernike初级球差系数kfringesphere。

利用图3所示的标准眼模型角膜改变量与标准眼模型全眼像差的定量关系，按照图4所示，反向求解需要产生患者全眼像差wfringe的患者眼模型角膜改变量，及对应的zernike系数ki，i代表第i项zernike系数，计算结果如图6所示。本实施例以在标准眼模型中心视场(零视场)处测量得到的彗差为例，说明对角膜改变量对应zernike系数ki的反向求解过程。根据图3可知，对应彗差的角膜改变量对应zernike系数在零视场产生的彗差系数为，

而通过临床测量手段已经获得患者眼的彗差系数同时，缩放因子a可通过光线追迹获得，因此可通过反向求解式(10)，得到角膜改变量对应zernike系数

从而可得对应水平彗差和垂直彗差的角膜改变量zernike系数k7和k8分别为，

s3.利用步骤s2求解得到的zernike系数ki表征标准眼模型的角膜前表面，即将ki带入式(1)的zernike多项式系数中，从而获得含有图5中实际波前像差分布wfringe的术后患者眼模型；

s4.使用调制传递函数、点扩散函数和二维图像仿真，通过所述术后患者眼模型完成患者手术后视觉质量综合的客观评价。

在本实施例中，通过codev光学软件，利用经典的像质评价功能对患者眼模型进行调制传递函数mtf分析，患者眼模型的mtf分析结果如图7所示。利用软件导出mtf曲线的离散坐标值，采用梯形积分法求取mtf曲线下归一化面积，梯形积分法原理如图8所示，将mtf曲线f(x)沿着离散等间隔空间频率xm(m＝1,2,...,n)进行积分，得到mtf曲线下归一化面积为：

其中，δh为采样空间频率间隔，smax＝nδh，为mtf曲线最大归一化面积。本实施例以s作为患者视觉质量评价指标，并规定mtf曲线下面积较大的患者眼具有较好的光学成像质量。mtf曲线参数以及求得的mtf下归一化面积如图9所示。同时对患者眼模型进行点扩散函数psf和二维图像仿真分析，患者眼模型的psf和二维图像仿真分析结果分别如图10和图11所示。根据归一化面积s、点扩散函数psf分析图、二维图像仿真分析图，对屈光手术术后患者视觉质量进行综合的客观评价。

如图12所示，本发明基于矢量像差理论的屈光手术后视觉质量评价系统包括标准眼模块、患者临床眼模块、患者眼模块、指标模块；

所述标准眼模块用于使用zernike自由曲面表征所述标准眼模型的角膜前表面，对所述标准眼模型全眼像差进行定量表述，建立标准眼模型角膜改变量与标准眼模型全眼像差的定量关系；

所述患者临床眼模块用于获取采用临床标准zernike多项式表示的临床患者全眼像差系数，根据所述标准眼模块建立的定量关系，得到患者眼模型的角膜改变量及对应的zernike系数；

所述患者眼模块用于根据所述患者临床眼模块得到的zernike系数，建立术后患者眼模型；

所述指标模块用于根据所述患者眼模块建立的所述术后患者眼模型，得到调制传递函数曲线下的归一化面积、点扩散函数分析图和二维图像仿真分析图。

以上所述的实施例仅是对本发明优选方式进行的描述，并非对本发明的范围进行限定，在不脱离本发明设计精神的前提下，本领域普通技术人员对本发明的技术方案做出的各种变形和改进，均应落入本发明权利要求书确定的保护范围内。