邓文波,陈 华,聂 雄

(广西大学 计算机与电子信息学院,南宁 530004)

引 言

在数字共焦显微技术的研究中，显微镜光学系统的3维点扩散函数(three-dimensional point spread function，3-D-PSF)的研究对3维显微图像复原有着重要的意义。3-D-PSF的准确估计对3维显微图像的复原效果有很大的关系[1-4]。在3维显微图像复原方法的研究以及实际工程应用中，采用的是理论点扩散函数或者实验点扩散函数。理论点扩散函数可以通过光学系统传输函数模型进行构建[5-7]，也可以通过高斯型近似理论模型[3]构建。实验点扩散函数是通过实验获取，主要是使用一个直径尽可能小的荧光小球作为成像物体，在实验室条件下通过获取荧光小球的成像信息构建显微成像系统的点扩散函数[8-11]。也可以以黑白跳变的图像作为观察物，利用CCD像平面上采集图像进行拟合，获得光学系统的点扩散函数[12]。

由于实验获取的3维点扩散函数针对特定显微镜光学系统，因此,采用以正确方法获取图像复原处理用的实验点扩散函数，比采用理论点扩散函数进行图像复原具有更为准确的复原效果。

本文中以荧光微珠模拟点光源，通过数字共焦显微镜采集一系列荧光微珠的图像，构造显微镜光学系统的实验3维点扩散函数，采用图像平均法降低噪声对点扩散函数的影响，获取图像复原处理用的实验点扩散函数。

1 光学显微镜成像系统

1.1 3维点扩散函数理论模型

3-D-PSF在离散空域里可视为一个3维矩阵。是由沿着3-D-PSF双锥体中心光轴(z轴)不同的径向横截面(x-y面)组成，每一个横截面对应不同散焦量下的2维点扩散函数(two-dimensional point spread function，2-D-PSF)，其中z=0处中间截面为焦平面2-D-PSF。

圆形孔径小散焦光学传递函数(optical transfer function,OTF)的数学公式为[13]：

(1)

式中，α=4πws/λ，βcos-1(q/fc),s=2q/fc,w=δzdNA2/2,J1(x) 为第1类1阶贝塞尔函数，w为最大光程差，λ为发射光波长，q为频率，fc为系统截止频率，dNA为数值孔径，δz为样本空间的沿z轴的散焦距离，即散焦量。当散焦量δz取不同值时，得到不同的光程差w，进而通过(1)式得到不同散焦量的光学传递函数，通过傅里叶逆变换可以得到对应散焦量的一系列2-D-PSF，进而构成3-D-PSF。

1.2 显微镜成像原理

将厚度为T的3维样本用平面模型表示为:

(2)

式中,N=T/Δz,Δz为物体平面沿z轴方向的微小间距。由3维显微镜成像系统可知，样本在焦平面z′处所成像是每个平面与对应2-D-PSF卷积之和，即：

(3)

式中，g为显微镜成像，f为原始3维样本，h为对应显微镜系统3-D-PSF，h(x,y,z′-iΔz),(i=1,2,…,N)表示散焦量δz为z′-iΔz的光学系统的脉冲响应函数，即2-D-PSF。当f为点光源且位于该焦平面时，即只有当x=0,y=0,iΔz=z′时，f才有不为0的值，其余情况下都为0，由以上公式可知，点光源通过显微镜系统后所得的显微图像即为显微镜焦平面2-D-PSF。通过改变点光源的位置，即改变(3)式中iΔz的大小，便可以得到一系列不同散焦量的2-D-PSF，进而可得出显微系统的实验3-D-PSF。

2 3-D-PSF实验制作

2.1 实验点光源制作

点光源实验样品的制备：本实验中采用嵌于聚丙乙烯酰胺凝胶中的荧光微珠(F8803，Thermo Fisher公司)作为模拟点光源实验样品，荧光微珠直径为0.1μm，激发光波长为460nm～485nm，发射光波长为510nm～525nm。将荧光微珠溶液稀释1000倍，快速混匀后滴加至载玻片，待凝固后使用。

本实验是利用尼康Nikon Ti-E显微镜在物镜放大倍数M=40、数值孔径dNA=0.55的条件下对荧光微珠进行序列图像采集。

在使用实验测量点扩散函数时，作为成像物体的荧光微珠其直径应尽可能小。根据显微镜成像分辨率的判断标准，两个点光源刚好被分辨出来的条件是两个点光源在物平面的距离满足瑞利距离或者阿贝距离。瑞利距离计算公式为[14-15]：

δ=0.61λ/dNA

(4)

取发射波长λ=520nm，则δ=0.57μm，即只要两个光源的距离大于0.57μm时，便可通过显微镜进行分辨。本实验中采用的荧光微珠直径为0.1μm，小于该瑞利距离，表明该荧光微珠可作为空间的一点。

2.2 图像采集

2.2.1 光学切片图像采集 实验中使用尼康Nikon Ti-E显微镜，其搭载的高速电动控制系统，可以缩短图像获取的曝光时间，减少光毒对生物细胞的影响。该型号通过手柄或控制器可以控制电动XY载物台和z轴，实现光轴最小步进为0.025μm的3维显微图像的采集。

在显微镜物镜放大倍数M=40、数值孔径dNA=0.55条件下，利用型号为Nikon DS-Qi1的非致冷单色CCD对荧光微珠图像进行采集，其中光轴步进间隔Δ=0.2μm。图1为在显微镜所采集的荧光微珠图像。其中图1b为图1a中所框荧光微珠对应的局部放大图(图中灰度值是经过归一化处理，以便显示和观察，并非反映实际强度值,下同)。可以发现，即使稀释了1000倍，载玻片上的荧光微珠依然较多。在选取微珠研究对象过程中，一是注意要间距,由于散焦的影响，相邻微珠相互干扰，在选取过程中尽量加大相邻微珠的间距，以减少相互之间的干扰，提高采集图像的真实性;二是要确保所选微珠均有较高能量(亮度),在显微镜系统中，系统噪声不可忽视，随着散焦的不断增加，信号越来越弱，当信号大小低于噪声大小时，主要信息将被噪声所淹没。

Fig.1 a—microscope picture of fluorescent microsphere b—partial enlargement

图2为利用显微镜采集的一组不同散焦量的荧光微珠图像，即不同散焦量的2-D-PSF (图中数值为对应散焦量δz的值,下同)。通过一系列不同散焦量的2-D-PSF进而可以构造得到实验3-D-PSF。由图2可以明显看出，随着散焦量向两侧的增加，散焦像斑越来越大，图像越来越模糊，同时两侧的散焦像斑存在差异。3维点扩散函数的3维空间坐标系z轴与显微镜光轴重合，原点为焦平面中心点。图3为对应实验3-D-PSF能量分布曲线(纵坐标表示图像进行归一化后像素大小，图例数值为对应的散焦量δz)。其中图3a为对应实验3-D-PSF沿着z轴(x=y=0)的能量分布，图3b为对应实验3-D-PSF沿着径向(x轴)方向的能量分布。从图3中曲线可以看出：(1)在光轴方向上，3-D-PSF在z=0的焦面处能量最大，证实了3-D-PSF的能量主要集中在双锥体中部的锥顶附近区域，离开焦面后能量迅速衰减；(2)在径向方向上，3-D-PSF越远离焦面其能量扩散范围越大。

Fig.2 Experimental 3-D-PSF

Fig.3 Energy distribution of experimental 3-D-PSFa—along z axis b—along radial direction

2.2.2 多次采集平均 图像采集采用的是非致冷CCD传感器，在采集过程中产生的噪声影响不可忽略。随着散焦量的逐步增大，散焦像的亮度会迅速降低，散焦2-D-PSF的分布状况受影响的程度迅速加大，进而影响所构建的3-D-PSF的准确性。为此，本实验中采用多图像平均法，对荧光微珠多次图像采集以进行图像平均运算，降低噪声对所构建的3-D-PSF的影响。

本实验中对荧光微珠进行4次图像采集。图像采集时，在同一实验条件下，以荧光微珠中心区域设置为固定的截面定位点，每次图像采集均以此截面定位点为基准，沿光轴从距离荧光微珠的上部12μm至下部12μm以0.2μm的步进距离进行微珠图像采集，获取微珠不同散焦量序列切片图像。之后对4次采集的同一散焦量的4幅切片图像进行图像平均运算，平均后的图像作为构建3-D-PSF的序列2D-PSF。

图4即为4次采集实验中3-D-PSF在每个平面中心点的能量分布情况(纵坐标表示图像实际像素值，图例T1～T4分别表示连续进行4次采集对应的实验3-D-PSF)。图4b为图4a中曲线顶点部分局部放大图，y轴为对应采集图像像素大小。通过曲线可知，在微珠图像采集过程中，随着光照时间增加，荧光微珠受激产生的能量逐渐降低。为保证实验结果的准确性，在挑选荧光微珠时，尽可能挑选能量较高、其处于焦面时的中心能量尽量与CCD最大量化值相接近，同时尽量缩短单次图像采集曝光时间，以降低荧光微珠能量的衰减程度。

Fig.4 Energy distribution of PSF made by 4 experimentsa—energy distribution of center point b—partial enlargement

图5是利用以上4次采集图像通过多图像平均法得出的不同散焦量的微珠序列切片图像，即2-D-PSF。对比图2单次采集的3-D-PSF切片图像可以发现，利用多次采样叠加平均法能够在一定程度上降低显微镜系统噪声对点扩散函数的影响。

Fig.5 Experimental 3-D-PSF after several images average

3 图像复原

3.1 理论点扩散函数制作

(1)式的光学系统光学传输函数是理论点扩散函数的计算公式，结合显微镜系统相关参量、物镜大小、CCD大小，以及显微镜成像特性及CCD图像的数字化原理，可以通过该式计算获得一组任意散焦量间隔的序列理论2-D-PSF。图6所示为通过理论计算得到的一系列不同散焦量的2-D-PSF图像。图7为由该序列理论2-D-PSF构建的理论3-D-PSF能量分布图(纵坐标表示图像进行归一化后像素大小，图例数值为对应的散焦量δz)。其中图7a为理论3-D-PSF能量沿轴向分布，图7b为理论3-D-PSF能量沿径向的分布。由图7与图3中曲线比较可看出，实验3-D-PSF与理论3-D-PSF能量都是主要集中在3-D-PSF双椎体结构的椎顶部分，能量分布沿轴向和径向方向基本一致，而在局部区域有所不同。

Fig.6 Theoretical 3-D-PSF

Fig.7 Energy distribution of theoretical 3-D-PSFa—along z axis b—along radial direction

3.2 图像复原分析

实验中荧光微珠和生物荧光组织3维切片间距均为0.2μm；采用3-D-PSF空间大小均为101pixel×101pixel×101pixel，3-D-PSF在z轴间距亦为0.2μm；图像复原方法采用迭代盲去卷积(iterative blind deconvolution，IBD)算法[16]，初始3-D-PSF分别为理论3-D-PSF和实验3-D-PSF，迭代次数50次。

(5)

(6)

实验1:对荧光微珠3维切片图像的复原。在相同的条件下，分别利用实验3-D-PSF和理论3-D-PSF对荧光微珠3维显微镜图像进行去卷积图像复原处理。图8为对荧光微珠3维显微图像复原的效果。其中图8a为采集的荧光微珠成像的原始3维重构光轴剖面(y-z面)图，图8b、图8c分别为利用理论3-D-PSF和实验3-D-PSF复原的微珠图像3维重构光轴剖面图。

Fig.8 Image restoration of microsphere

a—microsphere image b—theoretical 3-D-PSF restoration result c—experimental 3-D-PSF restoration result

由图8可知，由于显微镜系统散焦的作用，荧光微珠通过显微镜成像后在3维空间上形成双椎体结构，能量集中在椎体附近。采用两种点扩散函数对荧光微珠3维图像的复原均获得良好的效果，在3维空间复原为一个微珠(亮点)。同时复原图像的上下两端均存在一些杂散信息，主要是由于在去卷积过程中，边界部分无法确定造成的。但是对比发现，两者的复原效果不完全相同，利用实验3-D-PSF复原效果要比理论3-D-PSF复原效果要好。利用理论3-D-PSF复原的图像虽然在空间上也形成了一个亮点，但是在旁边还存在干扰，且亮度较小，而利用实验3-D-PSF复原的图像在空间中形成的亮点亮度较大，与原始亮度较为接近，且周围干扰较小。

实验2:对生物荧光组织3维切片显微镜图像的复原。在相同的条件下，利用显微镜对荧光生物组织进行3维切片的采集，然后利用实验3-D-PSF和理论3-D-PSF对荧光生物组织3维切片图像进行去卷积图像复原处理。

图9中显示了4组不同切片图像的复原效果。图9a为采集的荧光组织原始切片图像，图9b为对应利用理论3-D-PSF复原效果，图9c为利用实验3-D-PSF复原效果。图10是荧光组织的图像复原3维重构效果对比图。其中图10a为荧光组织原始3维重构显示图，图10b、图10c分别为利用理论3-D-PSF和实验3-D-PSF复原的3维重构显示图。

Fig.9 Image restoration of fluorescent tissue

a—microscope images collected by microscope b—theoretical 3-D-PSF image restoration c—experimental 3-D-PSF image restoration

Fig.10 Image restoration of 3-D fluorescent tissue

a—original fluorescent tissue b—theoretical 3-D-PSF restoration result c—experimental 3-D-PSF restoration result

由图9和图10可知，在原始样本中，由于相邻层散焦层对焦平面层之间的相互干扰，采集的切片图像即包含该焦平面图像信息，又含有相邻散焦面的信息，使得采集图像模糊，3维重构图像中无法看清内部结构。而通过实验3-D-PSF和理论3-D-PSF对图像进行复原时，能较好去除相邻散焦面对图像的影响，处理后图像在每层中保留了该焦平面的主要信息，均获得良好的复原效果，并很好地实现了生物组织的3维重构。另外，实验3-D-PSF和理论3-D-PSF两者的复原效果不完全相同，仔细观察，实验3-D-PSF的复原效果更好。

4 结 论

利用荧光微珠实际测量系统的3-D-PSF，可以弥补理论3-D-PSF无法表征实际光学系统在生产制造过程中存在不可避免的偏差问题和噪声因素的缺陷，使得实验3-D-PSF更符合该光学系统的光学传递特性。本文中以荧光微珠模拟点光源，通过数字共焦显微镜对荧光微珠的不同截面和不同散焦量的切片图像进行采集，通过4次采集得到相同的序列荧光微珠图像，采用多图像叠加平均法进行处理，降低噪声的影响，构建显微镜光学系统的实验3-D-PSF。以该实验3-D-PSF与理论3-D-PSF对荧光微珠切片图像以及荧光生物组织3维切片图像进行去卷积复原处理和对比分析。实验结果表明，实验3-D-PSF与理论3-D-PSF均获得良好的复原效果，而实验3-D-PSF具有更好的复原精度。