CN112168132B - 一种使用oct信号进行眼底屈光补偿判定与成像优化的方法 - Google Patents
一种使用oct信号进行眼底屈光补偿判定与成像优化的方法 Download PDFInfo
- Publication number
- CN112168132B CN112168132B CN202011236869.8A CN202011236869A CN112168132B CN 112168132 B CN112168132 B CN 112168132B CN 202011236869 A CN202011236869 A CN 202011236869A CN 112168132 B CN112168132 B CN 112168132B
- Authority
- CN
- China
- Prior art keywords
- image
- optical path
- imaging
- retina
- oct
- 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/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/12—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
- G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30041—Eye; Retina; Ophthalmic
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Theoretical Computer Science (AREA)
- Ophthalmology & Optometry (AREA)
- Molecular Biology (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Data Mining & Analysis (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Discrete Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Eye Examination Apparatus (AREA)
Abstract
本申请实施例公开一种使用OCT信号进行眼底屈光补偿判定与成像优化的方法,该方法包括:S1、控制参考臂光程:通过计算视网膜结构的凹凸性以及成像位置,在图像中找到视网膜结构并调整到图像中合适位置;S2、自动屈光补偿:控制参考臂光程与样品臂扫描光束光程的光程差为固定值,通过爬山法找到信号最强的图像进行屈光补偿;S3、视网膜结构左右微调:判定视网膜结构的左右偏移,根据偏移控制Y轴电机进行水平调整;S4、调整视网膜成像位置:根据图像强度再次调节参考臂光程以调整视网膜在图像中的高低位置。本申请的使用OCT信号进行眼底屈光补偿判定与成像优化的方法,能够降低成本并实现快速精准的屈光补偿,同时对视网膜成像进行优化。
Description
技术领域
本申请涉及视网膜图像成像技术领域,具体是一种使用OCT信号进行眼底屈光补偿判定与成像优化的方法。
背景技术
眼底视网膜是人体一种微观结构组织,很多眼部疾病以及部分全身性疾病都可以在眼底视网膜上有所体现。光学相干断层扫描成像设备(OCT)是一种眼底视网膜成像系统,可以获取三维眼底图像。OCT设备使用低相干光干涉成像技术,需要扫描探测光束与参考光束的光程差在一定范围内才可以清晰成像。为了使用OCT设备获取清晰的眼底图像,通常的OCT设备需要进行扫描光路的屈光补偿和成像高度调整两步。屈光补偿是指对于不同屈光度的人眼通过调节光路使扫描光束能聚焦到眼底视网膜,以保证成像的横向分辨率和反射信号强度。当屈光补偿调节好后,再调节OCT设备的参考臂光程,以使扫描光束与参考光束光束的光程差在成像范围内。
专利号为WO2012145882A1的发明专利《一种眼科OCT系统和眼科OCT成像方法》,提出了一种屈光补偿镜位于扫描光路和固视光路的屈光补偿设计,可以通过被测人员对固视靶标清晰度的反馈进行得知是否已补偿好屈光再进行OCT成像。但该方法需要被测人员观测靶标并进行反馈,需要对被测人员进行先行培训且反馈过程中可能会有头部移动而影响后续OCT测量。专利号为CN102824159B的发明专利《眼底OCT自动跳去系统、设备和方法》公开了一种使用微透镜阵列并用探测器探测的光束聚焦点相邻两点间距离进行判定的自动调屈设备,但该方法需要额外的探测器而增加硬件成本。
发明内容
本申请旨在解决上述技术问题,提供一种使用OCT信号进行眼底屈光补偿判定与成像优化的方法,能够降低成本并实现快速精准的屈光补偿,同时对视网膜成像进行优化。
为实现上述目的,本申请公开了一种使用OCT信号进行眼底屈光补偿判定与成像优化的方法,该方法包括:
S1、控制参考臂光程:通过计算视网膜结构的凹凸性以及成像位置,在图像中找到视网膜结构,并将视网膜结构调整到图像中合适位置;在OCT图像显示中,将图像底部作为扫描光束与参考光束的光程差OPD为0的位置;在OCT成像时,根据需要选择视网膜层靠近光程差OPD=0的位置或脉络膜层靠近光程差OPD=0的位置;
S2、自动屈光补偿:控制参考臂光程与样品臂扫描光束光程,使参考臂光程与样品臂扫描光束光程的光程差为固定值,通过爬山法找到信号最强的图像,完成自动屈光补偿;
S3、视网膜结构左右微调:通过比较黄斑位置或图像中间位置左右区域强度值的大小,判定视网膜结构的左右偏移,根据偏移控制Y轴电机进行水平调整;
S4、调整视网膜成像位置:根据图像强度再次调节参考臂光程以调整视网膜在图像中的高低位置。
作为优选,所述S1控制参考臂光程具体包括:
S1-1、将光程电机与屈光电机移动至预设的初始位置,完成初始化;
S1-2、步进控制光程电机,步长为step1,每前进一步,从OCT设备中获取光谱数据并通过重建算法重建为二维眼底图像;计算视网膜结构的凹凸性与成像的大致位置,进入步骤S1-3;
S1-3、设定条件一判断图像凹凸性与大致成像位置,条件一:前一位置的视网膜方向为凹,前一位置的成像位置在图像较高处,且当前位置的视网膜方向为凹,当前位置的成像位置在图像较低处;重复步骤S1-2,并判断条件一,若条件一成立,则说明已在图像中找到视网膜结构,并且成像位置处于合适状态,进入步骤S2;若条件一未成立,进入步骤S1-4;
S1-4、改变光程电机移动方向,移动至预设的初始位置,步进光程电机,获取光谱数据及重建图像,并计算凹凸性与成像大致位置;
S1-5、设定条件二判断图像凹凸性与大致成像位置,条件二:前一位置的视网膜方向为凹,前一位置的成像位置在图像较低处,且当前位置的视网膜方向为凹,当前位置的成像位置在图像较高处;重复操作S1-4,判断条件二是否成立;若条件二成立,则前一位置处的图像即为合适状态,将光程电机移至前一位置处,进入步骤S2;若未成立,则自动聚焦失败,转为手动聚焦。
作为优选,所述S1-3和所述S1-5中判断图像凹凸性具体包括:
定义重建图像列数为cols,提取图像第列作为左侧数据,提取图像第列作为中间数据,提取图像第列作为右侧数据,计算左侧数据、中间数据、右侧数据中最大像素值所在行的行数和,定义左侧数据行数和为maxLeftSum,中间数据行数和为maxMediumSum,右侧数据行数和为maxRightSum;若则表示视网膜结构为凹;若则表示视网膜结构为凸。
作为优选,所述S1-3和所述S1-5中判断大致成像位置具体包括:
作为优选,所述S2自动屈光补偿中的爬山法具体包括如下过程:
T1、定义所述S2自动屈光补偿中的移动步长为step2,同时步进光程电机与屈光电机,计算干涉谱幅度;
T3、依次重复步骤T1、步骤T2,计算总步进次数,设定条件step’<0.1,若达到条件,则认定完成对焦,进入步骤S3;若总步进次数超过设定的次数阈值还未达到该条件,则认定对焦失败,转为手动对焦。
作为优选,所述T1中计算干涉谱幅值的方法具体包括:
将光谱信号减去其平均值得到去直流噪声后的干涉谱信号,将OCT设备得到的图像以黄斑或视盘为中心,设定感兴趣区域ROI为干涉谱中间200列数据,将数据叠加除以列数,即可得到平均每列幅值,该数值即为干涉谱幅度值。
作为优选,所述S3视网膜结构左右微调具体包括:
S3-1、设定图像中间200列为感兴趣区域ROI,通过计算粗略NFL层,确定中心位置;
S3-2、将中心位置左边100列的最大值所在行求和得到sumLeft;将中心位置右边100列的最大值所在行求和得到sumRight,定义sumLeft与sumRight之差为sumDiff;若sumDiff的绝对值小于所设定的阈值,则认为视网膜结构左右位置合适;若其大于设定的阈值,则需控制Y轴电机往对应方向进行移动,直至视网膜结构左右位置合适。
作为优选,所述S3-1具体包括:
在感兴趣区域ROI中每一列Ascan的像素中,检测到第一个灰度值大于设定值的像素位置,该位置即为粗略的NFL层位置;判断NFL层的凹凸性,若为凹,则将NFL层数值最大的位置定义为黄斑位置,将中心位置设为黄斑位置;若为凸,则认定黄斑处有病变,将中心位置设定为图像中心位置,进入所述S3-2。
作为优选,该种使用OCT信号进行眼底屈光补偿判定与成像优化的方法的系统装置包括:低相干光源SLD、高速光谱仪、参考臂光路、屈服补偿模块、目镜、瞳孔定位模块;所述参考臂光路包第一光纤准直器、反射镜,所述屈服补偿模块包括第二光纤准直器、二维扫描振镜、内固视靶标、透镜、二色反射镜,所述内固视靶标、所述透镜、所述二色反射镜、所述目镜、所述瞳孔定位模块的中心位于同一光轴连线上;所述瞳孔定位模块包括右侧相机、左侧相机,所述右侧相机、所述左侧相机对称设置与所述目镜的两侧,用于瞳孔定位;所述低相干光源SLD、所述高速光谱仪通过2×2光纤耦合器与所述第一光纤准直器、所述第二光纤准直器相连;所述反射镜安装于电动平移台上,用于沿光传播方向改变光程;所述屈服补偿模块安装于电动平移台上,用于沿所述二色反射镜、所述目镜的光轴连线方向移动,以同步改变所述二维扫描振镜和所述目镜、所述内固视靶标和所述目镜之间的距离。
作为优选,所述光程电机为安装所述反射镜的电动平移台上的控制电机,所述屈光电机为安装所述屈光补偿模块的电动平移台上的控制电机。
有益效果:本申请的使用OCT信号进行眼底屈光补偿判定与成像优化的方法,首先将扫描光束的屈光度预置在正常眼位置,再调节参考光束的光程以找到眼底图像,再固定图像的成像高度并通过优化方法自动找到具有最强干涉信号的位置,则该位置就是屈光补偿的最佳位置。屈光补偿过程完成后,可能仍有患者因眼动漂移而导致与OCT设备的光轴对准发生变化,使拍摄到的视网膜结构在图像中的位置会有较大倾斜而影响后续图像的判读。通过分析OCT信号,比较OCT图像在两侧的高度差,并控制电机调节OCT设备的光轴可使重新与瞳孔中心对准。不需要被测人员进行反馈,且不需要额外的探测与强度成像光路,有利于降低成本并实现快速精准的屈光补偿。综上,本申请不需要被测人员进行反馈,且不需要额外的探测与强度成像光路,有利于降低成本并实现快速精准的屈光补偿。同时克服了屈光补偿过程中患者因眼动漂移而导致视网膜结构在图像中的位置会有较大倾斜的问题,对视网膜成像进行了优化。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅是本申请的一些实施例,对于本领域技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本申请实施例的系统装置图。
图2是本申请实施例的算法流程图。
图3是本申请实施例中光程电机步进在图像中找到视网膜结构的部分中间过程图。
图4是本申请实施例中基于模拟眼的在不用屈光电机位置处的干涉频谱归一化幅值折线图。
图5是本申请实施例中光程电机与屈光电机步进完成对焦的部分中间过程图。
图6是本申请实施例中瞳孔定位微调前后图。
附图标记:1、第一光纤准直器;2、反射镜;3、第二光纤准直器;4、二维扫描振镜;5、内固视靶标;6、透镜;7、二色反射镜;8、目镜;9、右侧相机;10、左侧相机;11、人眼。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
实施例:本发明的系统装置图如图1所示,系统装置包括:低相干光源SLD、高速光谱仪、参考臂光路、屈服补偿模块、目镜8、瞳孔定位模块;参考臂光路包第一光纤准直器1、反射镜2,屈服补偿模块包括第二光纤准直器3、二维扫描振镜4、内固视靶标5、透镜6、二色反射镜7,内固视靶标5、透镜6、二色反射镜7、目镜8、瞳孔定位模块的中心位于同一光轴连线上。瞳孔定位模块包括右侧相机9、左侧相机10,右侧相机9、左侧相机10对称设置与目镜8的两侧,用于人眼11的瞳孔定位,当目镜8的光轴与人眼11的瞳孔中心对准时,在右侧相机9、左侧相机10拍摄的眼表影像中,瞳孔都位于图像正中心。
反射镜2安装于电动平移台上,用于沿光传播方向改变光程;屈服补偿模块安装于电动平移台上,用于沿二色反射镜7、目镜8的光轴连线方向移动,以同步改变二维扫描振镜4和目镜8、内固视靶标5和目镜8之间的距离,从而实现对不同人眼11的不同屈光度补偿
第一光纤准直器1用于将光纤出射光准直为平行光。反射镜2,并置于电动平移台上,可以沿光传播方向改变光程。第二光纤准直器3,将光纤出射光准直为平行光。二维扫描振镜4将光沿垂直于光轴的两个正交方向偏折。内固视靶标5,可由彩色液晶屏组成,显示不同颜色和大小的靶标共被测者观察。二色反射镜7,对于扫描光束全反射,对于固视光束全透射。右侧相机9用于拍摄眼表。左侧相机10用于拍摄眼表。可以控制反射镜2的移动改变参考臂的光程,第二光纤准直器3、二维扫描振镜4、二色反射镜7、目镜8、人眼11组成样品臂扫描光路,可以控制二维扫描振镜4的扫描方式实现线扫描、区域扫描和环形扫描等不同的OCT成像扫描方式。内固视靶标5、透镜6、二色反射镜7、目镜8、人眼11组成内固视光路,使人眼在屈光补偿好时可以看到清晰的靶标。
低相干光源SLD、高速光谱仪通过2×2光纤耦合器与第一光纤准直器1、第二光纤准直器3相连。低相干光源SLD用于提供OCT成像所需的宽带激光。2X2光纤耦合器用于将光源出射光分光至参考臂光路和样品臂扫描光路,并将参考臂和样品臂眼底返回的光信号干涉后导入高速光谱仪并记录采集干涉条纹。
光程电机为安装反射镜2的电动平移台上的控制电机,屈光电机为安装屈光补偿模块的电动平移台上的控制电机。光程电机控制参考臂光路的反射镜2位置,屈光电机控制屈光补偿模块的透镜6位置。本实施例的等光程位置在图像最下方,因此成像位置在图像下方1/4左右处较为合适,另外本实施例所拍摄到的图像尺寸为1000×1000。
请参阅图2,图2为本方法的算法流程图包括:
S1、控制光程电机,通过计算视网膜结构的凹凸性以及成像的大概位置,在图像中找到视网膜结构,并将视网膜结构调整到图像中合适位置。不失一般性的,在OCT图像显示中,定义图像底部为扫描光束与参考光束的光程差(Optical Path Difference,OPD)为0的位置。由于光谱仪相机对干涉光谱的离散采样特性,OCT图像在OPD绝对数值越大信噪比会越低,因此在OCT成像过程中会要求在避免图像出现翻折的前提下尽量使组织结构在图像中的位置尽量靠近光程差OPD=0的位置。在OCT成像时,可根据实际需要选择视网膜层靠近光程差OPD=0的位置或是脉络膜层靠近光程差OPD=0的位置。由于眼睛近似于椭球的生理结构,在OCT成像中若视网膜层靠近光程差OPD=0的位置,可以观测到图像是凸出状的;若脉络膜层靠近OPD=0的位置,可以观测到图像时凹陷状的。为防止OCT图像中组织结构在光程差OPD=0处出现翻折,可以将组织结构固定在光程差OPD>0.5mm的地方,当然也可以固定在其它光程差OPD位置。
S2、同时控制屈光与光程电机,使用爬山法找到信号最强的图像,确定对焦位置,完成自动聚焦。
S3、通过比较黄斑位置或图像中间位置左右区域强度值的大小,判定视网膜结构的左右偏移,根据偏移控制Y轴电机进行水平调整。
S4、根据图像强度再次调整视网膜在图像中的高低位置,以克服对焦过程中眼动对成像位置的影响。
具体地,S1步骤包括如下过程:
S1-1、将光程电机与屈光电机移动至初始位置,本发明所设计的OCT设备光程电机初始位置在55mm处,屈光电机在45mm处。
S1-2、步进光程电机,步长设定为0.5mm,每步进一步,从OCT设备中获取光谱数据以及重建后的图像,计算视网膜的凹凸方向以及成像的大致位置。
具体地,视网膜凹凸方向的计算方法如下:找出图像100、200、300列中灰度值最大像素的所在行,将其相加得到maxLeftSum;找出图像450、500、550列中灰度值最大像素的所在行,相加得到中间数据行数和maxMediumSum;找出图像700、800、900列中灰度值最大像素的所在行,相加得到右侧数据行数和maxRightSum。若则认为视网膜结构为凹;若则认为视网膜结构为凸。
成像大致位置的计算方式如下:将图像中的像素按行求和,找到求和值最大的行rowPosition,该行即为视网膜成像的大致位置。若rowPosition>900,则认为成像位置太低,若rowPosition<750,则认定成像位置太高,其余状态则认定成像位置合适。
S1-3、设定条件一:前一位置的视网膜方向为凹,前一位置的成像位置在图像较高处,并且当前位置的视网膜方向为凹,当前位置的成像位置在图像较低处。重复S1-2步骤,判断条件一。若条件一成立,则说明已在图像中找到视网膜结构,并且成像位置处于合适状态,进入下一阶段步骤S2。若条件一未成立,进入步骤S1-4。在本实施例中,在S1-3阶段未达到条件,进入S1-4步骤往反方向寻找视网膜结构。
S1-4、改变光程电机移动方向,移动至初始位置,步进光程电机,获取光谱数据及重建图像,并计算凹凸性与成像大致位置。
S1-5、设定条件二:前一位置的视网膜方向为凹,前一位置的成像位置在图像较低处,并且当前位置的视网膜方向为凹,当前位置的成像位置在图像较高处。重复S1-4,判断条件二是否成立。若条件二成立,则前一位置处的图像即为合适状态,将光程电机移至前一位置处,进入下一阶段步骤S2;若条件而未成立,则自动聚焦失败,转为手动聚焦。
请参阅图3,在本实施例中,在S1-4步骤中找到了视网膜结构,(a)图为找到视网膜结构前;(b)图为已经找到视网膜结构,但视网膜方向为凸,且成像位置过高,光程电机继续往该方向移动设定步长;(c)经过图像底部视网膜结构进行了翻折,继续步进;(d)此时视网膜方向变为凹,但成像位置太低,故光程电机继续步进;(e)该图中视网膜结构为凹,且该成像位置合适。但此步骤仅仅在图像中找到了视网膜结构并定位在了合适位置,从人眼视觉来看,视网膜信号较弱,分层不清晰,故需要进入下一阶段(即S2步骤)完成对焦。
S2步骤运用爬山法找到对焦位置的具体步骤包括:
T2、重复T1步骤,计算总步进次数,设定条件step2’<0.1,若达到条件,则认定完成对焦,进入下一阶段步骤S3;若总步进次数超过设定的次数阈值还未达到该条件,则认定对焦失败,转为手动对焦。
具体地,计算干涉谱幅值的具体步骤包括:
将光谱信号减去其平均值得到去直流噪声后的干涉谱信号,为更看清图像中间部分的视网膜结构(所设计OCT设备得到的图像一般以黄斑/视盘为中心),设定感兴趣区域ROI为干涉谱中间200列数据,将数据叠加除以列数,即可得到平均每列幅值,该数值即为本发明所定义的干涉谱幅度值。
请参阅图4,图中为基于模拟眼的在不用屈光电机位置处的干涉频谱归一化幅值折线图。X轴表示屈光电机位置,Y轴代表干涉频谱幅值归一化后的数值。图中三种曲线代表三种干涉谱幅值的评价方法,分别是计算原干涉谱信号强度、计算干涉谱经快速傅立叶变换所得线性信号强度、计算干涉谱经快速傅立叶变换在做对数变换所得log域信号强度。本例采用计算原干涉谱信号强度的评价方法,应该强调的是,本发明不仅局限于此评价方法,任何对以上评价方法的替代与变型都应当视为属于本发明的保护范围。
请参阅图5,(a)图为未进行对焦的OCT重建图像,视网膜信号强度较弱,分层不明显;(b)图为对焦过程中的图像,可见信号强度有所增强,但仍能提高;(c)图为最终自动对焦结果,视网膜图像更加清晰,信号强度更强,特别是图像中间的黄斑区域。
S3步骤中判定视网膜结构的左右偏移的方法见下:
S3-1、设定图像中间200列为感兴趣区域(ROI),通过计算粗略NFL层,确定中心位置。具体地,遍历ROI中每一列Ascan的像素,检测到第一个灰度值大于设定值的像素位置,该位置即为粗略的NFL层位置。判断NFL层的凹凸性,若为凹,则将NFL层数值最大的位置定义为黄斑位置,将中心位置设为黄斑位置;若为凸,则认定黄斑处有病变,将中心位置设定为图像中心位置。进入下一动作(即步骤S3-2);此实施例中,NFL层为凹,所设中心位置为黄斑位置,及NFL层中最低像素的所在行。
S3-2、找到黄斑位置左边100列的最大值,将最大值所在行求和得到sumLeft;找到黄斑位置右边100列的最大值,将最大值所在行求和得到sumRight;定义sumLeft与sumRight之差为sumDiff,若sumDiff的绝对值小于所设定的阈值,则认为视网膜结构左右位置合适;若其大于设定的阈值,则需控制Y轴电机往对应方向进行移动,直至视网膜结构左右位置合适。
请参阅图6,(a)图中视网膜结构左侧偏高,通过S3步骤计算左右位置高低,控制Y轴电机移动,进行视网膜结构的左右微调。(b)图中黄斑大致在图像中间位置,视网膜结构向左的倾斜被较好地校正。
S4步骤中判定图像高低位置的方法,与S1步骤中判定大致成像位置的方法相同。本实施例在S3步骤后,判定视网膜结构已经在图像的合适位置,故不做上下位置的移动控制。
以上描述是为了进行图示说明而不是为了进行限制。通过阅读上述描述,在所提供的示例之外的许多实施方式和许多应用对本领域技术人员来说都将是显而易见的。因此,本教导的范围不应该参照上述描述来确定,而是应该参照所附权利要求以及这些权利要求所拥有的等价物的全部范围来确定。出于全面之目的,所有文章和参考包括专利申请和公告的公开都通过参考结合在本文中。在前述权利要求中省略这里公开的主题的任何方面并不是为了放弃该主体内容,也不应该认为申请人没有将该主题考虑为所公开的申请主题的一部分。
Claims (5)
1.一种使用OCT信号进行眼底屈光补偿判定与成像优化的方法,其特征在于,该方法包括:
步骤S1、控制参考臂光程:通过计算视网膜结构的凹凸性以及成像位置,在图像中找到视网膜结构,并将视网膜结构调整到图像中合适位置;在OCT图像显示中,将图像底部作为扫描光束与参考光束的光程差OPD为0的位置;在OCT成像时,根据需要选择视网膜层靠近光程差OPD=0的位置或脉络膜层靠近光程差OPD=0的位置;
步骤S2、自动屈光补偿:控制参考臂光程与样品臂扫描光束光程,使参考臂光程与样品臂扫描光束光程的光程差为固定值,通过爬山法找到信号最强的图像,完成自动屈光补偿;
步骤S3、视网膜结构左右微调:通过比较黄斑位置或图像中间位置左右区域强度值的大小,判定视网膜结构的左右偏移,根据偏移控制Y轴电机进行水平调整;
步骤S4、调整视网膜成像位置:根据图像强度再次调节参考臂光程以调整视网膜在图像中的高低位置;
所述步骤S2自动屈光补偿中的爬山法具体包括如下过程:
步骤T1、定义所述步骤S2自动屈光补偿中的移动步长为step2,同时步进光程电机与屈光电机,计算干涉谱幅值;
步骤T3、依次重复步骤T1、步骤T2,计算总步进次数,设定条件step2’<0.1,若达到条件,则认定完成对焦,进入步骤S3;若总步进次数超过设定的次数阈值还未达到该条件,则认定对焦失败,转为手动对焦;
所述步骤T1中计算干涉谱幅值的方法具体包括:
将光谱信号减去其平均值得到去直流噪声后的干涉谱信号,将OCT设备得到的图像以黄斑或视盘为中心,设定感兴趣区域ROI为干涉谱中间200列数据,将数据叠加除以列数,即可得到平均每列幅值,该平均每列幅值即为干涉谱幅值;
该种使用OCT信号进行眼底屈光补偿判定与成像优化的方法的系统装置包括:低相干光源SLD、高速光谱仪、参考臂光路、屈光补偿模块、目镜(8)、瞳孔定位模块;所述参考臂光路包第一光纤准直器(1)、反射镜(2),所述屈光补偿模块包括第二光纤准直器(3)、二维扫描振镜(4)、内固视靶标(5)、透镜(6)、二色反射镜(7),所述内固视靶标(5)、所述透镜(6)、所述二色反射镜(7)、所述目镜(8)、所述瞳孔定位模块的中心位于同一光轴连线上;所述瞳孔定位模块包括右侧相机(9)、左侧相机(10),所述右侧相机(9)、所述左侧相机(10)对称设置与所述目镜(8)的两侧,用于瞳孔定位;所述低相干光源SLD、所述高速光谱仪通过2×2光纤耦合器与所述第一光纤准直器(1)、所述第二光纤准直器(3)相连;所述反射镜(2)安装于电动平移台上,用于沿光传播方向改变光程;所述屈光补偿模块安装于电动平移台上,用于沿所述二色反射镜(7)、所述目镜(8)的光轴连线方向移动,以同步改变所述二维扫描振镜(4)和所述目镜(8)、所述内固视靶标(5)和所述目镜(8)之间的距离。
2.根据权利要求1所述的使用OCT信号进行眼底屈光补偿判定与成像优化的方法,其特征在于,所述步骤S1控制参考臂光程具体包括:
步骤S1-1、将光程电机与屈光电机移动至预设的初始位置,完成初始化;
步骤S1-2、步进控制光程电机,步长为step1,每前进一步,从OCT设备中获取光谱数据并通过重建算法重建为二维眼底图像;计算视网膜结构的凹凸性与成像的大致位置,进入步骤S1-3;
步骤S1-3、设定条件一判断图像凹凸性与大致成像位置,条件一:前一位置的视网膜方向为凹,前一位置的成像位置在图像较高处,且当前位置的视网膜方向为凹,当前位置的成像位置在图像较低处;重复步骤S1-2,并判断条件一,若条件一成立,则说明已在图像中找到视网膜结构,并且成像位置处于合适状态,进入步骤S2;若条件一未成立,进入步骤S1-4;
步骤S1-4、改变光程电机移动方向,移动至预设的初始位置,步进光程电机,获取光谱数据及重建图像,并计算凹凸性与成像大致位置;
步骤S1-5、设定条件二判断图像凹凸性与大致成像位置,条件二:前一位置的视网膜方向为凹,前一位置的成像位置在图像较低处,且当前位置的视网膜方向为凹,当前位置的成像位置在图像较高处;重复操作步骤S1-4,判断条件二是否成立;若条件二成立,则前一位置处的图像即为合适状态,将光程电机移至前一位置处,进入步骤S2;若未成立,则自动聚焦失败,转为手动聚焦;
所述步骤S1-3和所述步骤S1-5中判断大致成像位置具体包括:
3.根据权利要求2所述的使用OCT信号进行眼底屈光补偿判定与成像优化的方法,其特征在于,所述光程电机为安装所述反射镜(2)的电动平移台上的控制电机,所述屈光电机为安装所述屈光补偿模块的电动平移台上的控制电机。
5.根据权利要求1所述的使用OCT信号进行眼底屈光补偿判定与成像优化的方法,其特征在于,所述步骤S3视网膜结构左右微调具体包括:
步骤S3-1、设定图像中间200列像素为感兴趣区域ROI,通过计算粗略NFL层,确定中心位置;
步骤S3-2、将中心位置左边100列像素的最大值所在行求和得到sumLeft;将中心位置右边100列的最大值所在行求和得到sumRight,定义sumLeft与sumRight之差为sumDiff;若sumDiff的绝对值小于所设定的阈值,则认为视网膜结构左右位置合适;若其大于设定的阈值,则需控制Y轴电机往对应方向进行移动,直至视网膜结构左右位置合适;
所述步骤S3-1具体包括:
在感兴趣区域ROI中每一列Ascan的像素中,检测到第一个灰度值大于设定值的像素位置,该位置即为粗略的神经纤维层位置;判断神经纤维层的凹凸性,若为凹,则将神经纤维层数值最大的位置定义为黄斑位置,将中心位置设为黄斑位置;若为凸,则认定黄斑处有病变,将中心位置设定为图像中心位置,进入所述步骤S3-2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011236869.8A CN112168132B (zh) | 2020-11-09 | 2020-11-09 | 一种使用oct信号进行眼底屈光补偿判定与成像优化的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011236869.8A CN112168132B (zh) | 2020-11-09 | 2020-11-09 | 一种使用oct信号进行眼底屈光补偿判定与成像优化的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112168132A CN112168132A (zh) | 2021-01-05 |
CN112168132B true CN112168132B (zh) | 2022-11-08 |
Family
ID=73917200
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011236869.8A Active CN112168132B (zh) | 2020-11-09 | 2020-11-09 | 一种使用oct信号进行眼底屈光补偿判定与成像优化的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112168132B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113040701A (zh) * | 2021-03-11 | 2021-06-29 | 视微影像(河南)科技有限公司 | 一种三维眼动追踪系统及其追踪方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009291253A (ja) * | 2008-06-02 | 2009-12-17 | Nidek Co Ltd | 眼科撮影装置 |
JP2012213489A (ja) * | 2011-03-31 | 2012-11-08 | Nidek Co Ltd | 眼底撮影装置 |
CN102824159A (zh) * | 2012-09-13 | 2012-12-19 | 深圳市斯尔顿科技有限公司 | 眼底oct自动调屈系统、设备和方法 |
CN203643682U (zh) * | 2012-04-24 | 2014-06-11 | 深圳市斯尔顿科技有限公司 | 一种眼科光学相干断层成像系统 |
CN103989453A (zh) * | 2014-06-03 | 2014-08-20 | 深圳市莫廷影像技术有限公司 | 一种多功能眼科测量装置及测试人眼不同部位的方法 |
CN104224109A (zh) * | 2014-10-16 | 2014-12-24 | 深圳市斯尔顿科技有限公司 | 一种结合oct系统的眼底相机 |
CN107224268A (zh) * | 2017-05-31 | 2017-10-03 | 执鼎医疗科技(杭州)有限公司 | 一种多光路成像自动联调方法及其装置 |
CN109480762A (zh) * | 2018-12-14 | 2019-03-19 | 中国科学院光电技术研究所 | 一种相干门与焦面重合并动态调焦的视网膜全场oct成像仪 |
JP2019103746A (ja) * | 2017-12-14 | 2019-06-27 | キヤノン株式会社 | 眼底撮像装置 |
CN110755031A (zh) * | 2019-10-31 | 2020-02-07 | 天津迈达医学科技股份有限公司 | 一种眼前后节频域光相干层析成像系统 |
US10653311B1 (en) * | 2019-06-12 | 2020-05-19 | Notal Vision Ltd. | Home OCT with automatic focus adjustment |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008161406A (ja) * | 2006-12-28 | 2008-07-17 | Topcon Corp | 眼科装置 |
-
2020
- 2020-11-09 CN CN202011236869.8A patent/CN112168132B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009291253A (ja) * | 2008-06-02 | 2009-12-17 | Nidek Co Ltd | 眼科撮影装置 |
JP2012213489A (ja) * | 2011-03-31 | 2012-11-08 | Nidek Co Ltd | 眼底撮影装置 |
CN203643682U (zh) * | 2012-04-24 | 2014-06-11 | 深圳市斯尔顿科技有限公司 | 一种眼科光学相干断层成像系统 |
CN102824159A (zh) * | 2012-09-13 | 2012-12-19 | 深圳市斯尔顿科技有限公司 | 眼底oct自动调屈系统、设备和方法 |
CN103989453A (zh) * | 2014-06-03 | 2014-08-20 | 深圳市莫廷影像技术有限公司 | 一种多功能眼科测量装置及测试人眼不同部位的方法 |
CN104224109A (zh) * | 2014-10-16 | 2014-12-24 | 深圳市斯尔顿科技有限公司 | 一种结合oct系统的眼底相机 |
CN107224268A (zh) * | 2017-05-31 | 2017-10-03 | 执鼎医疗科技(杭州)有限公司 | 一种多光路成像自动联调方法及其装置 |
JP2019103746A (ja) * | 2017-12-14 | 2019-06-27 | キヤノン株式会社 | 眼底撮像装置 |
CN109480762A (zh) * | 2018-12-14 | 2019-03-19 | 中国科学院光电技术研究所 | 一种相干门与焦面重合并动态调焦的视网膜全场oct成像仪 |
US10653311B1 (en) * | 2019-06-12 | 2020-05-19 | Notal Vision Ltd. | Home OCT with automatic focus adjustment |
CN110755031A (zh) * | 2019-10-31 | 2020-02-07 | 天津迈达医学科技股份有限公司 | 一种眼前后节频域光相干层析成像系统 |
Also Published As
Publication number | Publication date |
---|---|
CN112168132A (zh) | 2021-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11058293B2 (en) | Image processing apparatus, ophthalmologic imaging apparatus, image processing method, and storage medium | |
CN102421351B (zh) | 光学摄像设备和光学图像摄像方法 | |
US9931033B2 (en) | System and method for controlling a fundus imaging apparatus | |
KR101447010B1 (ko) | 안과장치 및 혈액 유속 산출방법 | |
US9875541B2 (en) | Enhanced algorithm for the detection of eye motion from fundus images | |
JP3706940B2 (ja) | 眼特性測定装置 | |
EP2460461A1 (en) | Determination of the retinal blood flow velocity | |
US11813022B2 (en) | Ophthalmic apparatus, controlling method thereof, and recording medium | |
EP2322080A1 (en) | Adaptive optics apparatus, adaptive optics method, and imaging apparatus | |
US11896305B2 (en) | Optical measurement systems and processes with fixation target having Bokeh compensation | |
US9082010B2 (en) | Apparatus and a method for processing an image of photoreceptor cells of a fundus of an eye | |
JP2017000469A (ja) | 前眼部光干渉断層撮影装置および前眼部光干渉断層撮影方法 | |
JP2019013362A (ja) | 眼科装置 | |
JP7286422B2 (ja) | 眼科情報処理装置、眼科装置、眼科情報処理方法、及びプログラム | |
CN112168132B (zh) | 一种使用oct信号进行眼底屈光补偿判定与成像优化的方法 | |
US10733735B2 (en) | Ophthalmic apparatus | |
JP2020156622A (ja) | 眼科装置 | |
JP7306978B2 (ja) | 眼科情報処理装置、眼科装置、眼科情報処理方法、及びプログラム | |
CN220124665U (zh) | 大视场验光装置及验光系统 | |
CN116616695B (zh) | 人眼光学特征获取方法、大视场验光装置及验光系统 | |
US11963722B2 (en) | Methods and systems for determining change in eye position between successive eye measurements | |
EP4209170A1 (en) | Ophthalmic apparatus | |
JP7306977B2 (ja) | 眼科装置、及びその制御方法 | |
EP4209171A1 (en) | Ophthalmic apparatus | |
CN116616695A (zh) | 人眼光学特征获取方法、大视场验光装置及验光系统 |
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 |