CN105467806B - 单像素全息相机 - Google Patents

Abstract

本发明涉及一种单像素全息相机，属于全息成像领域。该系统的工作过程包括记录阶段和重建阶段。记录阶段是利用空间振幅调制装置的变换矩阵来对物体的菲涅尔衍射全息图的傅里叶频谱信息进行编码；由单像素强度检测模块完成时变信号的分析和记录；通过多次改变空间振幅调制装置的变换矩阵并进行相应的测量，最终得到物体菲涅尔衍射全息图的傅里叶频谱的压缩感知成像结果。重建阶段则是根据所选择的编码矩阵利用最小化l₁范数算法进行频谱信息的复原，再经过逆傅里叶变换重建出物体的菲涅尔衍射全息图。

Description

单像素全息相机

技术领域

本发明涉及一种单像素全息相机，属于全息成像领域。

背景技术

传统的全息成像技术是利用分束镜，将单色激光源发出的光成为物光和参考光两路，物光照射物体以后的散射光与参考光合束发生干涉现象，由全息干板对干涉图样进行记录。重建过程则是利用与参考光束相同的光照射全息干板，在干板后方得到重现的物光波。这种方法能够完全记录物体的菲涅尔衍射全息图的复振幅，但是由于记录时直流项的存在，使得重建过程中叠加有较大的背景光强；同时如果是同轴全息，则在其重建过程中还会产生共轭像。这些因素都会影响所重建物光波的观测质量，而且使用全息干板也不便于信息的存储和传输。后来产生了数字全息术，使用探测器阵列来取代全息干板，使得全息信息能够被直接记录在计算机中，极大地提高了信息存储与传输的方便性。但这种方法也有其缺陷，就是大规模探测器阵列造价比较昂贵，而且在某些情况下(例如非可见光波段)，探测器阵列有可能无法制造或者成本过高；此外，在航天领域等对成像器件的体积和质量有严格要求的应用场合，大规模探测器阵列可能无法满足应用要求。

近年来，Candes、Romberg、Tao和Tonoho等人提出了压缩感知理论(CompressiveSensing)，它指出如果信号在某个变换域是稀疏的，那么就可以通过一个不相关的观测矩阵将高维信号投影到低维空间中，通过求解最优化算法，以高概率重构原始信号。

2006年，美国RICE大学的研究人员根据压缩感知理论设计出了单像素相机，通过数字微镜阵列(DMD)的调制矩阵对含有图像信息的光场强度进行编码，然后通过透镜将所有的光收集起来，由单像素光电检测器件完成光电转换，利用最优化算法根据记录结果重建出原始图像。此单像素相机只能对光场的强度信息进行记录，而直接忽略掉了光场的相位信息，换句话说，这样记录的图像不含有空间坐标信息的。

发明内容

为了能够在较低探测器成本下，快速获得不含有背景光强以及共轭像的物体菲涅尔衍射全息图，并将全息图直接记录在计算机中，本发明提出一种单像素全息相机。

本发明是利用空间振幅调制装置的变换矩阵，对物体菲涅尔衍射全息图的傅里叶频谱信息进行编码，由单像素强度检测模块完成时变信号的分析和记录；通过多次改变空间振幅调制装置的变换矩阵并进行相应的测量，最终得到物体菲涅尔衍射全息图的傅里叶频谱的压缩感知成像结果。

单像素全息相机整个装置包括激光器1、半波片2、第一宽带分光棱镜3、第一扩束镜4、偏振分光棱镜5、空间振幅调制装置6、第一透镜7、第二宽带分光棱镜8、声光移频器9、光阑10、第二扩束镜11、反射镜12、物体13、第二透镜14、单像素强度探测模块15、计算机16。

单像素全息相机主体结构采用马赫泽德干涉仪的形式，分成上下两臂光路；其中第一宽带分光棱镜3、声光移频器9、光阑10、第二扩束镜11、反射镜12、第二宽带分光棱镜8组成上臂光路；而第一宽带分光棱镜3、第一扩束镜4、偏振分光棱镜5、空间振幅调制装置6、第一透镜7、第二宽带分光棱镜8则组成下臂光路；上臂光路中放置的光阑10必须恰好让声光移频器9所输出的一级衍射光通过；在下臂光路中必须使空间振幅调制装置6和第二宽带分光棱镜8分别位于第一透镜7的前后焦面。

激光器1发出的激光束，经过一个半波片2入射到主体结构中，从主体结构出射的光则照射在物体13上；由第二透镜14完成对经过物体13的所有透射光的汇聚收集，然后由单像素强度探测模块15对总光强进行光电转换并进行信号的提取和分析，最后送往计算机16 完成记录。

单像素强度检测模块15由光敏器件、带通滤波器、正余弦信号发生器、乘法器、低通滤波器、AD转换器件和PCI数据采集卡构成。

单像素全息相机的工作过程包括记录和重建两个阶段，具体实施方法包括以下步骤：

激光器1发出的频率为ω₀的光束经第一宽带分光棱镜3分束后分成上下臂两路光；上臂一路平行光束经过声光移频器9调制后产生频率偏移，从而使上下两路光波有差频Ω；对于下臂一路，平行光束经过偏振分光棱镜5反射到空间振幅调制装置6后，把紧贴其表面的光场复振幅表示为m(x,y)。

因为空间振幅调制装置6和第二宽带分光棱镜8分别位于第一透镜7的前后焦面，则第二宽带分光第二宽带分光棱镜8处的光场与前焦面光场满足傅里叶变换：

其中x,y,z为空间坐标，k_x,k_y代表空间频域坐标，代表傅里叶变换。

该光波从第二宽带分光棱镜8处经过距离为z的菲涅尔衍射，由菲涅尔衍射公式可知到达物体13前表面的光场为：

w₂(x,y；z)＝w₁(x,y；z)*h(x,y；z) (2)

其中*代表卷积运算，h(x,y；z)代表菲涅尔衍射过程的点扩散函数：

上下两路光合束后投射到物体上，由第二透镜14完成对所有经过物体13的透射光的汇聚收集，并由光电二极管完成光电转换。经过一个中心频率为Ω的带通滤波器并且忽略掉一些常数因子后电流变为：

其中Re[·]表示复数的实部，v(x,y)代表物体的透过率分布，且：

(4)式中所表示的电流被分成两路，分别通过由乘法器和低通滤波器组成的两路正交锁相检测电路，得到同时包含振幅和相位信息的两个电流分量和将这两个值记录到计算机中之后，可以通过复数的加法得到：

通过这种方式可以将(4)式所表示的复信号进行完整记录。

考虑到物体函数相应的空间频谱范围是有限宽度的，再加上光学系统的有限孔径导致频谱宽度进一步受限，在此基础之上对(4)式进行离散化处理可得到：

其中W₂(k_x,k_y；z),V(k_x,k_y；z),H(k_x,k_y；z)分别为w₂(x,y；z),|v(x,y；z)|²,h(x,y；z)的傅里叶变换；Δ为积分区间，它代表了空间频谱范围；k_i,k_j是空间频域坐标k_x,k_y的离散化；则是空间振幅调制装置6出射复振幅场m(x,y)的离散化表示；需要重点说明的是L(k_i,k_j；z)，其值可表示为(7)式，它是物体强度分布的傅里叶变换和菲涅尔衍射光学传递函数的乘积，其物理意义为物体经过距离为z的菲涅尔衍射后复光场的空间频谱分布。

L(k_i,k_j；z)＝H(k_x,k_y；z)V(k_x,k_y；z) (7)

对二维数组L(k_i,k_j；z)进行向量化处理得到N²×1的列向量m,L，则(6)式可进一步表示为

把上式中的当做一次的观测结果，通过相继改变空间光调制器的振幅调制函数为 m₁,m₂,...m_M进行M(M＜N²)次测量，即可得到M×1维的压缩感知成像结果：

利用最小化l₁范数算法进行频谱信息的复原，再经过逆傅里叶变换重建出物体的菲涅尔衍射全息图。

与现有技术相比，本发明具有以下有益效果：

(1)本发明使用时变光场对物体进行照明，差频信号的振幅和相位就表征了透射物波的空间振幅和相位信息，即包含了物体的横向透过强度与纵向深度信息；通过解析该信号对复光场的振幅和相位都进行记录。(2)本发明是一种数字全息术，菲涅尔衍射全息图存储在计算机中，方便进行后期处理以及传输。(3)根据本发明所获得的全息图进行光场重建时不含有背景光强以及共轭像，能够有更好的观测效果。(4)本发明采用单像素探测器进行记录，大大节约了成本，或者可以使用灵敏度更高的探头来获得更高质量的全息图；同时由于其体积小质量轻，能够满足一些特殊应用场合的要求。(5)本发明得到的是压缩感知记录结果，稀疏变换基的选择关系到重建质量。现有傅里叶变换基、小波变换基等已经可以重建出较好的图像，如果能够找到更好的稀疏变换基，那么就能获得更好的重建图像；(6)本发明的重建图像质量和记录次数有关，并且各次记录之间相互独立。也就是说，在需要快速得到图样的情况下，可以记录较少次数重建出全息图大致场景，然后在空闲时继续进行记录，逐步获得具有更加清晰细节的全息图，这是一个渐进的过程。这样，前期记录的数据在后期可以作为感知记录结果的一部分直接使用，而不需要重新记录，节约了大量时间。

附图说明

图1是单像素全息相机基本结构图。

图2是本发明—实施例中待测物体示意图。

图3(a)、3(b)是本发明—实施例中待测物体的理想菲涅尔衍射图样的实、虚部。

图4(a)、4(b)是本发明—实施例中待测物体的单像素全息相机压缩感知成像结果的实、虚部。

图5是本发明—实施例中待测物体的单像素全息相机最终成像结果。(a)、(b)分别为采样数是完全采样的15％时重建的菲涅尔衍射全息图的实、虚部；(c)、(d)分别为采样数是完全采样的33％建的菲涅尔衍射全息图的实、虚部。

具体实施方式

下面结合附图与实施例对本发明作进一步说明。

单像素全息相机的工作过程包括记录阶段和重建阶段；记录阶段是利用空间振幅调制装置的变换矩阵，来对物体的菲涅尔衍射全息图的傅里叶频谱信息进行编码；由单像素强度检测模块完成时变信号的分析和记录；通过多次改变空间振幅调制装置的变换矩阵并进行相应的测量，最终得到物体菲涅尔衍射全息图的傅里叶频谱的压缩感知成像结果。重建阶段则是根据所选择的编码矩阵利用最小化l₁范数算法进行频谱信息的复原，再经过逆傅里叶变换重建出物体的菲涅尔衍射全息图。

如图1所示，单像素全息相机记录阶段装置，激光器1发出的激光束经过一个半波片 2后入射到第一宽带分光棱镜3中，被第一宽带分光棱镜3分为透射和反射两束光；其中反射光紧接着通过第一扩束镜4进行扩束后入射到偏振分光棱镜5；偏振分光棱镜5的反射线偏振光入射到空间振幅调制装置6；空间振幅调制装置6加载合适的随机变换矩阵完成对光场的编码；激光束经过空间振幅调制装置6的反射后通过第一透镜7并入射到第二宽带分光棱镜8上，这里的空间振幅调制装置6和第二宽带分光棱镜8分别位于第一透镜7的前后焦面；第一宽带分光棱镜3的透射光入射到声光移频器9中，声光移频器9的出射光含有多级衍射光，我们利用光阑10只选择其中的一级衍射光；光阑10的出射光接着通过第二扩束镜 11进行扩束，然后再经过反射镜12反射到第二宽带分光棱镜8上，由第二宽带分光棱镜8 完成两束光的合束；合束以后的光照射在物体13上，由第二透镜14完成对所有经过物体13 的透射光的汇聚收集，接着单像素强度探测模块15将对总光强进行光电转换，并进行信号的提取和分析，然后送往计算机16完成记录。

本实施例中使用的是波长为λ＝532nm的固体激光器，空间振幅调制装置为振幅型空间光调制器，变换矩阵m(x,y)为灰度值满足高斯概率分布的随机矩阵，声光移频器的调制频率为Ω＝100kHz，所使用的物体如图2所示，尺寸为1cm×1cm，像素数为128×128。

单像素全息相机，记录阶段按照如下步骤实现：

固体激光器1发出的频率为ω₀的光束经第一宽带分光棱镜3分束后分成上下两路光。上面一路平行光束经过声光移频器调制，导致上下两束光产生频差Ω＝100kHz。对于下面一路，平行光束经过振幅型空间光调制器6调制后，其出射光场复振幅为高斯随机矩阵m(x,y)，经由第一透镜7完成傅里叶变换。

上下两束光在第二宽带分光棱镜8处进行合束，然后合束光对物体进行照明并透射。由于两束光之间存在远小于光波频率的差频，所以该照明光场强度是时变的，变化频率为Ω＝100kHz。这样相应差频信号的振幅和相位就表征了透射物波的空间振幅和相位信息，即包含了物体的横向透过强度与纵向深度信息。

第二透镜14对所有经过物体13的透射光进行汇聚收集，采用硅光电池作为光敏器件，将总光强转换为电信号；输出电信号i(t)经过一个中心频率为Ω的带通滤波器滤除直流分量、 50Hz工频信号和高频噪声等主要干扰后被分成了两路；这两路电信号分别经过由乘法器和低通滤波器组成的锁相检测电路，但它们所乘信号的相位必须相差90度；从数学的角度讲，这样处理的结果可以分别得到复信号的实部和虚部，也就同时获得振幅和相位两个未知量；从观测效果来看，此举可以成功地消除共轭像，改善成像质量；然后，通过AD模数转换器件，将两路电信号都存储到计算机中，方便进行处理和传输。

两路输出信号分别存储了复信号的实部和虚部，即通过这样的方式对该复信号的振幅和相位同时进行了记录，通过复数的加法可以还原下式所示复信号。

其中w₂(x,y；z)代表下面一路到达物体前表面的光场分布，v(x,y；z)为物体的透过率分布。

由于物体函数相应的空间频谱范围是有限宽度的，再加上光学系统的有限孔径导致频谱宽度进一步受限，在此基础之上对(10)式进行离散化处理可得到：

其中k_i,k_j是空间频域坐标k_x,k_y的离散化；L(k_i,k_j；z)是物体强度分布的傅里叶变换和菲涅尔衍射光学传递函数的乘积，表示物体经过距离为z的菲涅尔衍射后复光场的频谱分布。 m(x,y)是满足高斯概率分布函数的随机矩阵。

对二维数组L(k_i,k_j；z)进行向量化处理得到N²×1的列向量m,L，那么所记录的复信号进一步可表示为：

把上式中的当做一次的观测结果，通过相继改变振幅型空间光调制器的振幅调制函数为m₁,m₂,...m_M进行M(M＜N²)次测量，即可得到M×1维的压缩感知成像结果：

如图3所示，(a)、(b)分别为压缩感知成像结果的实部和虚部。

利用最小化l₁范数算法进行频谱信息的复原，再经过逆傅里叶变换重建出物体的菲涅尔衍射全息图；如图5所示是待测物体的单像素全息相机最终成像结果，其中(a)、(b)分别为采样数是完全采样的15％时重建的菲涅尔衍射全息图的实部和虚部，(c)、d)分别为采样数是完全采样的33％时重建的菲涅尔衍射全息图的实部和虚部。

Claims (2)

1.单像素全息相机，其特征在于：

所述单像素全息相机的组成结构：

整个装置包括激光器( 1) 、半波片( 2) 、第一宽带分光棱镜( 3) 、第一扩束镜( 4)、偏振分光棱镜( 5) 、空间振幅调制装置( 6) 、第一透镜( 7) 、第二宽带分光棱镜( 8) 、声光移频器( 9) 、光阑( 10) 、第二扩束镜( 11) 、反射镜( 12) 、物体( 13) 、第二透镜( 14) 、单像素强度探测模块( 15) 、计算机( 16) ；

主体结构采用马赫泽德干涉仪的形式，分成上下两臂光路；其中第一宽带分光棱镜(3) 、声光移频器( 9) 、光阑( 10) 、第二扩束镜( 11) 、反射镜( 12) 、第二宽带分光棱镜( 8) 组成上臂光路；而第一宽带分光棱镜( 3) 、第一扩束镜( 4) 、偏振分光棱镜( 5) 、空间振幅调制装置( 6) 、第一透镜( 7) 、第二宽带分光棱镜( 8) 则组成下臂光路；上臂光路中放置的光阑( 10) 必须恰好让声光移频器( 9) 所输出的一级衍射光通过；在下臂光路中必须使空间振幅调制装置( 6) 和第二宽带分光棱镜( 8) 分别位于第一透镜( 7)的前后焦面；

激光器( 1) 发出的激光束，经过一个半波片( 2) 入射到主体结构中，从主体结构出射的光则照射在物体( 13) 上；由第二透镜( 14) 完成对经过物体( 13) 的所有透射光的汇聚收集，然后由单像素强度探测模块( 15) 对总光强进行光电转换并进行信号的提取和分析，最后送往计算机( 16) 完成记录；

单像素强度检测模块( 15) 由光敏器件、带通滤波器、正余弦信号发生器、乘法器、低通滤波器、AD转换器件和PCI数据采集卡构成；

所述单像素全息相机的工作过程：

包括记录阶段和重建阶段；

记录阶段是利用空间振幅调制装置( 6) 的变换矩阵来对物体的菲涅尔衍射全息图的傅里叶频谱信息进行编码；由单像素强度检测模块( 15) 完成时变信号的分析和记录；通过多次改变空间振幅调制装置( 6) 的变换矩阵并进行相应的测量，最终得到物体菲涅尔衍射全息图的傅里叶频谱的压缩感知成像结果；

重建阶段则是根据所选择的变换矩阵利用最小化l₁范数算法进行频谱信息的复原，再经过逆傅里叶变换重建出物体的菲涅尔衍射全息图；

具体实施方法包括以下步骤，其中步骤1)至7)属于记录阶段，步骤8)属于重建阶段：

1)激光器( 1) 发出的频率为ω₀的光束经第一宽带分光棱镜( 3) 分束后分成上下臂两路光；上臂一路平行光束经过声光移频器( 9) 调制后产生频率偏移，从而使上下两路光波有差频Ω；对于下臂一路，平行光束经过空间振幅调制装置( 6) 调制后，其出射光场m(x,y)，经由第一透镜( 7) 完成傅里叶变换；

2)上下两臂光束在第二宽带分光棱镜( 8) 处进行合束，然后合束光对物体( 13) 进行照明并透射出去；由于两束光之间存在远小于光波频率的差频Ω，所以透射出来的照明光场强度是时变的，变化频率也为Ω；这样相应差频信号的振幅和相位就表征了透射物波的空间振幅和相位信息，即包含了物体( 13) 的横向透过强度与纵向深度信息；

3)第二透镜( 14) 对经过物体( 13) 的所有透射光进行汇聚收集，由光敏器件响应完成光电转换；其输出电信号i(t)经过一个中心频率为Ω的带通滤波器后分成了两路；两路电信号分别经过由乘法器和低通滤波器组成的锁相检测电路，但其中两个乘法器所乘信号的相位必须相差90度，以消除共轭像，改善成像质量；之后通过AD模数转换器件，将两路电信号都存储到计算机中，方便进行处理和传输；

4)实际上，记录的这两路输出信号分别代表了复信号的实部和虚部，通过复数加法还原下式所表示的复信号：

其中w₂(x，y；z)代表下臂光路到达物体( 13) 前表面的光场；v(x，y；z)为物体( 13) 的透过率分布；

5)由于物体函数相应的空间频谱范围是有限宽度的，再加上光学系统的有限孔径导致频谱宽度进一步受限，在此基础之上对上式进行离散化处理可得到：

其中W₂(k_x，k_y；z)，V(k_x，k_y；z)分别为w₂(x，y；z)和v(x，y；z)的傅里叶变换，H(k_x，k_y；z)表示菲涅尔衍射点扩散函数的傅里叶变换；k_x，k_y代表空间频域的频域坐标，k_i,k_j是k_x,k_y的离散化，N则为离散化点数；表示空间振幅调制装置( 6) 出射复振幅场m(x,y)的离散化表示；需要重点说明的是L(k_i,k_j；z)，其值可表示为H(k_x,k_y；z)V(k_x,k_y；z)，它是物体强度分布的傅里叶变换和菲涅尔衍射光学传递函数的乘积，其物理意义为物体经过距离为z的菲涅尔衍射后复光场的空间频谱分布；

6)对二维数组和L(k_i，k_j；z)进行向量化处理，分别得到N²×1维的列向量m和L，那么所记录的复信号进一步可表示为m与L的内积，即

7)把上式中的当做一次的观测结果，通过相继改变空间振幅调制装置( 6) 的振幅调制函数，即通过不同的变换矩阵m₁,m₂,...m_M进行M(M＜N²)次测量，从而得到M×1维的压缩感知成像结果：

其中上标表示维数；需要注意的是，此处的m表示的是M个N²×1维列向量所构成的M×N²矩阵；

8)利用最小化l₁范数算法复原频谱信息，再由傅里叶逆变换重建物体的菲涅尔全息图。

2.根据权利要求1所述的单像素全息相机，其中

空间振幅调制装置是对入射光振幅的空间分布进行调制，常用有振幅型液晶空间光调制器、数字微镜器件；

光敏器件是用来实现对入射光强的响应，常用有光电二极管、硅光电池。