CN116919336A - 一种视网膜成像方法 - Google Patents
一种视网膜成像方法 Download PDFInfo
- Publication number
- CN116919336A CN116919336A CN202210366681.8A CN202210366681A CN116919336A CN 116919336 A CN116919336 A CN 116919336A CN 202210366681 A CN202210366681 A CN 202210366681A CN 116919336 A CN116919336 A CN 116919336A
- Authority
- CN
- China
- Prior art keywords
- image
- compensation
- pupil center
- wavefront
- compensation value
- 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.)
- Pending
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 67
- 210000001525 retina Anatomy 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 claims abstract description 37
- 210000001747 pupil Anatomy 0.000 claims abstract description 28
- 239000011159 matrix material Substances 0.000 claims abstract description 23
- 230000003287 optical effect Effects 0.000 claims description 38
- 230000002207 retinal effect Effects 0.000 claims description 32
- 230000004075 alteration Effects 0.000 claims description 28
- 230000003044 adaptive effect Effects 0.000 claims description 16
- 210000005252 bulbus oculi Anatomy 0.000 claims description 15
- 201000009310 astigmatism Diseases 0.000 claims description 14
- 210000001508 eye Anatomy 0.000 claims description 14
- 230000033001 locomotion Effects 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 11
- 230000004424 eye movement Effects 0.000 claims description 7
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 13
- 238000012014 optical coherence tomography Methods 0.000 description 14
- 238000001514 detection method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 238000002073 fluorescence micrograph Methods 0.000 description 5
- 230000001276 controlling effect Effects 0.000 description 4
- 201000010099 disease Diseases 0.000 description 4
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 4
- 238000005457 optimization Methods 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000004397 blinking Effects 0.000 description 2
- 230000009191 jumping Effects 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 238000009990 desizing Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 210000003128 head Anatomy 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- 230000004256 retinal image Effects 0.000 description 1
- 238000004621 scanning probe microscopy Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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/103—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for determining refraction, e.g. refractometers, skiascopes
-
- 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/102—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for optical coherence tomography [OCT]
-
- 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/11—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for measuring interpupillary distance or diameter of pupils
- A61B3/112—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for measuring interpupillary distance or diameter of pupils for measuring diameter of pupils
-
- 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/113—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for determining or recording eye movement
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Ophthalmology & Optometry (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Human Computer Interaction (AREA)
- Eye Examination Apparatus (AREA)
Abstract
本发明提供了一种视网膜成像方法,先对瞳孔位置进行调整,以获得更佳的成像位置,从而提高成像的质量;首先获取瞳孔中心的位置,调节瞳孔中心初次定位于预设的目标瞳孔中心,然后根据波前图像的点阵分布情况进一步调节瞳孔中心的位置,使瞳孔中心二次定位于预设的目标瞳孔中心,使用两种方法先后进行定位,使得瞳孔定位的结果更加精确,从而达到更优的视网膜成像效果。
Description
技术领域
本发明涉及视网膜成像技术领域,更具体地说,涉及一种视网膜成像方法。
背景技术
视网膜是人眼重要的组成部分,目前全球大约超过十亿人患有视网膜相关疾病,为了更有效的实现对视网膜相关疾病的治疗,对其治疗装置的优化是必不可免的,也就是说视网膜的高分辨成像装置对于视网膜相关疾病的诊断和疗效评估有着重大意义。
AO-OCT(Adaptive Optics-OpticalCoherence Tomography,自适应光学-光学相干断层成像技术)技术应用于眼科成像,即应用前沿自适应光学技术将低相干干涉仪和共焦扫描显微术加上现代计算机图像处理技术显示细胞级别视网膜成像。
该系统可以消除人眼像差,获得视网膜具有接近理论极限的高分辨率图像;通过提取及对比视网膜细胞形态、密度分布的图像特征,能够判断早期病变,系统无标记,无需打入特定造影剂,在观查、检测过程,对人体无任何伤害。
但是,目前AO-OCT系统的高复杂度会导致操作方法比较困难,使采集的图像质量较差以及采集效率较低。并且成像时人眼瞳孔中心的位置与成像装置的光轴难以精确对准,位置偏差会影响成像效果。
发明内容
有鉴于此,为解决上述问题,本发明提供一种视网膜成像方法,技术方案如下:
一种视网膜成像方法,先对瞳孔位置进行调整,包括:
S1,获取瞳孔中心的位置,调节所述瞳孔中心初次定位于预设的目标瞳孔中心;
S2,获取波前图像,所述波前图像为入射光线经视网膜反射回来后在波前探测器中汇聚而成的点阵图像;
S3,基于所述波前图像中的点阵分布情况,进一步调节所述瞳孔中心的位置,以使瞳孔中心二次定位于预设的目标瞳孔中心。
优选的,在上述一种视网膜成像方法中,所述步骤S1包括:
采用至少两个光探测器采集眼部区域图像,每个光探测器的采集角度各不相同;
对所述眼部区域图像进行分析处理,识别所述瞳孔中心;
比较所述瞳孔中心与预设的目标瞳孔中心之间的偏移量;
基于所述偏移量,调节所述瞳孔中心向预设的目标瞳孔中心移动。
优选的,在上述一种视网膜成像方法中,所述步骤S3的具体方法为:当点阵分布情况不满足预设的条件时,调节所述瞳孔中心的位置,然后重复步骤S2和S3。
优选的,在上述一种视网膜成像方法中,所述预设的条件为:与目标瞳孔点阵分布的区域面积或轮廓的相似程度;调节所述瞳孔中心的位置时,根据目标瞳孔的点阵分布与实时的点阵分布之间的差别来确定瞳孔位置的调整方向。
优选的,在上述一种视网膜成像方法中,在所述步骤S3之后自动进行屈光度补偿,方法为:
选取至少两个屈光度,对光路进行补偿,并获得对应的波前图像;
评价各波前图像的质量;
选取质量满足预设要求的波前图像,以其对应的屈光度对光路进行自动补偿。
优选的,在上述一种视网膜成像方法中,所述自动进行屈光度补偿包括自动进行球镜和/或散光补偿,所述自动进行球镜或散光补偿具体为:
选取初始补偿值和步长;
以所述初始补偿值加上和减去所述步长分别得到第一补偿值和第二补偿值;
对光路分别补偿所述第一补偿值、初始补偿值和第二补偿值,并获得对应的波前图像;
评价各波前图像的质量;
选取质量最优的波前图像,以其对应的补偿值作为新的初始补偿值,重新确定步长,并返回执行以所述初始补偿值加上和减去步长分别得到第一补偿值和第二补偿值这一步骤,直至步长或波前图像质量满足预设要求。
优选的,在上述一种视网膜成像方法中,所述自动进行屈光度补偿包括自动进行球镜和/或散光补偿,所述自动进行球镜或散光补偿具体为:
在补偿值域内按固定间隔选取多个样本补偿值;
对光路分别补偿各样本补偿值,并获得对应的波前图像;
评价各波前图像的质量;
选取质量最优的波前图像,以其对应的补偿值加上和减去固定间隔作为新的补偿值域端点,并重新确定间隔,返回执行在补偿值域内按固定间隔选取多个样本补偿值这一步骤,直至间隔或波前图像质量满足预设要求。
优选的,在上述一种视网膜成像方法中,在所述自动进行屈光度补偿后,进行自适应光学像差补偿,包括:
获取波前图像;
计算所述波前图像中的点阵相对于参考平面波的点阵的偏移量,根据点阵的偏移量确定像差补偿反馈值;
以所述像差补偿反馈值对光路进行补偿,实现自适应光学像差补偿。
优选的,在上述一种视网膜成像方法中,所述视网膜成像为视网膜多模态成像,包括视网膜的第一模态图像和第二模态图像同步成像,其方法为:
在进行眼球运动跟踪时,获取第一模态子图像,将其与预设的第一模态参考图像进行配准,求得眼球的偏移量;所述第一模态子图像是一帧第一模态图像中的一段;
基于所述偏移量生成第一控制信号,根据第一控制信号调节光路状态进行眼球运动跟踪;并且,基于所述偏移量确定第二模态子图像在第二模态图像中的位置;所述第二模态子图像是获取第一模态子图像时同步获取到的部分第二模态图像;
按所述位置累积第二模态子图像,组合成第二模态图像输出。
优选的,在上述一种视网膜成像方法中,所述视网膜成像方法还包括:
在屈光度补偿完成之后,返回执行步骤S1重复进行瞳孔定位和屈光度补偿。
相较于现有技术,本发明实现的有益效果为:
本发明提供的一种视网膜成像方法先对瞳孔位置进行调整,以获得更佳的成像位置,从而提高成像的质量。具体是先获取瞳孔中心的位置,调节瞳孔中心初次定位于预设的目标瞳孔中心,然后根据波前图像的点阵分布情况进一步调节瞳孔中心的位置,使瞳孔中心二次定位于预设的目标瞳孔中心,使用两种方法先后进行定位,使得瞳孔定位的结果更加精确,从而达到更优的视网膜成像效果。
并且,本发明还在瞳孔中心定位后进行像差补偿,补偿像差后能够在视网膜成像时,消除人眼像差的影响,进一步提升成像质量;并且进行像差补偿时,先自动进行屈光度的补偿,再进行自适应光学像差补偿,简化了操作流程、降低了操作难度,避免了手动输入屈光度参数的繁琐过程及外部获取的屈光度数据不准确问题,更适用于大规模推广应用;且屈光度补偿为自适应光学像差补偿提供了一个较好的初始条件以获得更稳定的收敛效果。
此外,本发明在多模态中可进行实时的眼动跟踪,能消除眼球漂移、跳动和眨眼等情况对成像的影响,从而获得所需模态目标视野的稳定图像,且将一帧图像分为多段处理,能够加快处理速度,提高扫描式的光学成像系统的眼动跟踪的实时性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为本发明实施例提供的一种视网膜成像方法的流程示意图;
图2为本发明实施例提供的一种瞳孔初次定位的示意图;
图3为本发明实施例提供的一种瞳孔二次定位的示意图;
图4为本发明实施例提供的一种数据处理的示意图;
图5为本发明实施例提供的一种基于跟踪的A-line数据拼接成B-Scan的示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
参考图1,图1为本发明实施例提供的一种视网膜成像方法的流程示意图。
具体地,先对瞳孔位置进行调整,包括:
S1,获取瞳孔中心的位置,调节所述瞳孔中心初次定位于预设的目标瞳孔中心。
S2,获取波前图像,所述波前图像为入射光线经视网膜反射回来后在波前探测器中汇聚而成的点阵图像。
S3,基于所述波前图像中的点阵分布情况,进一步调节所述瞳孔中心的位置,以使瞳孔中心二次定位于预设的目标瞳孔中心。
在该实施例中,通过至少两次瞳孔定位,使瞳孔中心的位置精确定位于预设的目标瞳孔中心,该目标瞳孔中心包括但不限定于通过标定的方法确定,在瞳孔中心精确定位后有利于提高后续视网膜成像的质量,获得所需的高质量的视网膜成像图像。
可选的,在本发明另一实施例中,步骤S1获取瞳孔中心的位置,调节所述瞳孔中心初次定位于预设的目标瞳孔中心,具体为:
采用至少两个光探测器采集眼部区域图像,每个光探测器的采集角度各不相同。
对所述眼部区域图像进行分析处理,识别所述瞳孔中心。
比较所述瞳孔中心与预设的目标瞳孔中心之间的偏移量。
基于所述偏移量,调节所述瞳孔中心向预设的目标瞳孔中心移动。
参考图2,图2为本发明实施例提供的一种瞳孔初次定位的效果示意图,计算机获取至少两个由光探测器拍摄的不同拍摄角度的眼部区域图像,从各个眼部区域图像中识别出瞳孔中心的位置,之后与预先设定的目标瞳孔中心所在位置进行比较,分析得到二者瞳孔中心之间的偏移量。
之后,计算瞳孔中心相对目标瞳孔中心的三维距离,发送控制信号给三维运动结构做上、下、前、后、左、右的运动,调整瞳孔中心的位置,重复以上操作直至偏移量小于设定的阈值条件即可,此时认为瞳孔初次定位成功。
需要说明的是,该三维运动结构可通过控制信号自动进行调节,也可根据计算得到的三维距离进行手动调节。
可选的,在本发明另一实施例中,为了最大程度的提高最终视网膜成像的成像效果,在瞳孔初次定位成功后,并不能保证成像位置处于最佳位置,还需要进行进一步的精准定位,即获取波前图像,基于所述波前图像中的点阵分布情况,进一步调节所述瞳孔中心的位置,以使瞳孔中心二次定位于预设的目标瞳孔中心。
参考图3,根据自适应光学原理可知,在光学系统中通过分光镜等装置将由被视网膜反射回来的光线部分收集到波前探测器,以使得波前探测器获得的波前图像可以表征出反射光线的波前、强度分布等信息。
基于此,获取实时的波前图像,并检测波前图像中的光线强度分布;波前探测器获得的波前图像是入射光线经视网膜反射回来后在波前探测器中汇聚而成的点阵图像,这些点阵的分布区域越大说明进入眼球的光线越多,越有利于视网膜成像质量的提高。
因此,利用波前探测器进行定位的原则就是使得波前图像中的点阵越多越好;基于此,可以采用简易的试探法以较小的步长移动三维运动结构来迭代优化,也可以根据目标的点阵分布与实时的点阵分布之间的差别,以及光学系统的空间几何关系来确定三维运动结构的移动方向,保证波前图像中的点阵分布情况满足预设的条件,所述预设的条件为:与目标瞳孔点阵分布的区域面积或轮廓的相似程度。
如图3所示,左边的实时点阵图相对于右边的目标点阵图缺失了左下角的部分点阵,若三维运动结构与图像表示的方向是一致的,则计算机发送指令会使三维运动结构向左下角移动合适的距离。多步迭代优化之后需保证点阵图像的点阵分布情况满足预设的分布条件,例如实时点阵分布的区域达到目标点阵分布的区域面积的90%或接近圆形分布等。
若点阵分布情况不满足预设的条件时,调节所述瞳孔中心的位置,然后重复步骤S2和步骤S3,直至点阵分布情况满足预设的条件。
可选的,在本发明另一实施例中,在所述步骤S3之后所述视网膜成像方法还包括:自动进行屈光度补偿,具体方法为:
选取至少两个屈光度,对光路进行补偿,并获得对应的波前图像。
评价各波前图像的质量。
选取质量满足预设要求的波前图像,以其对应的屈光度对光路进行自动补偿。
在该实施例中,为了实现最优的视网膜成像效果,在进行视网膜成像之前,在对瞳孔中心定位完成之后需要进行实时像差补偿,该实时像差补偿具体包括自动屈光度补偿和自适应光学像差补偿;其中屈光度属于像差中的低阶像差,先进行屈光度补偿的目的是给后续的自适应光学像差补偿提供了一个较好的初始条件以获得更稳定的收敛效果。
自动屈光度补偿具体是通过控制补偿镜的补偿值来实现对眼球屈光度的补偿;并且自动的方法可以避免手动输入的过程繁琐以及外部获取的屈光度数据不准确的问题。
可选的,在本发明的另一实施例中,所述自动进行屈光度补偿包括自动进行球镜和/或散光补偿,所述自动进行球镜或散光补偿具体为:
选取初始补偿值和步长。
以所述初始补偿值加上和减去所述步长分别得到第一补偿值和第二补偿值。
对光路分别补偿所述第一补偿值、初始补偿值和第二补偿值,并获得对应的波前图像。
评价各波前图像的质量。
选取质量最优的波前图像,以其对应的补偿值作为新的初始补偿值,重新确定步长,并返回执行以所述初始补偿值加上和减去步长分别得到第一补偿值和第二补偿值这一步骤,直至步长或波前图像质量满足预设要求。
在该实施例中,假设选择初始补偿值为0度、初始步长为100度,以该初始步长选择初始补偿值的正负方向(即加上初始步长和减去初始步长)各一个值,即0度(相当于初始补偿值)、+100度(相当于第一补偿值)、-100度(相当于第二补偿值),0度即为中间值。
通过控制补偿镜给光路分别补偿所述第一补偿值、初始补偿值和第二补偿值,进行三次补偿,每次补偿同时获取波前探测器的波前探测信号,依据所述波前探测信号生成各自补偿时的波前图像,选定一种评价波前质量的度量方法,如波前图像中所有像素点的平方和,平方和越大表示补偿后的光线越接近平面波,然后在这三个值中确定下一步搜索的方向:
若是中间值的波前质量最优,则缩小步长后再重复以上步骤进行搜索,直到收敛至预设的阈值以下。例如0度的波前质量最优,则将初始步长缩小为50度,即选择0度、+50度、-50度,通过控制补偿镜给光路分别补偿这三个补偿值,重复以上步骤,直到收敛至预设的阈值以下,收敛至预设的阈值以下具体是指步长满足预设的精度要求(如步长小于25°)或者波前图像质量满足预设要求。
若中间值的波前质量不是最优,则将最优值作为新的中间值,步长保持不变。例如+100度的波前质量最优,则选择+100度、0度、+200度,通过控制补偿镜给光路分别补偿这三个补偿值,重复以上步骤,直到收敛至预设的阈值以下。
需要说明的是,假设选择初始补偿值为0度、初始步长为100度,以该初始步长还可以选择初始补偿值的正负方向(即加上初始步长和减去初始步长)各多个值(以两个值为例进行说明),即0度、+100度、+200度、-100度、-200度,0度为中间值。
需要说明的是,自动进行球镜的补偿和自动进行散光的补偿的补偿思路相同,可以只补偿其中一个,若要对二者都进行补偿,可以在其中一个补偿完成之后再进行另一个的补偿,例如先对球镜补偿后再进行散光的补偿。
可选的,在本发明另一实施例中,所述自动进行屈光度补偿包括自动进行球镜和/或散光补偿,所述自动进行球镜或散光补偿具体为:
在补偿值域内按固定间隔选取多个样本补偿值。
对光路分别补偿各样本补偿值,并获得对应的波前图像。
评价各波前图像的质量。
选取质量最优的波前图像,以其对应的补偿值加上和减去固定间隔作为新的补偿值域端点,并重新确定间隔,返回执行在补偿值域内按固定间隔选取多个样本补偿值这一步骤,直至间隔或波前图像质量满足预设要求。
在该实施例中,假设球镜的补偿值域为-400度到+400度,固定间隔为200度,则基于200度的固定间隔对-400度到+400度的补偿值域进行分割处理,获得-400度、-200度、0度、+200度、+400度的多个样本补偿值。
通过控制补偿镜给光路分别补偿-400度、-200度、0度、+200度、+400度的样本补偿值,进行多次补偿,每次补偿同时获取波前探测器的波前探测信号,依据所述波前探测信号生成各自补偿时的波前图像,选定一种评价波前质量的度量方法,如波前图像中所有像素点的平方和,平方和越大表示补偿后的光线越接近平面波,然后在这些样本补偿值中确定下一步搜索的方向:
若0度的波前质量最优,则确定-200度和+200度为新的补偿值域端点,则新的补偿值域为-200度到+200度,假设重新确定的间隔为100度,则基于100度的固定间隔对-200度到+200度的补偿阈值进行分割处理,获得-200度、-100度、0度、+100度、+200度的多个样本补偿值。
通过控制补偿镜给光路分别补偿-400度、-200度、0度、+200度、+400度的样本补偿值,重复以上步骤,直到间隔满足预设的精度要求(如间隔小于10度)或波前图像质量满足预设要求。
需要说明的是,自动进行球镜的补偿和自动进行散光的补偿的补偿思路相同,可以只补偿其中一个,若要对二者都进行补偿,可以在其中一个补偿完成之后再进行另一个的补偿,例如先对球镜补偿后再进行散光的补偿。
需要说明的是,在屈光度补偿完成之后,返回执行步骤S1重复进行瞳孔定位和屈光度补偿,以优化瞳孔定位和屈光度补偿的效果,之后再进行自适应光学像差补偿。
可选的,在本发明另一实施例中,在屈光度补偿之后方可直接执行自适应光学像差补偿,并且屈光度补偿为自适应光学像差补偿提供了一个较好的初始条件以获得更稳定的收敛效果。
进行自适应光学像差补偿,包括:
获取波前图像。
计算所述波前图像中的点阵相对于参考平面波的点阵的偏移量,根据点阵的偏移量确定像差补偿反馈值。
以所述像差补偿反馈值对光路进行补偿,实现自适应光学像差补偿。
在该实施例中,根据自适应光学原理,自适应光学像差补偿的方式是通过实时获取波前探测器的波前探测信号,依据所述波前探测信号生成波前图像,检测图像中的点阵位置并计算相对于参考平面波的点阵的偏移量,进而与校正获得的控制矩阵相乘,得到应当发送至补偿镜的反馈。
需要说明的是,该过程是一个闭环控制过程,可以通过调节增益的大小来获得较为稳定的收敛过程。
可选的,在本发明另一实施例中,在自适应光学系统稳定成像之后,此时在头托、固视点等装置的协助下,眼球的运动被限制在一个较小的范围内,主要分为漂移、跳动和眨眼等情况。
在一些应用场景下,为了获得目标视野的稳定图像,还需要进行实时的眼球跟踪,利用获得的经视网膜反射回来的实时视网膜图像进行眼球运动跟踪,具体的:
所述视网膜成像为视网膜多模态成像,包括视网膜的第一模态图像和第二模态图像同步成像,其方法为:
在进行眼球运动跟踪时,获取第一模态子图像,将其与预设的第一模态参考图像进行配准,求得眼球的偏移量;所述第一模态子图像是一帧第一模态图像中的一段。
基于所述偏移量生成第一控制信号,以调节光路状态进行眼球运动跟踪;并且,基于所述偏移量确定第二模态子图像在第二模态图像中的位置;所述第二模态子图像是获取第一模态子图像时同步获取到的部分第二模态图像。
按位置关系累积第二模态子图像,组合成第二模态图像输出。
在该实施例中,所述第一模态图像为二维反射图像,第二模态图像为OCT图像或荧光图像,参考图4,图4为本发明实施例提供的一种数据处理的示意图,为了提高眼动跟踪的实时性,可以将一帧图像数据划分为多段依次发送给计算单元处理,假设系统是按行扫描,每帧包含512行数据,可每采集到32或16行的数据后立即发送给计算单元,这样相对于等到一帧数据完全采集完再发送的方式提高到16或32倍的跟踪反馈频率。
具体的跟踪过程如下:
数据采集模块每获得预定好的一小段图像数据后立即发送至计算机单元,计算单元获得数据后经过插值、去正弦化等预处理,开始跟踪前会从视频流中选择合适的一帧或多帧图像处理后作为参考帧,开始跟踪后实时将预处理后的每一小段图像数据与选定的参考帧图像数据进行配准求得眼球的偏移量,然后根据偏移量生成第一控制信号,根据第一控制信号调节光路的状态,实现眼球运动跟踪,以达到稳定成像区域的目的。调节光路的状态可通过调整光路中振镜的偏转角度实现。
针对点扫描的视网膜成像结构而言,所述第一模态图像为二维反射图像,第二模态图像为OCT图像,光探测器和OCT探测器的采集频率不一致的情况下,具体处理方法为:依据跟踪选定的参考帧图像确定待成像的视野区域,光探测器和OCT探测器同步采集图像数据,即光探测器采集的每一行反射图像数据对应同时间段内的一条OCT A-line。在获得每16或32行反射图像数据及相应的A-line图像数据后,根据反射图像数据计算眼球的偏移量并最终反馈到相应振镜上。同时将采集的相应A-line图像数据根据跟踪的偏移量将各A-line数据放置到基于跟踪参考帧所确定的三维参考系的对应位置中。每完成一帧反射图像的跟踪,对应完成一个B-scan的采集,然后再进入下一个位置的扫描。若跟踪失败则丢弃相应数据并保持位置不变,重复扫描直到跟踪成功,最终按位置关系拼接成整个视野的三维OCT数据,参考图5,图5为本发明实施例提供的一种基于跟踪的A-line数据拼接成B-Scan的示意图。
针对光探测器和OCT探测器的采集频率相同的情况,处理方法不同之处在于每完成一行反射图像的跟踪,对应完成OCT一个B-scan的采集,在获得每16或32行反射数据及相应B-scan数据后即可进行数据传输和处理。
针对线扫描的视网膜成像结构,所述第一模态图像为二维反射图像,第二模态图像为OCT图像,光探测器和OCT探测器的采集频率相同,与点扫描处理方法不同之处在于每完成一条反射线的跟踪,对应完成OCT一个B-scan的采集。
针对荧光扫描的视网膜成像结构,所述第一模态图像为二维反射图像,第二模态图像为荧光图像,每完成一行反射图像的跟踪,对应完成一行荧光图像的采集,每一行反射图像数据对应同时间段内的一行荧光图像数据,在获得每段(如16或32行)反射图像数据及相应的荧光图像数据后,根据反射图像数据计算眼球的偏移量并最终反馈到相应振镜上。同时根据跟踪的偏移量将各段荧光图像数据放置到基于跟踪参考帧所确定的参考系的对应位置中,多段荧光图像数据按位置关系拼接成一张荧光图像数据,之后还可以再按位置关系累积多张荧光图像数据,最终叠加成一个总的荧光图像输出,达到图像增强效果。
采用实时跟踪技术的成像系统可以解决眼球运动导致的图像扭曲、错位等问题,对于较慢的OCT成像眼球运动会导致各个B-scan的位置不准确等问题,往往需要复杂而缓慢的后处理来减小这些影响。因此如果采用上述方法的实时跟踪可以在采集过程中实时展示有价值的图像数据,而不用等到后处理完才能看到想要的结果。
以上对本发明所提供的一种视网膜成像方法进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种视网膜成像方法,其特征在于,先对瞳孔位置进行调整,包括:
S1,获取瞳孔中心的位置,调节所述瞳孔中心初次定位于预设的目标瞳孔中心;
S2,获取波前图像,所述波前图像为入射光线经视网膜反射回来后在波前探测器中汇聚而成的点阵图像;
S3,基于所述波前图像中的点阵分布情况,进一步调节所述瞳孔中心的位置,以使瞳孔中心二次定位于预设的目标瞳孔中心。
2.根据权利要求1所述的一种视网膜成像方法,其特征在于,所述步骤S1包括:
采用至少两个光探测器采集眼部区域图像,每个光探测器的采集角度各不相同;
对所述眼部区域图像进行分析处理,识别所述瞳孔中心;
比较所述瞳孔中心与预设的目标瞳孔中心之间的偏移量;
基于所述偏移量,调节所述瞳孔中心向预设的目标瞳孔中心移动。
3.根据权利要求1所述的一种视网膜成像方法,其特征在于,所述步骤S3的具体方法为:当点阵分布情况不满足预设的条件时,调节所述瞳孔中心的位置,然后重复步骤S2和S3。
4.根据权利要求3所述的一种视网膜成像方法,其特征在于,所述预设的条件为:与目标瞳孔点阵分布的区域面积或轮廓的相似程度;调节所述瞳孔中心的位置时,根据目标瞳孔的点阵分布与实时的点阵分布之间的差别来确定瞳孔位置的调整方向。
5.根据权利要求1所述的一种视网膜成像方法,其特征在于,在所述步骤S3之后自动进行屈光度补偿,方法为:
选取至少两个屈光度,对光路进行补偿,并获得对应的波前图像;
评价各波前图像的质量;
选取质量满足预设要求的波前图像,以其对应的屈光度对光路进行自动补偿。
6.根据权利要求5所述的一种视网膜成像方法,其特征在于,所述自动进行屈光度补偿包括自动进行球镜和/或散光补偿,所述自动进行球镜或散光补偿具体为:
选取初始补偿值和步长;
以所述初始补偿值加上和减去所述步长分别得到第一补偿值和第二补偿值;
对光路分别补偿所述第一补偿值、初始补偿值和第二补偿值,并获得对应的波前图像;
评价各波前图像的质量;
选取质量最优的波前图像,以其对应的补偿值作为新的初始补偿值,重新确定步长,并返回执行以所述初始补偿值加上和减去步长分别得到第一补偿值和第二补偿值这一步骤,直至步长或波前图像质量满足预设要求。
7.根据权利要求5所述的一种视网膜成像方法,其特征在于,所述自动进行屈光度补偿包括自动进行球镜和/或散光补偿,所述自动进行球镜或散光补偿具体为:
在补偿值域内按固定间隔选取多个样本补偿值;
对光路分别补偿各样本补偿值,并获得对应的波前图像;
评价各波前图像的质量;
选取质量最优的波前图像,以其对应的补偿值加上和减去固定间隔作为新的补偿值域端点,并重新确定间隔,返回执行在补偿值域内按固定间隔选取多个样本补偿值这一步骤,直至间隔或波前图像质量满足预设要求。
8.根据权利要求5所述的一种视网膜成像方法,其特征在于,在所述自动进行屈光度补偿后,进行自适应光学像差补偿,包括:
获取波前图像;
计算所述波前图像中的点阵相对于参考平面波的点阵的偏移量,根据点阵的偏移量确定像差补偿反馈值;
以所述像差补偿反馈值对光路进行补偿,实现自适应光学像差补偿。
9.根据权利要求1所述的一种视网膜成像方法,其特征在于,所述视网膜成像为视网膜多模态成像,包括视网膜的第一模态图像和第二模态图像同步成像,其方法为:
在进行眼球运动跟踪时,获取第一模态子图像,将其与预设的第一模态参考图像进行配准,求得眼球的偏移量;所述第一模态子图像是一帧第一模态图像中的一段;
基于所述偏移量生成第一控制信号,根据第一控制信号调节光路状态进行眼球运动跟踪;并且,基于所述偏移量确定第二模态子图像在第二模态图像中的位置;所述第二模态子图像是获取第一模态子图像时同步获取到的部分第二模态图像;
按所述位置累积第二模态子图像,组合成第二模态图像输出。
10.根据权利要求5-7任一项所述的一种视网膜成像方法,其特征在于,所述视网膜成像方法还包括:
在屈光度补偿完成之后,返回执行步骤S1重复进行瞳孔定位和屈光度补偿。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210366681.8A CN116919336A (zh) | 2022-04-08 | 2022-04-08 | 一种视网膜成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210366681.8A CN116919336A (zh) | 2022-04-08 | 2022-04-08 | 一种视网膜成像方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116919336A true CN116919336A (zh) | 2023-10-24 |
Family
ID=88374382
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210366681.8A Pending CN116919336A (zh) | 2022-04-08 | 2022-04-08 | 一种视网膜成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116919336A (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050140984A1 (en) * | 2003-12-31 | 2005-06-30 | Hitzenberger Christoph K. | Efficient optical coherence tomography (OCT) system and method for rapid imaging in three dimensions |
CN101862178A (zh) * | 2010-06-02 | 2010-10-20 | 中国科学院光电技术研究所 | 一种基于自适应光学的反射式共焦扫描视网膜成像系统 |
US20120274783A1 (en) * | 2011-04-29 | 2012-11-01 | Optovue, Inc. | Imaging with real-time tracking using optical coherence tomography |
JP2015031812A (ja) * | 2013-08-02 | 2015-02-16 | 株式会社ニコン | 補償光学素子の設定方法及び顕微鏡 |
US20160317031A1 (en) * | 2015-04-29 | 2016-11-03 | University Of Rochester | Method of imaging multiple retinal structures |
US20210121063A1 (en) * | 2019-09-09 | 2021-04-29 | Suzhou Institute Of Biomedical Engineering And Technology, Chinese Academy Of Sciences | Large Field -of -View Adaptive Optics Retinal Imaging System and Method with Common Optical Path Beam Scanning |
CN112754420A (zh) * | 2021-01-08 | 2021-05-07 | 瑞尔明康(杭州)视光科技有限公司 | 波前像差测量系统、相关联的视知觉学习训练系统和方法 |
-
2022
- 2022-04-08 CN CN202210366681.8A patent/CN116919336A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050140984A1 (en) * | 2003-12-31 | 2005-06-30 | Hitzenberger Christoph K. | Efficient optical coherence tomography (OCT) system and method for rapid imaging in three dimensions |
CN101862178A (zh) * | 2010-06-02 | 2010-10-20 | 中国科学院光电技术研究所 | 一种基于自适应光学的反射式共焦扫描视网膜成像系统 |
US20120274783A1 (en) * | 2011-04-29 | 2012-11-01 | Optovue, Inc. | Imaging with real-time tracking using optical coherence tomography |
JP2015031812A (ja) * | 2013-08-02 | 2015-02-16 | 株式会社ニコン | 補償光学素子の設定方法及び顕微鏡 |
US20160317031A1 (en) * | 2015-04-29 | 2016-11-03 | University Of Rochester | Method of imaging multiple retinal structures |
US20210121063A1 (en) * | 2019-09-09 | 2021-04-29 | Suzhou Institute Of Biomedical Engineering And Technology, Chinese Academy Of Sciences | Large Field -of -View Adaptive Optics Retinal Imaging System and Method with Common Optical Path Beam Scanning |
CN112754420A (zh) * | 2021-01-08 | 2021-05-07 | 瑞尔明康(杭州)视光科技有限公司 | 波前像差测量系统、相关联的视知觉学习训练系统和方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9320424B2 (en) | Image display apparatus, image display method and imaging system | |
JP5297415B2 (ja) | 眼科装置及び眼科方法 | |
JP6167518B2 (ja) | 画像処理装置および方法、並びにプログラム | |
US9089280B2 (en) | Image processing apparatus, image processing method, and program storage medium | |
US20210373308A1 (en) | Microscope system, smart medical device, automatic focusing method and storage medium | |
US11826102B2 (en) | Ophthalmic device, control method therefor, and recording medium | |
CN113040701B (zh) | 一种三维眼动追踪系统及其追踪方法 | |
CN113827181A (zh) | 角膜图像生成方法及装置 | |
CN116919336A (zh) | 一种视网膜成像方法 | |
CN112155511B (zh) | 一种基于深度学习补偿oct采集过程中人眼抖动的方法 | |
US10674902B2 (en) | Information processing apparatus, operation method thereof, and computer program | |
JP2016209133A (ja) | 眼科撮像装置、その作動方法、及びコンピュータプログラム | |
US10123688B2 (en) | Information processing apparatus, operation method thereof, and computer program | |
US11857255B2 (en) | Ophthalmic apparatus | |
US11903646B2 (en) | Ophthalmic apparatus, method of controlling the same, method of ophthalmic optical coherence tomography, and recording medium | |
JP6587430B2 (ja) | 画像生成装置及び画像生成方法 | |
JP6602043B2 (ja) | 情報処理装置、その作動方法、及びコンピュータプログラム | |
JP2016086842A (ja) | 画像処理装置及び画像処理方法 | |
CN113040700B (zh) | 一种眼动追踪系统及其追踪方法 | |
JP5726229B2 (ja) | 眼科装置及び眼科方法 | |
CN116327112B (zh) | 一种动态机器视觉指导的全自动眼科oct系统 | |
CN117414103B (zh) | 一种用于眼部的激光成像方法与系统 | |
JP2017086807A (ja) | 画像取得装置及びその駆動方法、並びに、プログラム | |
WO2023008385A1 (ja) | 前眼部解析装置、前眼部解析方法、およびプログラム | |
US20230380680A1 (en) | Ophthalmic apparatus and method of controlling the same |
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 |