对典型伪装材料的高光谱特征及其识别方法进行了研究

引言

随着材料技术的高速发展，目标伪装手段和策略也不断的提高，尤其是先进涂层和伪装网的使用，在某些背景环境下能达到目标和背景“异物同谱”的效果，增加了伪装识别的难度。基于传统的可见光或多光谱遥感技术，由于其单个波段覆盖光谱范围广，很难识别出伪装目标，给伪装识别带来了很大的挑战。

高光谱遥感是近些年来发展迅速的一种多维信息获取技术。它能够同时获取目标区域的2维几何空间信息和1维光谱信息，具有光谱范围宽、波段多、光谱分辨率高（纳米级）等特点，为了解地物提供了极其丰富的信息，也为地物无损快速的理化性质检测提供了手段。目标与背景之间存在的差异是目标暴露的根本原因，地物的光谱信息对目标识别具有指纹效应，通过分析地物的光谱特征并据此提取的特征指数，与物质的特性相关，可为精细识别地物提供依据，尤其在地表资源环境探测中，高光谱技术有效的提高了地物识别和分类的精度。

目前，有许多学者探究高光谱技术在伪装材料特性分析及识别中的应用。刘凯龙等以瑞典、瑞士的森林伪装网和伪装服及不同种类、季节和环境的绿色植被光谱特性为基本识别特征，提出数学分析模型及统计计算和判别效果检验相结合的光谱特征检测新方法，伪装目标识别正确率达到了99%以上，该研究主要在实验室的可控环境下，通过实验求得判别函数及相关系数，在室内应用时，都能得到较为满意的结果，但其不能完全适用于不同环境下采样的数据。刘志明等对比绿色伪装涂料和植被的可见光-近红外反射光谱及一阶微分谱的差异，提出了550和780~1300nm波段反射率波动性是星载和机载高光谱遥感伪装探测中的有效特征，为绿色伪装材料识别提供了一定理论方法；卢云龙等分析比较伪装材料光谱曲线，利用高光谱图像实现验证了基于光谱微分、光谱积分、光谱指数和特征波段的伪装目标检测算法，具有很好的检测效果；付严宇等通过高光谱实验及分析得到了可见光图像难以区分的三组伪装材料（真假草坪、绿色木质和钢制材料、丛林迷彩服和丛林雨衣）的光谱特征差异，对伪装识别的波段选取提供了依据。虽然许多学者针对高光谱技术在伪装目标识别中的应用进行了较为深入的研究，但目前基于光谱特征指数识别伪装目标的研究还较少，已有的特征指数和识别方法不适用于北方环境中的丛林迷彩伪装网的识别，无法有效区分伪装目标与背景目标。

针对北方地区大型目标的伪装识别问题，以北方典型丛林迷彩伪装网为研究对象，利用地物光谱仪测量其在300~2500nm的光谱曲线，实验分析不同湿度条件下伪装网的光谱变化规律和特征，进一步揭示伪装网与地带性典型植被之间的光谱差异和可探测波段，提出基于特征光谱的伪装迷彩识别策略，并通过仿真伪装环境下的高光谱图像，对识别方法进行验证。本研究成果对北方环境中丛林迷彩伪装网的识别具有一定的现实指导意义，同时为基于高光谱图像快速、准确识别伪装以及伪装材料的改进提供了理论和技术基础。

实验部分

2.1 室内可见光-近红外光谱实验

实验对象为我国北方地区常用丛林迷彩伪装网，如图1所示。伪装网光谱采集可以使用莱森光学地物光谱仪，该仪器的波段范围为300~2500nm，有1024个通道，测试环境为密闭黑暗的，采用功率为50W标准直流锡丝石英卤素灯作为光源。为降低样本光谱各向异性的影响，测量时转动样本3次，每次转动角度约90°，获取样本4个方向的光谱曲线，每个方向采集5次，每个样品共采集20条样本线，算术平均后得到该伪装材料的实际光谱反射率。对于传统绿色伪装材料而言，因不具有类似植被含水的组织结构，其在近红外波段的光谱特征会区别于植被，不能在近红外波段与植被“同谱”，水的存在与否直接影响着伪装材料在近红外波段的光谱响应。为更好的分析伪装材料的光谱特征及响应机理，模拟多环境对光谱的影响，因此对伪装材料部分样品进行浸水处理。实验设计以浸水时间5min为间隔，测量不同浸水时间下的伪装网的光谱曲线。伪装网表层于25~35min左右出现水分饱和现象，此时水分附着在伪装网表层不再被吸收，水分含量对光谱吸收特征的影响达到极限，所测得伪装网光谱曲线逐渐趋于一致。因此剔除了浸水时间超过30min的光谱曲线，最终选择具有显著光谱差异的浸水时间（0，5，10，20和30min）进行伪装网高光谱特征分析。

图1丛林迷彩伪装网

2.2成像光谱伪装模拟

为了对丛林迷彩伪装网的有效识别波段进行验证，利用包含植被、砂砾、岩石、车辆等地物的仿真模型，以植被为背景目标，通过丛林迷彩伪装网对车辆进行伪装，并使用高光谱相机对伪装前后的模拟环境进行高光谱成像实验，如图2所示。

图2成像光谱实验和模拟伪装环境

左：无伪装环境；右：模拟伪装环境

2.3光谱预处理方法

2.3.1光谱相似性测度

光谱角余弦（SAC）将光谱数据看作一个n维的特征向量，通过计算参考光谱和目标光谱之间的“角度余弦”来确定两者之间的相似性，它能够同时反映光谱向量在数值和形状之间的差异。光谱角为两个光谱向量的夹角，夹角越小，SAC就越接近于1，两个光谱向量相似度越高。光谱角余弦的计算方法如式（1）所示

式（1）中：cos（x，y）为光谱角余弦值，x和y分别为参考光谱向量和目标光谱向量。

2.3.2连续统去除法

连续统去除法是将光谱曲线吸收强烈部分的波段特征进行转换，使光谱吸收特征归一化到共同基线基础上，即将反射率归一化为0~1.0，从而有利于进行光谱吸收特征的分析，并与其他光谱曲线进行特征比较，提取特征波段。计算方法如式（2）

式（2）中：λj是第j波段；RCj为波段j的包络线去除值，Rj是波段j的原始光谱反射率。λend和λstart是在吸收曲线的起始点和末端点的原始光谱反射率。λend和λstart是吸收曲线的起始点波长和末端点波长；K是吸收曲线起始点波段和末端点波段之间的斜率。

2.4特征提取

光谱指数法是将特定波段或特定波段处的反射率值等进行某种形式的组合。选取波段时应根据一定的分析基础作为参考，以包含地物波谱特性机理，可分性度量大为原则。如常见的归一化植被指数NDVI、比值植被指数RVI和差值植被指数DVI均能够与植被建立相关性，是植被生长状态及植被覆盖度的最佳指示因子。考虑到伪装网一般存在于植被覆盖度较高的区域，在此环境下比值指数具有更好的敏感性，因此选择比值光谱指数RI建立识别模型。如式（3）所示，其中R1，R2为波段1和波段2的反射率值。

结果与讨论

3.1 伪装材料及目标背景光谱特征

根据研究背景环境植被的类型，从ENVI标准波谱库中选取了针叶树（Conifer）、落叶树（Decidous）和草本（Grass）这三种植被类型的光谱曲线，以及北方地区针阔叶混交林中典型植被山杨、松树、冷杉光谱曲线。图3（a）是浸水前后迷彩伪装网和6种植被的反射光谱对比。可以发现不同浸水时间伪装网的光谱曲线的基本形态相似，可见光部分反射率较低，在400~480和660nm处的蓝、红光呈低谷，在530nm处有一个反射峰，在680~760nm波段反射率急剧增大，几乎增加到最大值，随后在760~1000nm基本保持不变。进入近红外波段后，反射率开始下降，在1200~1750和1850~2100nm附近处有不同程度的吸收谷特征，主要为水分子吸收带。浸水前后伪装网的返射光谱仅在1900nm附近存在差异性：当浸水时间较短时，伪装网反射光谱即呈现出明显的吸收特征；随着浸水时间的增加，反射率缓慢降低，光谱曲线整体呈下降趋势。

图3干燥/浸水状态迷彩伪装网与6种植被的平均光谱曲线（a）和1800~2250nm处干燥/浸水状态迷彩伪装网与6种植被的光谱曲线（b）

对比伪装网和各植被的光谱曲线，在可见光波段干燥和浸水状态下的伪装网都与植被光谱具有较好的相似性；在900~1600nm处，伪装网的吸收特征与典型绿色植被的吸收特征之间存在明显差异；对比不同状态下伪装网和各植被在1800~2250nm处的反射光谱，如图3（b）所示，可以发现伪装网在该波段的吸收特征受水分影响较大，分别计算干湿状态下伪装网和各植被的光谱角余弦值，干燥状态下的伪装网和各植被的相似度均值为0.895，而浸水状态下的伪装网与各植被的相似度均值上升到了0.939，表明伪装网在该波段处的吸收特征会随着浸水处理而相似于植被的光谱吸收特征。值得注意的是，在植被光谱中，1400nm波段附近是典型的水分子中O—H振动所导致的吸收特征波段。但在本工作中，由于伪装网材料本身特性与植被存在差异，干燥状态下的伪装网在1250~1650nm范围内出现了较宽的反射率低谷，随着浸水时长的改变，没有因水分含量变化产生明显的吸收特征。该特征在其他伪装材料的相关研究中也有所体现。

“绿色反射峰”、“红边”、“近红外高原”和“水分吸收带”是绿色植被反射光谱的共同特征，为探测绿色植被环境中的光谱异常提供了基本的光谱参量。针对上述4个特征参量波段，分别计算迷彩伪装网（干燥）和各植被的光谱角余弦，如表1所示。可以发现，在550nm波段的“绿色反射峰”和680~780nm波段的“红边”处，迷彩伪装网（干燥）和植被的反射光谱具有较好的匹配性，计算相似度平均值分别为0.9905和0.9628。在1300~2500nm波段的“水分吸收带”处，迷彩伪装网（干燥）和植被的光谱匹配较差，计算相似度平均值仅有0.7357，差异性较为明显。然而在780~1300nm波段的“近红外高原”处，迷彩伪装网（干燥）和植被的反射光谱相似度均在0.98以上。这是因为光谱角余弦是基于光谱曲线之间的几何测度，更多的是从方向上区分差异，而对绝对的数值不敏感，虽然在这个波段范围两者相似度较高，但是仍存在差异性。

表1典型特征参量波段处植被与伪装网（干燥）的光谱角余弦值

为有效识别迷彩伪装网，进一步分析差异波段范围的识别特征，对反射光谱进行包络线去除处理，处理后的光谱曲线如图4所示。相较于植被，伪装网（干燥）在970和1190nm处无明显吸收特征，在1440nm处具有较宽的吸收谷，两者区分度较大。在900~1300nm波段，伪装网光谱曲线急速下降，在1300~1600nm波段，伪装网的光谱曲线处于较宽的吸收谷底，波动平稳，斜率偏小；而植被的光谱曲线波动较大，在1150~1300nm和1300~1440nm波段分别存在陡升、陡降现象，斜率较大。

图4迷彩伪装网和各植被的包络线去除结果

通过以上分析可知，迷彩伪装网（干燥）和植被在可见光波段反射光谱相似，二者不易区分；在近红外波段，970，1190和1440nm附近处，迷彩伪装网（干燥）和植被反射光谱差异明显，可用于绿色植被环境中迷彩伪装网的识别。

3.2伪装材料光谱特征指数及识别策略

通过对迷彩伪装网和绿色植被的波谱特征分析，迷彩伪装网和各植被的光谱在900~1900nm波段范围内存在几处反射率斜率的明显差异。据此为了有效识别迷彩伪装网，选取吸收特征差异较为明显的1150~1440nm波段，尝试建立比值伪装指数（RCI）

式（4）中：R1190，R1270和R1440分别是1190，1270和1440nm波段处的反射率。其中，1190和1270nm波段的选取分别对应1150~1300nm范围内与伪装网有差异的植被吸收谷的波谷及右肩波峰处的特征波段；1270和1440nm波段的选取分别对应1300~1440nm范围内与伪装网有差异的植被吸收谷的左肩波峰及波谷处的特征波段。分别计算迷彩伪装网和各植被的RCI值，如表2所示。在1150~1300nm波段伪装网的反射率斜率均大于1.4，在1300~1440nm 波段伪装网的反射率斜率均小于1.4，且都与植被有所区别，如果能同时满足这两个条件可把伪装网从绿色环境中识别出来，故此设置识别阈值条件为：RCI1>1.4；RCI2<1.4。

表2迷彩伪装网和各植被的RCI值

为了进一步验证指数的实用性，将其应用到高光谱图像，对图像内的迷彩伪装网进行识别划分，结果如图5所示。在识别提取结果中白色代表迷彩伪装网，黑色部分代表其他背景地物。从结果可以看出，所提取的伪装网在形态和大小上基本与原始影像一致，利用交并比计算实验结果的识别检测精度，可达到0.95，说明识别效果较为良好，分类准确率较高。

图5基于RCI指数构建识别策略及识别检测结果

讨论

针对北方地区典型丛林迷彩伪装网，设计了室内高光谱观测实验以及高光谱成像观测实验，对伪装网和北方地区地带性典型植被的反射光谱特征进行了深入分析，并提取伪装网与植被的差异性光谱特征波段，构建光谱比值指数RCI，用于绿色环境中伪装目标的识别，得到以下结论：

（1）不同浸水时间的丛林迷彩伪装网的光谱曲线基本形态相似，且反射率随浸水时间的增加而整体呈下降趋势；1900nm波段是伪装网反射光谱对含水量响应最为明显的波段。

（2）丛林迷彩伪装网和植被在可见光波段的波动情况相似，相似度均高于0.95，但在近红外波段光谱特征差异明显，尤其在1300~2500nm波段，相似度平均值只有0.7357。实验分析得出970，1190和1440nm波段附近处是典型伪装网成像识别的特征波段，伪装网在1300~1600nm波段光谱曲线波动平稳，斜率偏小；而植被的光谱曲线波动较大，在1150~1300和1300~1440nm波段分别存在陡升、陡降现象，斜率较大。

（3）基于丛林迷彩伪装网和植被的光谱特征差异构建了RCI指数，用于基于高光谱图像的绿色环境中伪装目标的识别。室内成像实验结果表明，该指数识别伪装目标的准确率达到95%，可以很好的将伪装网识别出来。

推荐：

便携式地物光谱仪iSpecField-NIR/WNIR

专门用于野外遥感测量、土壤环境、矿物地质勘探等领域的最新明星产品，由于其操作灵活、便携方便、光谱测试速度快、光谱数据准确是一款真正意义上便携式地物光谱仪。

无人机机载高光谱成像系统iSpecHyper-VM100

一款基于小型多旋翼无人机机载高光谱成像系统，该系统由高光谱成像相机、稳定云台、机载控制与数据采集模块、机载供电模块等部分组成。无人机机载高光谱成像系统通过独特的内置式或外部扫描和稳定控制，有效地解决了在微型无人机搭载推扫式高光谱照相机时，由于振动引起的图像质量较差的问题，并具备较高的光谱分辨率和良好的成像性能。

便携式高光谱成像系统iSpecHyper-VS1000

专门用于公安刑侦、物证鉴定、医学医疗、精准农业、矿物地质勘探等领域的最新产品，主要优势具有体积小、帧率高、高光谱分辨率高、高像质等性价比特点采用了透射光栅内推扫原理高光谱成像，系统集成高性能数据采集与分析处理系统，高速USB3.0接口传输，全靶面高成像质量光学设计，物镜接口为标准C-Mount，可根据用户需求更换物镜。

审核编辑 黄宇