光 学 学 报
0517001-
3 高速成像方法
FF-OCT 是基于多步移相干涉测量的方法,对干涉条纹(相位)的调制与提取是关键。当样品臂和参考臂
之间的光程差小于光源的相干长度时满足干涉条件,CCD 上每个像素点
p(x,y)
接收的干涉信号为
I(x,y)
:
I
( )
x,y = I
ˉ
( )
x,y + A
( )
x,y cos[ϕ(x,y)]
, (1)
式中
I
ˉ
( )
x,y
为由非相干光形成的背景图像的光强,
ϕ
( )
x,y
为参考光与样品光的相位差,
A
( )
x,y
为层析图像的振幅。
为了 解 调 出一幅 正 面 光学断 层 图 像的 振 幅 ,需要对 干 涉 信号进 行 调 制。 移 相 法为调 制 相 位的方 法 之
一,通过在参考光束和样品光束间引入移相量,采集几幅干涉图像后,对其光强信号进行加减运算滤除干涉
项与背景图像信息来获得样品的层析图,该方法具有操作简单、信息存储和处理量小的优点,常采用三步移
相法、四步移相法。在本课题组之前的研究中,系统使用四步移相法,每步移相
φ = π/2
,那么 CCD 上每个像
素点采集到的信号
I
( )
x,y
表示为
I
( )
x,y = I
ˉ
( )
x,y + A
i
( )
x,y cos
[ ]
ϕ
( )
x,y +
( )
n - 3 φ ,n = 1, 2, 3, … , (2)
式 中 cos
[ ]
ϕ
( )
x,y +
( )
n - 3 φ 表 示 移 相 后 的 干 涉 项 ,随 样 品 和 参 考 镜 上 对 应 点
( )
x,y
的 相 位 差 变 化 而 变 化 ;
A
i
( )
x,y
表示第
i
层组织层析图像的振 幅,应用移相公式 滤除背景光强
I
ˉ
( )
x,y
和干涉项即可求出样 品图像信
息
[18]
。操作 中,需要手动输 入对应 的电压值来 控制 PZT 的 步进,输入 4 次电压 值解调 出一幅断层 图,重复性
较差,不能减小单次移相引入的误差。此外,采用手动方式控制 PZT 步进量的方法导致采集速度较慢,每次
移相操作需要 15 s 左右的时间,一幅断层图从调制到采集时间约为 60 s,无法对活体组织进行成像。虽然也
可对 PZT 持续施加控制电压以获得连续步进的方法,但光源的相干长度为 l
c
= λ
2
/Δλ = 1.5125 μm ,四步总移
相量
φ = 2π
对应 来 回光 程 为
Δl = 0.55 μm
,连续 采集三 幅 图后 将 引入
Δ = 3Δl = 1.65 μm
的光 程 差。 由 于光
程差
Δ
超出相干范围
l
c
,干涉条纹消失,无法发生干涉。
为了克服上述缺点,其他研究小组采用了不同的解决方案。例如日本山形大学采用基于移相法的双通
道外差检 测法
[19]
,由两臂的液晶移相器(LCS)与脉冲配合对相干光引入
π 2
的相位差 ,且由两台 CCD 分别接
收,同 时采集 两幅干涉图,使得成 像时间大为缩短,但系统使用的迈克 尔孙干 涉结构两臂不能独立调 节,且
双 CCD 的同步与采集信号的计算量很大。法国物理化学研究院采用偏振干涉
[20-21]
,由光弹 性双折射调 制器
产生频 率(50 kHz)极高的 正弦信 号调制 两臂,参考光与样品光通过 CCD 前的 检偏器 后发生 干涉。该方法采
集速 度 快,但 所 使用的 仪 器非 常 复杂 ,由 于利用 的 是偏 振 特性,该 研究院 采 用的 样 品都被 认 为没 有 二向色
性。而实际应用中,胶原蛋白、髓 鞘、弹 性纤维 都是具有二向色性的生 物组织 ,采用 该方法 对样品的选择有
局限性。与前两者克服采集速度慢的方法不同,本文采用了基于载频法实现信号采集的方案。载频法对单
片机的控制 电路有严格 要求
[22]
,但 可以实现很 高的 信噪比 和更 快的采集速 度,比双 CCD 与光弹性 双折射器
件的方案容易处理,具有更好的生物医学应用潜力。
图 2 同步 CCD 与 PZT 的触发信号
Fig.2 Synchronic signal of CCD and PZT
载频法移相过程中 CCD 上每个像素点
p(x,y)
在不同时刻 t 时采集到的信号
I
( )
x,y, t
可表示为
I
( )
x,y, t = I
ˉ
( )
x,y + A
( )
x,y cos
[ ]
ϕ
( )
x,y + Ψ sin
( )
2πft + θ , (3)
式中
Ψ sin
( )
2πft + θ
代表单片机调制 PZT 的 正弦信号。为了得到最好的信噪比,消去初始相位差的干扰,需
要计算 单片机调制信号的幅值与相位。每个周期内 CCD 采集 4 幅 干涉图,那么 CCD 每像 素采集到的 1/4 周
期的振幅
E
p
(x, y)
可分别表示为
[5]
E
p
(x, y) =
∫
( p - 1)T/4
pT/4
I(x,y, t)dt,p = 1,2,3,4. (4)
将信号强度表达式按第一类 n 阶贝塞尔函数展开 并对一 幅图进 行运算 将正弦 与余弦部分分离得到
∑
s
3