一种基于特征统计的动态伪装效果评估方法

杨鑫, 许卫东, 贾其

(陆军工程大学 野战工程学院， 江苏 南京 210007)

摘要： 现有的伪装效果评估主要针对静止的单幅图像，不能很好地模拟侦察人员对目标的判读过程。结合Mean shift目标跟踪技术，提出一种基于特征统计的动态伪装效果评估方法。该方法通过统计目标与背景8联通域的相关性特征数据，建立归一化联合高斯分布，利用概率密度的分布范围评估目标伪装效果。计算联合分布的概率密度时，提出对数放大概率，解决了高维联合分布概率密度数值敏感度低、不便于阈值设定的问题。引入样本更新策略，使样本库按照一定的概率随机更新，从而较好地适应了由于季节交替等因素引起的背景大范围变化。实验过程分别对某一指挥车实施1级伪装、2级伪装和3级伪装。采集数据后计算其对数放大概率并对曲线作出统计，结果表明：实际中划分的3种伪装状态与依据3σ 准则预先设定的3种伪装状态完全对应；该模型能够有效反映出目标的伪装效果。

关键词： 伪装效果； 效果评估； 特征统计； 概率密度； 高斯分布； 目标跟踪

0 引言

现代战争中，随着侦察监视技术的不断发展以及精确制导武器的广泛使用,国防工程等固定式目标面临越来越大的生存威胁。伪装作为一种重要作战保障措施，所面临的侦察监视和目标捕获等威胁正变得愈来愈严重。如何科学、全面、准确地评估伪装效果变得非常重要。伪装效果评估体系内涵丰富，其中伪装目标在单台侦察器材探测下的评估方法——设计级伪装效果评估，长期以来为业界所重视^[1]。设计级伪装效果评估目前主要有两种思路：一种是基于人员对目标发现概率的判读统计，该方法受观察人员的经验、生理条件的影响比较大；另一种是基于数字图像处理的方法，通过数字图像中伪装目标与周围背景的特征差别，分析伪装效果。基于这种思路，人们以图像区域的亮度、纹理、形状和结构等为特征进行计算处理；也有人基于心理学建立相应的刺激函数，或者基于反向传播(BP)神经网络进行特征组合^[2-5]。这些方法均属于单幅静止图像的特征评估。

实际工程中，伪装效果评估结果受到侦察器材角度、侦察环境以及目标和背景特性变化等因素的影响。例如在不同条件下侦察，目标和背景亮度、颜色、斑点形状等特性均有不同，但是符合一定的变化规律。因此伪装效果评估应具备动态性、统计性等特点。本文正是基于目标跟踪和特征统计技术，提出一种动态伪装效果评估方法。

第一，有组织犯罪的成员应当达到3人或3人以上。有组织犯罪作为一种组织性犯罪，其成员自然不可能是一个，而两个人也很难被称为组织。实际上，有组织犯罪的成员至少有四、五人，一般都能达到十几人、几十人，只有3个人的是极个别情况。而像组织遍布世界各地的意大利“黑手党”这种黑社会组织，其党徒有几十万之多。

1 动态评估过程分析

1 .1 连续侦察过程分析

实际侦察过程中，敌方侦察人员往往会从不同角度、在不同条件下对目标实施连续侦察。侦察过程如图1所示。相比于静态的单角度侦察过程，一方面，侦察员能够联系前后不同帧之间的图像关系，分析目标的特性；另一方面，目标在不同角度、不同条件下的斑点形状、颜色亮度等特征也不同。因此评估过程也应具备对连续图像分析的特性。

1 .2 特征提取

目标发现过程的实质可抽象为通过特征的提取和比对，分割目标与背景区域的过程。人眼对图像底层特征的提取主要包括亮度、颜色、纹理和斑点形状等。亮度和颜色的分布特性表征不同频率的光谱反射特性和背景空间分布及其组合规律。斑点形状特性是视觉上很重要的特征，可以通过矩形度、圆形度、形状描述子等方法描述。纹理是一种区域特征，可通过定量方式对纹理结构进行描述，如灰度共生矩阵等。工程实践中，应结合目标具体的特性选择相应的特征提取算法。

图1 连续侦察流程示意图
Fig.1 Flow chart of continuous reconnaissance

从评估的角度分析，划分目标与背景区域较合理的方法是文献[6]所提出的8联通域法，如图2所示。该方法只考虑目标区域附近的背景，并非将全部图像区域划为背景区域，从而符合人眼判读的规律。尤其对于斑驳背景而言，目标仅融于临近背景，与远区域背景无关。考虑到运算的简洁性，本文使用矩形对目标区域进行分割且采用图2中对背景区域的编号顺序编号。

赵之谦：田绍登[9]对赵之谦的研究角度则相对特别，通过对赵之谦和何绍基进行比较认为，赵之谦早年取法于何，行书基本是一派“颜家样”，全无碑趣。赵弱化了何书法中的北碑化因素，强化了北碑的笔法风格特征，并使之与帖学笔法加以有效的结合，化解了碑学危机，为晚期碑学实现碑帖融合的范式转换奠定了基础。在书法史上，就碑行书的原创性而言，赵之谦无疑达到了清代碑学的最高峰。他的卓荦之处在于在固守北碑生朴拙风格的同时，将北碑与帖学的融合推向完美的境地。在强调碑的骨力的同时，帖的韵致及帖学书法的精到表现，这在清代碑学家中是绝无仅有的。

图2 目标与背景区域的8联通域划分图
Fig.2 Eight-way domain partition map of target and background area

1 .3 统计与评估

同一目标的特征数据会受到光照、天气、观察角度等随机因素的影响。由中心极限定律可知，各种随机因素的叠加可近似认为其符合高斯分布。因此，统计伪装良好状态下的特征数据可以得到概率密度函数的参数，并由此建立多特征联合概率密度函数。评估过程即计算联合概率密度值，并通过设定的阈值比较，确定伪装效果等级。不同于以往的伪装效果评估方式，基于特征统计的伪装效果评估方法利用了目标特征的先验知识，能够更加客观地体现目标状态。

(1) 该地堆积体主要由粗、巨颗粒占优势，其占比约95%，细粒占比不到5%，堆积体密实度高。根据Shepard三角沉积分类，其样品点2 mm以下颗粒以粗颗粒砂占优势。

建立目标模型时，采用Epanechikov核函数：

(1)

(2)

菌种活化和菌悬液的配置：分别挑取菌种在相应的培养基上划线，然后置于恒温培养箱中培养（37℃，36 h）。用接种环挑取少量活化后的菌种接入装有无菌生理盐水的试管内，震荡均匀，制成菌悬液，然后用无菌水进行稀释，使活菌数为107-108cfu/ml，备用。

式中：a _ijk 表示第k 帧、第j 个背景、第i 个特征的值。由于标准正态分布的公式为

(3)

对于待检测帧o ，可以计算反映目标与背景特征关系变化的归一化联合概率密度：

裂缝外侧顺路线走向具有临空条件，为右岸坡面(倾向WWN)，因砂泥岩差异风化形成了高约2～4 m、深约2.5～3.0 m的凹岩腔(图8、9)，局部存在坍落、掉块现象，凹岩腔前缘发现一处新近相对滑移迹象，局部渗水，无泉点出露，滑移矢量方向约205°±5°，相对位移约3～5cm(图9、10)。

(4)

2 目标跟踪模型

设在k -1帧时目标所在的中心区域为f ₀、候选目标中心区域为f ，核函数窗口大小为h ，则候选区域的概率为

隐匿阴茎是一种先天性外生殖器畸形，指原本正常的阴茎被埋藏于皮下，包皮口与阴茎根距离短。病因是由于胚胎发育期间，正常延伸至生殖结节的尿生殖窦远端发育不全所致[3]。外观呈“鸟嘴样”或“山丘样”(图1)。按压阴茎周围皮肤可暴露正常的阴茎体，放手后恢复原状。儿童隐匿性阴茎的诊断成立，需要具备5个条件：①阴茎外观似宝塔状，②具有发育正常的阴茎体，③下按阴茎周围组织可显示阴茎全貌，松开后即恢复如初，④需排除其他先天性尿道疾患及海绵体发育不良的阴茎疾患，⑤排除肥胖病因。

2 .1 模型表示

以n 表示目标区域的像素点个数，(l _i ,s _i )代表第i 个像素点处的绝对坐标值，(l ₀,s ₀)为目标的中心点坐标，用z _i 表示其绝对位置，表示以目标中心点(l ₀,s ₀)为原点的相对位置，将灰度颜色空间均匀划分为m 个相等的区间构成灰度直方图，用u 代表直方图的灰度索引。则目标区域模型q _u 为

事实上，无论EQ或AQ，皆属于情绪教育的一环。情绪教育之核心，乃是在探讨如何培养个体体验他人情绪、控制和表达自身情绪的能力，目的在于培养个体正确感受自己与他人情绪，且适当表达自己情绪的能力，最终增进对自己与他人的信念与价值，促进个体情绪之正向发展[6]。借由情绪教育，可以提升个体的EQ与AQ。因此，情绪教育的重要性，已受到教育界的认同，并着力于思考如何协助学生舒缓情绪压力，以期能促进学生的心理健康。

(5)

式中：φ 用于确定z _i 处的灰度值在直方图中的索引；δ 与φ 函数结合用于判断z _i 处的颜色值是否属于u ；C 为归一化参数；K _E为Epanechikov核函数。和δ 函数的计算方法如下：

(6)

(7)

(8)

由于有8个背景区域，若收集到M 帧样本数据，需计算R 个特征的分布参数，联合分布的维度为8×R 维，其均值和方差估计公式如下：

(9)

在目标与背景分割过程中，为保证前后两帧目标区域位置准确，必须使用目标跟踪算法。现有目标跟踪算法可以分为4类:基于主动轮廓的跟踪、基于特征的跟踪、基于区域的跟踪和基于模型的跟踪^[7]。其中，基于特征的跟踪算法优点在于对运动目标尺度、形变和亮度等变化不敏感^[8]。本文选用基于灰度直方图特征的Mean shift均值漂移跟踪算法^[9]，应用密度梯度的非参数估计原理实现对移动目标区域的快速跟踪。

(10)

2 .2 相似性表达

在相似性度量方面，采用Bhattacharyya系数建立概率密度函数如下：

(11)

2 .3 迭代过程

购买蔬菜、水果时，不买异常气味的蔬菜和水果。有些商贩为了蔬菜更好看，用化学药剂进行浸泡，这些物质有异味，不容易被冲洗掉。不买颜色异常的蔬菜和水果，不买形状异常的蔬菜和水果。有的蔬菜使用了激素物质，会长成畸形。尽量不要食用“新奇”“野味”等不常食用或从未食用过的食物。

(12)

式中：

(13)

对(12)式求导并整理后，可以得到迭代后的坐标位置f _t+1 ：

对概率密度函数进行泰勒展开后化简，可以得到其近似表达：

(14)

通过限制迭代次数或设定迭代变化阈值结束迭代过程。设迭代结束后的位置为f _e，可以计算出新一帧目标的移动Δf 为

Δf =f _e-f ₀.

(15)

3 样本更新策略

为了适应气候更替、太阳角度变化等自然因素引起的背景区域缓慢变化，采用一种随机更新的策略更新样本库。具体过程如下：当某一帧的前v 帧被分类为背景点时，就有τ 的概率来更新样本库，在样本集中随机选择一个样本填充更新。τ 为样本库的更新概率，值越大、更新越及时。在选择要替换样本集中的样本时，随机选取一个样本值进行更新，可以保证样本值的平滑生命周期。设样本库大小为N ，由于是随机更新，一个样本值在时刻t 不被更新的概率是(N -1)/N . 设时间是连续的，则在dt 时间过去后样本值仍然保留的概率为

(16)

(16)式表明,在该策略下样本是否被替换的概率与时间因素无关。

4 模型建立与实验分析

4 .1 模型建立

首先选取目标与背景的相关性特征。目前没有解析方法可以指导特征的选取，通常是凭直觉和经验初步确定。一般情况下，人类视觉系统的目标识别过程主要有3类目标暴露特征，即统计特征、形状特征和纹理特征。考虑到以上因素，选取以下6种特征作为目标与背景区域相关性的描述，计算方法见(17)式～(26)式。其中：令x 、y 为两幅待计算的图像,和分别表示灰度均值；K 、W 分别为图像的长和宽。前3种为HSV空间中H通道、S通道和V通道的欧氏距离，如(17)式～(19)式所示，其中H 、S 、V 为目标图像的通道，H ′、S ′、V ′为背景图像的通道。

(17)

(18)

(19)

第4种为目标与背景区域的峰值信噪比^[10]：

(20)

(21)

第5种为感知哈希相似度，将图像缩小为8×8阶灰度矩阵后得到x ′_ij 和y ′_ij ，分别计算x ′_ij 和y ′_ij 的均值并获得感知哈希相似度的计算公式：

(22)

(23)

(24)

式中：⊕表示异或运算；b 函数的计算公式为

(25)

第6种为图像相关性运算^[11-12]：

(26)

根据上述模型所述，编制相应的MATLAB程序并进行计算验证。算法的特征参数提取流程图如图3所示。读入训练数据后，选出目标区域，得到背景8联通域(若背景处于边缘地区导致无法获取完整的8联通域，则丢掉此帧)。一边跟踪一边计算特征，得到均值方差并建立归一化联合分布。效果评估过程的算法流程图如图4所示。读入待检测数据并获取8联通域，将其特征代入归一化联合分布函数，计算联合概率密度。由于其值波动范围为0～6.94×10^-20，不利于阈值的设定。为增强密度值的灵敏度，对结果先放大10倍，再取对数，具体过程如下：

(27)

图3 联合分布建立过程流程图
Fig.3 Flow chart of joint distribution process

图4 伪装效果评估过程流程图
Fig.4 Flow chart of camouflage effect evaluation process

通过对数放大概率限定数值变化范围为-∞～66.41. 由此可以看出，数据的变化范围明显变大，能够有效减少计算误差。为了有效区分不同的伪装效果状态，并考虑到高斯分布的3σ 准则，将结果划分为3个伪装等级，分别为1级伪装(r ≥0)、2级伪装(0>r ≥-1 000)、3级伪装(r <-1 000)。

4 .2 实验过程与结果分析

实验过程于2018年3月中旬在南京市南郊某地进行，对某一指挥车实施连续的空中成像采集。飞行器高度50 m左右，飞行条件分别选取早晨、中午、下午、晴天、阴天等时间段或天候。将伪装完成后采集的数据作为训练样本(见图5)，设置1级伪装、2级伪装和3级伪装3种伪装状态测试模型的区分能力。随机地从不同天候、时间段等实验条件下抽取25帧3类数据作以展示。

计算并绘制的样本特征范围统计直方图如表1所示。表1中横向对应8个背景区域、纵向对应6种特征，共计48幅直方图。表2所示为特征的均值与方差，即联合高斯分布的均值与方差。综合表1和表2可见，由于数据量较少，其分布规律不是很明显，但是基本符合高斯分布的曲线特征：特征数据总体上分布得比较紧凑，没有分散现象，表明目标与区域背景的特征数据分布规律能够用高斯分布近似地表达。

图5 实验数据中的3种伪装状态
Fig.5 Three camouflage states in the experimental data

表1 样本特征统计直方图表

Tab.1 Feature statistical histogram of the samples

表2 样本特征的均值与方差统计表

Tab.2 Mean and variance table of sample features

图6 3种伪装状态的概率密度曲线图
Fig.6 Probability density curves of three camouflage states

通过分别计算图5所示的3类待检测数据，得到图6所示的概率密度曲线图。对这3条曲线进行数学统计，得到表3所示的结果。从图6和表3中可以看出：1级伪装的状态下曲线均值为54.475 8；2级伪装的状态下曲线均值为-707.033 6；3级伪装的状态下曲线均值为-1 111.226 2. 曲线的均值数据很好地反映了其伪装状态。曲线的波动较平稳，方差和极差值较稳定，表明该模型能够在动态检测过程中平稳运行。从1级伪装到2级伪装，改变了伪装网布置方位并去掉了上面布置的饰物，但是曲线已经接近3级伪装的阈值。从而表明样本数据对于数据空间的表达不够充分，下一步应合理补充样本数据。总之，基于特征统计的伪装效果评估模型对于目标伪装状态的变化能够给予准确的反映。

耳鸣程度分级：参照世界中医药学联合会中医耳鼻咽喉口腔科专业标准审定委员会及中华中医药学会耳鼻咽喉科分会推荐的《耳鸣严重程度评估与疗效评定参考标准》（2007，青岛）[7]根据以上耳鸣严重程度分级评定疗效。在治疗前和治疗后分别对患者的耳鸣严重程度进行评分，根据总分值对耳鸣严重程度进行分级（I～V级），再根据治疗前后耳鸣程度的级别改变判断为痊愈、显效、有效、无效。

表3 3种伪装状态概率密度曲线统计数据表

Tab.3 Statistical data of probability density curves of three camouflage states

5 结论

本文基于目标跟踪和特征统计技术，提出了一种动态伪装效果评估方法。通过统计伪装良好状态下目标与8联通域背景特征值的均值和方差，首先建立多维特征联合分布模型，然后计算对数放大概率密度得出目标的伪装效果。实验中对某一指挥车分别实施1级、2级和3级伪装后采集图像进行评估，根据高斯分布的3σ 准则，模型能够准确对所采集的伪装图像进行分级。相比于单幅图像的静态评估方法，基于特征统计的伪装效果评估模型能够从多帧图像、多个角度对目标实施评估，能够在一定范围内适应拍摄时的光照、天气等因素的影响。模型已初步应用于有关于军事巡查任务的研究课题之中。下一步需要分析不同特征算法对伪装分级的影响，使其更加符合视觉规律。

参考文献

[1] 许卫东, 王向伟. 伪装检测与效果评估理论和技术[M]. 北京:国防大学出版社, 2015:78-80．

XU W D, WANG X W. Camouflage detection and evaluation theory and technology[M]. Beijing: National Defense University Press, 2015:78-80. (in Chinese)

[2] 吕绪良, 林伟, 许卫东, 等. 应用迭代自组织数据分析技术模糊聚类法划分伪装效能等级[J]. 兵工学报, 2005, 26(5):681-684.

LÜ X L, LIN W, XU W D, et al. Applying the fuzzy clustering analysis of ISODATA to the classification of camouflage effectiveness[J]. Acta Armamentarii, 2005, 26(5):681-684. (in Chinese)

[3] 林伟, 陈玉华, 王吉远，等. 基于图像特征与心理感知量的伪装效果评价方法[J]. 兵工学报, 2013, 34(4):412-417.

LIN W, CHEN Y H, WANG J Y, et al. Camouflage assessment method based on image features and psychological perception quantity[J]. Acta Armamentarii, 2013, 34(4):412-417. (in Chinese)

[4] 于金, 朱立凡, 杜海璐. 基于BP神经网络的光学伪装效果评估模型[J]. 舰船电子对抗, 2009, 32(6):55-57.

YU J, ZHU L F, DU H L. Evaluation model of optical camouflage effect based on BP neural network[J]. Shipboard Electronic Countermeasure, 2009, 32(6):55-57. (in Chinese)

[5] QIN J F, QU L Y, ZHU L Y, et al. Optical camouflage effect objective evaluation method research under the condition of complex backgrounds[J]. MATEC Web of Conference, 2016, 61(10):1-4.

[6] 崔宝生, 薛士强, 姬艳，等. 图像特征的伪装效果评估技术[J]. 红外与激光工程, 2010, 39(6):1178-1183.

CUI B S, XUE S Q, JI Y, et al. Camouflage effectiveness evaluation based on image feature[J]. Infrared and Laser Engineering, 2010, 39(6):1178-1183.(in Chinese)

[7] LU X L, JIA Q, XU W D, et al. Research on camouflage assessment psychological impact of different pattern paintings based on the identification probability model[J]. Basic & Clinical Pharmacology & Toxicology, 2016, 118(36):30-40.

[8] HARITAOGLU I, HARWOOD D, DAVID L S. Real-time surveillance of people and their activities[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2003, 22(8):809-830.

[9] FUKUNAGE K, HOSTETLER L D. The estimation of the gradient of a density function with application in pattern recognition[J]. IEEE Transactions on Information Theory, 1975, 21(1):32-40.

[10] LIN C J, CHANG C C, LIU B S, et al. Developing and evaluating a target-background similarity metric for camouflage detection[J]. PloS ONE, 2014, 9(2):1-11.

[11] XUE F, WU F, WANG J W, et al. Camouflage texture design based on its camouflage performance evaluation[J]. Neurocomputing, 2018, 274(24):106-114.

[12] 杨骏堂, 许卫东, 曲杨， 等. 一种改进的Surendra运动目标伪装效果检测方法[J]. 兵工学报, 2017, 38(1):190-194.

YANG J T, XU W D, QU Y, et al. A Surendra-based improved detection method of moving target camouflage effect[J]. Acta Armamentarii, 2017, 38(1):190-194.(in Chinese)

A Dynamic Camouflage Effect Evaluation Method Based on Feature Statistics

YANG Xin, XU Weidong, JIA Qi

(Field Engineering College, Army Engineering University of PLA, Nanjing 210007, Jiangsu, China)

Abstract : The existing camouflage effect evaluation is mainly for a single still image, which can not simulate the process of the reconnaissance personnel’s interpretation of a target. A feature statistics-based dynamic camouflage effect evaluation method is proposed based on mean shift target tracking algorithm. The proposed method is to establish a normalized joint Gaussian distribution by using the data of correlation features between the target and the background of eight-way domain, and use the distribution range of probability density to evaluate the camouflage effect of target. A logarithmic amplification probability is proposed for calculating the probability density of joint distribution, which solves the problem that the high-dimensional joint distribution probability density has low numerical sensitivity and is inconvenient to set a threshold. At the same time, a sample update strategy is introduced to make the sample library update randomly according to a certain probability, so as to better adapt to the change in large-scale background caused by the turn of seasons and other factors. In the experimental process, the first-level, second-level and third-level camouflages are applied to a certain command vehicle. After collecting the data, the logarithmic amplification probability is calculated, and the statistics on the curves is made. The results show that 3 camouflages divided in reality completely correspond to the first-level，second-level and third-level camouflages pre-set by the 3σ criterion. The experimental results show that the model can effectively reflect the camouflage effect of target.

Keywords : camouflage effect； effect evaluation; feature statistics; probability density; Gaussian distribution; target tracking

中图分类号： E951.4

文献标志码： A

文章编号： 1000-1093(2019)08-1693-07

DOI :10.3969/j.issn.1000-1093.2019.08.019

收稿日期： 2018-09-14

基金项目： 江苏省自然科学基金项目(BK20180579)

作者简介： 杨鑫(1996—)，男，硕士研究生。E-mail:1435227062@qq.com

通信作者： 许卫东(1966—)，男，教授，博士生导师。E-mail: xweibing1968@aliyun.com

标签：伪装效果论文;  效果评估论文;  特征统计论文;  概率密度论文;  高斯分布论文;  目标跟踪论文;  陆军工程大学野战工程学院论文;