主成分分析用正交方向压缩遥感波段的相关变化

说明协方差、特征向量和主成分怎样重组多波段信息,并分析标准化、异常值与物理解释的限制。

主成分分析用正交方向压缩遥感波段的相关变化

多光谱和高光谱影像的波段之间通常存在较强相关性。主成分分析通过线性变换,把原始波段转换为一组互不相关的主成分,并按照样本方差从大到小排列。它可以压缩冗余、辅助显示和降低部分模型的输入维度,但主成分代表的是统计方差方向,并不天然对应明确的地物或物理过程。

变换如何形成

设每个像元由多个波段组成一个向量。首先对各波段减去均值,再计算协方差矩阵;对协方差矩阵进行特征分解后,特征向量给出新的坐标轴,特征值表示各主成分承载的方差。将中心化后的像元投影到这些轴上,就得到主成分影像。

如果使用相关系数矩阵,等价于在计算前按标准差对波段进行缩放。协方差矩阵保留原始量纲和方差差异,方差较大的波段会占更大权重;相关系数矩阵让各波段在起点上权重更接近。两种方式没有普遍优劣,应根据数据单位、噪声水平和分析目标选择,并在结果中明确说明。

数据准备决定结果

主成分轴由输入样本估计,因此云、无数据、饱和像元和异常条带会影响协方差结构。计算前应应用质量掩膜,并把填充值排除。若影像覆盖海洋、沙漠、城市和森林等差异很大的区域,最大方差可能主要来自场景构成,而不是研究对象内部的细微差异。

波段也应具有可比的物理含义。把未经缩放的整数产品、温度波段和反射率波段直接放入同一协方差矩阵,结果容易被量纲支配。多日期分析还需统一辐射处理、季节和有效像元规则。若每期影像分别拟合主成分,第一主成分在不同日期的载荷可能不同,不能直接把它们当成同一个指标比较。

大型影像通常不必用所有像元估计协方差矩阵,但抽样必须具有代表性。分层抽取不同地表类型和亮度区间,比简单取某一局部窗口更可靠。拟合完成后,应将同一组均值、缩放参数和载荷应用到待比较数据上。

如何解释主成分

每个主成分的含义应结合载荷和影像空间分布解释。若多个可见光和近红外波段载荷同号且数值接近,第一主成分往往反映总体亮度;载荷正负相反的成分可能表达波段间对比,例如植被与土壤差异。但这种解释依赖输入数据,不能仅凭“第几主成分”命名。

主成分符号本身是不确定的。一个特征向量乘以负一仍表示同一方向,因此不同软件或重新拟合后可能出现整幅影像正负颠倒。这不是计算错误。要比较不同运行结果,应比较载荷方向及其绝对关系,而不是假设亮值总表示同一种地物。

后部主成分承载的方差较小,但不等于只包含噪声。面积很小却重要的目标可能只贡献少量总体方差;相反,扫描条带、配准误差和残余云也可能形成清晰成分。保留多少成分不能只看累计解释方差,还要检查下游任务表现和重要目标是否被保留。

降维与分类

在分类或回归前使用主成分可以降低共线性和计算成本,但也会牺牲原始波段的可解释性。监督任务关注的是与标签相关的信息,而主成分只优化输入方差。一个方差很小的波段组合可能恰好对稀有类别最有区分力,因此应通过独立验证比较原始波段、选定主成分和其他特征方案。

数据泄漏也是常见问题。均值、标准差和主成分载荷只能用训练数据拟合,再应用于验证集和测试集。如果先在全部样本上拟合,验证数据的分布信息已经进入特征变换,评估会偏乐观。空间分块验证时,主成分拟合也应遵循同样的训练边界。

对高光谱数据,标准主成分可能过度关注高方差背景。最小噪声分离等方法会显式考虑噪声协方差,但需要可靠的噪声估计。无论使用哪种方法,降维都不能替代辐射校正、坏波段剔除和质量控制。

结果检查与交付

结果应记录输入波段顺序、有效像元规则、抽样方法、中心化或标准化参数、特征值、解释方差和载荷矩阵。主成分影像的存储精度也要谨慎;过早转换为低位整数可能截断负值和尾部差异。

可视化时常把前三个主成分映射到彩色通道,这有助于观察结构,却不意味着它们组成了真实颜色。分析报告应把统计变换与地学解释分开,并通过原始波段、载荷及地面知识验证所见模式。主成分分析最有价值的地方,是清楚揭示数据中主要的协同变化,而不是自动给出地物成因。

跨区域应用

在一个区域拟合的主成分轴应用到另一区域时,必须检查新数据是否超出原有均值和方差范围。若地表类型、季节或传感器发生明显变化,原轴可能不再承载相同含义。重新拟合可以提高局部方差表达,却会失去跨区域可比性;固定变换便于比较,但可能降低压缩效率。选择取决于任务目标,并应通过载荷、重建误差和下游指标共同验证。

参考资料