ENVI系列教程（十九）——目标探测与识别
目录
1 概述
2 去伪装目标探测
2.1 数据浏览与准备目标波谱
2.2 目标探测流程化工具
2.2.1 选择输入 / 输出文件
2.2.2 大气校正
2.2.3 输入目标波谱
2.2.4 输入背景波谱
2.2.5 执行 MNF 变换
2.2.6 选择及运行分析方法
2.2.7 浏览结果以及提取目标
2.2.8 结果后处理
2.2.9 输出结果
3 基于 BandMax 向导的 SAM 目标探测
3.1 工具功能介绍
3.2 详细操作流程
3.2.1 数据打开与目视解译
3.2.2 收集目标与背景波谱
3.2.3 启动 SAM Target Finder with BandMax
1 概述
高光谱图像分类方法与传统的多光谱分类有本质的区别，从高光谱图像的每个像元均可以获取一条连续的波谱曲线，就可以考虑用已知的波谱曲线和图上每个像元获取的波谱曲线进行对比，理想情况下，如果两条波谱曲线一样，就能说明这个像元是哪种物质。
高光谱图像的这个特性，使它除了应用于一般的图像分类，还应用于物质识别、目标探测等。图像分类更多关注的是地物覆盖和物质成分，目标识别和探测是对特定对象的搜索，其结果是 “有” 或者 “没有”。因此我们把高光谱图像分类、物质识别、探测等称为波谱识别。
本节以两个专题为例（去伪装目标探测与基于 BandMax 向导工具的 SAM 目标探测），学习 ENVI 中提供的高光谱目标探测与识别功能。
2 去伪装目标探测
去伪装目标探测也是利用高光谱图像的地物识别能力，从图像上探测遮掩或者伪装的目标，比如一种特殊物质、矿物甚至军事目标等。下面以从图像上探测一个目标为例（练习数据存放在 “..\18. 目标探测与识别 \ 数据 \1. 目标探测” 中），介绍 ENVI 的 Target Detection Wizard 工具的操作流程。
在本例中，示例数据是包含 384 个波段，波段覆盖 382nm~2500nm 的高光谱数据。主要过程如下：从图上目视解译一个目标（可以是多个目标），以这个目标的平均波谱作为参考，搜索整个图像，识别具有类似或者相同波谱的目标。
2.1 数据浏览与准备目标波谱
（1）选择菜单 File > Open，打开文件 nvis_sub1_hsi.img；注：可以采用 RGB 彩色显示模式，在 Data Manager 中选择波段 54、34、14 分别对应 RGB 三个通道，然后点击 Load Data 按钮可以近似真彩色的显示原始图像，如下图。
（2）在右边图层管理器中nvis_sub1_hsi.img右键选择New Region Of Interest打开ROI Tool，在 ROI Tool 中修改 ROI Name：target 1。
（3）在 ENVI 主界面中的工具栏 Goto 输入：130,650，回车定位到一辆装甲车目标。绘制这个目标的一个多边形 ROI。
注：练习数据中提供了目标 ROI 文件。其中 target1_sub1_roi.roi 是 ENVI Classic 或者 ENVI4.8 及更早版本的 ROI 文件，target1_sub1_roi.xml 是 ENVI5.1 或者更高版本的 ROI 文件。选择菜单 File > Open，可以打开 ROI 文件。
（4）浏览目标波谱。在 Layer Manager 中，右键点击 target 1，选择 Statistics，会对 ROI 的波谱进行统计，结果如下图所示。
2.2 目标探测流程化工具
2.2.1 选择输入 / 输出文件
（1）启动 Target Detection Wizard 工具，位于 Toolbox/Target Detection/Target Detection Wizard，在弹出的面板左侧为整个流程的描述信息，单击 Next 按钮。
（2）单击 Select Input File，选择高光谱数据 nvis_sub1_hsi.img。
（3）单击 Select Output Root Name，选择输出结果的根目录，如 D:\Temp\nvis_sub1_hsi。
（4）单击 Next 按钮进入 “Atmospheric Correction” 面板。
2.2.2 大气校正
此数据是经过大气校正的，选择 None / Already Corrected 选项，单击 Next 进入 “Select Target Spectra” 面板。
注：这个面板提供了快速大气校正（QUAC）、几个基于统计学的大气校正方法（IAR Reflectance、Log Residuals、Flat Field、Empirical Line）和简化的黑暗像元法（Dark Subtraction），如果输入数据未经过大气校正，用户可根据需要进行方法的选择。
2.2.3 输入目标波谱
在 “Select Target Spectra” 面板选择 Import > From ROI/EVF from input file，选择列表中的 target 1，点击 OK，选择后效果如下图所示。单击 Select All，选中目标波谱，然后单击 Next 按钮进入 “Select Non-Target Spectra” 面板。
注：① 如果需要探测多个目标，则输入多个目标的波谱，单击 Select All 选择列表中所有目标波，单击 Next 执行下一步操作。② 当需要应用 Orthogonal Subspace Projection（OSP）、Target-Constrained Interference-Minimized Filter（TCIMF）和 Mixture Tuned Target-Constrained Interference-Minimized Filter（MTTCIMF）三种波谱分析方法时，需要至少 2 个目标波谱，或提供背景波谱。
2.2.4 输入背景波谱
在 “Select Non - Target Spectra” 面板中可以选择易与目标波谱混淆的波谱作为背景波谱，有助于提高探测精度。这里选择 No（不选择背景波谱），单击 Next 进入 “Apply MNF Transform” 面板。
注：如果要使用 Orthogonal Subspace Projection（OSP）或 Target - Constrained Interference - Minimized Filter（TCIMF）方法，必须在第三步骤输入大于一种的目标波谱，或在第四步骤输入背景波谱。如果要使用 Mixture Tuned Target - Constrained Interference - Minimized Filter（MTTCIMF）分析方法，必须在第三步骤输入大于一种的目标波谱，或在第四步骤输入背景波谱，而且在下面的第五步要执行 MNF 变换。
2.2.5 执行 MNF 变换
（1）在 Apply MNF Transform 面板中，Apply MNF Transform？选择 Yes。
（2）单击 Show Advanced Options 按钮，默认选择全部的 MNF 波段。
（3）单击 Noise Stats Shift Diff Spatial Subset，默认选择全部图像区域用于统计噪声。
（4）单击 Next 执行 MNF 变换，计算完成后自动进入 “Target Detection Methods” 面板。
注：MNF 变换可以分离噪声，对数据降维以减少计算量。如果选择 No，那么将不能选择Mixture Tuned Matched Filter （MTMF）和 Mixture Tuned Target-Constrained Interference-Minimized Filter （MTTCIMF）识别方法。
2.2.6 选择及运行分析方法
在“Target Detection Methods” 面板选择 CEM、ACE 和 MTMF 三种方法（如图所示），单击 Next 按钮执行分析，之后自动进入“Load Rule Images and Preview Result”面板。
2.2.7 浏览结果以及提取目标
在 Load Rule Images and Preview Result 面板中，“Target” 列表中显示所有探测目标参考波谱，在 “Method” 列表中选择相应分析方法，其得到的规则自动显示在 Display 窗口中，选择规则阈值或者散点图上选择点云将目标分离。
（1）在 “Method” 列表中选择 CEM，规则阈值（Rule Threshold）：0.2。
（2）在 “Method” 列表中选择 ACE，规则阈值（Rule Threshold）：0.1。
（3）在 “Method” 列表中选择 MTMF，自动会生成一个 MF scores 和 infeasibility values 的散点图，选择高 MF scores 和低 infeasibility values 的点云，也就是散点图右下角部分点云，就是探测到的目标。
（4）单击 Next 执行从规则图像中分离目标，进入 “Filter Targets” 面板。
注：① 对于 MF, CEM, ACE, SAM, OSP, 和 TCIMF，ENVI 自动生成默认值阈值。当手动修改阈值时，调整阈值越小，得到的目标点越多，“假目标” 也随之增多，重新获取的目标点会在图上高亮显示。SAM 刚好相反。② 对于 MTTCIMF 和 MTMF，ENVI 自动生成 MF scores 和 Infeasibility values 的整个图像的散点图。用鼠标左键绘制多边形区域选择点云，鼠标右键结束选择。鼠标中建拉框可放大点云，单击中键回到上一个试图，同时之前选择的点云被取消，当选择错误时候用这个功能重新选择点云。散点图勾选区域如图所示，需要选择 Infeasibility values 较低，且右侧散落的点。
2.2.8 结果后处理
在 Filter Targets 面板中（如下图），按照默认设置单击 Next 按钮，进入 Export Results 面板。
注：这一步提供分类后处理的方法（Clumping 和 Sieving）用于去除结果中的小斑点。Clumping 是用卷积的方法定性去除小斑点；Sieving 是用定性的方法去除小斑点，通过设置最小聚类像素个数（Group Min Threshold）移除小斑点。
2.2.9 输出结果
在 Export Results 面板（如下图），可以将探测结果输出为感兴趣区（ROI）和矢量（Shapefile），按照默认设置（ROI）单击 Next 按钮，进入最后一个面板 “View Statistics and Report”，自动统计探测的结果，且探测的所有结果自动加载到 ROI Tool 中并显示在图上。定位探测结果（169，163）是一个树林掩盖的目标，点击 Finish 完成。
注：整个流程都可以使用 Prev 按钮回到上一步操作。
3 基于 BandMax 向导的 SAM 目标探测
3.1 工具功能介绍
基于 BandMax 向导的 SAM 目标探测工具可以引导我们完成高光谱影像的目标探测。向导的 BandMax 工具能帮助找到最佳的波谱子集从而区分背景和目标，并节省处理的时间。
向导主要有以下几个步骤：
（1）选择输入 / 输出文件。
（2）选择目标 — 选择目标波谱。
（3）选择背景 — 选择需要抑制的背景信息。
（4）利用 BandMax 计算有效波段 — 识别对 SAM 分析中有效的波段。
（5）选择最大角阈值 — 定义 SAM 最大角。
（6）检验制图成果 — SAM 分析以及成果检验。
如果处理结果理想，那么可以点击 Finish 退出向导。如果检验结果显示波段子集不充分，向导会返回到第二步，需要重新输入目标和背景波谱，然后利用 BandMax 选择最合适的波段子集以及利用 SAM 对输入数据重新分类。处理分析完成后，向导中将显示分析报告。如果需要，我们可以保存这个文件以备后用。向导特定步骤中得出的影像结果会显示在 ENVI 的 Data Manager 中。如果输入数据的质量不好或者设置了不合适的参数，很有可能得到不理想的结果。
向导的第一个面板主要介绍该向导的工作流程。向导中其他的面板主要介绍完成这个过程的 6 个步骤。这些处理步骤有时会因为分析方法的不同而反复出现。
3.2 详细操作流程
下面介绍详细操作流程，以 AVIRIS 高光谱数据为例（位于 “...\ 数据 \2. Target Finding With SAM and BandMax\ sandiego_reflectance.img”），此数据经过了大气校正，为地表反射率图像。
3.2.1 数据打开与目视解译
（1）启动 ENVI 5.1，打开 File > Open，选择 “...\ 数据 \2. Target Finding With SAM and BandMax\sandiego_reflectance.img” 文件，并加载到视图中，如图所示。
（2）从图中可以看出，影像中包含了一些飞机。绿色圆圈选中的三架飞机是我们的探测目标，像元位置为 [88, 10] 和 [142, 337] 等。另外最明显的红色圆圈选中的两架飞机将是我们需要抑制的背景，像元位置为 [231, 135] 和 [243, 143]。
注：像元坐标以 [0, 0] 为原点。
3.2.2 收集目标与背景波谱
（1）启动光谱曲线查看器，菜单为 Spectral > Profiles > Spectral，如下图所示。
（2）收集目标波谱。在 ENVI 主界面第一行工具栏最右侧的 “Go To” 框内手动输入 “88, 10”，然后按 Enter 键，将自动定位到目标像元，并在波谱查看器中显示目标波谱。
注：输入坐标不包含双引号，X 与 Y 坐标中间为英文状态下的逗号。
（3）在 Spectral Profile 面板中，选择 Options > New Window with Plots，打开一个新的 Plot 窗口，并保留了目标波谱，可以点击右侧中间三角形箭头显示出波谱列表。
（4）以同样的方法收集背景波谱。在第一行工具栏最右侧的 “Go To” 框内手动输入 “243, 143”，然后按 Enter 键，将自动定位到背景像元，波谱将自动显示在 Spectral Profile 面板中，点击右侧中间三角形箭头，可以显示波谱列表。
注：如果想收集多条光谱曲线，可以使用鼠标左键将波谱拖入到新建 Plot 窗口的波谱列表中，如下图所示，本例中目标和背景波谱均为一条曲线，可以不进行拖拽。
3.2.3 启动 SAM Target Finder with BandMax
（1）启动工具 Toolbox/Target Detection/SAM Target Finder with BandMax，工具第一步左侧文字介绍了工具的功能与操作步骤，点击 Next 进入下一步。
（2）点击 Select Input File，选择输入文件，即 “sandiego_reflectance.img”，ENVI 会自动设置输出文件夹路径，可以点击 Select Output Root Name 按钮修改输出路径。点击 Next 进入下一步。
（3）在 Target Spectra 面板，需要输入目标波谱。在表格空白处右键选择 [88, 10] 波谱曲线，将自动加载到波谱列表中，如下图所示。点击 Next 进入下一步。
（4）选择 Select Backgrounds to Reject 为 Yes，在下面用同样的方法将 [243 143] 输入作为背景波谱，点击 Next 进入下一步。
（5）在此步骤，将 Band Significance Threshold 设置为 0.2775，Number of Significant Bands 将变为 33，如下图所示，点击 Next 进入下一步。
注 1：ENVI 利用 BandMax 算法为输入文件的每一个波段计算特征值（Band Significance），范围 0~1，值越大表示目标与背景区分越好。BandMax 会选出特征值大于阈值（Band Significance Threshold）的波段进行目标探测。
注 2：默认的阈值为 0.2475，是根据大约 25% 的比例计算得出的。选出的波段数过大时，将不利于区分目标与背景。所以增加阈值到 0.2775，选出波段数为 33。用户可根据需要修改阈值或波段数，以获取不同的结果。
注 3：可以点击 Save Significant Bands to File 按钮，将保存波段子集到 txt 文本文件中。
（6）设置 SAM Maximum Angle 为 0.08，点击 Next 执行 SAM 分析。
注：降低阈值可以减少匹配出的像元数，可以更好的与目标波谱匹配。增大阈值，可能得到空间相干更好的结果。然后，阈值越低，整体像元匹配效果不一定会越好。
（7）在此步骤，可以点击按钮 Load SAM Class Result 加载 SAM 分类结果，如下图所示，红色区域为目标探测结果。
注 1：可以点击按钮 Load SAM Rule Image 查看规则图像，DN 值为像元波谱与目标波谱之间的波谱角。
注 2：如果目标没有被很好的区分，可以点击 Next 按钮返回到步骤 3。
注 3：点击 Show Report 按钮可以浏览整个过程的参数设置情况。
（8）如果结果理想，可以点击 Finish 完成整个流程。