一种对激光窃听器定位的识别方法及识别装置

摘要

本发明公开一种对激光窃听器定位的识别方法及识别装置，所述识别方法包括：待定位的激光窃听器向玻璃发射激光，且所述激光通过所述玻璃形成散射光和透射光；通过声光可调谐光谱成像AOTF装置接收所述AOTF装置与玻璃中间的凝视检测区域内的透射光，获得散射光斑图像；根据所述散射光斑图像确定所述待定位的激光窃听器的方位。本发明对激光窃听器定位的识别方法利用AOTF技术对近红外波段进行连续窄带光谱滤波，通过对窗户玻璃进行光谱成像，并且结合玻璃的后向散射光斑的图像特征，通过控制装置确定所述待定位的激光窃听器的方位，提高识别的准确度。

说明书

技术领域

本发明涉及定位追踪技术领域，特别是涉及一种对激光窃听器定位的识别方法及识别装置。

背景技术

激光窃听是指通过激光技术探测到远处人谈话的内容，这种窃听方式作用距离较远，而且不容易受到干扰，最重要的一个特点是不需要在窃听目标周围安装任何设备。

随着激光窃听行为的日益严重，各种商业、科技情报信息愈来愈容易被泄露，使得研究激光窃听告警技术是非常重要的问题，而且目前尚未有比较完整、系统的激光窃听告警系统的分析方法。

然而对于目前的激光告警技术，窃听激光束照射在冕牌玻璃、钢化玻璃、宝石蓝玻璃等任何一种玻璃材料上时，都会使光束在玻璃上发生Rayleigh散射现象，从而无法追踪到窃听激光器的方位信息。其次，激光窃听技术采用的激光器波段均为红外波段，人眼以及普通相机无法直接探测到散射光斑，进而给反窃听预警带来了极大困难。此外，由于现有的散射探测式激光告警系统主要是针对战场环境上的激光制导武器发出的导引光束，而激光窃听技术通常采用固体激光器，两者光束在激光波长、散射模式、照射光斑直径、照射模式等方面都存在着很大差异，因此，散射探测式激光告警定位研究方法不能直接运用于激光窃听领域，需重新进行系统结构设计以及方位反演方法的研究。

发明内容

本发明的目的是提供一种对激光窃听器定位的识别方法，可准确确定待定位的激光窃听器的方位。

为实现上述目的，本发明提供了如下方案：

一种对激光窃听器定位的识别方法，其特征在于，所述识别方法包括：

待定位的激光窃听器向玻璃发射激光，且所述激光通过所述玻璃形成散射光和透射光；

通过声光可调谐光谱成像AOTF装置接收所述AOTF装置与玻璃中间的凝视检测区域内的透射光，获得散射光斑图像；

根据所述散射光斑图像确定所述待定位的激光窃听器的方位。

可选的，所述确定所述待定位的激光窃听器的方位的方法包括：

对所述散射光斑图像进行二值化预处理，得到预处理图像；

对所述预处理图像进行边缘检测，得到灰度近似值对应的二值化图像；

提取所述二值化图像中对象的轮廓；

对所述轮廓进行椭圆拟合，获得椭圆参数；

根据所述椭圆参数构建激光窃听器与玻璃之间的数学模型，以确定待定位的激光窃听器的方位。

可选的，所述椭圆参数包括椭圆倾斜角度、椭圆长轴、短轴及中心点坐标。

可选的，所述待定位的激光窃听器的方位包括待定位的激光窃听器的方位角和俯仰角。

根据本发明提供的具体实施例，本发明公开了以下技术效果：

本发明对激光窃听器定位的识别方法利用AOTF技术对近红外波段进行连续窄带光谱滤波，通过对窗户玻璃进行光谱成像，并且结合玻璃的后向散射光斑的图像特征，通过控制装置确定所述待定位的激光窃听器的方位，提高识别的准确度。

本发明的目的是提供一种对激光窃听器定位的识别装置，可准确确定待定位的激光窃听器的方位。

为实现上述目的，本发明提供了如下方案：

一种对激光窃听器定位的识别装置，所述识别装置包括：

玻璃，所述玻璃在受到激光窃听器的发射激光照射后，使所述发射激光分为散射光和透射光；

AOTF装置，所述AOTF装置对应所述玻璃设置，且与所述玻璃之间形成凝视检测区域，所述AOTF装置用于获取所述凝视检测区域内经过散射后的透射光，形成散射光斑图像；

控制装置，与所述AOTF装置连接，用于根据所述散射光斑图像确定所述激光窃听器的位置。

可选的，所述控制装置包括：

预处理单元，用于对所述散射光斑图像进行二值化预处理，得到预处理图像；

边缘检测单元，用于对所述预处理图像进行边缘检测，得到灰度近似值对应的二值化图像；

提取单元，用于提取所述二值化图像中对象的轮廓；

拟合单元，用于对所述轮廓进行椭圆拟合，获得椭圆参数；

建模单元，用于根据所述椭圆参数构建激光窃听器与玻璃之间的数学模型，以确定待定位的激光窃听器的方位。

可选的，所述椭圆参数包括椭圆倾斜角度、椭圆长轴、短轴及中心点坐标。

可选的，所述待定位的激光窃听器的方位包括待定位的激光窃听器的方位角和俯仰角。

相对于现有技术，本发明对激光窃听器定位的识别装置与上述对激光窃听器定位的识别方法的有效效果相同，在此不再赘述。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1为本发明对激光窃听器定位的识别方法的流程图；

图2为本发明对激光窃听器定位的识别方法的倾斜椭圆的计算原理图一；

图3为本发明对激光窃听器定位的识别方法的倾斜椭圆的计算原理图二；

图4为待定位激光窃听器的位置信息图；

图5为待定位激光窃听器的俯仰角图；

图6为本发明对激光窃听器定位的识别装置的模块结构示意图；

图7为本发明对激光窃听器定位的识别装置的工作原理图。

符号说明：

激光窃听器1 玻璃2

AOTF装置3 控制装置4

预处理单元41边缘检测单元42

提取单元43拟合单元44

建模单元44散射光5

透射光6 凝视检测区域7。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明的目的是提供一种对激光窃听器定位的识别方法，通过声光可调谐光谱成像(Acousto-Optic Turnable Filter，AOTF)装置可以任意选择系统的光谱，选择具有目标特性不同波段的光谱波长，从而可以形成同一目标在不同光谱波长下的不同图像。利用AOTF技术对近红外波段进行连续窄带光谱滤波，通过对窗户玻璃进行光谱成像，并且结合玻璃的后向散射光斑的图像特征，通过控制装置确定所述待定位的激光窃听器的方位，提高识别的准确度。

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。

如图1所示，本发明对激光窃听器定位的识别方法包括：

步骤100：待定位的激光窃听器向玻璃发射激光，且所述激光通过所述玻璃形成散射光和透射光。

步骤200：通过声光可调谐光谱成像AOTF装置接收所述AOTF装置与玻璃中间的凝视检测区域内的透射光，获得散射光斑图像。

步骤300：根据所述散射光斑图像确定所述待定位的激光窃听器的方位。

其中，在步骤200中，通过多次调整探测角度或位置，对同一块玻璃进行照射，从而获取对应的散射光斑。

在步骤300中，所述确定所述待定位的激光窃听器的方位的方法包括：

步骤310：对所述散射光斑图像进行二值化预处理，得到预处理图像。

步骤320：对所述预处理图像进行边缘检测，得到灰度近似值对应的二值化图像。在本实施例中，采用Sobel算子法做边缘检测。

步骤330：提取所述二值化图像中对象的轮廓。一般利用bwboundaries函数获取二值化图像中对应的轮廓。

步骤340：对所述轮廓进行椭圆拟合，获得椭圆参数。

其中，所得的轮廓利用最小二乘法进行椭圆拟合。所述椭圆参数包括椭圆倾斜角度、椭圆长轴、短轴及中心点坐标。

步骤350：根据所述椭圆参数构建激光窃听器与玻璃之间的数学模型，以确定待定位的激光窃听器的方位。

所述对所述轮廓进行椭圆拟合的方法具体包括：

(1)计算倾斜椭圆的系数

设倾斜椭圆方程式为ax²+bxy+cy²+dx+ey-1＝0，利用最小二乘法即求的最小值，所以两边求导得矩阵表示为：

$>\left(\begin{array}{ccccc}\underset{i=1}{\overset{N}{\Sigma }}{x}_i^4& \underset{i=1}{\overset{N}{\Sigma }}{x}_i^3{y}_i& \underset{i=1}{\overset{N}{\Sigma }}{x}_i^2{y}_i^2& \underset{i=1}{\overset{N}{\Sigma }}{x}_i^3& \underset{i=1}{\overset{N}{\Sigma }}{x}_i^2{y}_i\\ \underset{i=1}{\overset{N}{\Sigma }}{x}_i^3{y}_i& \underset{i=1}{\overset{N}{\Sigma }}{x}_i^2{y}_i^2& \underset{i=1}{\overset{N}{\Sigma }}{x}_i{y}_i^3& \underset{i=1}{\overset{N}{\Sigma }}{x}_i^2{y}_i& \underset{i=1}{\overset{N}{\Sigma }}{x}_i{y}_i^2\\ \underset{i=1}{\overset{N}{\Sigma }}{x}_i^2{y}_i^2& \underset{i=1}{\overset{N}{\Sigma }}{x}_i{y}_i^3& \underset{i=1}{\overset{N}{\Sigma }}{y}_i^4& \underset{i=1}{\overset{N}{\Sigma }}{x}_i{y}_i^3& \underset{i=1}{\overset{N}{\Sigma }}{y}_i^3\\ \underset{i=1}{\overset{N}{\Sigma }}{x}_i^3& \underset{i=1}{\overset{N}{\Sigma }}{x}_i^2{y}_i& \underset{i=1}{\overset{N}{\Sigma }}{x}_i{y}_i^2& \underset{i=1}{\overset{N}{\Sigma }}{x}_i^2& \underset{i=1}{\overset{N}{\Sigma }}{x}_i{y}_i\\ \underset{i=1}{\overset{N}{\Sigma }}{x}_i^2{y}_i& \underset{i=1}{\overset{N}{\Sigma }}{x}_i{y}_i^2& \underset{i=1}{\overset{N}{\Sigma }}{y}_i^3& \underset{i=1}{\overset{N}{\Sigma }}{x}_i{y}_i& \underset{i=1}{\overset{N}{\Sigma }}{y}_i^2\end{array}\right)·\left(\begin{array}{c}a\\ b\\ c\\ d\\ e\end{array}\right)=\left(\begin{array}{c}\underset{i=1}{\overset{N}{\Sigma }}{x}_i^2\\ \underset{i=1}{\overset{N}{\Sigma }}{x}_i{y}_i\\ \underset{i=1}{\overset{N}{\Sigma }}{y}_i^2\\ \underset{i=1}{\overset{N}{\Sigma }}{x}_i\\ \underset{i=1}{\overset{N}{\Sigma }}{y}_i\end{array}\right)$>

若令X_i＝(x_i²,x_iy_i,y_i²,x_i,y_i)，上述矩阵可写成

(2)计算标准椭圆与倾斜椭圆的关系

由于倾斜椭圆与标准椭圆在坐标上存在着一定的关系，所以利用倾斜椭圆计算标准椭圆，存在两种情况，分别如下：

第一种情况下，设(x,y)是标准椭圆上一点坐标，(x′,y′)是倾斜椭圆上一点坐标，则由图2中(a)和(b)得两坐标关系：

x′＝r·cos(α-θ)＝r cosθcosα+r sinθsinα＝x cosα+y sinα；

y′＝r·sin(α-θ)＝r cosθsinα-r sinθcosα＝x sinα-y cosα；

其中，r表示(x,y)或(x′,y′)到原点的距离，θ表示(x,y)点到原点的直线与x轴的夹角，α表示(x′,y′)点到原点的直线与(x,y)点到原点的直线的夹角。

第二种情况下，设(x,y)是标准椭圆上一点坐标，(x′,y′)是倾斜椭圆上一点坐标，则由图3中(a)和(b)得两坐标关系：

x′＝r·cos(α+θ)＝r cosθcosα-r sinθsinα＝x cosα-y sinα

y′＝r·sin(α+θ)＝r cosθsinα+r sinθcosα＝x sinα+y cosα

(3)计算标准椭圆的系数

利用步骤(2)中倾斜和标准椭圆的关系分别计算标准椭圆的系数。设倾斜椭圆方程式为ax′²+bx′y′+cy′²+dx′+ey′-1＝0，标准椭圆方程式为a′x²+c′y²+d′x+e′y-1＝0。

第一种情况下，将x′＝x cosα+y sinα，y′＝x sinα-y cosα代入倾斜椭圆方程式，即