1. 课程导入与任务总览
地质遥感波段比值异常图,是利用多光谱或高光谱影像中不同波段反射率的比值运算,压制地形阴影、反照率差异和大气影响,突出与矿化蚀变、岩性差异相关的微弱光谱特征,最终形成的一张可供地质解译与找矿靶区圈定的专题图件。本课程要完成的核心任务,是把一景原始遥感影像,经过辐射与几何预处理、波段选择与比值计算、异常增强与阈值分割、噪声抑制与地质筛选,输出为一张带有异常等级分级、坐标信息与图例说明的波段比值异常图产品。
需要先建立三个基本认识。第一,比值图不是“算完就完”的中间结果,它必须是一张可交付的产品:有明确的投影与坐标、有可复现的处理参数、有分级标准与图例、有异常区编号与属性表。第二,比值异常的本质是“相对高值”,它反映的是像元之间的光谱差异,而不是绝对含量,因此异常图必须与地质背景、已知矿点、构造线等资料叠加判读,才能成为找矿信息。第三,整条工作流是有顺序依赖的:预处理不到位,比值结果会被阴影和大气差异污染;波段选择不合理,异常会被淹没在噪声里;阈值定得过松或过紧,异常图要么满图红点、要么漏掉目标。
成果形式与交付标准
本课程的最终成果是一张波段比值异常图,建议以 GeoTIFF 栅格加 PDF/PNG 制图版两种形式同时提交。栅格保留分级编码值(例如 1、2、3 分别代表一、二、三级异常,0 为背景),便于后续在 GIS 中统计面积、叠加分析;制图版则用于汇报与野外对照。异常分级建议采用三级制:以比值影像的均值 μ 与标准差 σ 为基准,μ+1σ 至 μ+2σ 为一级异常,μ+2σ 至 μ+3σ 为二级异常,大于 μ+3σ 为三级异常。这一分级方式在示例数据中便于操作,但实际项目应根据影像质量与已知矿点命中情况调整倍数,不能机械套用。
异常区属性表至少应包含:异常编号、中心坐标、面积(km²)、比值均值、最大值、所处地层单元、距最近断裂距离、是否与已知矿点重合。处理记录则应写清影像标识、波段编号与中心波长、比值公式、大气校正方法、滤波窗口大小、阈值参数与软件版本。这两项是保证成果可复现、可审查的关键,也是很多初学者容易遗漏的部分。
完整工作流的顺序依赖
整条工作流可以概括为五段:数据输入与检查、预处理、波段选择与比值计算、异常增强与分级、成果制图与判读。它们之间存在明确的顺序依赖:大气校正不充分时,比值影像会出现大范围渐变,阈值分割会把整片区域判为异常;几何校正不到位时,异常区与地质图叠加会错位,属性表中的地层单元统计就失去意义;波段选择不当,例如把两个相关性极高的相邻波段做比值,结果接近常数,异常信息几乎为零。
因此在实际操作中,建议在比值计算之后先做一次快速质检:查看比值影像的直方图是否近似单峰、是否存在明显的条带或云影残留、异常高值是否集中在特定地物上。若直方图出现双峰且低值峰对应水体或阴影,说明掩膜不完整;若高值呈规则条带状,多为传感器条带或地形阴影残留,需要回退到预处理阶段处理。这一质检环节不产生新成果,但决定了后续分级是否有意义。
本课程的组织方式
后续内容将按上述五段工作流依次展开:先讲清多光谱波段与蚀变矿物光谱特征的对应关系,说明为什么某些波段比值能突出蚀变;再进入预处理与比值计算的具体操作,给出可复现的公式与参数;然后讲异常增强、阈值分割与噪声抑制;最后讲制图输出、属性统计与地质判读。每一节都会围绕一个可操作的具体环节展开,并在关键处给出示例参数与判断依据。学习时建议对照本节的工作流图,明确当前所处阶段、输入是什么、输出是什么、以及输出不合格时应回退到哪一步。
需要强调的是,波段比值异常图是一种“找矿线索图”,它的价值在于缩小野外验证范围,而不是直接给出矿体位置。判读时必须结合区域成矿规律、地层岩性、构造展布和已知矿化信息综合解释。把异常图当作唯一依据去布设工程,是这类方法最常见的误用。本课程在成果输出部分会专门讨论异常筛选与靶区圈定的判读逻辑,帮助大家把图件转化为可执行的验证建议。