1. 课程导入与任务总览
在矿产勘查与地质填图的早期阶段,遥感蚀变信息提取是一种成本低、覆盖广、见效快的技术手段。其核心逻辑并不复杂:热液活动在改造围岩的过程中,会生成一批特定的蚀变矿物,而这些矿物在可见光—短波红外波段具有可识别的诊断性吸收特征。铁染蚀变以褐铁矿、赤铁矿、针铁矿等含铁矿物为代表,通常在蓝波段(约0.49 μm)附近出现强吸收,在红波段相对反射较高,整体呈现"蓝低红高"的光谱形态;羟基蚀变以高岭石、绢云母、明矾石、绿泥石等含羟基或碳酸根矿物为代表,其诊断性吸收位于短波红外约2.20 μm附近,并常在2.10 μm与2.30 μm两侧形成相对反射肩。把这两类光谱特征从影像中稳定地"翻译"成空间分布图,就是本课程要完成的核心任务。
从地质意义上看,铁染与羟基异常并不是矿本身,而是热液系统留下的"化学脚印"。铁染异常往往对应氧化带、黄铁矿风化帽或酸性淋滤带,可指示浅成低温热液金矿、斑岩铜矿的氧化部位;羟基异常则更多反映中低温热液蚀变带、泥化带与绢英岩化带,是寻找斑岩型、浅成低温热液型及部分VMS型矿床的重要线索。二者在空间上的套合关系尤其有价值:当羟基异常呈环带状分布、铁染异常位于其外侧或叠加于核部时,往往暗示一个较完整的蚀变分带系统。因此,本课程不满足于"提取出一张异常图",而是要求学员理解异常的空间结构,并能结合地质背景给出合理解释。
选择Sentinel-2作为数据源,是基于其波段配置与工程可用性的综合权衡。Sentinel-2 MSI在可见光—近红外—短波红外共设置13个波段,其中B2(0.49 μm)、B4(0.665 μm)、B6(0.74 μm)、B11(1.61 μm)、B12(2.19 μm)对铁染与羟基识别最为关键;B11与B12的20 m空间分辨率虽低于10 m可见光波段,但对区域尺度蚀变填图已足够。更实际的优势在于:数据免费开放、重访周期约5天(双星)、全球覆盖、且已进行过辐射定标与正射校正,配合大气校正产品(如Sen2Cor输出的L2A级地表反射率),可以大幅压缩预处理工作量。这意味着学员可以把主要精力放在光谱机理理解与异常提取方法上,而不是消耗在数据获取与几何校正上。
本课程的实战目标可归纳为三个层次。第一层是"会处理":能够独立完成L2A数据的波段选择、重采样、掩膜(云、水体、植被、阴影)与反射率归一化,得到可用于计算的地表反射率立方体。第二层是"会提取":掌握波段比值法、主成分分析(PCA)与光谱指数法三类主流思路,并能针对铁染与羟基分别构建合理的指数。例如铁染常用B4/B2或B4/B3比值增强"红高蓝低"特征;羟基常用B11/B12比值突出2.20 μm吸收造成的B12相对低值。第三层是"会解释":能对提取结果进行阈值分割、密度统计与地质叠合分析,判断异常是否具有地质合理性,而不是把噪声或云边缘误判为蚀变。
为了让目标可检验,课程设定了明确的成果交付要求。最终需提交一套完整的成果包,包括:①研究区L2A反射率数据与预处理记录(含掩膜规则与参数);②铁染与羟基异常分级图(建议分为三级:强、中、弱),并附阈值选取依据;③异常面积统计表与空间分布描述;④一份简短的结果解释说明,指出主要异常带的位置、走向及其与已知地质体或矿化点的关系。所有中间参数(如PCA特征向量、比值阈值、掩膜阈值)都应记录在案,保证结果可复现。
需要特别提醒的是,蚀变提取结果对阈值极为敏感。同一幅影像,阈值放宽一点,异常面积可能翻倍;收紧一点,又可能把真实弱异常滤掉。因此课程要求每位学员在提交成果时,必须说明阈值是如何确定的——是基于直方图拐点、均值加若干倍标准差,还是参照已知矿点的异常强度反推。这种"参数留痕"的习惯,比最终图件本身更能体现专业素养。