CN106949966B - 扫频光学相干层析成像系统的光谱标定方法 - Google Patents
扫频光学相干层析成像系统的光谱标定方法 Download PDFInfo
- Publication number
- CN106949966B CN106949966B CN201710183419.9A CN201710183419A CN106949966B CN 106949966 B CN106949966 B CN 106949966B CN 201710183419 A CN201710183419 A CN 201710183419A CN 106949966 B CN106949966 B CN 106949966B
- Authority
- CN
- China
- Prior art keywords
- array
- interference signal
- wave number
- sample
- transparent sheet
- 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
- 238000012014 optical coherence tomography Methods 0.000 title claims abstract description 44
- 238000001228 spectrum Methods 0.000 title claims abstract description 24
- 238000000034 method Methods 0.000 title claims abstract description 19
- 230000003287 optical effect Effects 0.000 claims abstract description 26
- 238000005070 sampling Methods 0.000 claims abstract description 13
- 238000012545 processing Methods 0.000 claims abstract description 7
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 6
- 230000009466 transformation Effects 0.000 claims description 15
- 230000008859 change Effects 0.000 claims description 6
- 230000001419 dependent effect Effects 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 5
- 239000006185 dispersion Substances 0.000 claims description 5
- 230000003595 spectral effect Effects 0.000 claims description 5
- 230000008878 coupling Effects 0.000 claims description 4
- 238000010168 coupling process Methods 0.000 claims description 4
- 238000005859 coupling reaction Methods 0.000 claims description 4
- 238000000844 transformation Methods 0.000 claims description 3
- 238000004587 chromatography analysis Methods 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 12
- 238000003384 imaging method Methods 0.000 description 8
- 230000001427 coherent effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 239000013307 optical fiber Substances 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 210000002294 anterior eye segment Anatomy 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 210000004351 coronary vessel Anatomy 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 210000001525 retina Anatomy 0.000 description 1
- 238000012421 spiking Methods 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/2823—Imaging spectrometer
Abstract
一种扫频光学相干层析成像系统的光谱标定方法，在扫频光学相干层析成像系统的参考臂中放置一个透明薄片，利用该透明薄片可以获得光程差为2nd的固定光程差干涉信号。通过密集采样采集系统的干涉信号，该信号包含薄片的干涉信号和样品干涉信号，利用薄片的干涉信号进行光谱标定，获得波数等间隔的样品干涉信号，再进行逆傅里叶变换的数据处理和图像重建算法，即可获得样品的层析结构图像。本发明利用一个透明薄片实现系统干涉信号的实时标定，结构简单，使用方便，避免了使用马赫‑曾德(MZI)干涉仪等额外的器件，降低了系统成本和复杂性。
Description
技术领域
本发明涉及扫频光学相干层析成像(Swept source Optical CoherenceTomography，简称SS-OCT)，特别是一种用于SS-OCT的光谱标定方法。
背景技术
光学相干层析成像(Optical Coherence Tomography，以下简称为OCT)是一种通过探测样品背向散射光的强度来获取样品不同深度结构信息的生物医学光学成像技术，具有非侵入、分辨率高、可在体检测生物组织内部微结构信息等特点。1991年，美国麻省理工学院的J.G.Fujimoto和D.Huang等人首先提出了此概念，并对视网膜和冠状动脉进行了离体成像。目前OCT技术已被广泛应用于眼科、皮肤科、心血管等领域的临床诊断和研究。
OCT可分为时域OCT(Time Domain OCT，TD-OCT)和频域OCT(Fourier Domain OCT，FD-OCT)，而频域OCT又可分为基于光谱仪的频域OCT(Spectral Domain OCT，SD-OCT)和基于扫频光源的频域OCT(又称扫频光学相干层析成像，Swept source OCT，SS-OCT)。扫频光学相干层析成像系统主要包括扫频光源1、迈克尔逊干涉仪2、平衡探测器3和数据采集卡4。扫频光源发出的光经迈克尔逊干涉仪产生干涉信号，由平衡探测器3探测，数据采集卡采集获得干涉信号随波长变化的强度分布(干涉谱)，然后对其做倒数变换后得到干涉信号在频域的强度分布，即频域干涉条纹。基于待测物体内各层光反射或背向散射界面的深度对应频域干涉条纹的不同频域的原理，SS-OCT对频域干涉条纹作逆傅里叶变换得到待测物体沿照明光光轴方向的深度分辨的反射率或背向散射率分布，即层析图。在SS-OCT中，样品深度和光源波数是一对傅里叶变换对，将采集的干涉信号沿波数K进行傅里叶变换(简称为FFT)，获得样品的层析结构图像。对非波数等间隔的干涉信号直接进行FFT会产生一个展宽因子，导致系统的轴向分辨率和测量精度降低，影响图像质量。通过光谱标定，可以获得波数等间隔的干涉信号。因此，光谱标定在频域OCT系统中具有重要的作用。
扫频光学相干层析成像(SS-OCT)系统中，数据采集卡密集采样时采集到的干涉信号是时间等间隔分布，而扫频光源输出光的波数通常并不是随时间线性变化的，因此需要在傅里叶变换之前获得波数等间隔的干涉信号。目前主要有硬件和软件两种方式获得波数等间隔的干涉信号。
硬件的方式主要有：
1)采用一个窄带固定波长滤波器(如光纤布拉格光栅)来标定扫频输出光的起始点，并且认为波长随时间线性分布，对干涉光谱信号进行插值，得到波数等间隔分布的干涉光谱采样点(参见在先技术[1]，B.J.Vakoc,S.H.Yun,J.F.de Boer,G.J.Tearney,andB.E.Bouma，“Phase-resolved optical frequency domain imaging”，Opt.Express，Vol.13,No.14，5483–5493，2005)；
2)在扫频OCT系统中，将光源的一部分输出光引入到光纤法布里-珀罗干涉仪(FFPI)，当扫频光源的输出光谱覆盖FFPI的自由光谱范围时，产生一个尖峰信号，触发采集卡采集一个采样点，获得波数等间隔的干涉信号(参见在先技术[2]，Michael A.Choma,Kevin Hsu and Joseph A.Izatt，“Swept source optical coherence tomography usingan all-fiber 1300-nm ring laser source”，Journal of Biomedical Optics，Vol.10，No.4，044009，2005)；
3)J.G.Fujimoto等人用马赫曾德干涉仪(MZI)代替FFPI，产生具有固定光程差的正弦干涉信号。通过最近邻域比较算法、遗传算法和插值相结合等方法得到等波数间隔采样点，同步触发数据采集卡，获得波数等间隔的干涉信号(参见在先技术[3]，R.Huber,M.Wojtkowski and J.G.Fujimoto，“Three-dimensional and C-mode OCT imaging witha compact,frequency swept laser source at 1300nm”，Opt.Express，Vol.13,No.25，10523–10538，2005)。
4)另外，XingdeLi等人采用马赫曾德干涉仪结合过零点探测电路得到等波数间隔采样时钟，同步触发数据采集卡，获得波数等间隔的干涉信号(参见在先技术[4]，JiefengXi,Li Huo,Jiasong Li and Xingde Li，“Generic real-time uniform K-spacesampling method for high-speed swept-Source optical coherence tomography”,Opt.Express Vol.18,No.9，9511-9517，2010)。
上述基于硬件方法的精度与采样点数受到触发同步硬件性能的限制，如采用窄带固定波长滤波器只能将标定起始点的波长固定于一点，FFPI只能产生有限多的采样点，而MZI干涉仪的光程差限制了系统的波数等间隔采样点数。此外，硬件的方法需要引入额外器件，增加了系统的复杂性。
软件的方式有两种：
1)Zhihua Ding等提出的基于非均匀离散傅里叶变换的方法，直接对采集的波长等间隔的干涉信号进行非均匀离散傅里叶变换，但该方法数据处理量大、耗时，不适用于实时显示(参见在先技术[5]，Tong Wu，Zhihua Ding，Kai Wang，Chuan Wang，“Swept sourceoptical coherence tomography based on non-uniform discrete fouriertransform”,中国光学快报(英文版),Vol.7,No.10，941-944，2009)；
2)Y.Yasuno等提出的基于反射镜干涉信号的校正方法，该方法先用反射镜作为样品获得干涉信号，再利用该干涉信号的相位信息进行光谱校准，但是需要已知光源波长随时间的变化关系，同时认为光源至少在数个月内保持稳定，且不能实时校准采集的样品干涉信号(参见在先技术[6]，Y.Yasuno,V.D.Madjarova,S.Makita,M.Akiba,A.Morosawa,C.Chong,T.Sakai,K.P.Chan,M.Itoh,and T.Yatagai,“Three-dimensional and high-speed swept-source optical coherence tomography for in vivo investigation ofhuman anterior eye segments,”Opt.Express Vol.13,No.25 10652–10664 2005)。
发明内容
本发明的目的是为了克服上述在先技术的不足，提供一种扫频光学相干层析成像系统的光谱标定方法。在扫频光学相干层析成像系统的参考臂中放置一个折射率为n、厚度为d的透明薄片，利用该透明薄片的可以获得具有固定光程差的干涉信号，再获得该干涉信号的相位信息，从而实现对样品干涉信号的光谱标定。
本发明的技术解决方案如下：
一种扫频光学相干层析成像系统的光谱标定方法，该扫频光学相干层析成像系统中扫频光源的输出端口与第一耦合器的第一端口相连，第二端口连接至第二耦合器的第一端口，第一耦合器将光分为两路，分别连接至样品光路和参考光路，其中样品光路中，第一耦合器的第三端口连接至第一准直透镜，准直后的光出射方向顺序放置二维扫描振镜、聚焦透镜和待测样品；参考光路中，输出光的出射方向顺序放置第二光束准直器、色散补偿块和第三准直透镜，输出端口连接至第二耦合器的第二端口，第二耦合器的两个输出端口和连接至平衡探测器，再通过数据采集卡与计算机相连，其特点在于，该光谱标定法包括以下步骤：
①在所述的参考光路中第二光束准直器和第三准直透镜之间放置一个透明薄片，该透明薄片的前后表面能够产生光程差为2nd的干涉信号：
当参考光路是透射式结构的参考臂，则所述的透明薄片的折射率为n，厚度为d；
当参考光路是反射式结构的参考臂，则所述的透明薄片的折射率为n，厚度为d/2；
②启动扫频光学相干层析成像系统，所述的透明薄片产生光程差为2nd的透明薄片干涉信号阵列Is，同时参考臂的光和样品臂的光在第二个耦合器发生干涉，得到样品干涉信号阵列Irs，透明薄片干涉信号阵列Is和样品干涉信号阵列Irs对应点的数值叠加得到信号阵列I，经探测器探测后通过数据采集卡传输到计算机；
所述的透明薄片干涉信号阵列Is、样品干涉信号阵列Irs和信号阵列I均是含有m*p个采样点数的阵列，m表示深度(轴向)方向的采样点数，p表示每帧(B-scan)包含的横向扫描点数，其中干涉信号阵列Is沿横向扫描点方向的数值相同；
③信号阵列I沿横向扫描点方向的对应点求和，再除以p获得透明薄片干涉信号阵列Is，或者通过傅里叶变换后找到透明薄片干涉信号阵列Is所对应的频率再滤波该信号的方法获得透明薄片干涉信号阵列Is，然后再将信号阵列I的每一列减去透明薄片干涉信号阵列Is的对应点的值得到样品干涉信号阵列Irs；
④对透明薄片干涉信号阵列Is进行傅里叶变换，分别求得傅里叶变换后的左右信号并滤出，接着对这两个信号逆傅里叶变换后求取其相位阵列并解包裹，获得两个解包裹的相位阵列，然后选出与光源波数变化趋势相同的相位阵列
⑤将光源的光谱密度函数S(k)的波长扫描范围等分成m-1份，得到含有m个点的初始波长阵列λ，并通过公式K0＝2π/λ得到与λ对应的初始波数阵列K0；
⑥设定最大循环次数N，将初始波数阵列K0作为自变量、相位阵列作为因变量进行三次幂以上的多项式曲线拟合，获得一组拟合系数，再根据求出的多项式中二次幂及更高次幂多项式曲线拟合系数，计算出二次幂及更高次幂对应的非线性相位由公式求出相应的波数变化量ΔK0，其中，Z_peak0＝2*π*P_Idx0./(kmax0-kmin0),P_Idx0表示干涉信号阵列Is傅里叶变换后相应的峰值位置，kmax0和kmin0分别表示初始波数阵列K0的最大值和最小值；然后，从初始波数阵列K0的每一项减去对应项的波数变化量ΔK0，获得新的波数阵列K1；
⑦令Ki＝K1，将波数阵列Ki作为自变量、相位阵列作为因变量进行三次幂以上的多项式曲线拟合，获得一组新的拟合系数，再根据求出的二次幂及更高次幂多项式曲线拟合系数计算出二次幂及更高次幂对应的非线性相位由公式其中Z_peaki＝2*π*P_Idxi./(kmaxi-kmini),P_Idxi表示透明薄片干涉信号阵列Is傅里叶变换后相应的峰值位置，kmaxi和kmini分别表示波数阵列Ki的最大值和最小值；求出相应的波数变化量ΔKi，直至满足最大循环次数N或者波数变化量ΔKi的均值最小，得到最终的波数阵列KN；
⑧在最终的波数阵列KN的起始点和终止点之间等间隔均分m-1份，获得含有m个点的等间隔的波数阵列KN＇，再通过公式λN＝2π/KN＇获得与KN＇对应的波数等间隔的波长阵列λN，利用初始波长阵列λ和波长阵列λN对样品干涉信号阵列Irs进行插值，得到波数等间隔分布的样品干涉信号阵列Irs(λN)，对其进行逆傅里叶变换的数据处理和图像重建算法，得到样品的层析图像。
实施上述方法的扫频光学相干层析成像系统主要包括扫频光源、迈克尔逊干涉仪、平衡探测器和数据采集卡。所述的迈克尔逊干涉仪具有两个接近等光程的干涉光路，分别为参考光路和样品光路。
所述的样品光路将光束聚焦到样品并进行扫描成像，可以通过二维扫描振镜扫描样品，也可以通过移动精密机械位移台扫描样品。
所述的平衡探测器是具有光电信号转换功能和放大差分信号的探测器。
该系统的工作情况如下：
扫频光源发出的光经迈克尔逊干涉仪后分为参考光路和样品光路，样品光路中的光经一系列光学器件后聚焦到待测样品；参考光路中须通过透明薄片，用于获得具有固定光程差的干涉信号。样品内部不同深度的背向散射光被重新收集并沿原光路返回，该返回光和参考臂的耦合光发生干涉后接入平衡探测器，探测器探测的干涉信号通过数据采集卡送入计算机进行数据处理，得到待测样品沿深度方向的层析图。通过二维振镜沿光轴垂直方向扫描获得样品上不同位置的层析图，得到待测样品的二维或三维层析图。
本发明与现有技术相比有益的效果是：
1.与在先技术[1]、[2]、[3]和[4]相比，本发明不需要MZI等额外复杂的器件，减少了系统的成本和复杂性。
2.与在先技术[5]相比，本发明的数据计算量减小，节省时间，可实时显示成像结果。
3.与在先技术[6]相比，本发明认为光源在一个B-scan内保持稳定，降低了光源不稳定性的影响，而且采用傅里叶变换后再滤波的方式可以实时校准。
附图说明
图1是光纤型扫频光学相干层析成像系统的结构示意图。
图2是采用精密位移平台对样品进行扫描的样品臂光路示意图。
图3是反射式结构的参考臂光路示意图。
图4是光谱标定的流程图。
图5是模拟的光谱标定前的波长等间隔的两层样品的A-line信号图。
图6是模拟的光谱标定后的波长等间隔的两层样品的A-line信号图。
具体实施方式
下面结合实施例和附图对本发明作进一步说明，但不应以此实施例限制本发明的保护范围。
请参阅图1，图1为本发明的光纤型扫频光学相干层析成像系统结构示意图。其中包括扫频光源1，扫频光源1的输出端口与第一耦合器21的第一端口211相连，第二端口212连接至第二耦合器24的第一端口241，第一耦合器21将光分为两路213和214，分别连接至样品光路22和参考光路23。样品光路22中，第一耦合器21的第三端口213连接至第一准直透镜221，准直后的光出射方向顺序放置二维扫描振镜222、聚焦透镜223和待测样品224；参考光路23中，第一耦合器21的另一个出射端口214，输出光的出射方向顺序放置第二光束准直器231、透明薄片232、色散补偿块233和第三准直透镜234，输出端口连接至第二耦合器24的第二端口242。第二耦合器24的两个输出端口243和244连接至平衡探测器3，再通过数据采集卡4与计算机5相连。
扫频光源1的输出扫频激光进入第一耦合器21中并分成两路213和214，一路213进入样品臂光路22，并经过第一准直透镜221、二维扫描振镜222和聚焦透镜223后聚焦到待测样品224上；另一路214进入参考臂光路23，经过第二准直透镜231、透明薄片232、色散补偿块233和第三准直透镜234后输出耦合光。从待测样品234内部不同深度处背向散射回来的光被重新收集回第一耦合器21，并由第二端口212输出。该第二端口212的输出光与参考光路23输出的耦合光，输入到第二耦合器24的第一端口241和第二端口242，在第二耦合器24产生干涉信号，该干涉信号经平衡探测器3将光信号转换成电信号，经数据采集卡4采集后送入计算机5中。该干涉信号在计算机5中进行数据处理后获得透明薄片232的干涉信号和样品的干涉信号。由透明薄片232的干涉信号获得波数等间隔采样点，利用这些采样点对系统中被测样品234的干涉信号插值获得波数等间隔的干涉信号，逆傅里叶变换后得到被测样品234沿深度方向的层析图。通过二维扫描振镜223对待测样品234沿光轴垂直方向做横向扫描获得被测样品234的二维或三维层析图。
本实施例中，扫频光源的中心波长是1050nm，带宽约为102nm。第一准直透镜、第二准直透镜和第三准直透镜具有相同的参数，色散补偿块233是由与聚焦透镜223相同的两个透镜组成的4f系统。
图4是本发明扫频光学相干层析成像系统光谱标定的流程图，如图所示，包括下列步骤：
①在所述的参考光路中第二光束准直器和第三准直透镜之间放置一个透明薄片232，该透明薄片232的前后表面能够产生光程差为2nd的干涉信号：
当参考光路是透射式结构的参考臂，则所述的透明薄片232的折射率为n，厚度为d；
当参考光路是反射式结构的参考臂，则所述的透明薄片232的折射率为n，厚度为d/2；
启动扫频光学相干层析成像系统，所述的透明薄片产生光程差为2nd的透明薄片干涉信号阵列Is，同时参考臂的光和样品臂的光在第二个耦合器24发生干涉，得到样品干涉信号阵列Irs，透明薄片干涉信号阵列Is和样品干涉信号阵列Irs对应点的数值叠加得到信号阵列I，经探测器3探测后通过数据采集卡4传输到计算机5；
所述的透明薄片干涉信号阵列Is、样品干涉信号阵列Irs和信号阵列I均是含有m*p个采样点数的阵列，m表示深度(轴向)方向的采样点数，p表示每帧(B-scan)包含的横向扫描点数，其中干涉信号阵列Is沿横向扫描点方向的数值相同；
②信号阵列I沿横向扫描点方向的对应点求和，再除以p获得透明薄片干涉信号阵列Is，或者通过傅里叶变换后找到透明薄片干涉信号阵列Is所对应的频率再滤波该信号的方法获得透明薄片干涉信号阵列Is，然后再将信号阵列I的每一列减去透明薄片干涉信号阵列Is的对应点的值得到样品干涉信号阵列Irs；
③对透明薄片干涉信号阵列Is进行傅里叶变换，分别求得傅里叶变换后的左右信号并滤出，接着对这两个信号逆傅里叶变换后求取其相位阵列并解包裹，获得两个解包裹的相位阵列，然后选出与光源波数变化趋势相同的相位阵列
④将光源的光谱密度函数S(k)的波长扫描范围等分成m-1份，得到含有m个点的初始波长阵列λ，并通过公式K0＝2π/λ得到与λ对应的初始波数阵列K0；
⑤设定最大循环次数N，将初始波数阵列K0作为自变量、相位阵列作为因变量进行三次幂以上的多项式曲线拟合，获得一组拟合系数，再根据求出的多项式中二次幂及更高次幂多项式曲线拟合系数，计算出二次幂及更高次幂对应的非线性相位由公式求出相应的波数变化量ΔK0，其中，Z_peak0＝2*π*P_Idx0./(kmax0-kmin0),P_Idx0表示干涉信号阵列Is傅里叶变换后相应的峰值位置，kmax0和kmin0分别表示初始波数阵列K0的最大值和最小值；然后，从初始波数阵列K0的每一项减去对应项的波数变化量ΔK0，获得新的波数阵列K1；
⑥令Ki＝K1，将波数阵列Ki作为自变量、相位阵列作为因变量进行三次幂以上的多项式曲线拟合，获得一组新的拟合系数，再根据求出的二次幂及更高次幂多项式曲线拟合系数计算出二次幂及更高次幂对应的非线性相位由公式其中Z_peaki＝2*π*P_Idxi./(kmaxi-kmini),P_Idxi表示透明薄片干涉信号阵列Is傅里叶变换后相应的峰值位置，kmaxi和kmini分别表示波数阵列Ki的最大值和最小值；求出相应的波数变化量ΔKi，直至满足最大循环次数N或者波数变化量ΔKi的均值最小，得到最终的波数阵列KN；
⑦在最终的波数阵列KN的起始点和终止点之间等间隔均分m-1份，获得含有m个点的等间隔的波数阵列KN＇，再通过公式λN＝2π/KN＇获得与KN＇对应的波数等间隔的波长阵列λN，利用初始波长阵列λ和波长阵列λN对样品干涉信号阵列Irs进行插值，得到波数等间隔分布的样品干涉信号阵列Irs(λN)，对其进行逆傅里叶变换的数据处理和图像重建算法，得到样品的层析图像。
图5是模拟的光谱标定前的波长等间隔的两层样品的A-line信号图。图6是模拟的光谱标定后的波长等间隔的两层样品的A-line信号图。从图中可以看出，当采集的样品干涉信号是非波数等间隔时，A-line信号发生了严重的展宽。而且当两层距离较近时，容易发生混叠，甚至无法分开。
本发明应用于扫频光学相干层析成像系统的光谱标定，利用透明薄片产生的干涉信号对波数进行校准，从而实现系统中样品干涉信号的光谱标定，获得波数等间隔的系统样品干涉信号。
Claims (1)
1.一种扫频光学相干层析成像系统的光谱标定方法，该扫频光学相干层析成像系统中扫频光源的输出端口与第一耦合器的第一端口相连，第二端口连接至第二耦合器的第一端口，第一耦合器将光分为两路，分别连接至样品光路和参考光路，其中样品光路中，第一耦合器的第三端口连接至第一准直透镜，准直后的光出射方向顺序放置二维扫描振镜、聚焦透镜和待测样品；参考光路中，输出光的出射方向顺序放置第二光束准直器、色散补偿块和第三准直透镜，输出端口连接至第二耦合器的第二端口，第二耦合器的两个输出端口连接至平衡探测器，再通过数据采集卡与计算机相连，其特征在于，该光谱标定法包括以下步骤：
①在所述的参考光路中第二光束准直器和第三准直透镜之间放置一个透明薄片(232)，该透明薄片(232)的前后表面能够产生光程差为2nd的干涉信号：
当参考光路是透射式结构的参考臂，则所述的透明薄片(232)的折射率为n，厚度为d；
当参考光路是反射式结构的参考臂，则所述的透明薄片(232)的折射率为n，厚度为d/2；
②启动扫频光学相干层析成像系统，所述的透明薄片(232)产生光程差为2nd的透明薄片干涉信号阵列Is，同时参考臂的光和样品臂的光在第二个耦合器(24)发生干涉，得到样品干涉信号阵列Irs，透明薄片干涉信号阵列Is和样品干涉信号阵列Irs对应点的数值叠加得到信号阵列I，经探测器(3)探测后通过数据采集卡(4)传输到计算机(5)；
所述的透明薄片干涉信号阵列Is、样品干涉信号阵列Irs和信号阵列I均是含有m*p个采样点数的阵列，m表示深度方向的采样点数，p表示每帧包含的横向扫描点数，其中干涉信号阵列Is沿横向扫描点方向的数值相同；
③信号阵列I沿横向扫描点方向的对应点求和，再除以p获得透明薄片干涉信号阵列Is，或者通过傅里叶变换后找到透明薄片干涉信号阵列Is所对应的频率再滤波该信号的方法获得透明薄片干涉信号阵列Is，然后再将信号阵列I的每一列减去透明薄片干涉信号阵列Is的对应点的值得到样品干涉信号阵列Irs；
④对透明薄片干涉信号阵列Is进行傅里叶变换，分别求得傅里叶变换后的左右信号并滤出，接着对这两个信号逆傅里叶变换后求取其相位阵列并解包裹，获得两个解包裹的相位阵列，然后选出与光源波数变化趋势相同的相位阵列
⑤将光源的光谱密度函数S(k)的波长扫描范围等分成m-1份，得到含有m个点的初始波长阵列λ，并通过公式K0＝2π/λ得到与λ对应的初始波数阵列K0；
⑥设定最大循环次数N，将初始波数阵列K0作为自变量、相位阵列作为因变量进行三次幂以上的多项式曲线拟合，获得一组拟合系数，再根据求出的多项式中二次幂及更高次幂多项式曲线拟合系数，计算出二次幂及更高次幂对应的非线性相位由公式求出相应的波数变化量ΔK0，其中，Z_peak0＝2*π*P_Idx0./(kmax0-kmin0),P_Idx0表示干涉信号阵列Is傅里叶变换后相应的峰值位置，kmax0和kmin0分别表示初始波数阵列K0的最大值和最小值；然后，从初始波数阵列K0的每一项减去对应项的波数变化量ΔK0，获得新的波数阵列K1；
⑦令Ki＝K1，将波数阵列Ki作为自变量、相位阵列作为因变量进行三次幂以上的多项式曲线拟合，获得一组新的拟合系数，再根据求出的二次幂及更高次幂多项式曲线拟合系数计算出二次幂及更高次幂对应的非线性相位由公式其中Z_peaki＝2*π*P_Idxi./(kmaxi-kmini),P_Idxi表示透明薄片干涉信号阵列Is傅里叶变换后相应的峰值位置，kmaxi和kmini分别表示波数阵列Ki的最大值和最小值；求出相应的波数变化量ΔKi，直至满足最大循环次数N或者波数变化量ΔKi的均值最小，得到最终的波数阵列KN；
⑧在最终的波数阵列KN的起始点和终止点之间等间隔均分m-1份，获得含有m个点的等间隔的波数阵列KN′，再通过公式λN＝2π/KN′获得与KN′对应的波数等间隔的波长阵列λN，利用初始波长阵列λ和波长阵列λN对样品干涉信号阵列Irs进行插值，得到波数等间隔分布的样品干涉信号阵列Irs(λN)，对其进行逆傅里叶变换的数据处理和图像重建算法，得到样品的层析图像。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710183419.9A CN106949966B (zh) | 2017-03-24 | 2017-03-24 | 扫频光学相干层析成像系统的光谱标定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710183419.9A CN106949966B (zh) | 2017-03-24 | 2017-03-24 | 扫频光学相干层析成像系统的光谱标定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106949966A CN106949966A (zh) | 2017-07-14 |
| CN106949966B true CN106949966B (zh) | 2018-05-08 |
Family
ID=59472526
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710183419.9A Active CN106949966B (zh) | 2017-03-24 | 2017-03-24 | 扫频光学相干层析成像系统的光谱标定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106949966B (zh) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107495921B (zh) * | 2017-09-29 | 2019-02-19 | 视微影像(河南)科技有限公司 | 一种光学相干断层扫描成像系统 |
| JP7252977B2 (ja) * | 2018-01-26 | 2023-04-05 | アキリス・ソシエテ・アノニム | 波長掃引式光干渉断層撮影システム用の収集装置 |
| CN108309245B (zh) * | 2018-02-23 | 2020-12-01 | 东北大学秦皇岛分校 | 一种谱域光学相干层析的相位补偿修正方法及系统 |
| CN109157187A (zh) * | 2018-09-06 | 2019-01-08 | 中国科学院上海光学精密机械研究所 | 增加扫频光学相干层析成像系统成像深度范围的方法 |
| CN109211415A (zh) * | 2018-11-20 | 2019-01-15 | 电子科技大学 | 一种基于光源光谱特征波长的波长标定方法 |
| CN109724695B (zh) * | 2018-12-14 | 2020-07-14 | 执鼎医疗科技（杭州）有限公司 | 一种光谱仪波长标定装置及方法 |
| CN111134612A (zh) * | 2020-01-15 | 2020-05-12 | 南京理工大学 | 基于超宽带光源谱域光学相干层析系统及光谱标定方法 |
| CN112683848B (zh) * | 2020-12-21 | 2022-09-02 | 中国科学院上海光学精密机械研究所 | 光学相干层析成像系统色散补偿方法 |
| CN112773335B (zh) * | 2021-02-07 | 2021-12-28 | 苏州大学 | 应用于光学相干弹性成像的传感器、成像系统及成像方法 |
| CN116861167B (zh) * | 2023-06-12 | 2023-12-15 | 河北工程大学 | 一种基于深度学习的fbg光谱循环去噪方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002517710A (ja) * | 1998-06-02 | 2002-06-18 | ゼテティック・インスティチュート | 波数ドメイン反射率計および背景振幅の減少および補償を使用する共焦干渉顕微鏡のための方法及び装置 |
| WO2008148237A1 (en) * | 2007-06-06 | 2008-12-11 | Exalos Ag | Optical coherence tomography sensor |
| WO2014043799A1 (en) * | 2012-09-24 | 2014-03-27 | Tornado Medical Systems Inc. | Pixel-shifting spectrometer on chip |
| CN105651374B (zh) * | 2016-01-27 | 2019-04-05 | 山西大学 | 单管共轴光声光谱测声器及采用该测声器的气体探测装置 |
-
2017
- 2017-03-24 CN CN201710183419.9A patent/CN106949966B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106949966A (zh) | 2017-07-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106949966B (zh) | 扫频光学相干层析成像系统的光谱标定方法 | |
| Izatt et al. | Theory of optical coherence tomography | |
| De Boer et al. | Twenty-five years of optical coherence tomography: the paradigm shift in sensitivity and speed provided by Fourier domain OCT | |
| CN109157187A (zh) | 增加扫频光学相干层析成像系统成像深度范围的方法 | |
| Götzinger et al. | High speed spectral domain polarization sensitive optical coherence tomography of the human retina | |
| Ju et al. | Advanced multi-contrast Jones matrix optical coherence tomography for Doppler and polarization sensitive imaging | |
| Yun et al. | Motion artifacts in optical coherence tomography with frequency-domain ranging | |
| Bachmann et al. | Resonant Doppler flow imaging and optical vivisection of retinal blood vessels | |
| US7697145B2 (en) | System for fourier domain optical coherence tomography | |
| US10502544B2 (en) | Efficient sampling of optical coherence tomography data for explicit ranging over extended depth | |
| Subhash | Full-field and single-shot full-field optical coherence tomography: a novel technique for biomedical imaging applications | |
| CN102818786B (zh) | 正弦相位调制的并行复频域光学相干层析成像系统和方法 | |
| US9636011B2 (en) | Systems and methods for spectrally dispersed illumination optical coherence tomography | |
| US10485422B2 (en) | System and method for imaging subsurface of specimen | |
| CN103439295A (zh) | 全深度频域多普勒光学相干层析成像方法 | |
| WO2007059206A2 (en) | Spectral domain phase microscopy (sdpm) dual mode imaging systems and related methods | |
| EP1870030A1 (en) | Apparatus and method for frequency domain optical coherence tomography | |
| EP1870028A1 (en) | Apparatus and method for frequency domain optical coherence tomography | |
| Yasuno et al. | Three-dimensional line-field Fourier domain optical coherence tomography for in vivo dermatological investigation | |
| EP3180581B1 (en) | Instantaneous time domain optical coherence tomography | |
| Poddar et al. | In vitro 3D anterior segment imaging in lamb eye with extended depth range swept source optical coherence tomography | |
| EP1870029A1 (en) | Apparatus and method for frequency domain optical coherence tomography | |
| Batista et al. | Swept-source Phase-Stabilized Optical Coherence Tomography Setup for Elastography. | |
| Wang et al. | A simple system of swept source optical coherence tomography for a large imaging depth range | |
| Wang et al. | Full-field swept source master-slave optical coherence tomography |
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 |