专利名称：并行复频域光学相干层析成像方法与系统的制作方法
技术领域：
本发明涉及频域光学相干层析成像(Fourier-Domain Optical CoherenceTomography，简称 FD-OCT)，尤其涉及一种基于空间载波(spatial carrier)的 并行复频域光学相干层析成像方法与系统。
背景技术：
光学相干层析成像(Optical Coherence Tomography，简称OCT)是近年来发展起 来的一种光学层析成像技术，它能够对高散射介质如生物组织内部几个毫米深度范围内的 微小结构进行高分辨率非侵入成像，在生物组织活体成像和医疗成像诊断等领域具有广泛 的应用前景。
频域光学相干层析成像系统(FD-OCT)是一种新型OCT系统，它通过探测干涉 谱并对其作逆傅里叶变换得到物体的层析图，相对于早先的时域光学相干层析成像系统 (Time-Domain Optical Coherence Tomography,简称 TD—0CT)具有无需深度方向扫描、成 像速度快和探测灵敏度高的优势，能更好地满足生物组织活体成像以及医疗成像诊断的实 时性要求。频域光学相干层析成像系统主要由低相干光源(宽光谱光源)、迈克尔逊干涉仪 和光谱仪(核心元件为分光光栅、聚焦透镜和CCD探测器)三部分组成，通过将低相干光源 发出的宽光谱光经迈克尔逊干涉仪产生的干涉信号送入光谱仪，获取干涉信号随波长(入) 变化的强度分布(干涉谱)，然后对其做倒数变换后得到干涉信号在频域(ν域，ν=1/ λ )的强度分布，即频域干涉条纹。基于待测物体内各层光反射或背向散射界面的深度对应 频域干涉条纹的不同频率的原理，FD-OCT对频域干涉条纹作逆傅里叶变换得到待测物体沿 照明光光轴方向的深度分辨的光反射率或背向散射率分布，即层析图。但是，FD-OCT获得 的层析图中包含着若干寄生像，限制了 FD-OCT的应用。这些寄生像分别是直流背景，自相 干噪声和复共轭镜像。其中，直流背景和自相干噪声的存在降低了 FD-OCT的信噪比，影响 了成像质量；而复共轭镜像的存在，使FD-OCT无法区分正负光程差(探测光路相对参考光 路的光程差)，测量时待测物体只能置于零光程差位置的一侧，导致有效探测深度范围减少 一半。复频域光学相干层析成像通过重建频域干涉条纹的复解析信号，对该复解析频域 干涉条纹信号作逆傅里叶变换得到物体的层析图，可以消除传统FD-OCT重建的层析图中 存在的寄生像，特别是复共轭镜像，从而使FD-OCT可以区分正负光程差，探测深度范围扩 大为原来的2倍，实现全深度探测。目前，已提出的复频域OCT方法主要包括基于移相干涉 术和基于外差干涉术的复频域OCT。基于移相干涉术(phase-shiftinginterferometry)的复频域 OCT2002年，A. F. Fercher等人最早基于移相干涉术重建复频域干涉条纹，实现了复 频域 OCT (参见在先技术[1]，M. ffojtkowski, A. Kowalczyk, R. Leitgeb andA. F. Fercher, "Full range complex spectral optical coherence tomographytechnique in eyeimaging”，Optics Letters, Vol. 27，No. 16，1415-1417，2002)。然而，由于该方法需要连 续或步进采集至少3幅相互之间具有固定相移量的移相干涉图，降低了频域OCT的成像 速度，并对干涉仪和样品的稳定性提出了严格的要求，因此该方法不适用于生物组织的 活体成像。2005年，Jowph A. Izatt等人提出基于同步移相干涉术实现复频域OCT(参 见在先技术[2], Μ. V. Sarunic, M. A. Choma, C. Yang and J. A. Izatt, “ Instantaneous complex conjugate resolved spectraldomain and swept—source OCT using 3X3 fiber couplers”，OpticsExpress, Vol. 13，No. 3，957-967，2005)。该方法虽然可以实现多幅移相 干涉条纹的同时获取，但需要使用NXN(N>3)光纤耦合器作为同步移相装置，增加了系统 的复杂性和成本，且移相精度容易受到环境温度变化的影响，从而影响复共轭镜像的消除 效果。基于外差干涉术(heterodyne interferometry)的复频域 OCT基于外差干涉术的复频域OCT通过在频域干涉信号中引入时间或空间载波实现 复频域干涉条纹的探测，与基于移相干涉术的复频域OCT相比，具有不受移相精度限制的 优点。2006年，Bachmarm等人采用两个声光晶体在干涉仪的参考光和探测光中引入光频 差，产生一个含有时间载波的频域干涉条纹，然后通过锁相探测频域干涉条纹的正交分量 重建复频域干涉条纹(参见在先技 术[3]，A. H. Bachmann, R. A. Leitgeb and Τ. Lasser, "Heterodyne Fourier domainoptical coherence tomography for full range p robing with high axialresolution，” Optics Express，Vol. 14，No. 4，1487-1496，2006)。该方 法通过时间相位调制在对应待测物体同一横向位置的频域干涉条纹中引入时间载波，虽然 避免了移相误差的影响，但每个时间相位调制周期内只能得到一个横向位置的全深度层析 图像(full range A-line)，降低了频域OCT的成像速度。2007年，Wang等人在对待测物 体的横向扫描过程中，通过同步扫描参考镜的光程在对应待测物体不同横向位置的频域干 涉条纹中引入空间载波，实现了横向扫描方向上对应连续横向位置的全深度层析图像的实 时狈Ij量。(参见在先技术[4]，R. K. Wang，“In vivo full rangecomplex Fourier domain optical coherence tomography”，Applied PhysicsLetters，Vol. 90，No. 054103，2007)。该 方法基于一种空间载波外差干涉术实现复频域0CT，具有不降低频域OCT成像速度的优点， 但是，它需要横向机械扫描，对数据采集和扫描位置控制的同步要求比较高；而且在探测动 态物体时，横向机械扫描速率必须大于待测物体的运动速率，否则图像会出现运动模糊，即 系统存在一个允许的物体最大速率，这限制了该方法的应用范围。由以上分析，目前复频域OCT存在成像速度慢、系统结构复杂、需要复杂的同 步扫描控制等问题，这些问题可以通过采用并行频域光学相干层析成像技术(parallel FD-0CT)解决。并行频域OCT与传统基于单点照明的频域OCT的主要区别在于，它通过 采用线状光照明样品实现频域OCT 二维层析图(B-scan)的并行探测。(参见在先技术 [5], Branislav Grajciar, Michael Pircher, Adolf F. Fercherand Rainer A. Leitgeb, "Parallel Fourier domain optical coherence tomographyfor in vivo measurement of the human eye”，Optics Express, Vol. 13，No. 4，2005)。该方法一般通过在光路中添 加柱面镜实现对待测物体的线状光照明，并利用二维光电探测阵列并行记录对应线照明光 长度方向上待测物体多个横向位置的频域干涉条纹，重建得到一幅待测物体的二维层析图 (B-scan)。并行频域OCT由于避免了对待测物体的横向机械式扫描，成像速度快，对运动模糊不敏感。但是，并行频域OCT仍然存在复共轭镜像等寄生像问题。
发明内容
本发明的目的是为了克服上述在先技术的不足，通过空间载波外差干涉术与并行 频域光学相干层析成像结合，提供一种并行复频域光学相干层析成像的方法与系统。本发 明只需一次曝光即可实现全深度的复频域光学相干层析成像，具有结构简单、成像速度快、 对运动模糊不敏感的特点。本发明的技术解决方案如下一种并行复频域光学相干层析成像的方法，该方法是在并行频域光学相干层析成 像方法的基础上，通过用倾斜的平面反射式衍射光栅代替干涉参考臂的参考平面反射镜， 并使入射参考光的一级衍射光沿原入射光路逆向返回，从而在二维光电探测器阵列获得的 二维频域干涉条纹沿并行探测方向上引入线性空间相位调制，即在二维频域干涉条纹中引 入空间载波；然后对含有空间载波的二维频域干涉条纹沿并行探测方向作傅里叶变换，接 着依次通过频域滤波、坐标平移和频域的逆傅里叶变换的过程，得到二维复频域干涉条纹， 最后再通过以波数为变量的逆傅里叶变换获得待测物体层析图。本发明并行复频域光学相干层析成像的方法的具体步骤如下
①在并行频域光学相干层析成像方法的基础上，在干涉参考臂的参考平面反射镜 改为倾斜放置的平面反射式衍射光栅，使得入射参考光的一级衍射光沿着原光路逆向返 回，参考光入射到平面反射式衍射光栅上的入射角θ应该满足关系式(1)sin^ = ^2-,(1)
2d其中λ ^是低相干光源的中心波长，d是平面反射式衍射光栅的周期常数。光谱仪并行记录对应待测样品上线状光照明区域内每一点的频域干涉条纹，即 一幅二维频域干涉条纹。参考光以入射角θ入射到平面反射式衍射光栅上，其一级衍射 光沿原光路逆向返回，在二维频域干涉条纹中沿X轴方向引入线性空间相位调制ψ(χ)=
f fItgQ
2kx · tg θ / σ，即在二维频域干涉条纹中引入空间载波Λο =H °其中λ代表波长，k = 2 π / λ代表波数；二维频域干涉条纹的两个维度分别对 应着待测样品上沿线状照明光长度方向的横向点经一维成像系统成像在光谱仪中二维光 电探测器阵列上的横向位置(X轴)和光源波长(y轴)；一维成像系统分别由迈克尔逊干涉 仪中平面反射式衍射光栅前第一聚焦透镜和待测样品前第二聚焦透镜与光谱仪中二维光 电探测器阵列前第三透镜组成，σ =F2ZiF1代表一维成像系统的横向放大率，F1代表迈克尔 逊干涉仪中在平面反射式衍射光栅和待测样品前第一、二聚焦透镜的焦距，F2代表光谱仪 中第三聚焦透镜的焦距；X'代表待测样品沿线照明光长度方向的横向位置，χ' =χ/σ ；②系统工作后，所述的二维光电探测器阵列记录的含有空间载波的二维频域干涉 条纹信号如式(2)所示 gik, X) = S_Q + 2 s{k)an (χ)
η
<formula>formula see original document page 7</formula>
<formula>formula see original document page 7</formula>
其中S(k)代表低相干光源的功率谱密度，β 0代表平面反射式衍射光栅一级衍 射的等效反射率，α n(x)、Qm(X)代表二维光电探测器阵列上横向位置X对应的待测样品上 横向位置χ'处第n、m层反射或散射界面的反射率或背向散射率，zn(X)、zm(X)代表二维光 电探测器阵列上横向位置χ对应的待测样品上横向位置X'处第n、m层反射或散射界面的 纵向深度。式(2)中前两项分别是平面反射式衍射光栅的反射光的自谱密度函数和待测样 品内各层深度处反射或背向散射光的自谱密度函数叠加项，第三项为待测样品内不同深度 处反射或背向散射光的互谱密度函数叠加项，第四项为平面反射式衍射光栅反射光和待测 样品内各层深度处反射或背向散射光的互谱密度函数叠加项。式(2)可以简化为式(3)<formula>formula see original document page 7</formula>
其中<formula>formula see original document page 7</formula>流分量和自相干噪声，它不受平面反射式衍射光栅引入的空间载波的调制。式⑶可以用式(4)表示<formula>formula see original document page 7</formula>
其中b (k,χ)= S<Jc\Jcc (χ)β0 exp[说(k,χ)]，Φη(1 ，χ) = 2kzn(x)，* 表示复共轭运 算；③对二维频域干涉信号式(4)作以χ为变量的傅里叶变换得到式(5)<formula>formula see original document page 7</formula>
其中G和B分别对应g和b的傅里叶频谱，fx代表对应χ轴的空间频谱，fx0代表 平面反射式衍射光栅引入的空间载波频率；并且其中《8和(^分别对应Gtl
和Σ5 的频谱带宽；
η④对步骤③所得的式(5)信号先乘上一个以fx(l为中心，2 (Ob为区间长度的矩形窗
函数炉(人)=化:It=^j==进行频域滤波，得到f5 0，y;-y;。)；⑤将Σ代伏，/c -Λ。)在频域坐标fx轴上左移fx(l，使其处于频域坐标轴原点上，得
η
到Σβ α，Λ);再将它作以fx为变量的逆傅里叶变换，即得到一个二维复频域干涉条纹，如
η
式(6)所示<formula>formula see original document page 8</formula>⑥对步骤⑤所得的二维复频域干涉条纹信号(6)作以k为变量的逆傅里叶变换得 到式(7)
<formula>formula see original document page 8</formula>其中符号表示以k为变量的逆傅里叶变换；Γ代表低相干光源功率谱的逆傅 里叶变换，即低相干光源的自相关函数。将关系式χ' = χ/σ代入式(7)得到式(8)<formula>formula see original document page 8</formula>⑦取/(Χ’,ζ)的幅度信息得到待测样品的二维层析图。⑧通过精密平移台对待测样品沿与线状照明光长度方向和该线状照明光的光轴 构成的平面垂直的水平方向作一维扫描，重复以上步骤② ⑦得到待测样品的三维层析 图。本发明方法得到的层析图与没有引入空间相位调制的并行频域OCT层析图式(9) 相比，消除了复共轭镜像(12)、直流背景(Itl)和自相干噪声(I1)三种寄生像，提高了信噪 比，实现了全深度探测的并行复频域光学相干层析成像。
<formula>formula see original document page 8</formula>
其中1。=㈡Γ(ζ)代表直流背景分量，
<formula>formula see original document page 8</formula>
代表自相干噪声分量，L=ZVqU力+ 代表复共轭镜像分量。
η实施上述方法的并行复频域光学层析成像系统包括低相干光源，在该低相干光源 的照明方向上顺序放置准直扩束器、柱面镜、迈克尔逊干涉仪，该迈克尔逊干涉仪的分光器 将入射光分为探测臂光路和参考臂光路，参考臂光路的末端为一个聚焦透镜和倾斜摆放的 平面反射式衍射光栅，探测臂光路的末端为一个聚焦透镜和待测样品，待测样品放置在一个精密移动平台上；迈克尔逊干涉仪输出端连接一光谱仪，该光谱仪由分光光栅，聚焦透镜 和二维光电探测器阵列组成；二维光电探测器阵列通过图像数据采集卡和计算机连接。该 系统的特点是所述的平面反射式衍射光栅倾斜摆放，使得参考光入射到平面反射式衍射光
栅的入射角为e =，从而使得入射光的一级衍射光正好沿原入射光路逆向返回。
2d所述的柱面镜将一束入射平行光会聚为一个线状聚焦光；所述的迈克尔逊干涉仪 中在平面反射式衍射光栅和待测样品前的聚焦透镜的焦距相同；所述的柱面镜与迈克尔逊 干涉仪中在平面反射式衍射光栅和待测样品前的聚焦透镜是共焦关系；所述的迈克尔逊干 涉仪中在平面反射式衍射光栅和待测样品前的聚焦透镜分别与光谱仪中的聚焦透镜是共 焦关系；所述的待测样品和平面反射式衍射光栅分别与二维光电探测器阵列在系统光路上 是物像共轭关系。所述的低相干光源为宽带光源，其光谱典型半高全宽为几十纳米到几百纳米，如 发光二极管(LED)或超辐射发光二极管(SLD)或飞秒激光器或超连续谱光源等。所述的准直扩束器由物镜和若干透镜组成。所述的迈克尔逊干涉仪，其特征在于具有两个接近等光程的干涉光路，一路为参 考臂光路，另一路为探测臂光路。所述的二维光电探测器阵列是面阵CXD或面阵CMOS或面阵InGaAs或其它具有光 电信号转换功能的二维探测器阵列。所述的精密移动平台可以沿三个互相垂直方向做微米级精度的平移。该系统的工作情况如下低相干光源发出的光经准直器扩束后，又经柱面镜在其会聚作用平面内聚焦，产 生一个线状聚焦光，然后在迈克尔逊干涉仪中待分成两束，一束光经过参考臂入射到平面 反射式衍射光栅上，其一级衍射光沿原入射光路逆向返回，另一束光经探测臂入射到待测 样品内，从平面反射式衍射光栅回来的衍射光和从待测样品内不同深度处反射或背向散射 回来的光波待收集并沿参考臂和探测臂返回，在迈克尔逊干涉仪中会合发生干涉，再送入 光谱仪分光并记录，经图像数据采集卡数模转换后送入计算机进行数据处理，得到待测样 品沿线状照明光长度方向和照明光光轴方向的一个二维层析图。通过精密平移台对待测样 品沿与线状照明光长度方向和照明光光轴垂直的方向作一维横向扫描，得到待测样品的三 维层析图。本发明与现有技术相比具有的有益效果是本发明并行复频域光学相干层析成像的方法的特点是将空间载波外差干涉术用 于并行频域光学相干层析成像，通过倾斜的平面反射式衍射光栅在并行探测方向上引入空 间载波，利用空间傅里叶变换分析方法重建低相干光频域干涉复振幅信号，消除了 FD-OCT 成像中存在的复共轭镜像、直流背景和自相干噪声三种寄生像，实现全深度探测的并行复 频域光学相干层析成像。与在先技术1相比，本发明只需通过一次曝光即可获得一幅复频域干涉条纹，对干涉仪和样品的稳定性要求不高。与在先技术2相比，本发明系统结构简单，成本低，抗环境干扰能力强。与在先技术3和4相比，本发明无需深度方向和横向的机械扫描，通过一次曝光即可获得一幅全深度的二维层析图，成像速度快；本发明不需要复杂的同步扫描控制，系统结 构简单，而且具有对运动模糊不敏感的优点。与在先技术5相比，本发明结合空间载波外差干涉术与并行频域OCT，解决了复共 轭镜像等寄生像问题，实现了全深度的并行复频域OCT测量。
图1为本发明并行复频域光学相干层析成像系统的侧视光路和系统结构示意图。图2为本发明并行复频域光学相干层析成像系统的俯视光路和系统结构示意图。
具体实施例方式下面结合实施例和附图对本发明作进一步说明，但不应以此限制本 发明的保护范围。请参阅图1和2。图1为本发明并行复频域光学相干层析成像系统的侧视光路和 系统结构示意图。图2为本发明并行复频域光学相干层析成像系统的俯视光路和系统结构 示意图。由图1和2可见，本发明并行复频域光学相干层析成像系统包括低相干光源1，在 该低相干光源1的照明方向上顺序放置准直扩束器2、柱面镜3、迈克尔逊干涉仪4，该迈克 尔逊干涉仪4的分光器41将入射光分为探测臂光路44和参考臂光路42，参考臂光路42的 末端为第一聚焦透镜46和倾斜摆放的平面反射式衍射光栅43，探测臂光路的末端为第二 聚焦透镜47和待测样品45，待测样品放置在一个精密移动平台(图中未示)上；迈克尔逊 干涉仪4的输出端连接一光谱仪5，该光谱仪5由分光光栅51，第三聚焦透镜52和面阵CXD 探测器53组成；面阵CXD探测器53通过图像数据采集卡6和计算机7连接。该系统的特 点是所述的平面反射式衍射光栅43倾斜摆放，使得入射光的一级衍射光正好沿原入射光 路逆向返回。所述的柱面镜3，它将一束入射平行光会聚为一条线状聚焦光；所述的第一聚焦 透镜46、第二聚焦透镜47的焦距相同；所述的柱面镜3与迈克尔逊干涉仪4中的第一聚焦 透镜46、第二聚焦透镜47是共焦关系；所述的迈克尔逊干涉仪4中的第一聚焦透镜46、第 二聚焦透镜47分别与光谱仪5中的第三聚焦透镜52是共焦关系；所述的待测样品45和平 面反射式衍射光栅43分别与面阵CCD探测器53在系统光路上是物像共轭关系。低相干光源1发出的宽光谱光经准直器2扩束后，又经柱面镜3在侧视光路平面 内会聚(见图1)，在俯视光路平面内平行透射(见图2)，产生一个线状聚焦光，然后在迈克 尔逊干涉仪4中待分光棱镜41分成两束，一束透射光经过参考臂光路42入射到平面反射 式衍射光栅43上，其一级衍射光沿原入射光路逆向返回，另一束反射光经探测臂光路44入 射到放置在精密平移台(图中未示)上的待测样品45内，从平面反射式衍射光栅43衍射 回来的一级衍射光和从待测样品45内不同深度处反射或背向散射回来的光波待收集并分 别沿参考臂光路42和探测臂光路44返回，在迈克尔逊干涉仪4中汇合发生干涉，再送入光 谱仪5待分光光栅51分光，经第三聚焦透镜52，成像在面阵CCD探测器53，转换成电信号 后，经图像数据采集卡6数模转换后送入计算机7进行数据处理，得到待测样品45沿线状 照明光长度方向和照明光光轴方向的一个二维层析图。所述的平面反射式衍射光栅43倾斜放置，使参考光入射到所述的平面反射式衍射光栅43上的入射角θ满足sin^ = ^.,(10)
Id其中λ ^是低相干光源1的中心波长，d是平面反射式衍射光栅43的光栅周期常
数。所述的光谱仪5并行记录了对应待测样品45上线状光照明区域内每一点的频域 干涉条纹，即一幅二维频域干涉条纹。参考光以入射角θ入射到平面反射式衍射光栅43 上，其一级衍射光沿原光路逆向返回，在二维频域干涉条纹中沿χ轴方向引入线性空间相
位调制Ψ(χ) =2kx*tg0/o，即在二维频域干涉条纹中引入空间载波/』=，。其中λ
σλ
代表波长，k = 2 π / λ代表波数；二维频域干涉条纹的两个维度分别对应着待测样品45上 沿线照明光长度方向的横向点经一维成像系统成像在光谱仪5中面阵CCD探测器53上的 横向位置(χ轴)和光源波长方向(y轴)；一维成像系统分别由迈克尔逊干涉仪4中第一 聚焦透镜46和第二聚焦透镜47与光谱仪5中第三聚焦透镜52组成，σ = F2ZF1代表一维 成像系统的横向放大率，F1代表聚焦透镜46、47的焦距，F2代表第三聚焦透镜52的焦距； x'代表待测样品45沿线照明光长度方向的横向位置，χ' = χ/ο。所述的CXD探测器53记录的二维频域干涉条纹信号为S(k,x) = g0(k,x) + 2ΣS叫α (χ)β0 cos[2k(zn(χ) + χ·tg9/σ)]，(11)
η其中S(k)代表低相干光源1的功率谱密度，β ^为平面反射式衍射光栅43 —级 衍射的等效反射率，α η(χ)、α m(x)代表CXD探测器53上横向位置χ对应的待测样品45上 横向位置χ'处第n、m层反射或散射界面的反射率或背向散射率，Zn(X)、zm(x)代表C⑶探 测器53上横向位置χ对应的待测样品45上横向位置χ'处第n、m层反射或散射界面的纵 向深度。其中g0 (k, χ) = + Σ S(k)a (χ) + S叫cc (x)am (χ) cos [2k(z (χ) - zm (χ))]为直
ηn*m
流分量和自相干噪声，它不受平面反射式衍射光栅引入的空间载波的调制。式(11)可以用式(12)表示S(k, x) = g0 (k, + (k, χ) exp [ 2π/χ0χ] + ^ b: (k, χ) exp [~ 2π/χ0χ], (12)
ηη其中& (Α:，χ)= Sik^an(X)β0 exp[从(Α:，χ)], Φη(1 ，χ) = 2kzn(x)，* 表示复共轭运
笪弁。然后对式(12)作以χ为变量的傅里叶变换得到式(13)G(k,fx) = G0(k,fx) + ^B (k,fx-fx0) + YB;(k,fx+fx0),(13)
ηη其中G和B分别对应g和b的傅里叶频谱，fx代表对应χ轴的空间频谱，fx0代表 平面反射式衍射光栅43引入的空间载波频率；并且厶>·^^，其中《8和(^分别对应
G0和Σ5 的频谱带宽。
η对式(13)先乘上一个以fx(1为中心，2cob为区间长度的矩形窗函数<formula>formula see original document page 12</formula>进行频域滤波，得到;-/^0);接着将它在频域坐标fx轴
上左移fx(1，使其处于频域坐标轴原点上，得到Σ5 @，λ);最后将它作以fx为变量的逆傅里
<formula>formula see original document page 12</formula>
叶变换，得到一个二维复频域干涉条纹，如式(14)所示<formula>formula see original document page 12</formula>对式(14)作以k为变量的逆傅里叶变换得到式(15)
<formula>formula see original document page 12</formula>其中符号表示以k为变量的逆傅里叶变换；Γ代表低相干光源1功率谱的逆 傅里叶变换，即低相干光源1的自相关函数。将关系式χ' = χ/σ代入式(15)得到式(16)<formula>formula see original document page 12</formula>取/(χ’，ζ)的幅度信息得到待测样品45的一个二维层析图。通过精密平移台(图中未示)对待测样品45沿与线状照明光长度方向和该线状 照明光的光轴构成的平面垂直的水平方向作一维扫描，重复以上过程得到待测样品45的 三维层析图。
权利要求
一种并行复频域光学相干层析成像的方法，特征在于该方法是在并行频域光学相干层析成像方法的基础上，通过用倾斜的平面反射式衍射光栅代替干涉参考臂的参考平面反射镜，并使入射参考光的一级衍射光沿原入射光路逆向返回，从而在二维光电探测器阵列获得的二维频域干涉条纹沿并行探测方向上引入线性空间相位调制，即在二维频域干涉条纹中引入空间载波；然后对含有空间载波的二维频域干涉条纹沿并行探测方向作傅里叶变换，接着依次通过频域滤波、坐标平移和频域的逆傅里叶变换的过程，得到二维复频域干涉条纹，最后再通过以波数为变量的逆傅里叶变换获得待测物体层析图。
2.根据权利要求1所述的并行复频域光学相干层析成像方法，其特征在于该方法的具 体步骤如下①在并行频域光学相干层析成像方法的基础上，在干涉参考臂的参考平面反射镜改为 倾斜放置的平面反射式衍射光栅，其倾斜角，即参考光入射到平面反射式衍射光栅的入射角<formula>formula see original document page 2</formula>X ^代表低相干光源的中心波长，d代表平面反射式衍射光栅的光栅周期 2d常数，则入射参考光的一级衍射光沿原入射光路逆向返回，从而在二维光电探测器阵列获 得的二维频域干涉条纹沿并行探测方向上引入线性空间相位调制V (X) = 2kx tg e / o，即在二维频域干涉条纹中引入空间载波<formula>formula see original document page 2</formula>其中代表波长，<formula>formula see original document page 2</formula>代表波 数，X代表待测样品和干涉参考臂平面反射式衍射光栅沿线照明光长度方向的横向位置经 一维成像系统成像在光谱仪中二维光电探测器阵列上的横向位置；所述的一维成像系统分 别由迈克尔逊干涉仪中平面反射式衍射光栅前的第一聚焦透镜和待测样品前的第二聚焦 透镜与光谱仪中二维光电探测器阵列前第三聚焦透镜组成，o = F2/Fi代表一维成像系统 的横向放大率，代表迈克尔逊干涉仪中在平面反射式衍射光栅前第一聚焦透镜和待测样 品前第二聚焦透镜的焦距，&代表光谱仪中二维光电探测器阵列前第三聚焦透镜的焦距； x'代表待测样品沿线照明光长度方向的横向位置，=x/o ；②系统工作后，所述的二维光电探测器阵列记录了含有空间载波的二维频域干涉信号<formula>formula see original document page 2</formula>其中<formula>formula see original document page 2</formula>表低相干光源的功率谱密度，3 o代表平面反射式衍射光栅一级衍射的等效反射率，a n(x)、 affl(x)代表二维光电探测器阵列上横向位置x对应的待测样品上横向位置x'处第n、m层 反射或散射界面的反射率或背向散射率，zn(x) >zffl(x)代表二维光电探测器阵列上横向位置 x对应的待测样品上横向位置x'处第n、m层反射或散射界面的纵向深度； 上式二维频域干涉信号又可以表示为其中bn(k,x) = S(k)ylan(x)J30 exp[吆(k,x)],小n(k，x) = 2kzn(x)，* 表示复共轭运算；③对二维频域干涉信号g(k，x)作以x为变量的傅里叶变换，得到 <formula>formula see original document page 3</formula>其中G和B分别对应g和b的傅里叶频谱，fx代表对应x轴的空间频谱；④将G(k，fx)乘上一个矩形窗函数进行频域滤波，得到∑Bn(k,fx-fx0)其中矩形窗函数为<formula>formula see original document page 3</formula>；wb为∑Bn的频谱带宽；⑤将∑Bn(k,fx-fx0)在频域坐标fx轴上左移fx(l，得到再以fx为变量作逆傅里叶变换得到二维复频域干涉信号gcomp(k，x)<formula>formula see original document page 3</formula>⑥将二维复频域干涉信号g。。mp(k，X)以k为变量作逆傅里叶变换，并代入关系式X'= x/0，得到<formula>formula see original document page 3</formula>其中r代表低相干光源功率谱的逆傅里叶变换，即低相干光源的自相关函数；⑦取/0c’，Z)的幅度信息得到待测样品的二维层析图。⑧通过精密平移台对待测样品(45)沿与所述的线状照明光长度方向和该线状照明光 的光轴构成的平面垂直的水平方向作一维扫描，重复以上步骤② ⑦得到待测样品(45) 的三维层析图。
3. 一种实现权利要求1或2所述方法的并行复频域光学相干层析成像系统，包括低相 干光源(1)，在低相干光源(1)的光束前进方向上顺序放置准直扩束器(2)、柱面镜(3)、迈 克尔逊干涉仪(4)，该迈克尔逊干涉仪(4)的分光器(41)将入射光分为探测臂光路(44) 和参考臂光路(42)，参考臂光路(42)的末端为第一聚焦透镜(46)和平面反射式衍射光 栅(43)，探测臂光路(44)的末端为第二聚焦透镜(47)和待测样品(45)，待测样品(45) 放置在一个精密移动平台上；该迈克尔逊干涉仪(4)的输出端连接一光谱仪(5)；该光谱 仪(5)由分光光栅(51)、第三聚焦透镜(52)和二维光电探测器阵列(53)组成；二维光电 探测器阵列(53)通过图像数据采集卡(6)和计算机(7)连接；其特征在于所述的平面反 射式衍射光栅(43)，倾斜摆放，使得参考光入射到所述的平面反射式衍射光栅的入射角为θ=arcsinλ0/2d所述的柱面镜(3)将一束入射平行光会聚为一个线状聚焦光；所述的第一聚焦透镜 (46)、第二聚焦透镜(47)的焦距相同；所述的柱面镜(3)与迈克尔逊干涉仪(4)中的第一 聚焦透镜(46)、第二聚焦透镜(47)是共焦关系；所述的迈克尔逊干涉仪(4)中的第一聚焦 透镜(46)、第二聚焦透镜(47)分别与光谱仪(5)中的第三聚焦透镜(52)是共焦关系；所 述的待测样品(45)和平面反射式衍射光栅(43)分别与二维光电探测器阵列(53)在系统 光路上是物像共轭关系。全文摘要
一种并行复频域光学相干层析成像方法与系统，该方法是在并行频域光学相干层析成像方法的基础上，通过用倾斜的反射光栅代替干涉参考臂的参考平面反射镜，并使入射参考光的一级衍射光沿原入射光路逆向返回，从而在二维光电探测器阵列获得的二维频域干涉条纹沿并行探测方向上引入线性空间相位调制，即在二维频域干涉条纹中引入空间载波；然后对含有空间载波的二维频域干涉条纹沿并行探测方向进行傅里叶变换，接着依次通过频域滤波窗滤波、坐标平移和沿频谱方向的逆傅里叶变换的过程，得到二维复频域干涉条纹，最后再通过沿光频方向的逆傅里叶变换获得待测物体层析图。本发明具有结构简单，成像速度快，只需一次曝光即可获得待测物体层析图。
文档编号G01J3/45GK101832817SQ20101011662
公开日2010年9月15日 申请日期2010年3月3日 优先权日2010年3月3日
发明者南楠, 步鹏, 王向朝, 黄炳杰 申请人:中国科学院上海光学精密机械研究所