CN109620131B - 共光路微透镜阵列多光束光学相干弹性测量系统及方法 - Google Patents

共光路微透镜阵列多光束光学相干弹性测量系统及方法 Download PDF

Info

Publication number
CN109620131B
CN109620131B CN201811539190.9A CN201811539190A CN109620131B CN 109620131 B CN109620131 B CN 109620131B CN 201811539190 A CN201811539190 A CN 201811539190A CN 109620131 B CN109620131 B CN 109620131B
Authority
CN
China
Prior art keywords
cornea
path
common
subsystem
detection
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
Application number
CN201811539190.9A
Other languages
English (en)
Other versions
CN109620131A (zh
Inventor
蓝公仆
迈克尔·图
黄燕平
安林
许景江
秦嘉
谭海曙
陈国杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Foshan University
Original Assignee
Foshan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Foshan University filed Critical Foshan University
Priority to CN201811539190.9A priority Critical patent/CN109620131B/zh
Publication of CN109620131A publication Critical patent/CN109620131A/zh
Application granted granted Critical
Publication of CN109620131B publication Critical patent/CN109620131B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

Abstract

本发明提供了共光路微透镜阵列多光束光学相干弹性测量系统及方法，采用角膜定位子系统对角膜的激励和测量位置进行定位、采用角膜激励子系统使角膜产生微量的机械波和形变，采用基于共光路微透镜阵列的多光束OCT探测子系统对各采样点的机械波和形变进行高精度探测，实现角膜弹性参数的快速测量，并结合共光路设计，采用共光路微透镜阵列以及超快速光开关作为关键的控制元件，简化了系统结构、降低了制造成本，既能保证角膜的多方向、多位置的弹性响应的快速测量，也能保证高相位稳定性和高探测灵敏度，提高弹性测量的精度。

Description

共光路微透镜阵列多光束光学相干弹性测量系统及方法
技术领域
本发明涉及到生物医学弹性成像领域，尤其涉及基于共光路微透镜阵列多光束光学相干弹性测量系统及方法。
背景技术
人眼角膜是人眼主要的屈光介质，角膜的弹性力学特征对维持角膜的正常结构和功能具有重要作用、是研究角膜的生理和病理特性的重要基础。角膜疾病(例如圆锥角膜、角膜膨隆)及角膜手术(角膜屈光手术、紫外线交联术CXL等)会引起角膜弹性力学特征的改变。常规的临床检测方法依据角膜形态(地形图、厚度、曲率等)和眼压参数进行诊断；尽管对角膜疾病的检出率已经很高，但仍有部分角膜疾病被漏诊。而角膜结构的微小变化可引起弹性力学特征的明显改变，因此针对角膜弹性力学特征的量化研究在角膜疾病的诊断与治疗上具有重要意义。为了实现临床上角膜弹性力学特征的量化，研发非接触式、在体的人眼角膜弹性成像测量技术已成为眼科和视觉科学研究的一大需求和热点。
目前各种角膜的弹性测量技术仍不成熟，所测的角膜力学参数差别为几个量级。以兔子角膜的杨氏模量测量为例，其估算范围从大约1kPa(Thomasy etc.Acta Biomater10(2),785-791(2014))到大约11MPa(Wollensak etc.Acta Ophthalmol 87(1),48-51(2009))。Ruberti等提出了角膜弹性成像测量技术面临的若干个未解决问题(Openquestions)，其首要的三个问题为：“如何开发在体角膜弹性测量的新技术/仪器”、“如何区分在体角膜各区域的弹性力学特征”与“如何采用弹性成像测量的手段进行角膜手术的术前风险评估”(Ruberti etc.Annu Rev Biomed Eng 13,269-295(2011))。在角膜的在体弹性测量方案中，各种弹性成像技术均有其技术瓶颈：或无法达到高测量分辨率要求(Voorhees etc.Experimental Eye Research,160,85-95(2017))、或无法实现针对角膜病变区域的局部测量(Luce etc.J Cataract Refract Surg31(1),156-162(2005)；Hongetc.IOVS 54(1),659-665(2013))、或需要过长的测量时间(Scarcelli etc.NatPhotonics 2(1),39-43(2008)；Scarcelli etc.IOVS 53(1),185-190(2012))，难以达到角膜在体弹性测量的临床需求。目前角膜弹性参数量化研究所面临的一个重要科学问题是：如何更准确的量化在体人眼角膜的弹性力学特征，尤其是如何进行角膜局部区域弹性参数的高精度测量、实现临床上角膜病变或手术区域与正常区域的边界区分。这需要研究符合临床需求的角膜弹性测量的新方法，开发可实现在体人眼角膜局部区域和多方位弹性参数量化的高精度、快速测量技术。
而在OCE技术中，目前尚难以采用“扫描探测式”OCE实现人眼角膜机械波传播的追踪和杨氏模量的准确在体测量。机械波传播速度约为几米/秒，而每个测量点均需要数毫秒至数十毫秒的时间来获得该点的“位移-时间”曲线。OCE难以通过单帧成像实现机械波在某个传播方向的追踪。此外，眼动对活体角膜机械波传播速度的测量引入很大的测量误差。以往的SD-OCT系统采取多次激励和多次探测的方法对样品和离体角膜各测量点的弹性响应进行拼接，从而实现对机械波传播速度的估算。但由于眼动的存在，该方法难以用于在体人眼角膜的弹性测量。高速扫频OCT(Swept source OCT)系统可提高采集速度，例如150万A-line/秒(Song etc.Applied Physics Letters 108(19)(2016).；Singh etc.Opt.Lett.40(11),2588-2591(2015).)。然而扫频OCT目前仍有较大的相位误差，需要采用额外的稳相技术来获得稳定的相位。此外，高速扫频OCT价钱昂贵，尚难以在临床上推广。线视场(Linefield)扫描的OCE方案(Liu etc.Biomedical Optics Express 7(8),3021-3031(2016).)虽然可以提高机械波的探测速度，但是由于它需要更强的光强信号，目前也尚未能成功应用于人眼角膜的测量之中。
另外现有技术中也有一些能够提高采集速度的方案诸如：(1)高速扫频OCT(Sweptsource OCT)系统可提高采集速度，例如150万A-line/秒(Song etc.Applied PhysicsLetters 108(19)(2016).；Singh etc.Opt.Lett.40(11),2588-2591(2015).)。然而扫频OCT目前仍有较大的相位误差，需要采用额外的稳相技术来获得稳定的相位。此外，高速扫频OCT价钱昂贵，尚难以在临床上推广。(2)使用线视场(Line field)扫描的OCE方案(Liuetc.Biomedical Optics Express 7(8),3021-3031(2016).)虽然可以提高机械波的探测速度，但是由于它需要更强的光强信号，目前也尚未能成功应用于人眼角膜的测量之中。此外，通用的OCT方案参考臂和样品臂处于不同的空间位置上。OCT样品臂和参考臂之间振动态的差异会动态的改变两臂间的光程差、影响干涉信号稳定性。
发明内容
为此，本发明所要解决的技术问题是：提出了一种基于共光路微透镜阵列的多光束探测OCE系统，采用了关键光学元件——共光路微透镜阵列样品臂，既能保证角膜的多方向、多位置的弹性响应的快速测量，也能保证高相位稳定性和高探测灵敏度，提高弹性测量的精度。
本发明是通过以下技术方案实现的：一方面，提供共光路微透镜阵列多光束光学相干弹性测量系统，包括角膜定位子系统、载荷激励子系统和OCT探测子系统；所述角膜定位子系统对角膜测量位置进行定位、所述载荷激励子系统用于刺激角膜以使角膜产生机械波和形变，所述OCT探测子系统采用共光路微透镜阵列作为关键光学元件，对角膜的弹性响应进行测量；在所述OCT探测子系统中，宽光谱光源连接光纤耦合器的输入端，所述光纤耦合器的输出端连接光开关一端，所述光开关的另一端分别连接样品臂的对应光路，各光路分别对应相应的探测方位，所述样品臂的每条光路均设置有相互配合的准直镜和柱面镜，所述样品臂角膜前的位置还设置有共光路微透镜阵列，所述共光路微透镜阵列中间设置有激励孔，各方位探测光经过所述样品臂后分成多束探测光束，角膜多位置点反射回的探测光束和所述共光路微透镜阵列的参考平面反射回的参考光束发生干涉，其干涉光谱由光谱探测器进行探测，所述载荷激励子系统穿过所述共光路微透镜阵列的激励孔刺激角膜以使角膜产生机械波。
进一步，所述的多光束光学相干在体角膜弹性测量系统还包括角膜定位子系统；所述角膜定位子系统包括定位相机，其中，角膜沿相干测量系统轴向方向的动态位移由OCT进行追踪，角膜沿相干测量系统横向方向的动态位移由定位相机进行跟踪定位。
进一步，所述样品臂对角膜激励点周围不少于4个方位进行测量；所述共光路微透镜阵列将每个探测方位上的光分为若干个测量点，其采样测量点数不少于2个，以对角膜该方位上的弹性反应进行采样；所述共光路微透镜阵列靠近角膜的一侧提供了参考平面，以反射参考信号与角膜多采样点反射回的探测信号发射干涉。
进一步，所述各方位各采样点的探测光束与所述参考光束的光程差不同，其最大光程差控制在所述OCT探测子系统可测光程差范围之内。
进一步，所述的光开关是超快速的光开关，其开/关转换时间不超过1ms，实现各方位按顺序的交替重复采集。
另一方面，共光路微透镜阵列多光束光学相干弹性测量方法，基于上述技术方案中的共光路微透镜阵列多光束光学相干弹性测量系统，包括步骤：
S1、使用定位相机对人眼角膜及瞳孔成像，根据瞳孔的边缘及中心位置对角膜进行横向定位，并通过OCT信号中A-scan上标志点的深度变化进行轴向追踪，当角膜的轴向位置和横向位置在测量误差范围之内，开始通过载荷激励子系统对样品角膜表面进行激励，并通过OCT探测子系统对产生的机械波及角膜的微量形变进行探测，评估其弹性特性；
S2、在OCT探测子系统中，宽光谱光源发出的光通过光纤传入光纤耦合器，测量光在光纤耦合器内分光后，通向超快速光开关，超快速光开关控制了多方位的探测光路；每路探测光路对应了一个测量方向，各探测光路的光经过准直镜与柱镜的组合作用后，其光束被延展成了长条形光斑，并分别照射到角膜前的微透镜阵列的对应方位上，并由微透镜阵列分成若干个测量点，聚焦于角膜上；各路探测光在超快速光开关的控制下，对角膜激励点周围多方位按顺序进行交替重复采集；角膜多位置点反射回的探测光束和共光路微透镜阵列的参考平面反射回的参考光束发生干涉，并被光谱探测器探测。
进一步，获取各测量点的图像强度信号和随时间变化的相位信号，通过计算各测量点的相位信号可以得到该测量点的表面形变信息。
进一步，获取被测量角膜的恢复曲线，并通过将恢复曲线拟合成指数衰减曲线，得到角膜的固有频率。
进一步，获取角膜的受激阻尼振动信息，进行快速傅里叶变换，得到角膜的固有频率。
进一步，获取各测量点的形变信息及时间延迟信息，根据机械波传播模型计算出杨氏模量。
本发明采用角膜定位子系统对角膜的激励和测量位置进行定位、采用角膜激励子系统使角膜产生微量的机械波和形变，采用基于共光路微透镜阵列的多光束OCT探测子系统对各采样点的机械波和形变进行高精度探测，实现角膜弹性参数的快速测量，同时以超快速光开关作为控制元件，简化了系统结构、降低了制造成本，既能保证角膜的多方向、多位置的弹性响应的快速测量，也能保证高相位稳定性和高探测灵敏度，提高弹性测量的精度。
附图说明
图1是本发明所述共光路微透镜阵列的多光束光学相干弹性测量系统的结构示意图；
图2是共光路微透镜阵列的结构示意图；
图3是共光路微透镜阵列的截面结构示意图；
图4是样品臂的结构示意图；
图5是激励点和测量点在角膜上的放大位置示意图；
图6是光谱探测器对角膜四个方位的多点干涉叠加信号按顺序进行交替重复收集图；
图7是对各测量点获取的信号图，(a)为各测量点的强度图；(b)为各测量点的相位图。
具体实施方式
下面结合附图对本发明进行详细描述。
结合图1、图2、图3、图4、图5和图6，本实施例提供了一种共光路的多光束光学相干弹性测量系统，包括载荷激励子系统1和OCT探测子系统，所述OCT探测子系统包括宽光谱光源21，所述宽光谱光源21连接光纤耦合器22的输入端，所述光纤耦合器22的输出端连接超快速光开关23一端，所述超快速光开关23的另一端分别连接样品臂24对应的光路。超快速光开关的开/关转换时间不超过1ms，本实施例的光开关的开/关转换时间为50nm,实现了角膜激励点周围各方位按顺序的快速交替重复采集。
本实施例中，样品臂24的每条光路均设置有相互配合的准直镜241和柱面镜242，将光纤发出的光延展成长条形光斑，并分别照射到角膜前的多方位，本实施例设置为4方位的共光路微透镜阵列243的对应方位上。所述共光路微透镜阵列243靠近角膜61的一侧为参考平面2432，共光路微透镜阵列243的中间设有激励孔2433，在共光路微透镜阵列243上，各方向的长条型光斑有部分光被微透镜阵列2431聚焦在角膜61各测量点上并反射回系统，这部分光为探测光；另有部分光被参考平面2432反射回来，为参考光，角膜多位置点反射回的探测光束和所述样品臂内反射回的参考光束发生干涉，其干涉光谱由与所述光纤耦合器连接的光谱探测器25进行探测。
在本共光路设计的探测光路中，各方位的光经过共光路微透镜阵列243该方位上多个(本实施例为4个)焦距不同且光程差各异的微透镜单元2431后，聚焦在角膜该方位上横向距离3mm到4mm区域内，并对该区域进行采样。各透镜焦距针对人眼角膜模型设计，使各光束更有效的聚焦在角膜的各测量点上。本设计在共光路微透镜阵列243的靠近角膜61侧的平面镀反射膜，形成了参考平面2432。参考平面2432与角膜61顶端距离约为0.3mm到1mm。
所述角膜定位子系统用于对待测角膜61的位置进行定位，所述角膜定位子系统包括定位相机31，其中，角膜61沿相干测量系统轴向方向的动态位移由OCT进行追踪，角膜沿相干测量系统横向方向的动态位移由定位相机进行跟踪定位。鉴于横向眼动会影响角膜的测量位置精度，本发明采用定位相机31对人眼瞳孔成像，并根据瞳孔的边缘及中心点位置对角膜进行横向定位，减少由横向眼动引入的测量位置误差。
所述样品臂24对角膜61的激励点O周围不少于4个方位进行测量，每个方位上的测量点数不少于2个。本实施例采用了四个方位的探测。分别为方位1、方位2、方位3和方位4，为了方便描述，图1中只画出方位1和方位3的样品臂24的光路。在同一个方位中，具有4个测量点，以方位1为例，则为测量点A、测量点B、测量点C和测量点D。
测量的初始位置为：以角膜顶点为激励点O，对角膜的上方(Superiordirection)、下方(Inferior direction)、颞侧(Temporal direction)和鼻侧(Nasaldirection)进行探测。各透镜焦距针对人眼角膜模型设计，使各光束更有效的聚焦在角膜的各测量点上。并且各透镜设计均引入特征光程，便于区分各测量点信号。
所述载荷激励子系统1为适合人眼角膜激励的装置。载荷激励子系统1可为适合人眼角膜激励的装置，例如微量气体脉冲等。该装置应对眼睛尤其是角膜组织无毒无害，具有安全性。载荷激励子系统1的作用是对角膜61激励，使角膜61产生微量的受激形变和机械波。
所述探测光束与所述参考光束的光程差不同，其最大光程差控制在所述OCT探测子系统可测光程差之内。
相应的，使用以上装置，本实施例还提供了基于共光路微透镜阵列的多光束相干弹性测量方法，包括步骤：
S1、使用定位相机31对人眼角膜及瞳孔成像，根据瞳孔的边缘及中心位置对角膜进行横向定位，并通过OCT信号中A-scan上标志点的深度变化进行轴向追踪，当角膜61的轴向位置和横向位置在测量误差范围之内，开始通过载荷激励子系统对角膜61表面进行激励，并通过OCT探测子系统对产生的机械波及角膜61的微量形变进行探测，评估其弹性特性。
S2、在OCT探测子系统中，宽光谱光源21发出的光通过光纤传入光纤耦合器22，测量光在光纤耦合器22内分光后，通向超快速光开关23，超快速光开关23控制了多方位的探测光路；每路探测光路对应了一个测量方向，各探测光路的光经过准直镜241与柱镜242的组合作用后，其光束被延展成了长条形光斑，并分别照射到角膜61前的共光路微透镜阵列243的对应方位上，并由共光路微透镜阵列分成若干个测量点，聚焦于角膜61上；各路探测光在超快速光开关23的控制下，对角膜激励点O周围多方位按顺序进行交替重复采集；角膜多位置点反射回的探测光束和共光路微透镜阵列243的参考平面2432反射回的参考光束发生干涉，并被光谱探测器25探测。
光谱探测器25获得的信息解析之后可以得到：各测量点随时间变化的相位信号随着距离的增大其形变幅度递减，且各点位信号存在时间延迟，通过计算得到表面形变信息：
Figure BDA0001906216540000081
其中tJ和t0是在一系列A-scan信号的时间节点，t0是参考时间点，λ0是中心波长，在计算出形变信息后可以根据主要形变的幅值估算角膜的硬度。
结合图7(a)与(b)给出的针对同时采集四测量点信号的模拟示意图，由此可以根据各测量点与参考平面2432的光程差特征(包括多光束设计时引入的光程差特征)确定各点位信号的干涉特征，对包括伪信号在内的干扰噪声进行抑制或滤除。主要形变(图7(b)中曲线)幅值受激励载荷的直接影响。在同样的驱动力下，相对硬的样品形变更小；而同等幅值的主要形变在相对硬的样品中衰减更快。因此可以通过主要形变的幅值来判定样品(角膜)的软硬程度。本项目针对主要形变幅值的测量方案为：1、测量同等气体脉冲压强下，测量点A处样品(角膜)主要形变的幅值大小。2、测量一定激励压强范围下，测量点A处样品(角膜)主要形变幅值随着气体脉冲压强变化的曲线。3、测量各测量点(例如测量点A到D)主要形变幅值随着传播距离增大的衰减。由此可以估算样品的硬度。
光谱探测器25获得的信息解析之后可以得到：测量角膜的恢复曲线并通过将恢复曲线拟合成指数衰减曲线，实现固有频率的间接测量；也可以对角膜的受激阻尼振动进行高分辨探测和快速傅里叶变换(FFT)来得到角膜的固有频率。
具体的，根据恢复曲线测量固有频率的方法为：由于恢复曲线(图7(b)中曲线)与样品的黏弹性(viscoelasticity)有关，可根据动态模型(该动态模型为文献WuC.etc.IOVS.2015,56(2):1292-1300中记载的动态模型)将恢复曲线拟合成指数衰减曲线，实现固有频率的间接测量。
Figure BDA0001906216540000091
其中，ξ是衰减系数(Damping Ratio)，f是固有频率。对阻尼振动的微分方程(2)可根据三种情况进行求解，分别为：临界阻尼(critical-damping，ξ＝1)、欠阻尼(under-damping，0≤ξ<1)和过阻尼(over-damping，ξ>1)。
Figure BDA0001906216540000092
其中，幅值常数A和B由恢复曲线的指数拟合所得。
根据阻尼振动测量固有频率的方法为：FFT的频率分辨率(f0)表示的是其所能分辨的最小频率间隔。可由以下公式表示：
Figure BDA0001906216540000093
其中，fs为采样频率，ts为采样时间(即时间分辨率)，N为采样点数，T为采样时间长度。为了提高FFT的频率分辨率，可在FFT之前对数据进行如下两步处理。1、补零(Zeropadding)：在OCT的采样时间内(例如30ms)，角膜的阻尼振动幅度将逐渐趋于零。因此可对数据采用补零的手段，拓展采样时间。2、拼接：共光路OCT具有超稳定的相位，弹性样品受激前后相位基本保持不变。因此可将相邻激励周期内的数据(已补零)通过拼接进行周期拓展，提高总采样时间及采样点数。
根据机械波传播速度计算杨氏模量的方法为：光谱探测器25获得的信息解析之后可以得到：机械波传播模型，并通过机械波传播模型计算出该区域杨氏模量：
Figure BDA0001906216540000101
其中ci,j＝(di-dj)/(ti-tj)，i、j代表任意两个测量点，di和dj表示任意两个测量点与受激点沿角膜表面的距离，ti和tj表示机械波在任意两测量点间的传播时间，ci,j表示波在两点间的传播速度，ρ是密度，ν是泊松比(≈0.5)。
当测量系统的探测头激励进行旋转时，可以实现对角膜的各个方向的弹性参数进行测量，角膜表面各方向上弹性特性的差别体现了角膜的各向异性。
以上结合附图所描述的实施例仅仅是本发明的优选实施方式，并不是对本发明的保护范围的设定，任何基于本发明的思路所做的改进都理应在本发明的保护范围之内。

Claims (5)

1.共光路微透镜阵列多光束光学相干弹性测量方法，利用共光路微透镜阵列多光束光学相干弹性测量系统，所述共光路微透镜阵列多光束光学相干弹性测量系统包括：包括角膜定位子系统、载荷激励子系统和OCT探测子系统；所述角膜定位子系统对角膜测量位置进行定位、所述载荷激励子系统用于刺激角膜以使角膜产生机械波和形变，所述OCT探测子系统采用共光路微透镜阵列作为关键光学元件，对角膜的弹性响应进行测量；在所述OCT探测子系统中，宽光谱光源连接光纤耦合器的输入端，所述光纤耦合器的输出端连接光开关一端，所述光开关的另一端分别连接样品臂的对应光路，各光路分别对应相应的探测方位，所述样品臂的每条光路均设置有相互配合的准直镜和柱面镜，样品臂角膜前的位置还设置有共光路微透镜阵列，所述共光路微透镜阵列中间设置有激励孔，各方位探测光经过所述样品臂后分成多束探测光束，角膜多位置点反射回的探测光束和所述共光路微透镜阵列的参考平面反射回的参考光束发生干涉，其干涉光谱由光谱探测器进行探测，所述载荷激励子系统穿过所述共光路微透镜阵列的激励孔刺激角膜以使角膜产生机械波；
所述的多光束光学相干在体角膜弹性测量系统还包括角膜定位子系统；所述角膜定位子系统包括定位相机，其中，角膜沿相干测量系统轴向方向的动态位移由OCT进行追踪，角膜沿相干测量系统横向方向的动态位移由定位相机进行跟踪定位；
所述样品臂对角膜激励点周围不少于4个方位进行测量；所述共光路微透镜阵列将每个探测方位上的光分为若干个测量点，其采样测量点数不少于2个，以对角膜该方位上的弹性反应进行采样；所述共光路微透镜阵列靠近角膜的一侧提供了参考平面，以反射参考信号与角膜多采样点反射回的探测信号发射干涉；
所述各方位各采样点的探测光束与所述参考光束的光程差不同，其最大光程差控制在所述OCT探测子系统可测光程差范围之内；
所述的光开关是超快速的光开关，其开/关转换时间不超过1 ms，实现各方位按顺序的交替重复采集；
其特征在于，还包括步骤：
S1、使用定位相机对人眼角膜及瞳孔成像，根据瞳孔的边缘及中心位置对角膜进行横向定位，并通过OCT信号中A-scan上标志点的深度变化进行轴向追踪，当角膜的轴向位置和横向位置在测量误差范围之内，开始通过载荷激励子系统对样品角膜表面进行激励，并通过OCT探测子系统对产生的机械波及角膜的微量形变进行探测，评估其弹性特性；
S2、在OCT探测子系统中，宽光谱光源发出的光通过光纤传入光纤耦合器，测量光在光纤耦合器内分光后，通向超快速光开关，超快速光开关控制了多方位的探测光路；每路探测光路对应了一个测量方向，各探测光路的光经过准直镜与柱镜的组合作用后，其光束被延展成了长条形光斑，并分别照射到角膜前的微透镜阵列的对应方位上，并由微透镜阵列分成若干个测量点，聚焦于角膜上；各路探测光在超快速光开关的控制下，对角膜激励点周围多方位按顺序进行交替重复采集；角膜多位置点反射回的探测光束和共光路微透镜阵列的参考平面反射回的参考光束发生干涉，并被光谱探测器探测。
2.根据权利要求1所述的共光路微透镜阵列多光束光学相干弹性测量方法，其特征在于：获取各测量点的图像强度信号和随时间变化的相位信号，通过计算各测量点的相位信号可以得到该测量点的表面形变信息。
3.根据权利要求1所述的共光路微透镜阵列多光束光学相干弹性测量方法，其特征在于：获取被测量角膜的恢复曲线，并通过将恢复曲线拟合成指数衰减曲线，得到角膜的固有频率。
4.根据权利要求1所述的共光路微透镜阵列多光束光学相干弹性测量方法，其特征在于：获取角膜的受激阻尼振动信息，进行快速傅里叶变换，得到角膜的固有频率。
5.根据权利要求1所述的共光路微透镜阵列多光束光学相干弹性测量方法，其特征在于：获取各测量点的形变信息及时间延迟信息，根据机械波传播模型计算出杨氏模量。
CN201811539190.9A 2018-12-14 2018-12-14 共光路微透镜阵列多光束光学相干弹性测量系统及方法 Active CN109620131B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811539190.9A CN109620131B (zh) 2018-12-14 2018-12-14 共光路微透镜阵列多光束光学相干弹性测量系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811539190.9A CN109620131B (zh) 2018-12-14 2018-12-14 共光路微透镜阵列多光束光学相干弹性测量系统及方法

Publications (2)

Publication Number Publication Date
CN109620131A CN109620131A (zh) 2019-04-16
CN109620131B true CN109620131B (zh) 2021-08-03

Family

ID=66074463

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811539190.9A Active CN109620131B (zh) 2018-12-14 2018-12-14 共光路微透镜阵列多光束光学相干弹性测量系统及方法

Country Status (1)

Country Link
CN (1) CN109620131B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110270009B (zh) * 2019-07-18 2021-06-04 上海交通大学 用于视网膜局部电刺激的眼球外电极阵列设备
CN112683812B (zh) * 2020-12-09 2023-10-31 佛山科学技术学院 一种非线性采样的多光束光学相干弹性测量系统及方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1713020A (zh) * 2005-07-21 2005-12-28 中国科学院光电技术研究所 利用微透镜或微棱镜阵列进行扫描的光学相干层析系统
CN201481402U (zh) * 2009-08-06 2010-05-26 浙江大学 一种oct中基于时空分光的宽光谱高分辨探测系统
CN103565401A (zh) * 2012-07-27 2014-02-12 上海威景生物科技有限公司 一种全眼球光学相干层析自适应成像系统及方法
CN103688133A (zh) * 2011-05-20 2014-03-26 梅德路米克斯有限公司 用于低相干干涉测量的扫描设备
CN104198030A (zh) * 2014-08-29 2014-12-10 中国石油天然气股份有限公司 基于相干瑞利散射的多路振动检测方法及其检测系统
CN104434028A (zh) * 2014-11-15 2015-03-25 中国科学院光电技术研究所 角膜弹性成像与眼前节结构成像相结合的系统与方法
CN104854423A (zh) * 2012-12-06 2015-08-19 周超 空分复用光学相干断层扫描设备及方法
WO2016120401A1 (de) * 2015-01-28 2016-08-04 Carl Zeiss Meditec Ag Optische kohärenztomographie zur messung an der retina
WO2017137567A1 (en) * 2016-02-12 2017-08-17 Carl Zeiss Meditec, Inc. Systems and methods for improved oct measurements

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8348429B2 (en) * 2008-03-27 2013-01-08 Doheny Eye Institute Optical coherence tomography device, method, and system
CN201653399U (zh) * 2009-09-30 2010-11-24 深圳市斯尔顿科技有限公司 一种光学测距装置及其便携式验光仪
WO2013126653A1 (en) * 2012-02-22 2013-08-29 Amo Development, Llc Preformed lens systems and methods
US20150011895A1 (en) * 2012-03-28 2015-01-08 University Of Washington Through Its Center For Commercialization Methods and Systems for Determining Mechanical Properties of a Tissue
CN102914259A (zh) * 2012-09-29 2013-02-06 哈尔滨工程大学 基于分光同步移相的干涉检测装置及检测方法
KR102064190B1 (ko) * 2013-02-28 2020-01-09 주식회사 휴비츠 광 간섭 단층촬영 장치
DE102014014566A1 (de) * 2014-09-29 2016-03-31 Carl Zeiss Meditec Ag Augenchirurgisches Verfahren
US10548479B2 (en) * 2015-02-16 2020-02-04 The Regents Of The University Of California Optical coherence elastography (OCE) method under acoustic radiation force excitation using OCT doppler variance methods and OCT correlation-based methods
CN106963338B (zh) * 2017-03-31 2018-09-07 杭州电子科技大学 晶状体生物力学和光学特性无创在体成像系统及测量方法
CN107260120A (zh) * 2017-05-31 2017-10-20 执鼎医疗科技（杭州）有限公司 一种相机与oct共光路成像装置
CN107468214B (zh) * 2017-08-22 2024-04-16 广东唯仁医疗科技有限公司 一种光学相干断层成像系统及光学相干断层成像方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1713020A (zh) * 2005-07-21 2005-12-28 中国科学院光电技术研究所 利用微透镜或微棱镜阵列进行扫描的光学相干层析系统
CN201481402U (zh) * 2009-08-06 2010-05-26 浙江大学 一种oct中基于时空分光的宽光谱高分辨探测系统
CN103688133A (zh) * 2011-05-20 2014-03-26 梅德路米克斯有限公司 用于低相干干涉测量的扫描设备
CN103565401A (zh) * 2012-07-27 2014-02-12 上海威景生物科技有限公司 一种全眼球光学相干层析自适应成像系统及方法
CN104854423A (zh) * 2012-12-06 2015-08-19 周超 空分复用光学相干断层扫描设备及方法
CN104198030A (zh) * 2014-08-29 2014-12-10 中国石油天然气股份有限公司 基于相干瑞利散射的多路振动检测方法及其检测系统
CN104434028A (zh) * 2014-11-15 2015-03-25 中国科学院光电技术研究所 角膜弹性成像与眼前节结构成像相结合的系统与方法
WO2016120401A1 (de) * 2015-01-28 2016-08-04 Carl Zeiss Meditec Ag Optische kohärenztomographie zur messung an der retina
WO2017137567A1 (en) * 2016-02-12 2017-08-17 Carl Zeiss Meditec, Inc. Systems and methods for improved oct measurements

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Common-path phase-sensitive optical coherence tomography provides enhanced phase stability and detection sensitivity for dynamic elastography;Gongpu Lan，Manmohan Singh,Kirill V.Larin,Michael D.Twa;《Biomedical Optics Express》;20171101;第8卷(第11期);5253-5266 *
In-vivo human corneal elasticity imaging: a phase sensitive optical coherence elastography method;Twa, Michael D.,Lan, Gongpu,Singh, Manmohan,Larin, Kirill;《INVESTIGATIVE OPHTHALMOLOGY & VISUAL SCIENCE》;20170630;第58卷(第8期);4324 *
光折变晶体中共面多光束干涉研究;陈宝东，温静，张宝光;《红外与激光工程》;20120131;第41卷(第1期);133-138 *
基于补偿干涉仪的串联式全场光学相干层析系统;郭英呈，高万荣，朱越;《激光与光电子学进展》;20171231;121-128 *

Also Published As

Publication number Publication date
CN109620131A (zh) 2019-04-16

Similar Documents

Publication Publication Date Title
CN109645954B (zh) 基于微透镜阵列的多光束光学相干的弹性测量系统及方法
Ambroziński et al. Acoustic micro-tapping for non-contact 4D imaging of tissue elasticity
CN109620132B (zh) 一种多探测光束光学相干在体角膜弹性测量系统及方法
US10197379B2 (en) Optical sensing based on measurements of displacements induced by optical forces in viscoelastic media using phase-sensitive optical coherence tomography
JP4777362B2 (ja) 光学コヒーレンス・トモグラフィ撮像における動き補正方法
US8570527B2 (en) Polarization-sensitive spectral interferometry
JP5623028B2 (ja) 光干渉断層画像を撮る撮像方法及びその装置
JP4892719B2 (ja) 断層計測装置及び断層計測方法
US10072920B2 (en) Optical sensing based on measurements of displacements induced by optical scattering forces in viscoelastic media using phase-sensitive optical coherence tomography
US9492077B2 (en) Optimized device for swept source optical coherence domain reflectometry and tomography
JP6026406B2 (ja) 眼構造の画像を改善するための装置
US8740795B2 (en) Reflective non-contact ocular pulse analyzer for clinical diagnosis of eye and cerebrovascular disease
EP2801814A1 (en) Swept source optical coherence tomograph and method for stabilizing phase thereof
CN109620130B (zh) 一种共光路的多光束光学相干弹性测量系统及测量方法
CN106963338A (zh) 晶状体生物力学和光学特性无创在体成像系统及测量方法
CN104434028A (zh) 角膜弹性成像与眼前节结构成像相结合的系统与方法
JP2014505562A (ja) 光干渉断層撮影によって物体の内部寸法を測定するためのシステムおよび方法
CN109620131B (zh) 共光路微透镜阵列多光束光学相干弹性测量系统及方法
US20150011895A1 (en) Methods and Systems for Determining Mechanical Properties of a Tissue
CN109674441B (zh) 基于光开关和微透镜阵列的多光束弹性测量系统及方法
RU187692U1 (ru) Устройство эндоскопической оптической когерентной томографии с коррекцией волнового фронта
CN112826440B (zh) 基于光学相干法的角膜结构-眼压-多力学参数测量仪
CN112683812B (zh) 一种非线性采样的多光束光学相干弹性测量系统及方法
WO2018140703A1 (en) Common-path phase-sensitive optical coherence tomography
Drexler¹ et al. New technology for ultrahigh resolution optical coherence tomography of the retina

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