CN106983491A - 角膜地形诊断的角膜像散分析方法 - Google Patents
角膜地形诊断的角膜像散分析方法 Download PDFInfo
- Publication number
- CN106983491A CN106983491A CN201710204358.XA CN201710204358A CN106983491A CN 106983491 A CN106983491 A CN 106983491A CN 201710204358 A CN201710204358 A CN 201710204358A CN 106983491 A CN106983491 A CN 106983491A
- Authority
- CN
- China
- Prior art keywords
- cornea
- astigmatism
- corneal
- analysis method
- corneal topography
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
- A61B3/107—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for determining the shape or measuring the curvature of the cornea
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Ophthalmology & Optometry (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Eye Examination Apparatus (AREA)
Abstract
一种角膜地形诊断的角膜像散分析方法,根据眼前节生物测量设备获取的角膜离散点高度数据,利用高次多项式对测量数据进行拟合并精确重构角膜面型,通过求解角膜曲面的一阶、二阶偏导并结合微分几何原理,计算出角膜各点的两个主曲率及其正交方位,将主曲率的差异转化为角膜屈光力差异并作为角膜像散,以主曲率的方位表示像散的轴向。本发明涉及的分析方法适用于任意的连续的角膜地形,包括病变的或术后的角膜,既可以计算整体角膜的规则散光,又能定量评估局部角膜的非规则散光,有助于角膜疾病的临床诊断及诊疗后的角膜屈光评估。
Description
技术领域:
本发明涉及角膜地形与光学性能技术领域,具体涉及一种角膜地形诊断的角膜像散分析方法。
背景技术:
人眼的屈光不正包括近视、远视与散光,其中散光又根据角膜主子午线的最大和最小曲率经线是否正交,分为规则性散光与非规则性散光。角膜屈光力约占全眼屈光系统的70%,角膜屈光力主要取决于角膜前表面。角膜疾病或者角膜屈光手术都有可能引起角膜形态变化,在病变或者手术区域造成局部的非规则散光,严重的将会影响视觉质量。因此,眼科医生需要利用眼前节生物测量设备检测的角膜地形数据,对角膜形态与成像质量进行针对性的判断与分析。临床上的角膜地形图可以为眼科医生提供角膜屈光力分布和角膜整体的规则散光,无法提供角膜的像散分布及其对应的轴向信息,不适合进行非规则散光的定量评估。
角膜像散是由角膜曲面的两个主曲率之差来决定的,角膜表面任意一点的曲面位置上都存在两个主曲率,它们分别是过曲面该点法线的主平面与角膜相交形成的最大曲率与最小曲率,且两者方向相互垂直,主曲率的大小与方向通过微分几何的数学理论进行求解,因此,角膜曲面各点位置上的主曲率之差和曲率方向可以全面描述角膜表面的像散分布情况,像散分布的不规律性直观地反映出角膜的非规则散光特性。
发明内容:
为了解决现有技术的缺陷,本发明提出了一种角膜地形诊断的角膜像散分析方法,对角膜面型的表达式进行二阶微分运算,根据微分几何计算方法求解角膜表面各点的两个主曲率,转化为角膜的像散分布及其方位角信息,可以用于角膜规则散光与非规则散光分析。本发明涉及的方法适用于任意的连续的角膜形态,包括病变的以及术后的角膜,可以对整体角膜的规则散光和局部角膜的非规则散光进行评估,有助于角膜疾病的临床诊断及角膜屈光手术的效果评估。
本发明采用的技术解决方案是:一种角膜地形诊断的角膜像散分析方法,所述的角膜像散分析方法包括以下步骤:
(1)利用高次多项式对角膜离散点高度数据进行拟合并精确重构角膜面型;
(2)通过求解角膜曲面的一阶、二阶偏导并结合微分几何原理获得角膜各点的两个主曲率及其正交方位;
(3)将主曲率转化为对应的角膜屈光力与角膜像散,以角膜像散及其轴向表示角膜的非规则散光;
(4)将角膜屈光力转化为屈光力矢量;
(5)在孔径范围内,计算屈光力矢量中各个分量的加权均值得到角膜的规则散光数值。
所述的Zernike多项式的径向阶数要超过16阶,角膜面型的拟合误差RMS值要小于0.35μm。
所述的角膜面型各点处的两个主曲率是指该点的法线面与角膜曲面相交所形成的截面曲线中最大曲率κ1和最小曲率κ2,所述的最大曲率κ1和最小曲率κ2方向正交。
所述的主曲率是通过求解以下一元二次方程来实现:
(EG-F2)κ2-(EN+GL-2FM)κ+(LN-M2)=0;其中,Sx、Sy是角膜曲面的一阶偏导,Sxx、Syy和Sxy是角膜曲面的二阶偏导,
主曲率所在方位为
所述的主曲率对应的角膜屈光力为
Pi(x,y)=1000(nc-1)×κi(x,y),i=1,2;其nc是角膜折射率,角膜像散定义为A(x,y)=P2(x,y)-P1(x,y),轴向取α1。
所述的角膜屈光力转化为屈光力矢量[M,J0,J45],即
其中M(x,y)是等效球镜屈光力,J0(x,y)是0度角散光分量,J45(x,y)是45度角散光分量。
在孔径范围内,计算屈光力矢量[M,J0,J45]中各分量M(x,y)、J0(x,y)和J45(x,y)的加权均值,得到平均等效球镜屈光力Mave、平均0度角散光J0ave和平均45度角散光J45ave,其中J0ave和J45ave代表角膜的规则散光。
本发明得到的有益效果是:本发明提供一种角膜地形诊断的角膜像散分析方法,对角膜面型的表达式进行二阶微分运算,根据微分几何计算方法求解角膜表面各点的两个主曲率,转化为角膜的像散分布及其方位角信息,可以用于角膜规则散光与非规则散光分析。本发明涉及的方法适用于任意的连续的角膜形态,包括病变的以及术后的角膜,可以对整体角膜的规则散光和局部角膜的非规则散光进行评估,有助于角膜疾病的临床诊断及角膜屈光手术的效果评估。
附图说明
图1为本发明的工作流程图。
图2为本发明所述的多项式拟合并重构的角膜面型高度图。
图3为根据本发明所述的方法获得的普通人眼角膜像散及其轴向分布图。
图4为根据本发明所述的方法获得的角膜塑形术后角膜像散及其轴向分布图。
图5为基于角膜像散分布图提取角膜塑形术的中央治疗区域。
图6为根据本发明所述的方法获得的未经手术治疗的圆锥角膜像散及其轴向分布图。
具体实施方式:
实施例1
普通人眼角膜像散的计算与分析:
(1)采集普通人眼的眼前节生物测量数据,获得全角膜高度分布,由角膜地形定位测量中心与瞳孔测量中心坐标获得瞳孔偏心量;
(2)定义包含角膜所有有效测量点在内的最大孔径圆半径;
(3)采用最小二乘法,以Zernike多项式拟合孔径圆半径内所有的角膜测量点高度数据,Zernike径向阶数为20阶;
(4)根据曲面拟合获得的Zernike多项式系数重构出连续、完整的角膜面型(见图2),面型的采样间距为0.12mm*0.12mm;
(5)求解Zernike多项式的一阶、二阶偏导表达式,结合Zernike多项式系数计算出角膜面型的一阶、二阶偏导数;
(6)根据主曲率计算公式,求解角膜面型上任意点的两个主曲率(即最大曲率与最小曲率)及其所在方位;
(7)由主曲率之差以及角膜界面两侧折射率之差计算角膜的像散分布,在角膜地形测量的有效范围内绘制角膜像散及其方位角分布图(见图3),像散及其轴向的非均匀分布反映了普通人眼角膜的非规则散光。
(8)由主曲率对应的角膜屈光力以及角膜像散计算角膜的屈光力矢量分布,并在一定的孔径范围内计算屈光力矢量中各个分量的加权均值,得到角膜的规则散光数值。
实施例2
角膜塑形镜戴后角膜像散的计算与分析:
(1)采集配戴角膜塑形镜后患者的眼前节生物测量数据,获得全角膜高度分布,由角膜地形定位测量中心与瞳孔测量中心坐标获得瞳孔偏心量;
(2)定义包含角膜所有有效测量点在内的最大的孔径圆半径;
(3)采用最小二乘法,以Zernike多项式拟合孔径圆半径内所有的角膜测量点高度数据,Zernike径向阶数为20阶;
(4)根据曲面拟合获得的Zernike多项式系数重构出连续、完整的角膜面型,面型的采样间距为0.12mm*0.12mm;
(5)求解Zernike多项式的一阶、二阶偏导表达式,结合Zernike多项式系数计算出角膜面型的一阶、二阶偏导数;
(6)根据主曲率计算公式,求解角膜面型上任意点的两个主曲率(即最大曲率与最小曲率)及其所在方位;
(7)将主曲率转化为对应的角膜屈光力与角膜像散,以角膜像散及其轴向表示角膜的非规则散光,在角膜地形测量的有效范围内绘制角膜像散及其方位角分布图(见图4)。
(8)由主曲率对应的角膜屈光力以及角膜像散计算角膜的屈光力矢量分布,并在一定的孔径范围内计算屈光力矢量中的各个分量的加权均值。
(9)角膜塑形镜后表面是由多弧段组成不同工作区域,根据图4可知不同区域压迫角膜形成的分区情况以及像散分布,像散及其轴向的非均匀分布反映了塑形后的角膜的非规则散光,由像散极值划分出角膜塑形的中央治疗区域及其相对于瞳孔的偏心量即适配偏心(见图5)。
实施例3
未经治疗的圆锥角膜像散的计算与分析:
(1)选取圆锥角膜病例,采集圆锥角膜病例的眼前节生物测量数据,获得全角膜高度分布,由角膜地形定位测量中心与瞳孔测量中心坐标获得瞳孔偏心量;
(2)定义包含角膜所有有效测量点在内的最大的孔径圆半径;
(3)采用最小二乘法,以Zernike多项式拟合孔径圆半径内所有的角膜测量点高度数据,Zernike径向阶数为20阶;
(4)根据曲面拟合获得的Zernike多项式系数重构出连续、完整的角膜面型,面型的采样间距为0.12mm*0.12mm;
(5)求解Zernike多项式的一阶、二阶偏导表达式,结合Zernike多项式系数计算出角膜面型的一阶、二阶偏导数;
(6)根据主曲率计算公式,求解角膜面型上任意点的两个主曲率(即最大曲率与最小曲率)及其所在方位;
(7)由主曲率之差以及角膜界面两侧折射率之差计算角膜的像散分布,在角膜地形测量的有效范围内绘制角膜像散及其方位角分布图(见图6),像散及其轴向的非均匀分布反映了圆锥角膜的非规则散光。
(8)由主曲率对应的角膜屈光力以及角膜像散计算角膜的屈光力矢量分布,并在一定的孔径范围内计算屈光力矢量中各个分量的加权均值,得到圆锥角膜的规则散光数值。
具体实施方式只用于对本发明进行进一步说明,不能作为对本发明保护范围的限定,同时该领域的技术人员根据上述发明的内容对本发明作出一些非本质的改进和调整,都位于本发明的保护范围内,本发明的保护范围以权利要求书为准。
Claims (8)
1.一种角膜地形诊断的角膜像散分析方法,其特征在于:所述的角膜像散分析方法包括以下步骤:
(1)利用高次多项式对角膜离散点高度数据进行拟合并精确重构角膜面型;
(2)通过求解角膜曲面的一阶、二阶偏导并结合微分几何原理获得角膜各点的两个主曲率及其正交方位;
(3)将主曲率转化为对应的角膜屈光力与角膜像散,以角膜像散及其轴向表示角膜的非规则散光;
(4)将角膜屈光力转化为屈光力矢量;
(5)在孔径范围内,计算屈光力矢量中各个分量的加权均值得到角膜的规则散光数值。
2.根据权利要求1所述的角膜地形诊断的角膜像散分析方法,其特征在于:所述的高次多项式为Zernike多项式,多项式重构的角膜面型为:
其中是卡笛尔坐标系下的多项式表达式,为多项式各项的系数。
3.根据权利要求2所述的角膜地形诊断的角膜像散分析方法,其特征在于:所述的Zernike多项式的径向阶数要超过16阶,角膜面型的拟合误差RMS值要小于0.35μm。
4.根据权利要求1所述的角膜地形诊断的角膜像散分析方法,其特征在于:所述的角膜面型各点处的两个主曲率是指该点的法线面与角膜曲面相交所形成的截面曲线中最大曲率κ1和最小曲率κ2,所述的最大曲率κ1和最小曲率κ2方向正交。
5.根据权利要求4所述的角膜地形诊断的角膜像散分析方法,其特征在于:所述的主曲率是通过求解以下一元二次方程来实现:
(EG-F2)κ2-(EN+GL-2FM)κ+(LN-M2)=0;其中,Sx、Sy是角膜曲面的一阶偏导,Sxx、Syy和Sxy是角膜曲面的二阶偏导,
主曲率所在方位为
6.根据权利要求1所述的角膜地形诊断的角膜像散分析方法,具特征在于:所述的主曲率对应的角膜屈光力为
Pi(x,y)=1000(nc-1)×κi(x,y),i=1,2;其nc是角膜折射率,角膜像散定义为A(x,y)=P2(x,y)-P1(x,y),轴向取α1。
7.根据权利要求6所述的角膜地形诊断的角膜像散分析方法,其特征在于,所述的角膜屈光力转化为屈光力矢量[M,J0,J45],即
其中M(x,y)是等效球镜屈光力,J0(x,y)是0度角散光分量,J45(x,y)是45度角散光分量。
8.根据权利要求7所述的角膜地形诊断的角膜像散分析方法,其特征在于,在孔径范围内,计算屈光力矢量[M,J0,J45]中各分量M(x,y)、J0(x,y)和J45(x,y)的加权均值,得到平均等效球镜屈光力Mave、平均0度角散光J0ave和平均45度角散光J45ave,其中J0ave和J45ave代表角膜的规则散光。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710204358.XA CN106983491B (zh) | 2017-03-30 | 2017-03-30 | 角膜地形诊断的角膜像散分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710204358.XA CN106983491B (zh) | 2017-03-30 | 2017-03-30 | 角膜地形诊断的角膜像散分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106983491A true CN106983491A (zh) | 2017-07-28 |
CN106983491B CN106983491B (zh) | 2019-04-02 |
Family
ID=59412003
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710204358.XA Active CN106983491B (zh) | 2017-03-30 | 2017-03-30 | 角膜地形诊断的角膜像散分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106983491B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110797122A (zh) * | 2019-11-20 | 2020-02-14 | 杭州明视康眼科医院有限公司 | 一种用于角膜屈光矫正的定量调整纵向球差的方法 |
CN112587087A (zh) * | 2020-12-21 | 2021-04-02 | 上海美沃精密仪器股份有限公司 | 一种定位人眼角膜屈光地形图平陡k轴位的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080165324A1 (en) * | 2006-12-31 | 2008-07-10 | Lindacher Joseph M | Generalized presbyopic correction methodology |
CN101957502A (zh) * | 2010-08-31 | 2011-01-26 | 吉林大学 | 个性化角膜接触镜的设计方法 |
CN102508944A (zh) * | 2011-10-10 | 2012-06-20 | 沈阳理工大学 | 一种角膜面形拟合方法 |
CN102566085A (zh) * | 2012-03-20 | 2012-07-11 | 天津宇光光学有限公司 | 基于波前技术的非球面眼镜的设计方法 |
CN104321008A (zh) * | 2012-02-24 | 2015-01-28 | 诺埃尔·奥米·阿尔平斯 | 用于角膜散光分析和向量规划治疗的地形半子午线参数的评估 |
CN105608314A (zh) * | 2015-12-17 | 2016-05-25 | 深圳市斯尔顿科技有限公司 | 一种角膜屈光度的计算方法及系统 |
-
2017
- 2017-03-30 CN CN201710204358.XA patent/CN106983491B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080165324A1 (en) * | 2006-12-31 | 2008-07-10 | Lindacher Joseph M | Generalized presbyopic correction methodology |
CN101957502A (zh) * | 2010-08-31 | 2011-01-26 | 吉林大学 | 个性化角膜接触镜的设计方法 |
CN102508944A (zh) * | 2011-10-10 | 2012-06-20 | 沈阳理工大学 | 一种角膜面形拟合方法 |
CN104321008A (zh) * | 2012-02-24 | 2015-01-28 | 诺埃尔·奥米·阿尔平斯 | 用于角膜散光分析和向量规划治疗的地形半子午线参数的评估 |
CN102566085A (zh) * | 2012-03-20 | 2012-07-11 | 天津宇光光学有限公司 | 基于波前技术的非球面眼镜的设计方法 |
CN105608314A (zh) * | 2015-12-17 | 2016-05-25 | 深圳市斯尔顿科技有限公司 | 一种角膜屈光度的计算方法及系统 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110797122A (zh) * | 2019-11-20 | 2020-02-14 | 杭州明视康眼科医院有限公司 | 一种用于角膜屈光矫正的定量调整纵向球差的方法 |
CN110797122B (zh) * | 2019-11-20 | 2024-04-16 | 杭州明视康眼科医院有限公司 | 一种用于角膜屈光矫正的定量调整纵向球差的方法 |
CN112587087A (zh) * | 2020-12-21 | 2021-04-02 | 上海美沃精密仪器股份有限公司 | 一种定位人眼角膜屈光地形图平陡k轴位的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106983491B (zh) | 2019-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Guillon et al. | Corneal topography: a clinical model | |
US8708485B2 (en) | Method for providing a spectacle ophthalmic lens by calculating or selecting a design | |
Kosaki et al. | Magnitude and orientation of Zernike terms in patients with keratoconus | |
Atchison | Optical models for human myopic eyes | |
US10973407B2 (en) | Method for determining an ophthalmic lens having unwanted astigmatism | |
Budak et al. | Evaluation of relationships among refractive and topographic parameters | |
US11561413B2 (en) | Population of an eye model using measurement data in order to optimize spectacle lenses | |
McDonnell et al. | Topographic analysis and visual acuity after radial keratotomy | |
CN110785113B (zh) | 用于确定受试者的眼睛的眼睛转动中心的位置的方法以及相关装置 | |
Tang et al. | Characteristics of keratoconus and pellucid marginal degeneration in mean curvature maps | |
Tan et al. | How keratoconus influences optical performance of the eye | |
Suzaki et al. | Visual performance and optical quality of standardized asymmetric soft contact lenses in patients with keratoconus | |
Li et al. | Contribution of the retinal contour to the peripheral optics of human eye | |
CN106983491A (zh) | 角膜地形诊断的角膜像散分析方法 | |
Preussner et al. | Corneal model | |
El Hage et al. | Empirical advanced orthokeratology through corneal topography: the University of Houston clinical study | |
Asher et al. | Etiology of Keratoconus: proposed biomechanical pathogenesis | |
CN113197543B (zh) | 基于矢量像差理论的屈光手术后视觉质量评价方法和系统 | |
US20180210228A1 (en) | Population of an eye model for optimizing spectacle lenses with measurement data | |
Barcik et al. | Optical performance of the eye with progressive addition lens correction | |
Collins et al. | Tolerance to spherical aberration induced by rigid contact lenses | |
Rojo et al. | Generalized ray tracing method for the calculation of the peripheral refraction induced by an ophthalmic lens | |
Langenbucher et al. | Three-axis ellipsoidal fitting of videokeratoscopic height data after penetrating keratoplasty | |
CN112530596A (zh) | 一种近视发生风险评估装置 | |
Talu et al. | An overview on mathematical models of human corneal surface |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |