光场的光强与相位信息具有重要意义,光强信息可以确定物体的表面形态,而相位信息用来确定物体的表面轮廓甚至内部结构,光场的光强及相位信息在进行生物医学成像[1-2]、光学成像测量及矫正[3-4]、三维成像及重建[5]等领域至关重要。然而光的相位信息不能直接观测,1949年,Gabor提出了第一个定量技术,利用干涉信息恢复复杂的光场[6],奠定了现代全息术的基础。目前,相位成像的检测手段多是基于干涉理论,将参考光束与待测光束进行干涉,利用干涉条纹信息解调出待测场相位信息。虽然它们在进行精确相位测量方面非常强大,但它们对环境扰动(如机械振动和温度变化)的高度敏感性和对参考光束的需求阻碍了其在便携式和紧凑干涉成像系统中的发展[7]。
近年来,单像素成像不仅可以对光强成像,而且可以实现对相位成像[8-9]。单像素成像系统主要由光学掩模与光电探测器构成。物体的光场信息经过光学掩模调制,通过透镜汇聚于无分辨率的光电探测器中,每张掩模对应一个信号,光场信息经过一系列信号调制后,可根据每张掩模与其对应信号的相关性求解得出。单像素成像系统可用于宽光谱成像,只要改变光电探测器的光谱范围,就可以很容易地改变成像的光谱范围。这一特性对于通常需要昂贵相机的光谱范围内的成像非常方便,特别是红外成像[10-12]和太赫兹成像[13-14]。单像素探测器的灵敏度高于CCD阵列面,可适用于微弱信号的成像[12]。正是由于单像素系统这些不可代替的优点,引起了科研人员的广泛关注。
本文利用空间光调制器加载相位型光学掩模对光场信息进行调制,用单像素探测器对光场极大值进行光强探测,利用相位提升(phaselift)算法[15]对探测到的光强信息及其对应的光学掩模信息进行求解,从而得到待求光场的复振幅。实验中利用该单像素复振幅成像系统分别对衍射环和透明的薄膜物体进行成像,其中物体的刻线与非刻线的相位差为0.053,刻线宽度为220 μm,实验得到图像相位差为0.046和刻线宽度为256 μm,图像信息与传统的检测手段得到的信息非常接近,该成像系统无需参考光束,便于成像系统的集成,促进便携式成像系统发展。
1 成像原理单像素复振幅成像原理如图1所示,物体信息被一系列复振幅光学掩模调制,调制后的光束经透镜聚焦,针孔位于透镜焦点处,经过滤波后的光强信息被无分辨率的单像素探测器接收。
|
图 1 系统原理图 Figure 1 Schematic diagram of the system |
为简单起见,物体的光场信息用一维向量X={x1, x2, x3, …, xN}表示,向量总长度为N,第n个像素的光场信息为xn,
| ${x_n} ={a_n}\exp ( - \rm{i}{\varphi _n})$ | (1) |
式中:
| ${\varPhi _m}_n ={b_n}\exp ( - {\rm{i}}{\theta _n}){\rm{ }}$ | (2) |
式中:bn为Φmn的振幅;
| ${V_{mn}} ={{X}}{\varPhi _{mn}} ={\varPhi _{mn}}{a_n}\exp ( - {\rm{i}}{\varphi _n}){\rm{ }}$ | (3) |
当光场经过透镜聚焦,在聚焦平面的主极大强度为
| ${y_m} =| < {{X}}, {\varPhi _m} > {|^2}$ | (4) |
光电探测器可以测量出每张掩模对应的光场强度ym,所以相位恢复问题即已知ym和Φm,求解未知向量X,根据Candés等[15]提出的相位提升原理可将式(4)进行求解。相位提升的求解基本思想是将式(4)中待求解的复向量X问题根据Z=XX*提升为一个秩为1的矩阵Z的求解问题,使得一般的非凸的相位检索问题转化为满足二次约束的凸规划问题,然后利用凸规划问题的求解算出X最优解[16],从而求得光场的相位信息及光强信息。
2 实验结果和分析 2.1 实验平台及过程单像素复振幅成像系统的实验装置图如图2所示,该实验装置由扩束系统,调制系统与探测系统构成。He–Ne激光器发出波长632.9 nm的光束经过扩束系统扩束,扩束系统由焦距f=25 mm的物镜,针孔P1和焦距f=125 mm的透镜L1构成,对直径为0.5 mm左右的激光光束进行5倍扩束,然后光束再经过物体。物体的大小为2 cm×2 cm,其中有效成像区域为2.048 mm×2.048 mm。通过待成像物体后的物光光束经过4f系统(f=100 mm),4f系统由两个相同焦距的透镜L2和L3与针孔P2构成,P2的作用是去除光路中的杂散光束,提高光束的成像质量。实验中采用的是Jasper Display 公司的JD8554N空间光调制器SLM,4f系统将物件的光强信息与相位信息成像在空间光调制器的列阵面上,该空间光调制器通过晶体的变化反射出指定相位信息的光束。光束调制后经过分束镜BS反射进入聚焦透镜L4,该聚焦透镜的焦距f=200 mm,因此聚焦焦斑直径大小d=2fλ/D≈10 μm,其中D为光束直径。该透镜的焦距选择较长的,以便增大聚焦斑,便于针孔P3选择聚焦光斑的极大值。在焦面上的针孔大小为10 μm,单像素探测器(PDA10A-EC, THORLABS)位于针孔后方,探测经过针孔后的傅里叶零级光强。
|
图 2 实验装置图 Figure 2 Experimental setup |
光学掩模的选择对于单像素成像至关重要,一个适合的光学掩模可以缩短单像素成像系统的图像恢复时间得到更高的图像质量。随机高斯矩阵构成的掩模曾被证明只要对光场的调制次数大于NlogN次(N为待恢复光场矩阵的长度),就可以通过phaselift算法恢复光场的光强与相位信息[17]。本实验采用随机高斯矩阵作为光学掩模,实验中掩模的像素数目为32×32,其中每个像素的尺寸为64 μm。首先构造大小为32×32的矩阵Φ,矩阵Φ中每一个元素独立地服从均值为0,方差为1/32的高斯分布,即:
| $ {{{\varPhi}} _{i,j}}\sim N\left( {0,\frac{1}{{32}}} \right) $ | (5) |
| $ \varPhi '{\rm{ = }}\frac{{\varPhi {\rm{ - min}}\left( \varPhi \right)}}{{\max \left( \varPhi \right)}} \times {\rm{255}} $ | (6) |
然后将根据式(6)将矩阵Φ中的数值扩大到0~255,式(6)中min(Φ)表示矩阵Φ中值最小的元素,max(Φ)表示矩阵Φ中值最大的元素,首先将矩阵Φ归一化,然后将每个元素乘以255扩大到0~255范围,因此Φ' 矩阵就是一个光学调制掩模矩阵Φm。根据上述原理生成1024张光学调制掩模,将其加载到空间光调制器中之后,光学掩模中的0~255就对应于相位的0~2π。实验中采用的是相位调制的方式,光场的光强信息不进行调制,相位信息由空间光调制器加载的掩模进行调制。
2.2 实验数据及分析当待测光束为高斯光束经过针孔发生衍射时,利用单像素成像系统恢复的结果如图3所示,图3(a)与(c)为高斯光束衍射后的光强图,图3(b)与(d)为高斯光束衍射后的相位图,从图中可以明显看到中心原斑与第一衍射环。由于调制区域大小为2 mm左右,无法观测到更高级次的衍射环。图3(b)与(d)中得到的相位图像在第二衍射环外虽有噪点,但是图像轮廓可以清晰看出。在实验中所得图像均为满采样得到的图像,即测量次数为1024次,当进行欠采样重构时,得到的图像质量并不好,造成这种现象可能的原因是实际实验中,对傅里叶频谱零级的提取并不精确,并且由于环境光及光路中反射所造成的测量误差是无法避免的,导致单像素探测器测量的信号误差较大,得到的欠采样重构图像质量较差。
|
图 3 衍射环图像 Figure 3 The image of diffraction ring |
当目标物体换成相位物时,得到的复振幅图像如图4所示。此时的物体为聚二甲基硅氧烷(PDMS)薄膜,该薄膜具有很好地透光性,多用于纳米压印、软光刻技术等[18]。待观测的薄膜厚度约为1 mm,在薄膜上方为一个印有“上海”字的商标,如图4(a)所示。本次实验对此商标的部分进行成像,如图4(b)所示,成像范围为2.048 mm×2.048 mm。由于该物体的制作是胶体烘干,所以刻线并不绝对平整,利用Phaseview光束质量分析仪对该物体进行测量,刻线的线条宽度约为220 μm,刻线与非刻线部分的相位差平均约为0.053。
|
图 4 目标物体图像 Figure 4 The image of object |
将图4恢复的目标物体图像与光束分析仪得到的结果进行对比。图4(c)和图4(d)为恢复的图像,恢复图像的分辨率为32×32。从图4(c)中可以看出线条平均宽度为4个像素,尺寸为256 μm,与光束分析仪测量所得结果相差少于一个掩模的像素宽度。实验中采用的图像分辨率为32×32,测量误差为1个像素,对应的尺寸为64 μm,所以提高图像的测量精度,提高光学掩模的分辨率是至关重要的。对图4(d)物体的相位图中的四个区域的相位求取平均值,得到的数据如表1所示,从表1中得出的刻线与非刻线的相位差平均为0.046,与光束分析仪测量所得结果接近。我们的单像素复振幅成像系统为光场的复振幅成像提供了一种有效的测量方式,该方式无需参考光束,为便携式、可集成式成像系统提供了参考。
|
|
表 1 各区域相位的平均值结果 Table 1 Average results for the regions |
单像素复振幅成像系统可通过相位型光学掩模调制的方式实现光场的复振幅成像。实验中分别对衍射光场及相位物体进行成像,实验得到的衍射光场的光强与相位图像可以清晰的看到光场的光强信息与相位信息;对相位差为0.053,刻线宽度为220 μm的物体成像,实验得到图像相位差为0.046和刻线宽度为256 μm,图像信息与传统的检测手段得到的信息非常接近。单像素复振幅成像系统图像的精度受光学掩模分辨率的控制,虽然目前得到的图像质量无法达到基于干涉理论得到的图像质量,但是该系统的成像光路简单不需要参考光束,便于系统集成,促进便携式成像系统的发展,而且该成像方法可应用于宽光谱成像,理论上可减少采样次数,提高成像效率,在光场复振幅成像方面具有较大发展潜力。
| [1] | ZERNIKE F. How I discovered phase contrast[J]. Science, 1955, 121(3141): 345–349. DOI:10.1126/science.121.3141.345 |
| [2] | POPESCU G. The power of imaging with phase, not power[J]. Physics Today, 2017, 70(5): 34–40. DOI:10.1063/PT.3.3553 |
| [3] | LIANG J Z, GRIMM B, GOELZ S, et al. Objective measurement of wave aberrations of the human eye with the use of a Hartmann-Shack wave-front sensor[J]. Journal of the Optical Society of America A, 1994, 11(7): 1949–1957. DOI:10.1364/JOSAA.11.001949 |
| [4] | ARTAL P. Optics of the eye and its impact in vision: a tutorial[J]. Advances in Optics and Photonics, 2014, 6(3): 340–367. DOI:10.1364/AOP.6.000340 |
| [5] | NEHMETALLAH G, BANERJEE P P. Applications of digital and analog holography in three-dimensional imaging[J]. Advances in Optics and Photonics, 2012, 4(4): 472–553. DOI:10.1364/AOP.4.000472 |
| [6] | BRAGG W L. Microscopy by reconstructed wave-fronts[J]. Nature, 1950, 166(4218): 399–400. |
| [7] | SOLDEVILA F, DURÁN V, CLEMENTE P, et al. Phase imaging by spatial wavefront sampling[J]. Optica, 2018, 5(2): 164–174. DOI:10.1364/OPTICA.5.000164 |
| [8] | ZHANG Z B, MA X, ZHONG J G. Single-pixel imaging by means of Fourier spectrum acquisition[J]. Nature Communications, 2015, 6: 6225. DOI:10.1038/ncomms7225 |
| [9] | OTA K, HAYASAKI Y. Complex-amplitude single-pixel imaging[J]. Optics Letters, 2018, 43(15): 3682–3685. DOI:10.1364/OL.43.003682 |
| [10] | ASPDEN R S, GEMMELL N R, MORRIS P A, et al. Photon-sparse microscopy: visible light imaging using infrared illumination[J]. Optica, 2015, 2(12): 1049–1052. DOI:10.1364/OPTICA.2.001049 |
| [11] | RADWELL N, MITCHELL K J, GIBSON G M, et al. Single-pixel infrared and visible microscope[J]. Optica, 2014, 1(5): 285–289. DOI:10.1364/OPTICA.1.000285 |
| [12] | 霍娟, 李明飞, 杨然, 等. 基于单像素探测器的高灵敏度近红外成像系统[J]. 红外与激光工程, 2016, 45(S1): S104001. |
| [13] | CHAN W L, CHARAN K, TAKHAR D, et al. A single-pixel terahertz imaging system based on compressed sensing[J]. Applied Physics Letters, 2008, 93(12): 121105. DOI:10.1063/1.2989126 |
| [14] | STANTCHEV R I, SUN B Q, HORNETT S M, et al. Noninvasive, near-field terahertz imaging of hidden objects using a single-pixel detector[J]. Science Advances, 2016, 2(6): e1600190. DOI:10.1126/sciadv.1600190 |
| [15] | CANDÈS E J, STROHMER T, VORONINSKI V. PhaseLift: exact and stable signal recovery from magnitude measurements via convex programming[J]. Communications on Pure and Applied Mathematics, 2013, 66(8): 1241–1274. DOI:10.1002/cpa.21432 |
| [16] | 张芬. 基于低秩矩阵填充的相位检索方法研究[D]. 合肥: 安徽大学, 2015. |
| [17] | CANDÈS E J, LI X D, SOLTANOLKOTABI M. Phase retrieval from coded diffraction patterns[J]. Applied and Computational Harmonic Analysis, 2015, 39(2): 277–299. |
| [18] | 王晨光, 杨江涛, 康宁, 等. PDMS梯度光栅结构制备技术研究[J]. 光谱学与光谱分析, 2015, 35(12): 3529–3533. |
2020, Vol. 42
Issue (5): 57-62

