空间人造目标混合光谱材料数目确定方法仿真研究

李庆波，苗兴晋

北京航空航天大学仪器科学与光电工程学院，精密光机一体化技术教育部重点实验室，北京 100191

引 言

随着航天领域科学技术的不断发展，空间人造目标的观测与识别具有重要的研究意义和价值。空间人造目标主要包括大气外绕地运行的各种在轨工作或者停止工作的航天器、人造卫星等。空间人造目标光谱观测是利用光谱仪获取目标的光谱信息，可以根据目标光谱特征来识别分析目标的材料信息。光谱观测系统灵活、方便，对空间人造目标观测具有重要意义。由于观测距离和观测设备空间分辨率的限制，空间人造目标光谱通常是一个“混合光谱”，为空间人造目标的某个瞬时视场中多种纯物质材料的光谱特征的组合[2]。因此，在利用光谱信息分析空间人造目标材料组成和特性之前，需要确定混合光谱中包含的纯物质材料的数目，即，“端元数目”。尽管端元数目确定方法受到了广泛关注和研究，但其仍然是一个具有挑战性的研究领域。

目前，针对混合光谱中端元数目的确定，很多算法被开发出来。最著名的是一种基于“虚拟维度”(virtual dimension，VD)的方法[3]，称为噪声白化(noise whitening)的Harsanyi-Farrand-Chang(NWHFC)方法。其依赖于Neyman-Pearson检测理论，根据一个输入的虚警概率进行决策。该方法的缺点是需要人工输入一个错误警报参数，最终结果对该参数非常敏感。另一种方法是最小误差光谱信号子空间识别方法(hyperspectral signal subspace identification by minimum error，HySime)[4]。该方法不需要输入参数，在满足最小均方根误差的条件下，确定包含数据有效信息的最佳线性子空间的维数。但是该方法在信噪比较低、输入光谱数据维度较高时会有较大误差。此外，有学者提出了一种基于随机矩阵理论(random matrix theory，RMT)的方法[5]。该方法基于RMT的最新研究，根据一个理论阈值以区分信号协方差矩阵中由噪声引起的特征值和由信号引起的特征值。但是该方法需要一个由光谱数据统计分布信息所估计的阈值，当该阈值估计误差较大时，算法结果的误差也大。因此，本文提出了一种鲁棒特征值极大似然(robust eigenvalue maximum likelihood，REML)方法实现空间人造目标混合光谱端元数目的确定，该方法能够有效抑制高斯白噪声和频谱相关噪声的干扰，且不依赖于任何输入参数和固定阈值。

1 实验部分

1.1 原理与算法

1.1.1 数据预处理

传统的端元数目确定方法假定噪声为独立同分布的高斯噪声，这对于实际的空间人造目标光谱是不正确的，其往往会受到频谱相关噪声的干扰。因此，需要一个预处理步骤——噪声特性估计和噪声白化，来抑制频谱相关噪声的干扰。本文采用基于多元回归的方法来估计噪声协方差矩阵[6-7]，该方法当存在非高斯白噪声时，效果较好。接下来，使用一个噪声白化步骤将光谱数据转换到一个新的主成分空间来实现抑制频谱相关噪声的影响，具体方法如下：

F=PTKP

(1)

(2)

(3)

则经过上述变换后，原始光谱数据即被转换到一个新的主成分空间中，转换光谱数据由频带不相关的噪声组成，且该方法允许在噪声协方差矩阵估计过程中存在一定误差。

1.1.2 特征值极大似然方法

l=1,2,…,L

(4)

(5)

上述统计分布对应的离散的对数联合似然函数见式(6)

(6)

该对数联合似然函数将会在l=M+1时取得最大值，则混合光谱所包含的端元数目的估计值见式(7)

(7)

因此，该方法不需要任何输入参数和固定阈值即可完成空间人造目标混合光谱端元数目自动确定。

1.2 方法

1.2.1 仿真实验一

为了分析评价本文提出的空间人造目标混合光谱端元数目自动确定方法(简称REML)，使用实验室测得的五种空间人造目标常用的材料(白漆、太阳能电池板、聚酰亚胺、铝箔、外包覆)的可见/近红外光谱特征作为纯物质材料端元光谱进行仿真实验。实验装置示意图如图1所示，实验所使用的光源为Laser-1卤钨灯，光谱仪选用美国海洋光学公司的USB4000光纤光谱仪；使用QR400-7-VIS-BX型Y型光纤探头测量材料的反射光谱，测量波长范围是400～1 000 nm，光谱波长分辨力为0.2 nm。测量参数设置为：积分时间12ms，平滑度10，平均次数16，每条光谱数据包含3 091个波段，实验测得光谱结果如图2所示。

图2 五种空间人造目标材料的光谱特性Fig.2 Spectral signatures of five artificial space target materials

在获取到空间人造目标纯物质材料实验室光谱特性后，根据迪利克雷(Dirichlet)分布函数生成相应比例参数，将这五种材料光谱按该比例参数进行混合，从而形成空间人造目标混合光谱仿真数据集，共4 000条混合光谱。对该数据集分别添加20 dB，30 dB的高斯白噪声和频谱相关噪声，以模拟实际光谱获取过程中所受到的噪声干扰。将本文提出的REML方法和一些常用的空间人造目标端元数目确定方法：噪声白化的NWHFC方法[3]，HySime方法[4]，特征间隔法(noise-whitened eigengap approach，NWEGA)[5]，进行比较，为了探究算法的准确性和稳定性，实验中每种算法均重复运行50次，结果取50次运行的平均值，保留两位小数。

1.2.2 仿真实验二

由于未知空间人造目标的组成材料种类和数目多变，并且不完全已知，因此，需要进一步探究本文所提出的方法在空间人造目标混合光谱端元数目变化的鲁棒性。本实验使用美国地质勘测局(United States Geological Survey，USGS)光谱数据库中随机选择3，5，7和9条光谱数据作为纯物质材料光谱构建混合光谱仿真数据集，其包含224个波段(波长范围是370～2 500 nm)，共5 000条混合光谱。为了模仿数据采集过程，将具有不同信噪比的频谱相关噪声分别添加到仿真混合光谱数据中。同样，将本文提出的REML方法和一些常用的空间人造目标端元数目确定方法(NWHFC，HySime以及NWEGA)进行比较，实验中每个方法均重复运行50次，结果取50次运行的平均值，保留两位小数。

2 结果与讨论

2.1 第一组仿真实验的结果

表1和表2显示了NWHFC，HySime，NWEGA以及本文提出的REML算法在不同信噪比、不同类型噪声干扰下对由五种空间人造目标材料生成的混合光谱数据的端元数目的估计结果。其中NWHFC_1，NWHFC_2 和NWHFC_3分别表示输入虚预警参数分别为{10-3,10-4,10-5}的NWHFC算法。图3表示在信噪比为20dB时不同类型的噪声干扰下，REML算法的对数联合似然函数值关于特征值指数的变化。

图3 信噪比为20 dB时不同类型噪声干扰时REML的对数联合似然函数值Fig.3 Logarithmic joint likelihood function values of REML for different types of noise (SNR: 20 dB)

表1 高斯白噪声干扰下各个算法实验结果比较Table 1 Comparison of experimental results of various algorithms under the Gaussian white noise

表2 频谱相关噪声干扰下各个算法实验结果比较Table 2 Comparison of experimental results of various algorithms under the spectrally correlated noise

从表中实验结果可以看出，对于不同种类和信噪比的噪声，HySime算法的实验结果出现较大误差。REML算法和NWHFC算法对端元数目的确定均取得了较好的实验结果。但是，当信噪比较低时，NWHFC算法的效果有所下降；对于NWEGA算法而言，其运行结果的稳定性较差，这可能是由于其存在一个基于数据统计信息(易受噪声干扰)的硬阈值导致了计算结果的波动性。另外，从图3中可以看出，REML中的对数联合似然函数的值在不同类型噪声干扰下均能够在正确的特征值指数下取得最大值。因此，本文所采

表3 端元数目不同时频谱相关噪声干扰下各个算法实验结果比较Table 3 Comparison of experimental results of various algorithms under the spectrally correlated noise with different numbers of endmembers

用的REML算法则展现出对频谱相关噪声和白噪声干扰均有较强的鲁棒性。

2.2 第二组仿真实验的结果

第二组实验旨在探究不同算法对混合光谱端元数目变化的鲁棒性，实验结果如表3所示，为噪声干扰为频谱相关噪声的实验结果。

就端元数目变化对算法的影响而言，当混合光谱包含的纯物质材料数目变多时，NWHFC出现了低估现象，NWEGA存在高估和低估现象，表明这两种算法难以区分出不同信号源的差异性。当虚预警参数发生变化时，NWHFC算法的实验结果也发生变化，这也是该算法一个局限性。另外，当噪声信号存在频谱相关性时，HySime算法出现了严重的高估现象，完全无法使用。REML算法的实验结果具有很好的准确性和稳定性，展现出对端元数目变化较强的鲁棒性。

3 结 论

提出了一种用于确定空间人造目标混合光谱纯物质材料数目的全自动方法，该方法可以有效解决空间人造目标光谱观测过程中所存在的频谱相关噪声的干扰问题。其优势和特点如下：首先，该方法无输入参数，且不基于固定阈值；其次，当光谱信号受到白噪声和频谱相关噪声的干扰时，由于利用了一个有效且允许有一定估计误差的噪声特性估计和噪声白化方法，能够有效抑制噪声对计算结果的干扰。实验结果表明，该方法能够准确且稳定地确定空间人造目标混合光谱中纯物质材料的数目。因此，该方法在利用光谱信息进行空间人造目标材料分析领域具有广阔应用前景。