基于图像处理的目标活动分析方法_高飞.pdf

工程应用航天电子对抗2023年第 1期基于图像处理的目标活动分析方法高飞，徐才进，梁璟，倪辰辰，杨玲（中国电子科技集团公司第二十九研究所，四川 成都 610036）摘要：提出了基于图像处理技术的目标活动航路分析和区域分析方法，对于稀疏航迹点通过侦察数据图像化、二值化以及图像滤波实现了目标活动航路提取，对于密集航迹点通过侦察数据热度图生成、图像去噪、图像增强实现了目标活动区域提取。使用仿真数据验证了将图像处理技术引入目标活动分析的有效性和优点。关键词：图像处理；活动分析；情报中图分类号:TN97 文献标志码:AAnalysis method of target activity based on image processingGao Fei，Xu Caijin，Liang Jing，Ni Chenchen，Yang Ling(The 29th Research Institute of China Electronics Technology Group Corporation，Chengdu 610036，Sichuan,China)Abstract：The method of target activity route analysis and regional analysis based on imagine processing technology is presented.The target activity route is extracted through reconnaissance data imaged，binaryzation and image filtering for sparse track，and the target activity regional is extracted through reconnaissance data heat map generating，image denoising and image enhancement for dense track.The effectiveness and advantage of the method used to analyse target activity with image processing is proved by simulation data.Key words：image processing；activity analysis；intelligence0 引言随着现代战场环境日益复杂以及战场目标活动愈发频繁，传统目标情报中包括的基本属性、信号参数等静态信息已经难以满足联合作战对支援情报的需求。基于历史侦察数据的目标活动规律分析可以挖掘出重点目标的典型航路和频繁活动区域，不仅能为联合作战指挥控制、作战筹划等业务提供战场目标数据支撑，还能在战场态势感知中用于目标活动异常检测和预警等。受益于战场侦察能力的提升和情报信息系统能力的增强，战场目标的历史侦察数据在系统中持续积累，为目标活动分析提供了良好的条件。作为战场情报分析的主要组成部分，目标活动分析依赖战场侦察数据开展工作，传统的目标活动分析方法 1 以侦察数据统计、挖掘为主，辅以人工标注等干预手段，虽然能得到较为完整准备的目标活动分析结论，但是存在挖掘算法设计困难、大数据计算耗时较长以及人工标注效率较低等问题。蓬勃发展的图像处理技术为目标活动分析提供了新的技术思路，将侦察数据在地图上进行图像化，使用图像处理方法进行目标活动分析，可以有效减少算法设计难度、提升数据分析效率并增强分析结论可视化水平。本文创新性地将图像处理方法引入目标活动规律分析。由于环境和装备能力的影响，不同目标的历史侦察数据情况有所不同，根据数据量的数据分布情况可以分为历史稀疏航迹点和历史密集航迹点 2种情况：对于稀疏航迹点通过侦察数据图像化、二值化以及图像滤波实现了目标活动航路提取；对于密集航迹点通过侦察数据热度图生成、图像去噪、图像增强实现了目标活动区域提取。最后，使用仿真数据验证了将图像处理技术引入目标活动分析的有效性。1 稀疏航迹点目标活动航路分析针对稀疏型航迹点数据，主要研究基于历史稀疏型航迹点的航路分析，具体处理分析流程如图1所示。收稿日期：20220810；20221210修回。作者简介：高飞（1990-），男，工程师，硕士，主要研究方向为电子对抗数据情报处理。-41DOI:10.16328/j.htdz8511.2023.01.0012023，39（1）航天电子对抗1.1原始数据解析与航迹图平滑化通过战场情报侦察等手段获取的目标位置信息主要是经纬度信息，在进行基于数字图像处理的活动航路分析前，需要将这种地理坐标系统转换为平面坐标系统，形成平面坐标航迹图。考虑到侦察时的环境影响、装备状况和目标混叠等情况，原始侦察数据的经纬度信息中可能存在一定误差和噪声，可以采用数字图像处理中平滑模糊滤波方法 2，使航迹图像平滑模糊化并具有一定的容错性。图像空间滤波器一般由一个邻域和对邻域所包围图像像素执行的预定义操作组成。对每个像素滤波器产生一个新像素，其坐标等于邻域中心坐标，其像素值是滤波的响应结果。滤波器依次访问输入图像中的每个像素后，就生成了滤波后的图像。对图像中的任意一点(x，y)，其像素值为f(x，y)，滤波器的响应g(x，y)是滤波器系数与由该滤波器所包围的图像像素的乘积之和。滤波器邻域一般为奇数尺寸，最小为 33。对大小为m n的邻域，设奇数m=2a+1，奇 数n=2b+1，其 中a、b为 正 整 数。使 用 大 小 为m n的滤波器对大小为M N的图像进行线性空间滤波，可表示为：g(x，y)=s=-aat=-bbw(s，t)f(x+s，y+t)（1）例如，对于 33邻域的线性空间滤波器：g(x，y)=w(-1，-1)f(x-1，y-1)+w(-1，0)f(x-1，y)+w(0，0)f(x，y)+w(1，1)f(x+1，y+1)（2）1.2航迹图二值化进行目标活动航路分析时，图像中只存在航路区域和非航路区域，是一个典型的二值化问题。图像的二值化就是通过确定一个阈值，把灰度图像的每个像素按其灰度值分为 2类，使得原图像变为二值图像，即灰度图像的二值化 3。设灰度图像 P尺寸为 mn，经过上文滤波后处理的每一像素灰度值为 x（i，j）。图像像素灰度均值为：x =i=1mj=1nx(i，j)/(m n)（3）设整数 min、max，其中 min 为最小灰度像素灰度值向下取整，max为最大灰度像素灰度值向上取整。设整数 k min，max，对图像进行遍历，得到图像中灰度值近似为 k 的像素点个数Ak。遍历开始时Ak=0，对任一像素点，若 k-0.5 x（i，j）|n-x|的 n时，认为该 n为另一波峰灰度值。所要取的阈值 T 应在 n与 n之间：T=n+n-x（4）取得阈值后，对图像进行二值化运算得到边缘清晰的二值化图像：x(i，j)=|1，x(i，j)T0，x(i，j)T（5）1.3航路区域提取提取出二值化后的航迹图像后，研究采用数字图像处理中的数学形态学滤波方法，分析提取出航路区域。数学形态学可以很好地对图像中存在的一些空洞或断点进行处理，得到更加本质的图像形态，主要包括膨胀、腐蚀、细化、粗化、开运算和闭运算等。膨胀是将目标区域的背景点合并到该目标物中，使目标物边界向外部扩张的处理。选取结构元素B，将图像X中的每一个像素点扩大为B+x，膨胀可以使图像断点连合，成为一个整体。XB=X x+B|x X（6）腐蚀是将目标区域的边界收缩，直观上就是使得目标区域变小，与膨胀是互为对偶的过程，但是并非经过 2 种运算后图像可以还原。选取结构元素，将图像中的的B+x 收缩为点 x。腐蚀可以消除边界处的干扰点。XB=x|B+x X（7）图 1稀疏航迹点目标活动航路分析流程-42工程应用高飞，等：基于图像处理的目标活动分析方法2023，39（1）另外，开运算和闭运算是为了解决膨胀和腐蚀运算改变图像大小而提出来的。开运算是先腐蚀后膨胀的过程，可以消除图像中细小噪声，平滑物体边界。闭运算是先膨胀后腐蚀的过程，可以消除图像中存在的空洞，平滑物体边界。2 密集航迹点目标活动区域分析针对密集型航迹点数据，主要研究基于历史密集型航迹点的区域分析，具体处理分析流程如图2所示。2.1原始数据解析与航迹点热度图生成通过解析原始数据，将经纬度信息转换为平面坐标系下的航迹图，便于进行图像处理。统计每个像素点区域目标出现的热度（像素点位置出现一次航迹点，便将像素点灰度值加 1 或者其他固定值），分析处理后，目标出现越频繁的像素点区域灰度值越高，图像对应区域就越高亮显示。2.2航迹点热度图去噪通过图像学中值滤波去除热度图中的“非频繁噪声点”。中值滤波的基本原理是把数字图像或数字序列中一点的值用该点的一个邻域中各点的中值代替，从而消除孤立的噪声点 4。设二维矩 xij 表示一副数字化的图像（其中i、j表示 各 像 素 点 的 位 置，xij表 示 各 像 素 点 的 灰 度 值），W xij表 示 对 图 像xij中 的 点xij做 窗 口 操 作，Med xij表示对窗口 W xij 内的所有点取中值，若yij 为图像 xij 经过中值滤波后的输出图像，滤波窗口的大小为m n（其中m、n为奇整数），那么：yij=Med(W xij)（8）2.3航迹点热度图增强在数字图像处理中，直方图均衡化（HE）是一种最常用的空域直方图修正技术，其以概率理论为基础，运用灰度点运算来实现对图像直方图的变换，从而达到增强图像对比度的目的。HE 是把给定图像的直方图分布改造成均匀直方图分布，使输出像素灰度的概率密度均匀分布。基本原理是当图像中所有的灰度级出现的概率是一个均匀分布时，图像所暴露出来的信息量最大 5。空域图像增强的基本方法是对图像做灰度变换T：f(x，y)g(x，y)。直方图修正试图寻找灰度变换函数T：它能把f(x，y)映射为g(x，y)，使得g(x，y)的直方图Pg(r)具有指定的形式qr，因此它不是事先指定变换函数，而是从Pf(r)及qr中计算得到T。任一图像f(x，y)的直方图Pf(r)在变换s=T(r)下具有以下形式：s=T(r)=Pf(r)=01Pf(v)dv（9）式中，T(r)在区间 0，1 上是单调上升的单值函数，对于0 r 1，有0 T(r)1。实际应用当中，处理的图像都是数字图像，因此需要离散化，HE 的离散化表达式为：Pf(rk)=nk/(NL NS)，k=0，1，255（10）式中，NL、NS分别是图像行数和列数，从而可以得到：sk=255Pf(rk)=255i=0kPf(ri)=255i=0kni/(NL NS)（11）3 仿真验证基于上述目标活动航路和活动区域分析方法，构造对应的稀疏、密集仿真数据，验证其理论算法和分析流程的有效性。3.1稀疏航迹点目标活动航路分析验证构建稀疏的仿真目标活动位置信息，形成原始的平面坐标点迹图，如图 3所示。基于时间序列绘制其航路，中间缺失的航迹点采用直线连接表示，如图 4所示。图 2密集航迹点目标活动区域分析流程图 3原始点迹图-432023，39（1）航天电子对抗采用高斯滤波进行原始航迹图平滑化处理。根据航迹点在图像上的成像分辨率大小，选择一定半径的高斯核对航迹图像进行平滤波滑，针对仿真数据选择的高斯核大小为 55，滤波效果如图 5所示。航迹图平滑模糊滤波后，进行图像二值化处理，便于进行航路区域分析提取，试验效果如图 6所示。在二值化航迹图像的基础上，对二值化图像重复进行 33 或 55 等区域大小的图像学膨胀与腐蚀处理，重复进行 2至 3次，试验效果如图 7所示。为验证航路分析提取效果，将其与原始航迹数据叠加进行验证，试验效果如图 8所示。3.2密集航迹点目标活动区域分集验证构建密集的仿真目标活动位置信息，并将地球经纬度转换为平面像素坐标系下的绝对坐标值，形成如图 9所示的原始航迹点图像。在原始航迹点数据基础上，通过灰度图亮度值累加，生成航迹点热度图