应用卫星图像和数据处理增强技术进行Kerdous穹隆东南区域的地质填图

《The Egyptian Journal of Remote Sensing and Space Sciences》：Applying satellite imagery and data processing enhancements for geological mapping in the southeast area of the Kerdous inlier

　　

遥感已成为地质调查中不可或缺的工具，特别是在岩性填图和岩石单元区分方面。它对于研究难以到达的地形以及获取大范围区域的详细地质信息尤为有效。本研究采用Landsat 8 Operational Land Imager (OLI) 多光谱数据，利用一系列先进图像处理技术提取岩性和热液蚀变信息，包括波段比值 (BR)、主成分分析 (PCA)、相关系数 (CC) 和最优指数因子 (OIF)。研究人员应用了基于波段7、5、1的假彩色合成、选定的波段比值 (4/2、5/6、6/7) 以及主成分图像 (PC5、PC2、PC1) 来增强和分析研究区岩性单元的光谱特征。此外，研究人员使用热液蚀变指数来圈定富含金属元素的地带。遥感结果与地质野外观察相结合，并在地理信息系统 (GIS) 环境下进行分析，以对矿化区域进行多准则表征。研究结果表明，遥感技术显著改进并更新了该区域的现有地质图。利用Landsat 8 OLI多光谱影像，研究人员成功识别了关键的地质接触带以及具有经济价值的矿化岩脉和矿脉。ONHYM（国家碳氢化合物与矿产办公室）在填图任务中进行的野外验证和地球化学采样证实了遥感解译的准确性和可靠性。

**论文解读：应用卫星图像与数据处理增强技术进行Kerdous穹隆东南区域地质填图**

**研究背景与问题**

地质填图是矿产资源勘探和基础地质研究的基础工作。传统的野外地质填图方法在复杂地形、植被覆盖稀少但交通不便的干旱半干旱区域面临巨大挑战。遥感技术，特别是多光谱卫星影像，因其能够提供大范围、多光谱的地表信息，已成为岩性填图和热液蚀变识别的重要工具。然而，如何从海量遥感数据中有效提取地质信息，特别是区分不同岩性单元和识别与矿化相关的热液蚀变带，仍是研究的重点。在摩洛哥Anti-Atlas地区的Kerdous穹隆，该区域构造复杂，矿产资源丰富（如金、铜、黄铜矿、辉铜矿），但部分区域（如Pan-African链段）研究程度较低，传统地质图件存在更新需求。因此，开展本研究旨在利用Landsat 8 OLI多光谱数据，结合先进图像处理技术，提升该区域的地质填图精度，验证并优化此前ONHYM（国家碳氢化合物与矿产办公室）关于热液蚀变带的成果，并展示遥感在困难区域地质勘探中的辅助作用。该论文发表在《The Egyptian Journal of Remote Sensing and Space Sciences》。

**主要关键技术方法**

研究人员采用了以下关键遥感数据处理与分析方法：
1. **数据预处理**：对Landsat 8 OLI影像进行辐射定标和大气校正。
2. **最优指数因子 (OIF)**：计算波段组合的OIF值，以选取信息量最大、冗余最小的最佳假彩色合成方案（最终选定波段组合7,5,1）。
3. **主成分分析 (PCA)**：对波段1-7进行PCA，生成独立的主成分图像（PC1-PC7），利用特征向量统计量获取各岩石单元的光谱信息，并选择PC5、PC2、PC1等组合进行彩色合成，以增强岩性区分。
4. **波段比值 (BR)**：采用特定波段比值（如4/2、5/6、6/7）分别映射铁氧化物、铁镁矿物和粘土矿物（如高岭石、伊利石），以识别热液蚀变矿物。
5. **相关系数 (CC)**：分析波段间相关性，辅助理解光谱特征。
6. **野外验证**：将遥感解译结果与ONHYM在PNCG计划填图任务中进行的野外地质观测、地球化学采样及现有地质图进行对比验证。

**研究结果**

**4.1 岩性填图**
通过OIF计算，研究人员确定波段组合7-5-1能最佳表征研究区不同地质体。该组合（图3a）清晰区分了不同花岗岩体（如Afra花岗岩呈褐色蘑菇状、Amanouz花岗岩呈米色树叶状）、变质沉积岩（如Had-N-Tahala组呈深蓝色）以及基性和酸性岩脉（呈橙色至红色）。组合7-4-1（图3b）进一步补充了细节，区分了Tahala花岗岩（褐紫色）、Anmmar和Tasserhit花岗岩（黄色）以及流纹质熔结凝灰岩（深紫色）。组合7-2-1（图3c）以红色调反映大多数花岗岩，深绿色调反映变质沉积岩。主成分分析（PCA）的结果（图4、图5）在灰度以及彩色合成图像（PC5-2-1和PC5-4-2）中进一步突出了不同岩性单元，特别是Afra花岗岩与Phad变质沉积岩组之间的接触带，该接触带在组合PC5-2-1中呈粉红色，在PC5-4-2中呈黄橙色。通过多种波段比值组合（图6），研究人员实现了更清晰的岩性区分：例如组合6/4、5/3、4/2（图6a）检测到Afra花岗岩与Phad组之间的接触带（紫色）；组合4/2、5/6、6/7（图6b）以蓝色区分变质沉积岩、以绿色区分花岗岩；组合6/5、7/6、4/7（图6c）以黄色突出显示了不同岩脉和矿脉。

**4.1.1 热液蚀变填图**
研究人员利用波段比值4/2、5/6、6/7分别映射铁氧化物、铁镁矿物和粘土矿物（图7、图8）。波段比值6/7（图8c）以粉色识别出粘土矿物（如高岭石）富集区，与流纹质熔结凝灰岩及Awgmoute花岗岩中的裂隙吻合，其中含硫酸盐矿化，且由PIMA研究确定存在低温高岭石蚀变。波段比值4/2（图8a）以红色映射铁氧化物，与洪积锥、火山岩及Anmmar、Tasserhit、Tar?wat花岗岩分布区吻合，后者中可见含赤铁矿的裂隙。波段比值5/6（图8b）以蓝色映射铁镁矿物，主要对应于Had-N-Tahala组的变质沉积岩、Anzi砂岩（位于Afra花岗岩南侧）、流纹质熔结凝灰岩及下Adoudounien砾岩。这些结果与ONHYM的野外观察和分析一致。最终的综合热液蚀变图（图9）清晰地展示了蚀变区域的空间分布。

**讨论与结论**

研究表明，尽管遥感存在空间分辨率不足以识别小尺度特征、仅提供地表信息、受气象条件影响等局限性，但结合Landsat 8 OLI多光谱影像与OIF、PCA、波段比值等图像增强技术，能够显著提升对Kerdous穹隆东南区域岩性单元和热液蚀变带的识别能力。研究结果成功验证并更新了ONHYM此前的地质调查成果，精确界定了Afra花岗岩与变质沉积岩之间的主要变质接触带，并圈定了具有经济潜力的热液蚀变区域。研究结论强调，Landsat 8 OLI遥感影像在划分该区域岩相方面具有重要贡献，所采用的改进方法（如OIF、密度分割、MNF、PCA、彩色组合等）实现了对几乎所有露头的精确分析。遥感技术与GIS及野外验证的结合，被证明是复杂地形区域地质勘探和矿产资源评价的有效途径。