3. 多源遥感数据选择与获取

铁矿找矿遥感工作能否顺利推进,很大程度上取决于数据选得对不对、拿得全不全、预处理做得规不规范。不同卫星在空间分辨率、光谱配置、幅宽、重访周期和数据获取成本上差异明显,没有任何一种数据能同时满足“看得广、看得细、看得准、拿得到”的全部要求。实际工作中通常采用“多源协同”的策略:用中分辨率多光谱数据做区域蚀变异常提取,用高分辨率数据做矿区尺度构造解译与验证,用热红外或短波红外数据补充特定蚀变矿物识别。本节围绕铁矿找矿需求,梳理主流数据源的适用性差异,并给出可操作的下载与预处理流程。

一、各数据源在铁矿找矿中的适用性对比

Landsat 8/9 OLI 是区域尺度铁染与羟基蚀变提取的“基本盘”。它拥有 30 m 多光谱、15 m 全色、100 m 热红外,幅宽约 185 km,重访 16 天,且全球存档免费、时间序列长。OLI 的 Band 2(蓝)、Band 5(近红外)、Band 6/7(短波红外)组合可有效支撑铁染(Band 2/ Band 5)与羟基(Band 6/Band 7)比值或主成分分析。其不足是短波红外仅两个宽波段,难以细分具体蚀变矿物种类。

ASTER 的价值集中在短波红外与热红外。其 VNIR 15 m、SWIR 30 m、TIR 90 m,SWIR 有 6 个波段(Band 4–9),可区分高岭石、绢云母、绿泥石、碳酸盐等矿物组合,对矽卡岩型、热液型铁矿的围岩蚀变分带识别很有帮助;TIR 5 个波段可用于硅酸盐指数与石英含量估算。ASTER 幅宽 60 km,重访约 16 天,但 2008 年后 SWIR 探测器故障导致部分数据质量下降,选数据时需关注获取时间与质量标记。

Sentinel-2 MSI 的优势是“高重访 + 红边波段 + 免费”。其 10 m 可见光/近红外、20 m 红边与短波红外,幅宽 290 km,双星重访约 5 天,云量筛选后可用影像多,适合做时间序列变化检测和植被覆盖区的蚀变信息增强。红边波段(Band 5–7)对铁氧化物与植被胁迫响应敏感,可在植被覆盖度较高的地区辅助识别与矿化相关的异常。

WorldView-3 属于高分辨率商业数据,全色 0.31 m、多光谱 1.24 m、短波红外 3.7 m,且拥有 8 个 SWIR 波段,可做矿物填图级别的精细识别,适合矿区尺度构造解译、采坑与尾矿识别、蚀变带边界圈定。缺点是幅宽仅约 13 km、单景成本高、存档覆盖有限,通常只用于重点靶区验证。

国产卫星方面,高分一号/二号(GF-1/GF-2)提供 2 m/0.8 m 全色与 8 m/3.2 m 多光谱,适合矿区高分辨率底图;资源一号02D(ZY1-02D)搭载 166 谱段高光谱成像仪,5 m/10 m 分辨率,对蚀变矿物精细识别有潜力;高分五号(GF-5)AHSI 高光谱 30 m、330 个波段,可用于矿物填图研究。国产数据在获取时效、区域覆盖和成本上具有优势,但高光谱数据的辐射定标与大气校正处理链条相对复杂,需要更谨慎的质量控制。

多源遥感数据适用性对比与铁矿找矿协同分工 区域尺度 · 蚀变异常提取 Landsat 8/9 OLI 30 m 多光谱 / 185 km 幅宽 铁染:B2/B5;羟基:B6/B7 Sentinel-2 MSI 10/20 m · 290 km · 5 天重访 红边波段增强植被区异常 用途:圈定找矿远景区 矿物尺度 · 蚀变分带识别 ASTER SWIR 6 波段:高岭石/绢云母 TIR 5 波段:硅酸盐指数 ZY1-02D / GF-5 AHSI 高光谱 166/330 波段 矿物填图与端元识别 用途:判定蚀变类型与分带 矿区尺度 · 精细解译验证 WorldView-3 全色 0.31 m / SWIR 3.7 m 8 个 SWIR 波段精细填图 GF-1 / GF-2 0.8–2 m 全色底图 构造、采坑、尾矿识别 用途:靶区圈定与野外核查 缩小靶区 重点验证 下载与预处理主流程(以 Landsat 8 OLI 为例) ① 数据检索 USGS / 地理空间数据云 ② 云量筛选 云量 < 10% ③ 辐射定标 DN → 辐亮度 ④ 大气校正 FLAASH / 6S ⑤ 几何校正 正射 / 配准 ⑥ 裁剪 研究区掩膜 预处理输出:地表反射率影像 + 质量波段 + 研究区裁剪结果 后续用于铁染/羟基比值、主成分分析、蚀变异常分级与靶区圈定 注:Sentinel-2 用 Sen2Cor 做大气校正;ASTER 需关注 SWIR 质量标记

二、下载渠道与检索要点

Landsat 与 ASTER 数据可通过 USGS EarthExplorer 检索下载,也可使用地理空间数据云、国家对地观测科学数据中心等国内镜像。Sentinel-2 数据在 Copernicus Open Access Hub 或欧空局相关平台获取,国内可通过国家对地观测数据中心镜像下载。WorldView-3 需通过商业代理商按景订购,重点靶区可先申请历史存档样例评估质量。国产高分数据通过中国资源卫星应用中心或省级卫星应用中心申请,部分数据对科研用户免费。

检索时建议按“研究区范围 + 时间窗口 + 云量阈值 + 数据级别”四个条件组合。以铁矿区为例,时间窗口优先选择植被枯黄期或裸露基岩出露期(如我国北方 3–5 月、10–11 月),云量阈值一般设为 10% 以下;若研究区多云,可放宽到 20% 并配合多时相合成。数据级别上,Landsat 选 Collection 2 Level-1(已做几何与辐射校正)或 Level-2(地表反射率产品),Sentinel-2 选 L1C 或 L2A,ASTER 选 AST_L1T。下载前务必查看浏览图与质量波段,避免选到条带、云影或传感器异常的数据。

三、预处理流程与关键参数

以 Landsat 8 OLI 为例,标准预处理链为:辐射定标 → 大气校正 → 几何精校正/配准 → 研究区裁剪。辐射定标将 DN 值转为辐亮度或表观反射率,定标系数从元数据 MTL 文件中读取;大气校正推荐使用 FLAASH 或 6S 模型,输入需包含传感器类型、成像时间、太阳高度角、气溶胶模型等参数,输出为地表反射率。几何校正以已有正射影像或 DEM 为参考,控制点残差控制在 0.5 像元以内。裁剪时用研究区矢量掩膜,并保留质量波段用于后续云掩膜。

Sentinel-2 的 L1C 数据需用 Sen2Cor 做大气校正生成 L2A,再重采样到 10 m 或 20 m 统一分辨率;ASTER SWIR 数据需检查 2008 年后的质量标记,必要时用同期其他数据替代。所有数据预处理完成后,应统一投影(如 UTM/WGS84)、统一像元大小和辐射量纲,才能进入后续蚀变信息提取环节。

实操建议:首次处理某区域时,先选一景云量低、质量好的 Landsat 8 影像完整走通“定标—大气校正—裁剪”流程,检查地表反射率曲线是否合理(如水体在近红外应接近 0,植被在红边应抬升),确认无误后再批量处理其他数据源,可显著减少返工。
常见问题:直接对 DN 值做比值运算会放大地形阴影和传感器差异,导致铁染异常假阳性;未做大气校正的影像在短波红外波段受水汽影响明显,羟基异常提取结果不稳定。务必在反射率数据基础上开展蚀变信息提取。

总体而言,多源数据选择应坚持“区域用中分辨率、矿物识别用短波红外/高光谱、验证用高分辨率”的分工原则,下载时严把云量与质量关,预处理时统一辐射与几何基准。只有数据链条规范,后续的铁染与羟基异常提取、靶区圈定和野外验证才具备可靠基础。