一种互补单像素质心探测系统及方法

1.本发明属于计算成像和质心探测的技术领域，尤其涉及一种互补单像素质心探测系统及方法。


背景技术：

2.在质心探测领域，传统的方法均是通过面阵相机先获取物体图像，再由图像经过相应图像处理算法计算出目标物体质心，捕获的物体图像质量决定了探测质心的准确度。
3.在低信噪比的情况下，面阵相机拍摄到图像的质量会受到较大影响，并且一些不可见波段(红外，太赫兹等)面阵相机要么造价昂贵要么无法有效工作。传统探测质心的方法使用范围受到限制。单像素成像技术与传统的ccd(charge coupled device,电荷耦合元件)和cmos(complementary metal oxide semiconductor,互补金属氧化物半导体)等面阵感光元件为核心的数字成像技术相比，单像素成像技术使用时空变换的光照射目标物体或者调制物体图像，通过只有一个像素单元的探测器来进行物理信息采样，最后利用相应算法重建物体信息。
4.单像素成像的优势主要体现在两个方面，由于单像素探测器具有较宽的光谱响应范围和较高的光敏感度，单像素成像可以应用在阵列相机无法响应或者造价昂贵的波段以及弱光成像。将单像素成像方法与质心探测结合可以在弱光和宽波谱响应范围进行工作，同时摆脱成像帧频不高的限制。
5.现有技术中一种基于单像素探测器的质心探测方法(申请号：202010412816.0)在不成像的前提下使用三张二维调制信息直接目标物体的质心，这样速度较慢。


技术实现要素：

6.为了提高质心探测的速度，为此，本发明提出了一种互补单像素质心探测系统及方法，具体方案如下：
7.一种互补单像素质心探测系统，包括透镜组、单像素探测器组件、光调制器dmd、采集处理单元；
8.采集处理单元生成二维调制矩阵a和二维调制矩阵b并加载到光调制器dmd中,令二维调制矩阵a每列元素值等于列数，二维调制矩阵b每行元素值等于行数；
9.透镜组将目标物体反射或透射后的光进行处理，使目标物体图像成像在光调制器dmd上；
10.光调制器dmd对目标物体图像信号按照二维调制信息进行调制；
11.单像素探测器组件包括第一单像素探测器和第二单像素探测器,分别获取由光调制器dmd两个互补方向上反射光的强度值；
12.采集处理单元分别与第一单像素探测器和第二单像素探测器的数据输出端连接，计算出目标物体的质心。
13.具体地说，所述光调制器dmd由若干个微镜组成，每个微镜均可在正负角度之间来
回转动，对应着光调制器dmd微元的“0”和“1”状态；光调制器dmd在两个反射方向上的光强是互补。
14.具体地说，所述光调制器dmd的型号为dlp7000，由768x1024个微镜组成。
15.具体地说，还包括光源，所述光源用于照射目标物体。
16.具体地说，所述采集处理单元包括连接的计算机和数据采集器，所述计算机与光调制器dmd连接，所述数据采集器端口与第一单像素探测器和第二单像素探测器对应连接。
17.一种互补单像素质心探测方法，包括以下步骤：
18.s1、采集处理单元中的计算机生成二维调制矩阵a和二维调制矩阵b，并且分别使用二维函数表示；
19.二维函数所在的坐标系的横坐标和纵坐标方向分别对应二维调制矩阵的行和列方向，二维函数值对应二维调制矩阵的元素值；二维调制矩阵a和二维调制矩阵b所对应二维函数分别满足以下关系式：
20.s1(x,y)＝x
ꢀꢀꢀ
(1)
21.s2(x,y)＝y
ꢀꢀꢀ
(2)
22.式中s
n
(x,y)(n＝1,2)表示二维调制矩阵中对应着坐标(x,y)处的元素值；
23.s2、将生成调制信息加载至光调制器dmd中，对目标物体图像进行调制；
24.s3、数据采集器利用第一单像素探测器和第二单像素探测器获取由光调制器dmd两个互补方向上的反射光的强度值；
25.s4、根据质心求解算法求取质心。
26.优化的，当步骤s1还包括：采用空间抖动的方法将二维调制矩阵a和二维调制矩阵b转换为二值化调制矩阵，步骤s2中的调制信息为二值化调制矩阵。
27.具体地说，步骤s3计算第一单像素探测器(和第二单像素探测器所获取的反射光强度值表示为：
[0028][0029][0030]
式中f(x,y)为目标物体图像的二维分布函数，上标取值0和1，分别对应着光调制器dmd的正负角度的状态,下标取值1和2分别代表单像素探测器和单像素探测器，k取1和2分别表示光调制器dmd两个二维调制信息,和分别表示第k个调制信息下互补的两个调制信息。
[0031]
具体地说，所述质心求解算法求解质心(x
c
,y
c
)的计算公式如下：
[0032][0033][0034]
通过联立(3)、(4)式与(5)、(6)式，目标物体的质心位置参数可以表示为：
[0035][0036][0037]
本发明的有益效果在于：
[0038]
(1)本发明在不成像的前提下直接目标物体质心，并且系统利用数字微镜阵列光调制器dmd互补结构特点，进一步减少所用调制信息的数量，提高质心探测的速度。
[0039]
(2)光源的设置用于照射目标物体，可以使得目标物体更好的成像。
[0040]
(3)方法中当步骤s1采用空间抖动的方法转换成二值化调制矩阵满足了高速调制的需求，并充分利用光调制器dmd高速二值调制性能。
附图说明
[0041]
图1为互补单像素质心探测的流程图；
[0042]
图2为互补单像素质心探测系统结构图；
[0043]
图3是光调制器(dmd，digital micromirror device)中微镜结构示意图；
[0044]
图4是生成的两张调制信息图。
[0045]
图中：
[0046]
1、光源；2、目标物体；3、透镜组；41、第一单像素探测器；42、第二单像素探测器；5、光调制器dmd；6、采集处理单元。
具体实施方式
[0047]
如图2所示，一种互补单像素质心探测系统，包括透镜组3、单像素探测器组件4、光调制器dmd5、采集处理单元；优化的，还可以包括用于给目标物体照明的光源1。所述采集处理单元包括连接的计算机6和数据采集器。
[0048]
计算机生成二维调制矩阵a和二维调制矩阵b并加载到光调制器dmd5中,令二维调制矩阵a每列元素值等于列数，二维调制矩阵b每行元素值等于行数；
[0049]
透镜组3将目标物体2反射或透射后的光进行处理，使目标物体图像成像在光调制器dmd5上；
[0050]
光调制器dmd5对目标物体图像信号按照维调制信息进行调制。图4展示了其中部分微镜结构的示意图，图3中每个微镜均可在正负角度之间来回转动，对应着光调制器dmd5微元的“0”和“1”状态。光调制器dmd5在两个反射方向上的光强是互补的，也就是说两束反射光的强度之和就等于光调制器dmd5作为平面镜时反射光的强度。在该实施例中，所述光调制器dmd5的型号为dlp7000，由768x1024个微镜组成，正负角度为
±
12
°
。
[0051]
单像素探测器组件包括第一单像素探测器41和第二单像素探测器42,分别获取由光调制器dmd5两个互补方向上反射光的强度值，所述第一单像素探测器41和第二单像素探测器42均为单点光电探测器。
[0052]
采集处理单元中的数据采集器分别与第一单像素探测器41和第二单像素探测器42的数据输出端连接，数据采集器将采集的数据上传到计算机6中，计算机6计算出目标物
体2的质心。
[0053]
如图1所示，基于上述系统的互单像素质心探测方法，包括以下步骤：
[0054]
s1、采集处理单元中计算机成二维调制矩阵a和二维调制矩阵b，并且分别使用二维函数表示；
[0055]
二维函数所在的坐标系的横坐标和纵坐标方向分别对应二维调制矩阵的行和列方向，二维函数值对应二维调制矩阵的元素值；二维调制矩阵a和二维调制矩阵b所对应二维函数分别满足以下关系式：
[0056]
s1(x,y)＝x
ꢀꢀꢀ
(1)
[0057]
s2(x,y)＝y
ꢀꢀꢀ
(2)
[0058]
式中s
n
(x,y)(n＝1,2)表示二维调制矩阵中对应着坐标(x,y)处的元素值；
[0059]
若需要高速调制，为了充分利用光调制器dmd高速二值调制性能，可采用空间抖动的方法利用计算机6将二维调制矩阵a和二维调制矩阵b转换为二值化调制矩阵，如图4所示，空间抖动的方法的运用可以充分利用光调制器dmd5的高速二值化调制频率。
[0060]
s2、将图4中的调制信息加载至光调制器dmd5中，对由目标物体图像进行调制；当不需要高速调制时，所述调制信息为二维调制矩阵a和二维调制矩阵b对应的函数，当需要高速调制时，所述调制信息为二值化调制矩阵。
[0061]
s3、利用第一单像素探测器41和第二单像素探测器42获取由光调制器dmd5两个互补方向上的反射光的强度值。第一单像素探测器41和第二单像素探测器42所获取的反射光强度值表示为：
[0062][0063][0064]
式中f(x,y)为目标物体(2)图像的二维分布函数，上标取值0和1，分别对应着光调制器dmd(5)的正负角度的状态,下标取值1和2分别代表单像素探测器(41)和单像素探测器(42)，k取1和2分别表示光调制器dmd(5)两个二维调制信息,和分别表示第k个调制信息下互补的两个调制信息。
[0065]
s4、求取质心(x
c
,y
c
)，计算公式如下：
[0066][0067][0068]
通过联立(3)、(4)式与(5)、(6)式，目标物体2的质心位置参数可以表示为：
[0069]
[0070][0071]
本发明将单像素成像方法与质心探测方法相结合，充分利用光调制器dmd5的结构特点，在不成像的前提下仅用两张调制信息直接获取目标物体2的质心，在光调制器dmd5的最高调制频率22khz下，利用该技术可以实现11khz以上的质心探测帧频，可以应用在快速目标物体2追踪等领域。
[0072]
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。