论文部分内容阅读
目的:
分别使用构造正则方程、奇异值分解、Householder变换三种解法解决角膜像差波前Zernike重建线性模型的最小二乘解问题,比较重建精度、耗时以及可靠性,为应用于实时矫正眼球像差的自适应光学系统中波前重建最优算法的选择提供参考及理论依据。
方法:
使用Pentacam角膜地形图仪获取20只正常眼的角膜地形图前表面高度数据。数据范围为中央6mm直径的瞳孔区,根据Pentacam提供的参数生成最佳拟合球面(BFS,Best Fit Sphere),二者之差即为模拟的角膜前表面的像差数据。数据采样率为100μm(仪器设置),编制采样程序对其分别使用100μm,300μm,500μm的分辨率进行采样,对每个波前像差数据矩阵计算梯度,求得x方向及y方向的一阶偏导,此即为波前像差的斜率数据(即为模拟夏克-哈特曼传感器采样)。获得波前斜率数据后使用zernike多项式(分别使用1至130个模式项数)重建波前数据,对于上述三种采样率,都分别使用构造正则方程法、奇异值分解与Householder变换三种算法获得重建波面的模式系数。根据获得的模式系数获得重建波面,使用其与原始波面的差值的残差均方根(root mean square,rms)考察重建精度,同时使用MATLAB函数tic、tio考察重建耗时,并计算线性模型的条件数、观察模式系数矩阵解的合理性以考察算法的可靠性。
结果:
1.重建精度:三种采样率时,Householder变换与奇异值分解的表现是一致的,而构造正则方程法则在分辨率为500μm时表现较差,主要为高阶(zernike模式数K>88)时rms值明显的不稳定(表现为无规律的大幅度的波动)。
2.重建耗时:三种采样率时,随着阶数的增加,奇异值分解的耗时增加明显高于Householder变换与构造正则方程法,而另二者的差异性表现得并不直观。同时,采样率为500μm或300μm时,奇异值分解法在K>9时即与另二者出现统计学意义上的差异(P<0.05),而Householder变换与构造正则方程法在K>50才出现差异;采样率为100μm时,奇异值分解法在K>11时即与另二者出现统计学意义上的差异(P<0.05),而Householder变换与构造正则方程法在K>121才出现差异。
3.算法可靠性:采样率越高,使用的zernike多项式的阶数越低,结果越可靠。同时构造正则方程法在采样率较低时的计算不稳定性明显比另二者严重。
结论:
1.本研究模拟夏克-哈特曼传感器的斜率重建过程,同时选择角膜前表面像差代表人眼真实像差,“过程模拟”结合“数据可比性”,同时解决了仪器测人眼像差没有真实波面数据,仪器测角膜像差不符合像差测量原理的矛盾。
2.考察了三种经典解法在重建真实人眼波前像差数据时的优劣性,主要以重建精度和重建耗时两方面作为评价指标;而Householder变换在这两方面优于构造正则方程与奇异值分解法。
3.通过考察模式系数矩阵了解线性模型解的合理性,避免了只使用重建精度考察重建效能的不足。越高的采样率能得到较合理的解(即zernike模式系数),而低采样率则会在高阶处出现异常解,使得算法的适用范围下降。
4. Householder变换在精确以及高效两方面优于另二者,同时高采样率时,三种解法的可靠性几乎相同。而低采样率时,Householder变换仍有相对较稳定的表现,此结果为应用于眼科医学的自适应光学系统中波前重建最优算法的选择提供了参考及理论依据。