CN104434028B - 角膜弹性成像与眼前节结构成像相结合的系统与方法 - Google Patents
角膜弹性成像与眼前节结构成像相结合的系统与方法 Download PDFInfo
- Publication number
- CN104434028B CN104434028B CN201410654412.7A CN201410654412A CN104434028B CN 104434028 B CN104434028 B CN 104434028B CN 201410654412 A CN201410654412 A CN 201410654412A CN 104434028 B CN104434028 B CN 104434028B
- Authority
- CN
- China
- Prior art keywords
- imaging
- cornea
- anterior segment
- segment structure
- signal
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 87
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000002091 elastography Methods 0.000 title abstract 2
- 210000004087 cornea Anatomy 0.000 claims abstract description 65
- 230000005284 excitation Effects 0.000 claims abstract description 38
- 238000001228 spectrum Methods 0.000 claims abstract description 38
- 230000008569 process Effects 0.000 claims abstract description 21
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 13
- 230000008859 change Effects 0.000 claims abstract description 8
- 239000000523 sample Substances 0.000 claims description 46
- 230000003287 optical effect Effects 0.000 claims description 38
- 239000013307 optical fiber Substances 0.000 claims description 15
- 238000005070 sampling Methods 0.000 claims description 14
- 230000001360 synchronised effect Effects 0.000 claims description 10
- 238000005286 illumination Methods 0.000 claims description 9
- 230000000630 rising effect Effects 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 6
- 238000013519 translation Methods 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 4
- 230000001427 coherent effect Effects 0.000 claims description 4
- 238000013481 data capture Methods 0.000 claims description 4
- 239000011521 glass Substances 0.000 claims description 4
- 230000008054 signal transmission Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000003595 spectral effect Effects 0.000 claims description 2
- 238000003745 diagnosis Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 abstract description 6
- 238000001356 surgical procedure Methods 0.000 abstract description 6
- 201000010099 disease Diseases 0.000 abstract description 5
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 abstract description 5
- 238000013461 design Methods 0.000 abstract description 3
- 208000021921 corneal disease Diseases 0.000 abstract description 2
- 238000013399 early diagnosis Methods 0.000 abstract description 2
- 230000002980 postoperative effect Effects 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 abstract description 2
- 238000011156 evaluation Methods 0.000 abstract 1
- 239000000835 fiber Substances 0.000 description 20
- 238000005516 engineering process Methods 0.000 description 14
- 230000006870 function Effects 0.000 description 9
- 238000001514 detection method Methods 0.000 description 8
- 238000006073 displacement reaction Methods 0.000 description 7
- 230000010287 polarization Effects 0.000 description 6
- 239000006185 dispersion Substances 0.000 description 5
- 238000005481 NMR spectroscopy Methods 0.000 description 4
- 230000007812 deficiency Effects 0.000 description 4
- 230000007935 neutral effect Effects 0.000 description 4
- 210000004872 soft tissue Anatomy 0.000 description 4
- 210000001519 tissue Anatomy 0.000 description 4
- 238000012546 transfer Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000004587 chromatography analysis Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 206010023365 keratopathy Diseases 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012014 optical coherence tomography Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 206010006187 Breast cancer Diseases 0.000 description 1
- 208000026310 Breast neoplasm Diseases 0.000 description 1
- 206010010996 Corneal degeneration Diseases 0.000 description 1
- 201000002287 Keratoconus Diseases 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 210000002159 anterior chamber Anatomy 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 201000008275 breast carcinoma Diseases 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 208000030533 eye disease Diseases 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000013011 mating Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008450 motivation Effects 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000010897 surface acoustic wave method Methods 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Eye Examination Apparatus (AREA)
Abstract
本发明公开了一种角膜弹性成像与眼前节结构成像相结合的系统与方法。由喷气激励系统和可横向二维扫描的扫频OCT成像系统组成；喷气激励系统提供与成像光束同轴、沿角膜横向连续分布的气流，使角膜产生应变；同时由OCT系统对角膜成像，将成像系统获取的干涉光谱信息通过相位敏感算法来求取角膜内随深度变化的应变分布；然后可重建出角膜生物力学属性各参量的二维或三维分布图，实现了角膜的弹性成像。本发明还可对眼前节结构进行二维或三维成像。弹性信息对于角膜疾病的早期诊断、以及手术处理过程中的术前预估、手术设计、过程模拟、和术后评价均起着重要作用，它和眼前节结构信息的结合，能为眼科疾病的诊断和手术治疗提供更完整的依据。
Description
技术领域
本发明涉及角膜光学相干弹性成像和眼前节光学相干层析成像技术，尤其是涉及一种采用扫频光学相干层析成像技术、并结合相位敏感算法对角膜的生物力学信息和眼前节的结构信息进行在体测量的系统与方法。
背景技术
生命系统的力学属性是生物力学领域的重要研究内容之一，而力学属性的研究方法之一是基于弹性成像。弹性成像包括激励系统和应变探测成像系统两部分，激励系统对样品施加动态或静态的应力使之发生形变，样品形变区域内与深度相关的应变分布则由各种位移探测成像方法来获得，从而可得到样品形变区域内与深度相关的应力、应变、应力-应变、杨氏模量、泊松比等力学性能的分布。按照位移探测成像方法的不同，弹性成像的主要方法有：超声弹性成像、核磁共振弹性成像、和基于光学相干层析成像(OCT)技术的光学相干弹性成像(OCE)技术等。其中，超声弹性成像需借助耦合剂接触样品，核磁共振弹性成像的成本高昂、耗时、技术复杂，且它们的分辨率均只能达到0.1～1mm水平，根本不能满足临床上对高分辨率成像的要求，很难探测到微小或亚临床症状表现的组织弹性变化。而OCE技术运用OCT技术来探测形变，完全继承了OCT具有的非接触、无侵入、高分辨、和实时等特点，非常适合于组织的在体高分辨层析弹性成像。目前OCE的主要应用领域有：皮肤科中良性组织与肿瘤组织的区别、乳腺癌、角膜弹性成像等。
由于角膜的生物力学性能对角膜形状及正常生理功能的维持、以及在角膜疾病诊断、屈光手术设计、人工角膜研发、和角膜接触镜设计等方面具有的重要作用，使得利用弹性成像方法研究角膜的生物力学性能成为目前的热点。圆锥角膜、球形角膜、角膜变性等眼科疾病、以及屈光手术引起的角膜膨胀，都会导致角膜生物力学性能的改变，并进一步引起角膜形状的异常。角膜生物力学性能的异常表现较形状异常要早，属于亚临床特征，所以对角膜生物力学性能的监测成为临床上角膜疾病早期诊断的有效方法。它对角膜疾病手术治疗或屈光手术的作用则贯穿从术前预估、手术设计、过程模拟、术后受力分析、到效果评价的全过程。
基于OCT技术的角膜OCE技术是研究角膜生物力学性能的一种非常新颖的在体实时层析成像技术，主要由应力激励系统和应变探测成像系统(OCT系统)这两部分构成，激励源通常为声辐射力、热诱导、机械力、气流等。对眼科成像而言，气流激励为最简单和最安全的方式，如临床上广泛使用的非接触式眼压计及眼反应分析仪就属此种激励，因此是角膜OCE成像中激励源的最佳选择。
在DavidAlonso-Caneiro等人的文献(Assessmentofcornealdynamicswithhigh-speedsweptsourceopticalcoherencetomographycombinedwithanairpuffsystem.OpticsExpress,2011,19(15):14188-99)里，提出了基于扫频OCT技术的角膜OCE成像系统，采用气流对角膜施加激励，通过OCT成像系统来实时探测这一激励过程中角膜的位移(应变)量；通过对比角膜变形前后的结构图，来得到角膜某点处对应的整个轴向(A-scan)结构随时间的变化量。该系统存在着以下不足：1)角膜位移图基于OCT结构图，限制了位移分辨率(低于OCT的轴向分辨率，通常在μm量级)；2)为了探测到位移的变化，必须对角膜施加μm及以上量级的应变；3)到达平衡探测器两光路的光能不平衡，不能有效消除干涉光谱中直流项的干扰；4)样品臂中的喷气导管直径限制了系统的成像范围，不利于对角膜全局形变的观察；5)只对某一感兴趣的固定点进行观察。随后KarolKarnowski等人在文献(SweptsourceOCTwithairpuffchamberforcornealdynamicsmeasurement.Proc.ofSPIE,2012,8209:82090R)里对前述系统作了改进，弥补了3)和4)的不足。
在休斯顿大学JiasongLi等人的文献(Combiningafocusedair-puffsystemwithphase-sensitiveopticalcoherencetomographyforthedetectionofsoft-tissuetumorsbasedonelasticitymeasurement.Proc.ofSPIE,2013,8571:8571B-1)里，提出了采用扫频OCT技术及喷气激励系统的软组织弹性成像系统，在样品臂里利用喷气激励系统对软组织施加应力激励，产生表面弹性波，通过相位敏感型扫频OCT来获得表面弹性波的传播速度，从而得到软组织的杨氏模量分布。该系统的具体方法为：在离激励点2mm处每间隔0.5mm或1mm探测弹性波振幅及传播延迟时间，激励为聚焦的短脉冲气流。将激励点与测量点分离，样品臂中光束与激励气流方向不同轴，方便了激励系统的安装、和降低了对物镜焦距的限制。随后该组在文献(Dynamicopticalcoherencetomographymeasurementsofelasticwavepropagationintissue-mimickingphantomsandmousecorneainvivo.JournalofBiomedicalOptics,2013,18(12):121503)中，在样品臂及参考臂中加入了偏振控制器，减少了因偏振态失配导致的系统轴向分辨率下降。但在上述两系统中，也存在着以下不足：1)对实数干涉光谱信号进行反傅里叶变换产生的共轭镜像限制了系统的成像深度；2)聚焦的激励气流对软组织施加点激励时，需要在每个成像点探测随时间变化的激励信息，降低了系统的实时性；3)成像点为一些孤立的点，由这些点来重构表面弹性波的传播，因此获得的是离散弹性信息、而不是完整信息；4)需要保证激励与扫描的同步，且需要提供可重复控制的激励，增加了系统的复杂性。
此外，眼前节的结构信息，如：角膜厚度分布与曲率半径、前房深度与宽度、房角结构等，是目前临床上眼科疾病诊断的主要依据，而OCT技术是获取这些信息的有力工具。眼前节OCT成像仪已在临床上获得了广泛应用，成为了一种常规检查手段。
角膜的生物力学信息和眼前节的结构信息，在临床上表现出不同的形式，因此在疾病诊断和手术治疗上均起着重要作用，但各自均只表现了角膜的部分信息，只有把二者结合在一起，才能为眼科疾病的诊断和手术治疗提供更加多元和完善的信息。但到目前为止，还未见将此二者结合的任何报道。
发明内容
为了克服背景技术提到的不足，本发明的目的是提供一种角膜弹性成像与眼前节结构成像相结合的系统与方法。该系统采用喷气激励方法来对角膜施加生理范围内可承受的应力，激励为沿角膜横向连续分布的脉冲气流，气流方向与成像光束方向一致，作用区域与成像区域相一致；因激励导致的角膜应变，则由横向二维扫描的扫频OCT成像系统来探测，最后通过相位敏感算法来获得角膜的弹性信息。该系统还可对眼前节的结构进行OCT成像，最终获得包括角膜弹性信息和眼前节结构信息在内的完整信息。
本发明解决其技术问题所采用的技术方案是：一种角膜弹性成像与眼前节结构成像相结合的系统，其特征在于：由喷气激励系统和可横向二维扫描的扫频OCT成像系统组成；在喷气激励系统的作用下，角膜发生形变，同时通过扫频OCT成像系统对发生形变的角膜进行成像，将扫频OCT成像系统获取的干涉光谱信息通过相位敏感算法来求取角膜内随深度变化的应变分布，然后可重建出角膜的弹性分布图及眼前节结构图；
照明端中扫频光源提供的光束，被第一宽带光纤耦合器分成两束，分别进入样品臂和参考臂；调节参考臂中的平移台，使其带着参考镜移动，直至由样品臂中角膜或眼前节返回的光信号和由参考臂中参考镜返回的光信号之间形成干涉光谱信号；同时通过样品臂中的第一声光移频器和参考臂中的第二声光移频器引入的频率差来调制干涉光谱信号，被调制的干涉光谱信号经探测端中的第二宽带光纤耦合器分束后由平衡探测器接收；通过正交解调可获得复数干涉光谱信号，该复数干涉光谱信号随扫频光源的扫频时间而变化，对其进行反傅里叶变换可得到消除了镜像的眼前节结构信息，实现眼前节结构成像；另外，对该复干涉光谱信号通过相位敏感算法处理可得到角膜内与深度相关的弹性分布信息，实现角膜弹性成像；
扫频光源进行波长扫描的同时发出采样触发信号，去控制数据采集卡同步采集由平衡探测器接收到的干涉光谱信号；由函数发生卡的两个输出通道提供的扫描驱动信号与扫频光源发出的采样触发信号同步，分别控制垂直扫描器和水平扫描器进行扫描；数据采集卡采集到的信号传输至计算机进行处理。
进一步的，所述的喷气激励系统在面向扫描透镜的一侧为可使光束高效透过的玻璃窗口，面向角膜的一侧为气流出口；喷气激励系统提供沿角膜横向连续分布的气流，气流的方向与光束方向一致；角膜在气流的作用下发生变形。
进一步的，所述的可横向二维扫描的扫频OCT成像系统由照明端、样品臂、参考臂、探测端、函数发生卡、数据采集卡、和计算机组成。
进一步的，所述的扫频光源为输出波长随时间高速扫描的近红外波段宽光谱光源，同时也提供波数k空间的等间隔采样触发信号。
进一步的，所述的第一声光移频器和第二声光移频器，能在样品光及参考光之间引入频率差，正交解调后得到复数干涉光谱信号，对它进行反傅里叶变换来得到眼前节的结构图像，该图像不受共轭镜像、直流项、和自相干项的干扰。
进一步的，所述的第二宽带光纤耦合器为具有50:50分光比的2×2耦合器。
另外本发明还提供一种角膜弹性成像与眼前节结构成像相结合的方法，主要是将由可横向二维扫描的扫频OCT成像系统获得的复数干涉光谱信息传输至计算机，然后实施眼前节结构成像信息处理、和采用相位敏感算法的角膜弹性成像信息处理，具体分别如下：
对在每个横向扫描点采集到的复数干涉光谱信息，实施沿轴向的一维反傅里叶变换，即得到该扫描点处沿深度方向的一维OCT结构图；再结合横向一维或二维扫描，即可获得整个扫描区域内与深度相关的二维或三维眼前节结构图；
对每个横向扫描点采集到的复数干涉光谱信息，实施沿轴向的一维反傅里叶变换，由变换所得的实部和虚部经反余切运算，可获得包裹在-π～π之间的相位值，解包裹处理后得到所需的相位值，该相位值与角膜的形变量相关联；通过对比角膜变形前后的相位值，可求得角膜内的形变量δz，从而可计算出每个横向扫描点处沿深度方向的一维弹性分布；再结合横向一维或二维扫描，即可获得到整个扫描区域内与深度相关的二维或三维角膜弹性图。
本发明与现有技术相比的有益效果是：
(1)本发明能同时获取角膜的弹性信息和眼前节的结构信息，二者的结合能为眼科疾病的诊断和手术治疗提供更加多元和完善的信息。目前还未见同时针对此二者进行探测的系统与方法。
(2)本发明采用的激励方式为沿角膜横向连续分布的气流，气流的方向与成像光束方向一致，作用区域与成像区域相一致，使得弹性成像操作更简单、结果更准确，安全性也更高。而现有系统采用的点激励方式，需探测点激励传输至被成像点处的激励值，这一过程较为复杂；或者采用斜激励方式(激励方向与光束方向不一致)，则会同时引起沿角膜纵向(成像方向)与横向的应变，且此二者相互耦合，使得沿任一方向的计算结果均不准确。
(3)本发明成像的范围可覆盖整个角膜，也可只针对某一特定区域进行成像，而现有系统只针对某些点进行一维成像。由于OCT的扫频速率极高(101-102KHz量级)，成像光束可进行横向一维或二维扫描，从而可实现角膜弹性和眼前节结构的在体实时二维(帧频大于25fps)或三维成像。
(4)本发明采用相位敏感型算法，可得到nm量级的形变探测灵敏度。扫频OCT探测到某点的复干涉频谱为，其中：Rr为参考镜反射率、Rs(n)为对应采样点n处的样品反射率、ωn为与光程差nΔz相关的角频率、为与形变引起的光程差变化nΔz+δz相关的相位、n为轴向采样点数、Δz为OCT的轴向分辨率。I(t)经反傅里叶变换后为，得到的OCT结构信息为Rs(n)，其分辨率只有μm量级。而利用相位敏感型算法提取的相位值与形变引起的光程差变化nΔz+δz相关，由不同形变状态之间的相位差得到对应的形变量δz，其形变探测灵敏度由样品中反射信号的信噪比决定，即，所以其形变探测灵敏度可达nm量级。
附图说明
图1是本发明的系统结构示意图；
图2是本发明的控制系统示意图；
图3是本发明的信号处理方法及流程图；
图4是本发明的水平扫描光路及喷气激励系统结构图；
图5是本发明的扫描成像原理图。
图中：1.照明端、2.样品臂、3.参考臂、4、探测端、5.扫频光源、6.第一宽带光纤耦合器、7.第一声光移频器、8.第一光环形器、9.第一偏振控制器、10.第一准直透镜、11.垂直扫描器、12.水平扫描器、13.扫描透镜、14.喷气激励系统、15.角膜、16.眼前节、17.第二声光移频器、18.第二光环形器、19.第二偏振控制器、20.第二准直透镜、21.中性密度滤光片、22.色散补偿片、23.参考镜、24.平移台、25.第二宽带光纤耦合器、26.平衡探测器、27-37.第一至第十一单模光纤、38.函数发生卡、39.数据采集卡、40.计算机。
具体实施方式
下面结合附图以及具体实施例进一步说明本发明。
本发明的系统结构如图1、图2示，包括照明端1、样品臂2、参考臂3、探测端4、喷气激励系统14、函数发生卡38、数据采集卡39、和计算机40。
图1中照明端1由扫频光源5、第一单模光纤27、第一宽带光纤耦合器6组成。样品臂2由第一声光移频器7、第一光环形器8、第一偏振控制器9、第一准直透镜10、垂直和水平扫描器11和12、扫描透镜13、第二至第四单模光纤28-30、第八单模光纤34组成。参考臂3由第二声光移频器17、第二光环形器18、第二偏振控制器19、第二准直透镜20、中性密度滤光片21、色散补偿片22、参考镜23、平移台24、第五至第七单模光纤31-33、和第九单模光纤35组成。探测端4由第二宽带光纤耦合器25、平衡探测器26、第十和第十一单模光纤36和37组成。
由照明端1的扫频光源5发出的光信号经第一单模光纤27传输至第一宽带光纤耦合器6后分成两路：一路进入样品臂2、另一路进入参考臂3。扫频光源5输出波长随时间高速扫描的近红外波段宽光谱光源，同时也提供波数k空间的等间隔采样触发信号。
在样品臂2里，由第二单模光纤28传输至第一声光移频器7，经第三单模光纤29传输至第一光环形器8的端口a，再从第一光环形器8的端口b出射，出射光束经第四单模光纤30传输和第一准直透镜10准直后，依次通过垂直扫描器11和水平扫描器12后、被扫描透镜13聚焦在角膜15或眼前节16上。从角膜15或眼前节16返回的样品光，沿原路返回至第一光环形器8的端口b，再从第一光环形器8的端口c出射，并由第八单模光纤34传输至第二宽带光纤耦合器25。喷气激励系统14位于扫描透镜13和角膜15之间，在面对扫描透镜13的一侧为可使光束高效透过的玻璃窗口，面对角膜15的一侧为气流出口。由喷气激励系统14提供的激励气流作用在角膜15上，使之发生变形，激励方式为沿角膜横向连续分布的气流，激励气流的方向与成像光束方向一致。
在参考臂3里，由第五单模光纤31传输至第二声光移频器17，经第六单模光纤32传输至第二光环形器18的端口d，再从第二光环形器18的端口e出射，出射光束经第七单模光纤33传输和第二准直透镜20准直后，通过中性密度滤光片21和色散补偿片22后，垂直入射固定在平移台24上的参考镜23。从参考镜23返回的参考光，沿原路返回至第二光环形器18的端口e，再从第二光环形器18的端口f出射，并由第九单模光纤35传输至第二宽带光纤耦合器25。第一、二声光移频器分别在样品臂及参考臂中引入频率差调制干涉信号；中性密度滤光片21衰减参考臂中的光强；色散补偿片22用来匹配因样品臂中扫描透镜的引入造成的参考臂和样品臂之间的色散失配。
在探测端4，从第八单模光纤34传输的样品光和第九单模光纤35传输的参考光，通过第二宽带光纤耦合器25后，各自分成两部分，并分别由第十单模光纤36和第十一单模光纤37传输至平衡探测器26的负极和正极接收端。第二宽带光纤耦合器25为具有50:50分光比的2×2耦合器。
平移台24带着参考镜23作直线移动，直至由眼角膜15返回的样品光和从参考镜23返回的参考光形成干涉条纹。
照明端1、样品臂2、参考臂3、探测端4构成可横向二维扫描的扫频OCT成像系统，由可横向二维扫描的扫频OCT成像系统测量喷气激励系统14引起的角膜15内部形变分布。
本发明的控制系统如图2所示。扫频光源5进行波长扫描的同时发出采样触发信号，去控制数据采集卡39同步采集由平衡探测器26接收到的干涉光谱信号；由函数发生卡38的两个输出通道提供的扫描驱动信号与扫频光源5发出的采样触发信号同步，分别控制垂直扫描器11和水平扫描器12进行扫描；数据采集卡39采集到的信号传输至计算机40进行处理。
图3是本发明的信号处理方法及流程图。一种角膜弹性成像与眼前节结构成像相结合的方法，主要是将由可横向二维扫描的扫频OCT成像系统获得的复数干涉光谱信息传输至计算机40，然后实施眼前节结构成像信息处理、和采用相位敏感算法的角膜弹性成像信息处理，具体包括以下步骤：
步骤S1：系统参数计算与设置：计算第一和第二声光移频器7和17需在样品臂2及参考臂3之间引入的移频差，确定横向扫描范围、扫描间距及扫描方式等；
步骤S2：采集横向B-scan的干涉光谱信号：扫频光源5开始波长扫描，同时发出采样触发信号；函数发生卡38的一个输出通道提供锯齿波信号同步驱动水平扫描器12进行横向B-scan扫描，数据采集卡39同步采集横向B-scan的干涉光谱信号，它由一系列轴向、即A-scan信号构成；
步骤S3：接收端信号解调：平衡探测器26接收到的干涉光谱信号，经带通滤波器后，利用正交解调器得到干涉信号的实部及虚部，可重构出每一A-scan的复干涉光谱信号；
步骤S4：信号处理与二维眼前节结构图像重建：对复干涉光谱信号沿轴向进行一维反傅里叶变换，可得到消除镜像的A-scan结构图，通过水平扫描由B-scan包含的所有A-scan结构图可重建出眼前节16的二维OCT结构图；
步骤S5：采用相位敏感算法重建角膜二维弹性图：对复干涉光谱信号沿轴向进行一维反傅里叶变换，再由变换结果的实部和虚部经反余切处理，可获得包裹在-π～π之间的相位值，解包裹处理后得到所需的相位值，该相位值与角膜15的形变量相关联；通过对比角膜15变形前后的相位值，可求得角膜15的形变量δz，从而可计算出A-scan对应的弹性模量、形变位移等弹性信息，通过水平扫描由B-scan包含的所有A-scan弹性信息，可重构出角膜15的二维弹性图；
步骤S6：横向二维干涉光谱信号采集，及眼前节三维结构图和角膜三维弹性图的重建：在扫频光源5采样触发信号的触发下，函数发生卡38的另一输出通道提供的阶梯信号同步驱动垂直扫描器11进行扫描，重复步骤S2至步骤S5，可获得沿垂直扫描方向的多副B-scan二维图像，由它们可分别重建出眼前节16的三维结构图像和角膜15的三维弹性图像。
图4是本发明的水平扫描光路及喷气激励系统结构图。扫频光束由水平扫描器12偏转后透过扫描透镜13，并经喷气激励系统14的玻璃窗口14-1及喷气嘴14-2聚焦在位于焦平面处样品上的不同位置；扫描驱动信号控制水平扫描器12振动到b、a、b’位置，相应的扫频光束则聚焦到焦平面处样品面上的点1b、1a、1b’。
图5是本发明所用的扫描成像原理图。扫频光源发出的扫频光信号分成两束分别进入样品臂及参考臂，从样品臂及参考臂返回的光信号形成干涉光谱信号并被平衡探测器探测，该信号经正交解调后得到复数干涉光谱信号；扫频光源进行波长扫描的同时发出采样触发信号，去控制数据采集卡同步采集解调后的复数干涉光谱信号；由函数发生卡的两个输出通道提供的扫描驱动信号与扫频光源发出的采样触发信号同步，分别控制二维扫描系统进行扫描；数据采集卡采集到的信号传输至计算机进行眼前节结构成像和角膜弹性成像处理，完成成像的重构及显示。
本发明未详细阐述部分属于本领域技术人员的公知技术。
上述具体实施方式用来解释说明本发明，而不是对本发明进行限制。在本发明的精神和权利要求的保护范围内，对本发明作出的任何修改和改变，都落入本发明的保护范围。
Claims (7)
1.一种角膜弹性成像与眼前节结构成像相结合的系统，其特征在于：由喷气激励系统(14)和可横向二维扫描的扫频OCT成像系统组成；在喷气激励系统(14)的作用下，角膜(15)发生形变，同时通过扫频OCT成像系统对发生形变的角膜(15)进行成像，将扫频OCT成像系统获取的干涉光谱信息通过相位敏感算法来求取角膜内随深度变化的应变分布，然后可重建出角膜的弹性分布图及眼前节结构图；
照明端(1)中扫频光源(5)提供的光束，被第一宽带光纤耦合器(6)分成两束，分别进入样品臂(2)和参考臂(3)；调节参考臂(3)中的平移台(24)，使其带着参考镜(23)移动，直至由样品臂(2)中角膜(15)或眼前节(16)返回的光信号和由参考臂(3)中参考镜(23)返回的光信号之间形成干涉光谱信号；同时通过样品臂(2)中的第一声光移频器(7)和参考臂(3)中的第二声光移频器(17)引入的频率差来调制干涉光谱信号，被调制的干涉光谱信号经探测端(4)中的第二宽带光纤耦合器(25)分束后由平衡探测器(26)接收；通过正交解调可获得复数干涉光谱信号，该复数干涉光谱信号随扫频光源(5)的扫频时间而变化，对其进行反傅里叶变换可得到消除了镜像的眼前节结构信息，实现眼前节结构成像；另外，对该复数干涉光谱信号通过相位敏感算法处理可得到角膜(15)内与深度相关的弹性分布信息，实现角膜弹性成像；
扫频光源(5)进行波长扫描的同时发出采样触发信号，去控制数据采集卡(39)同步采集由平衡探测器(26)接收到的干涉光谱信号；由函数发生卡(38)的两个输出通道提供的扫描驱动信号与扫频光源(5)发出的采样触发信号同步，分别控制垂直扫描器(11)和水平扫描器(12)进行扫描；数据采集卡(39)采集到的信号传输至计算机(40)进行处理。
2.根据权利要求1所述的角膜弹性成像与眼前节结构成像相结合的系统，其特征在于：所述的喷气激励系统(14)在面向扫描透镜(13)的一侧为可使光束透过的玻璃窗口，面向角膜(15)的一侧为气流出口；喷气激励系统(14)提供沿角膜(15)横向连续分布的气流，气流的方向与光束方向一致；角膜(15)在气流的作用下发生变形。
3.根据权利要求1所述的角膜弹性成像与眼前节结构成像相结合的系统，其特征在于：所述的可横向二维扫描的扫频OCT成像系统由照明端(1)、样品臂(2)、参考臂(3)、探测端(4)、函数发生卡(38)、数据采集卡(39)和计算机(40)组成。
4.根据权利要求1所述的角膜弹性成像与眼前节结构成像相结合的系统，其特征在于：所述的扫频光源(5)为输出波长随时间高速扫描的近红外波段宽光谱光源，同时也提供波数k空间的等间隔采样触发信号。
5.根据权利要求1所述的角膜弹性成像与眼前节结构成像相结合的系统，其特征在于：所述的第一声光移频器(7)和第二声光移频器(17)，能在样品光及参考光之间引入频率差，正交解调后得到复数干涉光谱信号，对它进行反傅里叶变换来得到眼前节(16)的结构图像，该图像不受共轭镜像、直流项和自相干项的干扰。
6.根据权利要求1所述的角膜弹性成像与眼前节结构成像相结合的系统，其特征在于：所述的第二宽带光纤耦合器(25)为具有50:50分光比的2×2耦合器。
7.一种角膜弹性成像与眼前节结构成像相结合的方法，主要是将由可横向二维扫描的扫频OCT成像系统获得的复数干涉光谱信息传输至计算机(40)，然后实施眼前节结构成像信息处理和采用相位敏感算法的角膜弹性成像信息处理，具体分别如下：
对在每个横向扫描点采集到的复数干涉光谱信息，实施沿轴向的一维反傅里叶变换，即得到该扫描点处沿深度方向的一维OCT结构图；再结合横向一维或二维扫描，即可获得整个扫描区域内与深度相关的二维或三维眼前节结构图；
对每个横向扫描点采集到的复数干涉光谱信息，实施沿轴向的一维反傅里叶变换，由变换所得的实部和虚部经反余切运算，可获得包裹在-π～π之间的相位值，解包裹处理后得到所需的相位值，该相位值与角膜(15)的形变量相关联；通过对比角膜(15)变形前后的相位值，可求得角膜(15)内的形变量δz，从而可计算出每个横向扫描点处沿深度方向的一维弹性分布；再结合横向一维或二维扫描，即可获得到整个扫描区域内与深度相关的二维或三维角膜弹性图。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410654412.7A CN104434028B (zh) | 2014-11-15 | 2014-11-15 | 角膜弹性成像与眼前节结构成像相结合的系统与方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410654412.7A CN104434028B (zh) | 2014-11-15 | 2014-11-15 | 角膜弹性成像与眼前节结构成像相结合的系统与方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104434028A CN104434028A (zh) | 2015-03-25 |
| CN104434028B true CN104434028B (zh) | 2016-04-27 |
Family
ID=52881344
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410654412.7A Active CN104434028B (zh) | 2014-11-15 | 2014-11-15 | 角膜弹性成像与眼前节结构成像相结合的系统与方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104434028B (zh) |
Families Citing this family (26)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105167805A (zh) * | 2015-08-19 | 2015-12-23 | 深圳市亿领科技有限公司 | 一种测量角膜弹性的方法及装置 |
| CN105231990B (zh) * | 2015-11-17 | 2017-04-19 | 深圳市亿领科技有限公司 | 基于oct三维成像分析角膜生物力学性能的装置及方法 |
| CN105411530A (zh) * | 2015-12-17 | 2016-03-23 | 天津求实飞博科技有限公司 | 新型光学相干层析装置 |
| CN105615824A (zh) * | 2016-03-25 | 2016-06-01 | 中国科学院光电技术研究所 | 采用贝塞尔环带照明方式的长焦深暗场视网膜oct系统 |
| US11413007B2 (en) | 2016-05-26 | 2022-08-16 | University Of Washington | Non-contact acoustic radiation force based (ARF-based) generation of broad bandwidth mechanical waves using air-coupled ultrasound |
| CN105891434B (zh) * | 2016-06-03 | 2017-09-26 | 深圳职业技术学院 | 一种海水盐度的在线检测方法及其装置 |
| CN108572142B (zh) * | 2017-03-09 | 2020-09-22 | 清华大学 | 一种高信噪比的光计算光学相干成像系统 |
| CN106963338B (zh) * | 2017-03-31 | 2018-09-07 | 杭州电子科技大学 | 晶状体生物力学和光学特性无创在体成像系统及测量方法 |
| CN111225597A (zh) * | 2017-08-18 | 2020-06-02 | 通用医疗公司 | 用于布里渊光谱学和组织的成像的系统和方法 |
| CN108733913B (zh) * | 2018-05-17 | 2022-03-29 | 北京化工大学 | 一种基于dwpso算法的眼科oct设备横向分辨率检测方法 |
| CN109171639B (zh) * | 2018-09-04 | 2023-10-20 | 温州医科大学 | 一种基于光学相干层析成像技术的在体角膜参数的测量装置及测量方法 |
| CN109793496A (zh) * | 2018-12-04 | 2019-05-24 | 苏州大学 | 一种结合传感器和ss-oct技术测试眼压的装置及应用 |
| CN109620131B (zh) * | 2018-12-14 | 2021-08-03 | 佛山科学技术学院 | 共光路微透镜阵列多光束光学相干弹性测量系统及方法 |
| CN109645954B (zh) * | 2018-12-14 | 2021-06-29 | 佛山科学技术学院 | 基于微透镜阵列的多光束光学相干的弹性测量系统及方法 |
| CN109620132B (zh) * | 2018-12-14 | 2021-08-03 | 佛山科学技术学院 | 一种多探测光束光学相干在体角膜弹性测量系统及方法 |
| CN109674441B (zh) * | 2018-12-14 | 2021-06-29 | 佛山科学技术学院 | 基于光开关和微透镜阵列的多光束弹性测量系统及方法 |
| CN109691973B (zh) * | 2018-12-29 | 2023-06-13 | 佛山科学技术学院 | 一种用于测量眼球脉动的光学相干层析成像系统 |
| CN109875504B (zh) * | 2019-01-15 | 2021-07-30 | 温州医科大学 | 一种基于喷气式光学相干弹性成像技术无创测量角膜粘弹性的方法 |
| CN109965839A (zh) * | 2019-04-08 | 2019-07-05 | 广东唯仁医疗科技有限公司 | 眼前节oct成像装置 |
| CN109965841B (zh) * | 2019-04-23 | 2024-07-02 | 广东唯仁医疗科技有限公司 | 一种人工晶体植入术的弹性分析装置及方法 |
| CN110974148B (zh) * | 2019-12-03 | 2022-01-07 | 温州医科大学 | 一种基于喷气式光学相干弹性成像技术测量在体人眼角膜弹性模量的方法 |
| CN112826440B (zh) * | 2021-01-05 | 2022-12-20 | 佛山科学技术学院 | 基于光学相干法的角膜结构-眼压-多力学参数测量仪 |
| CN112587086A (zh) * | 2021-03-04 | 2021-04-02 | 季华实验室 | 一种双模式偏振光学相干成像系统及其成像方法 |
| CN115682988B (zh) * | 2022-12-13 | 2023-05-09 | 浙江大学 | 人体组织边缘qme检测系统 |
| CN116026682B (zh) * | 2023-03-30 | 2023-07-04 | 浙江大学 | 基于qme的快速弹性成像计算方法 |
| CN116725475B (zh) * | 2023-05-15 | 2024-03-12 | 北京大学 | 眼角膜生物力学性质的测量系统及方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202211675U (zh) * | 2011-07-13 | 2012-05-09 | 苏州新视野光电技术有限公司 | 双通道全眼光学相干层析成像系统 |
| CN202568206U (zh) * | 2012-02-14 | 2012-12-05 | 苏州微清医疗器械有限公司 | 一种视网膜三维成像装置 |
| CN202619628U (zh) * | 2012-05-07 | 2012-12-26 | 浙江大学 | 基于视见函数调控的全眼前节谱域oct成像系统 |
| EP2769666A1 (en) * | 2011-10-20 | 2014-08-27 | Asociación Industrial de Óptica, Color e Imagen - AIDO | Multi-view fundus camera |
| CN104013383A (zh) * | 2014-06-04 | 2014-09-03 | 南京航空航天大学 | 双焦点眼前后节同步成像系统及成像方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007195726A (ja) * | 2006-01-26 | 2007-08-09 | Mitaka Koki Co Ltd | 眼底検査装置 |
| JP5679630B2 (ja) * | 2008-12-26 | 2015-03-04 | キヤノン株式会社 | 光断層撮像装置及びその方法 |
-
2014
- 2014-11-15 CN CN201410654412.7A patent/CN104434028B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202211675U (zh) * | 2011-07-13 | 2012-05-09 | 苏州新视野光电技术有限公司 | 双通道全眼光学相干层析成像系统 |
| EP2769666A1 (en) * | 2011-10-20 | 2014-08-27 | Asociación Industrial de Óptica, Color e Imagen - AIDO | Multi-view fundus camera |
| CN202568206U (zh) * | 2012-02-14 | 2012-12-05 | 苏州微清医疗器械有限公司 | 一种视网膜三维成像装置 |
| CN202619628U (zh) * | 2012-05-07 | 2012-12-26 | 浙江大学 | 基于视见函数调控的全眼前节谱域oct成像系统 |
| CN104013383A (zh) * | 2014-06-04 | 2014-09-03 | 南京航空航天大学 | 双焦点眼前后节同步成像系统及成像方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104434028A (zh) | 2015-03-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104434028B (zh) | 角膜弹性成像与眼前节结构成像相结合的系统与方法 | |
| Ambroziński et al. | Acoustic micro-tapping for non-contact 4D imaging of tissue elasticity | |
| CN106963338B (zh) | 晶状体生物力学和光学特性无创在体成像系统及测量方法 | |
| CN104523239B (zh) | 全深度谱域光学相干层析成像装置及方法 | |
| CN101869466B (zh) | 基于自适应光学技术的共焦扫描与光学相干层析成像仪 | |
| US11413007B2 (en) | Non-contact acoustic radiation force based (ARF-based) generation of broad bandwidth mechanical waves using air-coupled ultrasound | |
| CN104958061A (zh) | 双目立体视觉三维成像的眼底oct成像方法及其系统 | |
| EP2683289B1 (en) | Multi-channel optical coherence tomography for imaging and temperature and force sensing | |
| JP6014026B2 (ja) | 強化された眼科測定方法および装置 | |
| CN102944521B (zh) | 非接触式光声和光学相干断层双成像装置及其检测方法 | |
| CN100488440C (zh) | 共路型内窥光学相干层析成像方法及系统 | |
| CN101884524B (zh) | 基于自适应光学技术的宽视场光学相干层析仪 | |
| EP2683288B1 (en) | Low cost high efficiency signal interrogation for multi-channel optical coherence tomography | |
| JP5679630B2 (ja) | 光断層撮像装置及びその方法 | |
| CN109645954B (zh) | 基于微透镜阵列的多光束光学相干的弹性测量系统及方法 | |
| CN102578993A (zh) | 一种内窥oct成像装置及成像方法 | |
| CN109171639A (zh) | 一种基于光学相干层析成像技术的在体角膜参数的测量装置及测量方法 | |
| CN111449629B (zh) | 一种光学相干弹性成像方法及装置 | |
| KR101011556B1 (ko) | 광학적 도플러 단층 촬영 장치 및 그의 도플러 주파수 검출방법 | |
| CN109620130A (zh) | 一种共光路的多光束光学相干弹性测量系统及测量方法 | |
| CN111466884B (zh) | 一种光学相干弹性成像装置 | |
| CN202681903U (zh) | 一种内窥oct成像装置 | |
| CN105167805A (zh) | 一种测量角膜弹性的方法及装置 | |
| CN106680214A (zh) | 一种同时检测弹性和粘性的光声腹腔镜装置及其方法 | |
| CN111436910A (zh) | 一种活体组织的光学相干层析多模态成像装置及方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |