本发明涉及彩色成像技术领域，特别涉及一种基于单像素探测器的彩色成像方法和系统。
背景技术：
彩色成像技术主要用于对同一观测目标的多个波段同时进行成像比对,以实现对其成分组成、分布状态等特性进行观测分析。相对于传统的单一波段成像系统,彩色成像系统能够获得更为全面的目标色彩信息。现有彩色成像方法大多采用面阵相机，且通过分光元件，如滤波片，光栅或者迈克尔逊双光束干涉光路等，将宽光谱光束分解成多个单波长的光束，成像装置比较复杂，而且需要多次测量，测量数据大，难以实现压缩采样。
单像素成像是近年来新兴的一种成像方式。该技术通过在空间光调制器(如dmd)上加载一组调制图案(如哈达码、傅里叶基底图案)，用于将物体的信息编码到变换域(如频域)，接着利用单像素探测器(如单个光电二极管)采集调制后的信号，最后通过相关算法(如傅里叶单像素算法)计算重建物体的图像。与传统的以ccd或者cmos相机为核心的成像器件相比，单像素探测器结构简单，且具有较宽的光谱响应范围和较强的光敏感度。此外，由于自然图像在变换域(如频域)中具有积聚特性，采用单像素成像方法可以实现压缩采样。
为了实现基于单像素探测器的彩色成像，文献“s.s.welsh,m.p.edgar,r.bowman,p.jonathan,b.sun,andm.j.padgett,"fastfull-colorcomputationalimagingwithsingle-pixeldetectors,"opt.express21,23068-23074(2013)”提出了一种基于三个单像素探测器的彩色成像系统，该方案通过在每个单像素探测器前分别加红、绿和蓝色滤光片，能够分别重建出被测物体的红、绿和蓝色图像，然后将三个重建结果组合，即可获得被测物体的彩色图像。然而，这种方法至少需要三个单像素探测器，系统复杂。文献“l.h.bian,j.l.suo,g.h.situ,z.w.li,j.t.fan,f.chen,andq.h.dai,"multispectralimagingusingasinglebucketdetector,"scientificreport6，24752(2016)”提出了一种基于一个单像素探测器的彩色成像系统，但该方案需要一个高速旋转的光栅将宽光谱光束分解成多个光谱频带的调制光，然后利用单像素探测器分别采集每个光谱的信号，测量次数比较多。此外，由于需要一个高速旋转的光栅，光栅、dmd和单像素探测器需要同步，系统比较复杂。文献“s.l.jin,w.w.wang,y.l.wang,k.c.huang,q.s.shi,c.f.ying,d.q.liu,q.ye,w.y.zhouandj.g.tian,hyperspectralimagingusingthesingle-pixelfouriertransformtechnique.scientificreports7,45209(2017)”提出了一种基于迈克尔逊双光束干涉分光的彩色单像素成像系统，该方案通过精确移动干涉光路中的反射镜，获得不同光谱反射的光强信号。由于需要采用干涉光路，测量系统对环境震动比较敏感。此外，该方案需要利用单像素探测器对每个光谱的反射光强进行测量，测量次数比较多。因此，如何提供一种高效、简单的彩色单像素成像方法和系统成为单像素成像领域中的一项技术难题、亟待解决。
技术实现要素：
本发明的目的在于克服现有技术中的缺点与不足，提供一种基于单像素探测器的彩色成像方法，能够使物体的轮廓信息和彩色信息分别积聚在傅里叶变换域的不同位置，实现压缩采用，并且图像重建过程无需迭代运算，重建算法简单，运算速度快。
本发明的另一目的在于提供基于单像素探测器的彩色成像方法的成像系统。
为实现以上目的，本发明采取如下技术方案：
一种基于单像素探测器的彩色成像方法，包括下述步骤：
s1、利用计算机产生一系列不同频率、不同相位的傅里叶基底灰度图案；
s2、将所述傅里叶基底灰度图案转换成基于拜尔滤波模板的傅里叶基底彩色图案；
s3、通过数字投影仪将所述傅里叶基底彩色图案依次投射到被测物体表面，同时利用一个单像素探测器采集经物体反射后的光强；
s4、利用傅里叶单像素成像算法处理所述的采集信号，得到物体的马赛克图像；
s5、采用空间色彩插值算法处理所述的马赛克图像，获得物体的彩色图像。
作为优选的技术方案，步骤s1具体为：
利用公式
作为优选的技术方案，步骤s2具体包括下述步骤：
s21、利用如下公式生成三个与傅里叶基底灰度图案像素数相等的拜尔滤波模板：
其中，mr、mb和mg分别表示红色、蓝色和绿色三种颜色的拜尔滤波模板；
s22、将步骤s21中的三个拜尔滤波模板分别乘以傅里叶基底灰度图案，得到三组傅里叶基底图案，分别对应红色、绿色和蓝色三种颜色：
其中，c表示颜色；
s23、将步骤s22得到的三组傅里叶基底图案组合成一个三维矩阵，即可得到一系列不同频率、不同相位的傅里叶基底彩色图案；由于采用拜尔滤波模板，傅里叶基底彩色图案中的每个像素只显示一种颜色。
作为优选的技术方案，步骤s3具体包括下述步骤：
s31、通过数字投影仪将一系列的傅里叶基底彩色图案依次投射到被测物体表面，同时利用一个单像素探测器采集经物体反射后的光场强度，并通过数据采集卡存储于计算机内；
s32、经物体反射后的光场总强度表示为：
s33、单像素探测器每次采集并存储的光强信号表示为：
作为优选的技术方案，步骤s4具体包括下述步骤：
s41、按照如下公式计算被测物体与频率为(fx,fy)对应的傅里叶频谱系数：
其中，j表示虚部单位，
s42、重复步骤s41的方法，计算被测物体与其他频率对应的傅里叶频谱系数，获得被测物体的傅里叶频谱分布图；
s43、对步骤s42得到的傅里叶频谱分布图求傅里叶反变换，然后对变换后的结果取模并进行归一化，即可得到被测物体的图像；由于投射到被测物体表面的图案为基于拜尔滤波模板的彩色傅里叶基底，重建结果中的每个像素只能获得红色、绿色和蓝色三种颜色中的一种颜色信息，其它两种颜色的信息被丢失，因此，该重建结果又称为被测物体的马赛克图像。
作为优选的技术方案，步骤s5具体为：
重建结果中任何一像素缺失颜色都用最临近的4个同颜色像素为基础，采用空间色彩插值算法获得；所述空间色彩插值算法包括：双线性插值算法、ronkimmel算法或optimalrecovery算法。
一种基于单像素探测器的彩色成像方法的成像系统，包括：傅里叶基底彩色图案生成模块、投影采集模块、以及数据处理模块；所述傅里叶基底彩色图案生成模块通过计算机产生傅里叶基底彩色图案，经数字投影仪将该傅里叶基底彩色图案投射到待测物体上，同时单像素探测器对经被测物体表面反射后的光场强度进行采集，并将采集到的光强信号通过数据采集与控制模块传输至数据处理模块，所述数据处理模块通过计算机算法处理光强信号，最后获得被测物体的彩色图像。
作为优选的技术方案，所述傅里叶基底彩色图案生成模块具体是利用计算机产生一系列不同频率、不同相位的傅里叶基底灰度图案，并将该基底图案转换成基于拜尔滤波模板的傅里叶基底彩色图案：
其中mr、mb和mg分别表示红色、蓝色和绿色三种颜色的拜尔滤波模板；d表示颜色；
作为优选的技术方案，所述单像素探测器采集经物体反射后的光场强度，并通过数据采集卡存储于计算机内；所述单像素探测器选用光敏二极管或光电倍增管。
作为优选的技术方案，所述数据处理模块通过傅里叶单像素成像算法处理所述的采集信号，得到被测物体的马赛克图像，再采用空间色彩插值算法处理所述的马赛克图像，获得被测物体的彩色图像；所述空间色彩插值算法包括：双线性插值算法、ronkimmel算法或optimalrecovery算法。
本发明相对于现有技术具有如下的优点和效果：
(1)本发明的彩色成像系统只需采用一个单像素探测器，系统简单；
(2)本发明的彩色成像方法采用物体的轮廓信息和彩色信息分别积聚在傅里叶变换域的不同位置，便于实现压缩采样；
(3)本发明的彩色成像方法重建过程不需要迭代运算，算法简单，运算速度快。
附图说明
图1为本发明基于单像素探测器的彩色成像方法流程图；
图2为本实施方式中，基于单像素探测器的彩色成像系统的软件组成示意图；
图3为本实施方式中，生成一张傅里叶基底彩色图案的局部放大图；
图4为本实施方式中，基于单像素探测器的彩色成像系统的示意图；附图标号：1、单像素探测器；2、数字投影仪；3、被测物体；
图5(a)为本实施方式中，利用基于单像素探测器的彩色成像方法和系统，采用压缩采样方式，获得被测物体的傅里叶频谱图；
图5(b)为本实施方式中，利用基于单像素探测器的彩色成像方法和系统，采用压缩采样方式，获得被测物体的彩色图像；
图5(c)为本实施方式中，利用基于单像素探测器的彩色成像方法和系统，采用全采样方式，获得被测物体的频谱分布图；
图5(d)为本实施方式中，利用基于单像素探测器的彩色成像方法和系统，采用全采样方式重建的彩色图像。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
实施例
如图1所示，一种基于单像素探测器的彩色成像方法，包括下述步骤：
s1、利用计算机产生一系列不同频率、不同相位的傅里叶基底灰度图案；
具体利用公式
s2、将所述傅里叶基底灰度图案转换成基于拜尔滤波模板的傅里叶基底彩色图案；具体包括下述步骤：
s21、利用如下公式生成三个与傅里叶基底灰度图案像素数相等的拜尔滤波模板：
其中，mr、mb和mg分别表示红色、蓝色和绿色三种颜色的拜尔滤波模板；
s22、将步骤s21中的三个拜尔滤波模板分别乘以傅里叶基底灰度图案，得到三组傅里叶基底图案，分别对应红色、绿色和蓝色三种颜色：
其中，c表示颜色；如图3所示，傅里叶基底灰度图案乘以一个拜尔滤波模板得到的傅里叶基底图案。
s23、将步骤s22得到的三组傅里叶基底图案组合成一个三维矩阵，即可得到一系列不同频率、不同相位的傅里叶基底彩色图案；由于采用拜尔滤波模板，傅里叶基底彩色图案中的每个像素只显示一种颜色。
s3、通过数字投影仪将所述傅里叶基底彩色图案依次投射到被测物体表面，同时利用一个单像素探测器采集经物体反射后的光强；具体包括下述步骤：
s31、通过数字投影仪将一系列的傅里叶基底彩色图案依次投射到被测物体表面，同时利用一个单像素探测器采集经物体反射后的光场强度，并通过数据采集卡存储于计算机内；
s32、经物体反射后的光场总强度表示为：
s33、单像素探测器每次采集存储的光强信号可表示为：
s4、利用傅里叶单像素成像算法处理所述的采集信号，得到物体的马赛克图像；具体包括下述步骤：
s41、按照如下公式计算被测物体与频率为(fx,fy)对应的傅里叶频谱系数：
其中，j表示虚部单位，
s42、重复步骤s41的方法，计算被测物体与其他频率对应的傅里叶频谱系数，获得被测物体的傅里叶频谱分布图；
s43、对步骤s42得到的傅里叶频谱分布图求傅里叶反变换，然后对变换后的结果取模并归一化，即可得到被测物体的图像；由于投射到被测物体表面的图案为基于拜尔滤波模板的彩色傅里叶基底，重建结果中的每个像素只能获得红色、绿色和蓝色三种颜色中的一种颜色信息，其它两种颜色的信息被丢失，因此，该重建结果又称为被测物体的马赛克图像。
s5、采用空间色彩插值算法处理所述的马赛克图像，获得物体的彩色图像；具体过程如下：
重建结果中任何一像素缺失颜色都用最临近的4个同颜色像素为基础，采用空间色彩插值算法获得；所述空间色彩插值算法包括：双线性插值算法、ronkimmel算法或optimalrecovery算法等。
如图2所示，本实施例的一种基于单像素探测器的彩色成像系统，包括傅里叶基底彩色图案生成模块、投影采集模块、以及数据处理模块；通过傅里叶基底彩色图案生成模块产生傅里叶基底彩色图案，经数字投影仪将该傅里叶基底彩色图案投射到待测物体上，同时单像素探测器对经被测物体表面反射后的光场强度进行采集，并将采集到的光强信号传输至数据处理模块，所述数据处理模块通过算法处理光强信号，最后获得被测物体的彩色图像。在本实施例中，所述单像素探测器选用光敏二极管或光电倍增管。
为了测试基于单像素探测器的彩色成像方法的有效性，对一彩色物体进行测量。搭建如图4所示的测量系统，包括单像素探测器1、数字投影仪2、以及被测物体3，首先利用计算机生成一组频率不同、相位不同的傅里叶基底图案，然后将该组图案转换成基于拜尔滤波模板的傅里叶基底彩色图案；接着通过数字投影仪向被测物体表面投射傅里叶基底彩色图案，同时利用一个单像素探测器采集反射后的光场强度，并存储于计算机内；随后采用傅里叶单像素算法处理采集的信号，获得被测物体的马赛克图像；最后利用空间色彩插值算法处理马赛克图像，获得被测物体的彩色图像。图5(c)为采用全采样方式，获得被测物体的完整傅里叶频谱分布图。从图5(c)可知，被测物体的轮廓信息和彩色信息分别积聚在傅里叶变换域的不同位置。图5(d)为对重建的频谱做反傅里叶变换，并通过空间色彩插值算法恢复的被测物体的彩色图像。
由于物体的轮廓信息和彩色信息分别积聚在傅里叶变换域内的不同位置，可以实现压缩采样。图5(a)为采用压缩方式，获得被测物体的部分傅里叶频谱分布图。图5(b)为对重建的频谱图做反傅里叶变换，并通过空间色彩插值算法恢复的被测物体的彩色图像。比较图5(b)和图5(d)的重建结果可知，采用压缩采样和全采样的方式对彩色物体测量，重建的结果基本一致，但采用压缩采样方式能够显著减少所需测量次数，缩短测量时间。
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以权利要求所述为准。