2. 上海理工大学 光电信息与计算机工程学院,上海 200093
2. School of Optical-Electrical and Computer Engineering, University of Shanghai for Science and Technology, Shanghai 200093,China
临床上,常用的眼轴长度测量方法主要有超声测量法和光学测量法[1-5]。虽然超声生物测量法被广泛使用,但仍有许多限制,如需要麻醉和接触式测量,会使受测者增加交叉感染的风险,同时对操作者有严格的要求,并且测量结果小于实际的眼轴长度[6-8]。光学生物测量所测得的眼轴长度为沿视轴方向从泪膜到视网膜色素上皮层的距离,相比于超声测量结果更接近真正意义上的视轴[3]。光学测量方法主要是基于光学相干技术的测量,目前应用较广的人眼测量光学仪器有IOL Master(德国,Carl Zeiss公司)、OA-1000(日本,Tomey)及Lenstar LS 900(瑞士,HAAG-ATREIT公司),它们分别基于部分相干干涉原理(PCI)和低相干反射测量原理(OLCR)。光学测量无需接触眼睛,安全方便,速度快且具有更高的分辨率和精度(可精确到10 μm)。
本文针对眼轴长度测量,设计了一套眼轴长度测量系统。该测量系统包括硬件设计和软件设计两部分。硬件设计中采用PIN光电二极管作为光电探测器,并设计一套信号处理电路,检测从眼球前后表面反射回来的两组干涉信号,并由磁栅尺位移传感器记录下两组干涉信号之间参考臂运动平台移动的距离,即眼球前后表面的光学距离。干涉信号包络信息和位置信息传给软件部分进行处理,然后根据眼球屈光介质的平均折射率可以计算得到眼轴长度。
1 测量系统硬件设计 1.1 硬件整体设计从眼睛各界面的反射光来分析,最外层界面(角膜前表面)的反射光的光强大于眼内其他界面的反射光,而眼轴长度测量中的一个重要干涉信号来自视网膜色素上皮层即眼球后表面,该表面的反射光相对较弱(有研究表明人眼角膜前表面反射的杂散光是视网膜有效反射率的10~100倍[9])。本系统采用光纤进行空间条纹提取,光纤芯径为62.8 μm,接收干涉图样中的某一级干涉条纹。这些因素都直接影响眼球前表面和后表面信号的强度,因此需要设计一个低噪声高性能的弱信号处理电路。
结合弱信号特点和测量系统要求,本系统的硬件设计框图如图1所示。
|
图 1 系统硬件设计框图 Figure 1 Block diagram of the hardware design |
该电路将PIN光电二极管探测到的弱光信号进行选频放大至百毫伏量级,便于控制系统量化采样。PIN光电二极管将光信号转换成便于处理的电流信号,经过前置放大电路将电流信号转化成电压信号并进行放大,以便于观察直流本底信号,协助光路调节。高通滤波器负责滤除直流本底信号以及有效干涉信号频率以外的信号。由于眼球前后表面返回的光强大小具有明显的差别,常规的线性放大电路极易使得前表面的信号满偏,干涉信号会失去其明显的幅值特征,影响最终的测量值。因此采用增益自动控制的非线性放大电路,即前表面的信号强,采用小放大倍数,后表面的信号弱,电路自动切换成高增益模式。带通滤波器用于提取干涉信号,干涉信号再通过同相比例放大电路输出至模数转换器。
1.2 前置放大电路设计弱信号处理电路的灵敏度依赖于光电探测器的选择,高灵敏度的光电探测器可以探测到更加精确的弱信号[10]。本系统采用的光电探测器是北京敏光科技生产的LSSPD-0.5的PIN光电二极管,该光电二极管具灵敏度高和暗电流小的特点,响应率为0.45 mA/mW,可用于400~1 100 nm波长的微弱光信号检测[11]。该PIN光电二极管可以将nW级的弱光信号转化为nA级的电流,但需要前置放大电路将该电流信号转化成电压信号并进行放大[12-13],以利于后续电路中信号的获取,因此该前置放大电路的设计尤为重要。
集成运算放大器AD8066具有低噪声和精密失调特性,特别适合作光电二极管的电流/电压(I/V)转换放大器。AD8626是一种JFET型输入放大器,其单位增益稳定,输入电压噪声为
|
图 2 前置放大电路图 Figure 2 Preamplifier circuit |
该眼轴长度测量系统中外差干涉信号作为载频进入电路,当直流电机与导轨相配合控制参考臂以45 mm/s的速度匀速移动时,对应的特征频率为115 kHz。利用ADI滤波器设计向导设计一个中心频率为115 kHz、通带为30 kHz的八阶巴特沃兹带通滤波电路,如图3所示。通过该部分电路,可以提取干涉信号并滤除带宽范围外的噪声。
|
图 3 带通滤波电路图 Figure 3 Bandpass filter circuit |
本系统选择 Max 10系列中的 10 M25SAE 作为主控芯片,其内部集成了模数转换器(ADC)、片上闪存(FLASH)、可编程逻辑阵列(LABs)及基本输入输出等。其中,内置的模数转换器不仅可以满足本系统的工程需求,而且能够减小印刷电路板(PCB)的设计面积,降低电路设计的难度。
弱信号处理电路的输出信号在传送至模数转换器的输入引脚时,如果输入信号线走线不合理,输入信号之间的串扰相对较大,进而影响测量精度。因此在PCB设计时要遵守两个设计原则:(1)电源线和地线尽可能粗短,输入信号线不平行于电源线、地线及数字IO 口走线;(2)外接电源串联磁珠接至电源引脚,并联一个10 μF电容接地,组成一个RC低通滤波器,以抑制电源线上的高频噪声,同时电源引脚使用100 nF去耦电容。
2 测量系统软件设计 2.1 软件整体设计测量系统的软件设计包括初始化、系统自检、干涉信号处理和串口通信等部分。系统上电后,初始化相关功能模块,根据串口接收的指令执行系统自检,自检通过则等待串口传送的上位机其他指令(如眼轴长度测量指令),内置模数转换器采集干涉信号并用峰值检波算法提取干涉信号的峰值,最后将测量数据发送给上位机。
2.2 干涉信号峰值检波程序设计眼轴长度测量依赖于对干涉信号的有效识别,干涉信号呈现高斯分布形状,记录其峰值信息。该干涉信号是由光干涉形成的载波调制信号,需要经过解调电路才能提取包络信号的峰值。常用的模拟型解调电路是由二极管、电容与电阻共同构成的RC充放电式的检波电路,由于二极管的存在使得该电路具有门限电压,降低了信号处理模块的动态输入范围,另外模拟电路容易受到外界干扰,也增加了后续处理的难度。Max 10具备丰富的可编程逻辑和内置的高采样率模数转换器,可直接使用FPGA实现数字电路解调,并将其包络信号存放至缓冲区,通过峰值检测算法计算出该波形的峰值,可以提高系统的抗干扰能力。
干涉信号的数字解调模块使用VHDL语言(一种用于电路设计的高级语言)进行编写,模数转换器输出的数字化后的干涉信号先去除直流分量,再将非正信号取反,实现波形取绝对值的功能,再通过数字低通滤波器实现解调功能,并将数据存放在缓冲区以供软核处理器使用。
在实验中采用示波器显示测量得到的干涉信号,展开后的信号宽度约1 ms,故配置20 μs定时器中断,并构造一维数组(包含50个元素)用来存放转换出的数字。通过窗口滑动比较数组内50个数的数值,找到波形的峰值并保存,同时记录对应位置磁栅尺的读数。峰值检测算法的流程如图4所示。
|
图 4 窗口滑动峰值检测算法 Figure 4 The flow chart of the peak detection algorithm |
为了方便快速识别有效数据包,提高数据传输效率,系统自定义一套用于PC机与下位机数据传输的通信协议,包括帧头、帧尾、数据长度、有效数据及求和校验等。帧头为两个字节的数组uart_buff[0]和uart_buff[1],当串口缓存区的数据uart_buff[0]=0×55且uart_buff[1]=0×AA时,认为该数据包为有效数据,第三个字节表明该数据包的类型,下位机根据命令进行相应的动作。
3 实验结果为验证本文设计的测量系统的有效性,采用蔡司公司的模拟眼作为被测样品,该标准模拟眼的轴长为20.80 mm,其折射率约为1.384 6,通过MSO400系列示波器直接显示各级输出波形。参考人眼激光安全规定,设定激光器输出电压为5 V,采用索雷博公司的光功率计(型号为PM100D)测得此时790 nm波长的红外光进入模拟眼的光功率为710 μW,测得经眼球内部反射回来的光功率为80 nW。通过系统中的弱信号处理电路,得到眼球前表面和后表面干涉信号分别为1 880 mV和900 mV,如图5所示。
|
图 5 模拟眼的弱信号处理电路结果图 Figure 5 Weak signal processing result of the simulation eye |
通过精度为5 μm的磁栅尺位移传感器得到两个峰值对应的位置信息,分别记为D1和D2,则被测器件的光学长度L表达式为:L=D2−D1。根据光学长度和折射率,可以得到被测物的物理长度。任意选取10组模拟眼眼轴长度测量结果进行分析,如表1所示。
|
|
表 1 模拟眼眼轴长度测量数据 Table 1 The measurement data of eye axis length |
测量结果显示,该系统测量误差为±0.01 mm。因此,该测量系统可以稳定地探测到模拟眼的眼轴前表面及后表面信号,并通过与其他机械结构配合得到模拟眼的眼轴长度。
4 结 论本文基于外差干涉原理,设计了一套用于测量眼轴长度的测量系统。硬件系统设计中选择带宽为30 kHz、中心频率为115 kHz的八阶带通滤波电路,该电路可以有效提高系统的信噪比。另外,系统选择数字型峰值检波电路替代传统的模拟型检波电路,提高了系统的抗干扰能力。实验结果表明,该测量系统较好地实现了光电信号的转换和信号提取,能有效且稳定地完成眼轴长度的测量。
| [1] | TAKEI K, SEKINE Y, OKAMOTO F, et al. Measurement of axial length of eyes with incomplete filling of silicone oil in the vitreous cavity using x ray computed tomography[J]. British Journal of Ophthalmology, 2002, 86(1): 47–50. DOI:10.1136/bjo.86.1.47 |
| [2] | 吴钦星. 磁共振在硅油眼眼轴测量中的应用[J]. 国际眼科杂志, 2010, 10(5): 1019–1020. DOI:10.3969/j.issn.1672-5123.2010.05.078 |
| [3] | 罗廷浩. IOL Master和A超眼轴测量比较研究[J]. 国际眼科杂志, 2012, 12(6): 1213. |
| [4] | TRIVEDI R H, WILSON M E. Prediction error after pediatric cataract surgery with intraocular lens implantation: contact versus immersion A-scan biometry[J]. Journal of Cataract & Refractive Surgery, 2011, 37(3): 501–505. |
| [5] | 杜敏晖. B型超声测量眼轴的临床研究[D]. 青岛: 青岛大学, 2008. |
| [6] | SUGIHARA K, TANITO M, TAKAI Y, et al. Causes of intracapsular cataract extraction, explantation of intraocular lenses and suture scleral fixation of intraocular lenses in the modern era[J]. Acta Ophthalmologica, 2018, 96(2): e262–e263. DOI:10.1111/aos.2018.96.issue-2 |
| [7] | REITBLAT O, LEVY A, KLEINMANN G, et al. Accuracy of intraocular lens power calculation using three optical biometry measurement devices: the OA-2000, Lenstar-LS900 and IOLMaster-500[J]. Eye, 2018, 32(7): 1244–1252. DOI:10.1038/s41433-018-0063-x |
| [8] | TRIVEDI R H, WILSON M E. Axial length measurements by contact and immersion techniques in pediatric eyes with cataract[J]. Ophthalmology, 2011, 118(3): 498–502. DOI:10.1016/j.ophtha.2010.06.042 |
| [9] | 程少园, 曹召良, 胡立发, 等. 消除角膜前表面反射杂散光方法的比较[J]. 中国光学, 2010, 3(3): 257–262. DOI:10.3969/j.issn.2095-1531.2010.03.008 |
| [10] | ENG P C, SONG S, PING B. State-of-the-art photodetectors for optoelectronic integration at telecommunication wavelength[J]. Nanophotonics, 2015, 4(3): 277–302. |
| [11] | TANZID M, MOHAMMEDY F M. Analysis of electrical and optical characteristics of InGaSb PIN photodiodes[J]. IET Optoelectronics, 2012, 6(5): 242–249. DOI:10.1049/iet-opt.2011.0101 |
| [12] | LU G N, PITTET P, SOU G, et al. A novel approach to implementing on-chip synchronous detection for CMOS optical detector systems[J]. Analog Integrated Circuits and Signal Processing, 2003, 37(1): 57–66. DOI:10.1023/A:1024828902878 |
| [13] | TORFS T, AARTS A A, ERISMIS M A, et al. Two-dimensional multi-channel neural probes with electronic depth control[J]. IEEE Transactions on Biomedical Circuits and Systems, 2011, 5(5): 403–412. DOI:10.1109/TBCAS.2011.2162840 |
| [14] | 余明, 陈锋, 李抄, 等. 微弱光强信号采样电路设计[J]. 光学仪器, 2014, 36(3): 253–257. DOI:10.3969/j.issn.1005-5630.2014.03.014 |
2019, Vol. 41
Issue (4): 74-79


