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 与 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)、统一像元大小和辐射量纲,才能进入后续蚀变信息提取环节。
总体而言,多源数据选择应坚持“区域用中分辨率、矿物识别用短波红外/高光谱、验证用高分辨率”的分工原则,下载时严把云量与质量关,预处理时统一辐射与几何基准。只有数据链条规范,后续的铁染与羟基异常提取、靶区圈定和野外验证才具备可靠基础。