CN108814544B - 一种眼前节形态分析方法 - Google Patents
一种眼前节形态分析方法 Download PDFInfo
- Publication number
- CN108814544B CN108814544B CN201810376739.0A CN201810376739A CN108814544B CN 108814544 B CN108814544 B CN 108814544B CN 201810376739 A CN201810376739 A CN 201810376739A CN 108814544 B CN108814544 B CN 108814544B
- Authority
- CN
- China
- Prior art keywords
- anterior segment
- scanning
- data
- crystal
- cornea
- 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.)
- Active
Links
Images
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/14—Arrangements specially adapted for eye photography
-
- 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/117—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for examining the anterior chamber or the anterior chamber angle, e.g. gonioscopes
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
技术领域
本发明属于眼科医疗设备技术领域，涉及的是一种眼前节形态分析方法。
背景技术
医学眼解剖概念范畴中，眼前节包括角膜至晶状体的眼组织部分。近年来眼前节患病率呈逐年上升趋势，对角膜情况的检查显得尤为重要，同时晶体病变情况在对青光眼、白内障以及屈光视力问题的诊断中，也是重要的参考依据。
传统眼前节影像技术一般可以观察角膜至晶体前表面范围内的眼前节情况，无法覆盖到全部晶体区域。在眼前节形态分析方面，基于Pacido盘原理的角膜地形图是目前应用最为广泛的技术之一。该技术在临床使用中容易受到光线或泪液的影响，使得测量结果不准确，同时该技术只能获取角膜前表面形态数据，无法满足对于整个眼前节的评估。
另一种较为常见的是基于Scheimpflug照相机技术的剖面成像仪器，如OCULUS公司生产的pentacam眼前节分析诊断系统。该技术利用光学成像原理对眼剖面成像，可以得到角膜至晶体范围内的二维照片，再经三维重建得到眼前节整体数据。这种仪器存在光学畸变相差，需要进行非常复杂的图像矫正，为了获得三维数据采用超高速的面阵CCD，价格非常昂贵。同时，照明光源采用蓝光(475nm)，对人眼形成非常大的刺激，患者体验不好，也无法获得放松状态下眼前节的形态。
光学相干层析(OCT)技术发展迅速，以其快速、非侵入等优点广泛用于医疗研究当中。在专利CN106963337中，提出了一种实现大景深眼前节分析系统。该技术采用离轴扫描振镜，扫描振镜沿旋转轴旋转将入射光以不同反射角度射出，光束旋转角度为a，旋转轴与入射光轴在x扫描方向上的偏移量为t，由于扫描振镜的旋转角度较小，根据小角度弧长公式，因扫描振镜在不同角度而形成的光程差b可以近似认为b＝at，即在任意扫描位置，入射光在样品臂所走过的光程都不同于相邻扫描位置处所走过的光程，且入射光光程随扫描振镜旋转角度改变，呈单调性。因此，对样品进行采样时，离轴扫描振镜在相邻的扫描点引入近似固定的光程差，产生时域多普勒效应，即沿x方向加入了一个调制频率。将采集数据在复数域中进行解调，可以消除在实数域做傅里叶变换形成的伪像。高分辨率光谱仪设计光谱分辨率大于或等于25皮米，对应相干长度大于14毫米，对角膜前表面到晶体后表面范围内的结构同时成像。该系统只能进行影像观察，不能定量分析眼前节形态参数。
发明内容
本发明的目的是为了克服现有技术的不足，提供一种眼前节形态分析方法。
本发明解决其技术问题是采取以下技术方案实现的：
一种眼前节形态分析方法，该方法采用的数据采集系统包括监视模块，分光模块，接目镜，光学相干层析模块及计算机，对人眼进行二维振镜非接触式光学扫描，采集眼前节三维数据块，该方法包括步骤如下：
(1)选择扫描模式，扫描模式采用旋转式扫描，所使用的二维振镜驱动信号为：两方向均采用锯齿波驱动，并且保证时序上的一一对应，通过改变驱动信号幅度，使振镜扫描步长保持一定比例，即可实现以不同角度的直线扫描，计算旋转式扫描的扫描线角度间隔对应的驱动信号比例，完成扫描；
(2)数据采集，光学相干断层扫描在每个扫描点上一次得到深度方向的Z个数据点，二维扫描决定了扫描点数为X*Y个，最终得到的采样数据点总数为X*Y*Z个，监视模块同步采集帧同步图像，使用阈值分割的方法提取虹膜边缘和瞳孔边缘，得到位置反馈信息以及瞳孔直径和白对白直径；
(3)边缘提取，对深度方向Z个数据点进行梯度运算，使用最短路径法，得到眼前节轮廓对应的数据点，在扫描方向上重复该步骤，数据重组后得到眼前节边缘轮廓数据集及全眼前节边缘轮廓图像；
(4)粗略对准，活体数据采集无法避免微动带来的图像偏移，使用监视图像提供位置反馈信号，用以矫正扫描方向上的位置偏移，根据眼前节的结构特性，在已经提取的眼前节边缘轮廓数据集中选择房角顶点和角膜顶点作为每个剖面的对准参考点，通过数据空间映射，将每个数据剖面中的对准参考点对齐到同一平面中，即消除空间上的位置偏差；
(5)精细对准，在数据集中引入Zernike多项式，其表达式为：
公式中各物理量的含义及单位为，C(r，θ)是角膜表面的数学表达式，P是多项式中序列索引，Zp(r，θ)代表Zernike多项式中第p项，ap是该项的模型系数，ε是实验模型误差，r代表径向位置，θ代表旋转角度；
使用Zernike多项式分别对数据集中角膜前表面数据子集，角膜后表面数据子集，晶体前表面数据子集，晶体后表面数据子集进行曲面拟合，再进行等间隔重采样，曲面上的任一点都对应计算出一个微分斜率，通过该微分斜率的集合得到对应曲面地形图数据，即角膜前表面地形图，角膜后表面地形图，晶体前表面地形图，晶体后表面地形图；计算相邻曲面位置差，得到角膜厚度分布及晶体厚度分布，同时，经过对准匹配的三维数据块，用以提供三维空间内任意方向的眼前节断层剖面图像。
而且，所述步骤(1)中的旋转式扫描由栅格式扫描代替，即一个维度上振镜以等幅度锯齿波驱动，其幅度决定了扫描线的范围，另一个维度的驱动信号为阶梯式波形，即将扫描线平移一定距离进行下一次扫描，阶梯信号的幅值变化对应扫描线的间隔距离。
而且，所述步骤(5)中任意方向的眼前节断层剖面图像包括重建后的角膜三维图像，房角三维图像，虹膜三维图像，晶体三维图像，利用这些图像进行角膜屈光力评估，角膜手术预评估，不规则角膜病变评估，角膜波前相差分析，白内障观察及人工晶体手术评估。
而且，所述步骤(5)中角膜前表面地形图，角膜后表面地形图，晶体前表面地形图，晶体后表面地形图，角膜厚度分布，晶体厚度分布作为从角膜前表面至晶体后表面全部眼前节形态基础生物参数数据集，使用该数据集中一种或多种数据子集进行联合分析，进一步计算角膜变形，晶体变形及前房变形的高级形态变化。
本发明的优点和积极效果是：
1、本发明可同时获得眼前节断层图像，角膜前表面地形图，角膜后表面地形图，角膜厚度分布，晶体前表面地形图，晶体后表面地形图，晶体厚度分布，角膜三维图像，房角三维图像，虹膜三维图像，晶体三维图像。
2、本发明不存在空间光学像差，降低了算法复杂程度，保证了测量精度。
3、本发明为非侵入方式，无刺激，易于被患者接受。
附图说明
图1是大景深光学相干层析扫描系统结构示意图；
图2是本发明眼前节形态分析方法流程图；
图3是本发明中旋转式扫描示意图；
图4是本发明中栅格式扫描示意图；
图5是本发明中角膜和虹膜边缘示意图；
图6是本发明中全眼前节数据错位示意图；
图7是本发明中全眼前节示意图。
具体实施方式
以下结合附图对本发明实施做进一步详述，以下实施例只是描述性的，不是限定性的，不能以此限定本发明的保护范围。
一种眼前节形态分析方法，如图1所示，该方法采用的数据采集系统包括监视模块1，分光模块2，接目镜3，光学相干层析模块5及计算机6，对人眼4进行二维振镜非接触式光学扫描，采集眼前节三维数据块，如图2所示，该方法包括步骤如下：
(1)选择扫描模式，扫描模式采用旋转式扫描，所使用的二维振镜驱动信号为：两方向均采用锯齿波驱动，并且保证时序上的一一对应，通过改变驱动信号幅度，使振镜扫描步长保持一定比例，即可实现以不同角度的直线扫描，计算旋转式扫描的扫描线角度间隔对应的驱动信号比例，可完成如图3所示的扫描形式；其中，所述旋转式扫描可以由栅格式扫描代替，如图4所示，一个维度上振镜以等幅度锯齿波驱动，其幅度决定了扫描线的范围，另一个维度的驱动信号为阶梯式波形，即将扫描线平移一定距离进行下一次扫描，阶梯信号的幅值变化对应扫描线的间隔距离。
(2)数据采集，光学相干断层扫描在每个扫描点上可以一次得到深度方向的数据点Z个，二维扫描决定了扫描点数为X*Y个，最终得到的采样数据点总数为X*Y*Z个。监视模块同步采集帧同步图像，使用阈值分割的方法提取虹膜边缘7和瞳孔边缘8，得到位置反馈信息以及瞳孔直径和白对白直径，如图5所示；
(3)边缘提取，对深度方向Z个数据点进行梯度运算，使用最短路径法，得到眼前节轮廓对应的数据点，在扫描方向上重复该步骤，数据重组后得到眼前节边缘轮廓数据集及全眼前节边缘轮廓图像，如图5所示；
(4)粗略对准，活体数据采集无法避免微动带来的图像偏移，如图6所示，使用监视图像提供位置反馈信号，用以矫正扫描方向上的位置偏移，根据眼前节的结构特性，在已经提取的眼前节边缘轮廓数据集中选择房角顶点9和角膜顶点10，如图7所示，作为每个剖面的对准参考点，通过数据空间映射，将每个数据剖面中的对准参考点对齐到同一平面中，即消除空间上的位置偏差；
(5)精细对准，为了进一步消除边缘提取误差带来的数据对准不平滑，在数据集中引入Zernike多项式，其表达式为：
公式中各物理量的含义及单位为，C(r，θ)是角膜表面的数学表达式，P是多项式中序列索引，Zp(r，θ)代表Zernike多项式中第p项，ap是该项的模型系数，ε是实验模型误差，r代表径向位置，θ代表旋转角度；
Zernike多项式多用于构建曲面表面模型，经过泽尼克多项式矫正，降低数据集中存在的离散对齐不准确程度，然后分别对数据集中角膜前表面数据子集，角膜后表面数据子集，晶体前表面数据子集，晶体后表面数据子集进行曲面拟合，再进行等间隔重采样，曲面上的任一点都对应计算出一个微分斜率，该微分斜率的集合可得到对应曲面地形图数据，即角膜前表面地形图，角膜后表面地形图，晶体前表面地形图，晶体后表面地形图；计算相邻曲面位置差，得到角膜厚度分布及晶体厚度分布。
同时，经过对准匹配的三维数据块可以提供三维空间内任意方向的眼前节断层剖面图像，并且重建后的角膜三维图像，房角三维图像，虹膜三维图像，晶体三维图像可用于评估角膜屈光力，角膜手术预评估，不规则角膜病变评估，角膜波前相差分析，白内障观察，人工晶体手术评估。
在本发明的具体实施中，角膜前表面地形图，角膜后表面地形图，晶体前表面地形图，晶体后表面地形图，角膜厚度分布，晶体厚度分布作为从角膜前表面至晶体后表面全部眼前节形态基础生物参数数据集，使用该数据集中一种或多种数据子集进行联合分析，可以进一步计算角膜变形，晶体变形，前房变形等高级形态变化。
上述实施例只为说明本发明的技术构思及特点，其目的在于简要的阐述本发明的内容并据以实施，并不能限制本发明的保护范围。凡根据本发明的精神实质所做的等效变化或修饰都应涵盖在本发明的保护范围之内。
Claims (3)
1.一种眼前节形态分析方法，该方法采用的数据采集系统包括监视模块，分光模块，接目镜，光学相干层析模块及计算机，对人眼进行二维振镜非接触式光学扫描，采集眼前节三维数据块，该方法包括步骤如下：
(1)选择扫描模式，扫描模式采用旋转式扫描，所使用的二维振镜驱动信号为：两方向均采用锯齿波驱动，并且保证时序上的一一对应，通过改变驱动信号幅度，使振镜扫描步长保持一定比例，即可实现以不同角度的直线扫描，计算旋转式扫描的扫描线角度间隔对应的驱动信号比例，完成扫描；
(2)数据采集，光学相干断层扫描在每个扫描点上一次得到深度方向的Z个数据点，二维扫描决定了扫描点数为X*Y个，最终得到的采样数据点总数为X*Y*Z个，监视模块同步采集帧同步图像，使用阈值分割的方法提取虹膜边缘和瞳孔边缘，得到位置反馈信息以及瞳孔直径和白对白直径；
(3)边缘提取，对深度方向Z个数据点进行梯度运算，使用最短路径法，得到眼前节轮廓对应的数据点，在扫描方向上重复该步骤，数据重组后得到眼前节边缘轮廓数据集及全眼前节边缘轮廓图像；
(4)粗略对准，活体数据采集无法避免微动带来的图像偏移，使用监视图像提供位置反馈信号，用以矫正扫描方向上的位置偏移，根据眼前节的结构特性，在已经提取的眼前节边缘轮廓数据集中选择房角顶点和角膜顶点作为每个剖面的对准参考点，通过数据空间映射，将每个数据剖面中的对准参考点对齐到同一平面中，即消除空间上的位置偏差；
(5)精细对准，在数据集中引入Zernike多项式，其表达式为：
公式中各物理量的含义及单位为，C(r，θ)是角膜表面的数学表达式，P是多项式中序列索引，Zp(r，θ)代表Zernike多项式中第p项，αp是该项的模型系数，ε是实验模型误差，r代表径向位置，θ代表旋转角度；
使用Zernike多项式分别对数据集中角膜前表面数据子集，角膜后表面数据子集，晶体前表面数据子集，晶体后表面数据子集进行曲面拟合，再进行等间隔重采样，曲面上的任一点都对应计算出一个微分斜率，通过该微分斜率的集合得到对应曲面地形图数据，即角膜前表面地形图，角膜后表面地形图，晶体前表面地形图，晶体后表面地形图；计算相邻曲面位置差，得到角膜厚度分布及晶体厚度分布，同时，经过对准匹配的三维数据块，用以提供三维空间内任意方向的眼前节断层剖面图像。
2.根据权利要求1所述的一种眼前节形态分析方法，其特征在于：所述步骤(1)中的旋转式扫描由栅格式扫描代替，即一个维度上振镜以等幅度锯齿波驱动，其幅度决定了扫描线的范围，另一个维度的驱动信号为阶梯式波形，即将扫描线平移一定距离进行下一次扫描，阶梯信号的幅值变化对应扫描线的间隔距离。
3.根据权利要求1所述的一种眼前节形态分析方法，其特征在于：所述步骤(5)中角膜前表面地形图，角膜后表面地形图，晶体前表面地形图，晶体后表面地形图，角膜厚度分布，晶体厚度分布作为从角膜前表面至晶体后表面全部眼前节形态基础生物参数数据集，使用该全部眼前节形态基础生物参数数据集集中一种或多种数据子集进行联合分析，进一步计算角膜变形，晶体变形及前房变形的高级形态变化。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810376739.0A CN108814544B (zh) | 2018-04-25 | 2018-04-25 | 一种眼前节形态分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810376739.0A CN108814544B (zh) | 2018-04-25 | 2018-04-25 | 一种眼前节形态分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108814544A CN108814544A (zh) | 2018-11-16 |
| CN108814544B true CN108814544B (zh) | 2020-10-16 |
Family
ID=64154807
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810376739.0A Active CN108814544B (zh) | 2018-04-25 | 2018-04-25 | 一种眼前节形态分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108814544B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109919946B (zh) * | 2019-02-19 | 2021-04-20 | 温州医科大学 | 一种基于光学相干层析成像技术预测巩膜透氧型接触镜镜后泪液形态变化止点的方法 |
| CN110517219B (zh) * | 2019-04-01 | 2023-04-07 | 中山大学中山眼科中心 | 一种基于深度学习的角膜地形图判别方法及系统 |
| CN109998477B (zh) * | 2019-04-12 | 2021-12-10 | 复旦大学附属眼耳鼻喉科医院 | 一种高度近视白内障手术预后智能预判系统 |
| CN110236485B (zh) * | 2019-07-16 | 2022-02-18 | 天津市索维电子技术有限公司 | 一种用于视网膜地形图测定的设备及方法 |
| CN110717884B (zh) * | 2019-08-30 | 2022-02-22 | 温州医科大学 | 一种基于眼表结构变化一致性表达角膜不规则改变的方法 |
| CN111260610B (zh) * | 2020-01-08 | 2023-08-01 | 上海美沃精密仪器股份有限公司 | 一种基于房角开放距离曲线的眼前节状态表征方法及系统 |
| CN113989198B (zh) * | 2021-10-09 | 2022-04-29 | 中国医学科学院生物医学工程研究所 | 一种房角开合角度获取方法 |
| CN114041747A (zh) * | 2021-11-15 | 2022-02-15 | 邓宏伟 | 高度近视儿童眼容积的三维重建系统 |
| CN116342636B (zh) * | 2023-05-23 | 2023-12-29 | 广东麦特维逊医学研究发展有限公司 | 眼前节oct图像轮廓拟合方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2818252A1 (en) * | 2010-12-03 | 2012-06-07 | Jay Wei | Scanning and processing using optical coherence tomography |
| CN102283635B (zh) * | 2011-07-13 | 2016-03-09 | 广东福地新视野光电技术有限公司 | 双通道全眼光学相干层析成像系统及成像方法 |
| CN103654721B (zh) * | 2013-12-27 | 2016-06-08 | 深圳市斯尔顿科技有限公司 | 一种角膜顶点精确对准的方法 |
| CN103989453B (zh) * | 2014-06-03 | 2017-01-25 | 深圳市莫廷影像技术有限公司 | 一种多功能眼科测量装置及测试人眼不同部位的方法 |
| CN105231990B (zh) * | 2015-11-17 | 2017-04-19 | 深圳市亿领科技有限公司 | 基于oct三维成像分析角膜生物力学性能的装置及方法 |
| CN106963337B (zh) * | 2017-03-29 | 2018-04-24 | 天津市索维电子技术有限公司 | 一种实现大景深眼前节分析系统 |
-
2018
- 2018-04-25 CN CN201810376739.0A patent/CN108814544B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108814544A (zh) | 2018-11-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108814544B (zh) | 一种眼前节形态分析方法 | |
| CN104958061B (zh) | 双目立体视觉三维成像的眼底oct成像方法及其系统 | |
| CN110934563B (zh) | 眼科信息处理装置、眼科装置及眼科信息处理方法 | |
| US9364144B2 (en) | Method and device for recording and displaying an OCT whole-eye scan | |
| EP1976424B1 (en) | A method of eye examination by optical coherence tomography | |
| US8770753B2 (en) | Scanning and processing using optical coherence tomography | |
| CN108618749B (zh) | 基于便携式数字化眼底照相机的视网膜血管三维重建方法 | |
| US10893799B2 (en) | Method for determining the topography of the cornea of an eye | |
| US11839510B2 (en) | Composite ultrasound images | |
| US9918680B2 (en) | Automated analysis of the optic nerve head via optical coherence tomography measurements, methods and representations | |
| CN103142210B (zh) | 一种基于oct技术的周边屈光测量方法 | |
| CN113100829B (zh) | 一种眼前节三维超声扫描成像装置以及方法 | |
| CN105608675A (zh) | 一种眼底组织oct图像运动伪差的矫正方法 | |
| JP2020137672A (ja) | 眼科情報処理装置、眼科装置、眼科情報処理方法、及びプログラム | |
| WO2022150294A1 (en) | Apparatus for calibrating retinal imaging systems and related methods | |
| CN111493832A (zh) | 一种基于En face-OCT的内窥成像方法 | |
| JP2023138773A (ja) | 眼科情報処理装置、眼科装置、眼科情報処理方法、及びプログラム | |
| CN103519781A (zh) | 眼前节光学相干层析图像角膜中央反光的消除方法及装置 | |
| CN118415577A (zh) | 基于眼前眼后光路切换的oct成像系统及眼轴测量方法 | |
| CN117064326A (zh) | 一种基于oct主体的双探头双目同步眼前节成像系统 | |
| JP2023123245A (ja) | 走査型イメージング装置、情報処理装置、走査型イメージング装置の制御方法、情報処理装置の制御方法、情報処理方法、プログラム、及び記録媒体 | |
| Khachikian et al. | Evolution of Topographic/Tomographic Imaging | |
| Restori | Ultrasound in opthalmic diagnosis |
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 |