[发明专利]一种适用于花岗岩体岩性划分的综合技术方法

说明书

本发明属于多元信息融合领域，具体涉及一种适用于花岗岩体岩性划分的综合技术方法。具体包括以下步骤：步骤一、对航空放射性数据进行空间插值,形成铀钍钾彩色合成图；步骤二、ETM遥感数据的波段选取，选取中红外ETM5波段的遥感数据ETM；步骤三、数据融合，将步骤二中的遥感数据ETM和步骤一中的铀钍钾彩色合成图数据融合，得到融合图；步骤四、遥感解译。本发明通过对花岗岩体内不同岩性进行划分，从而为地质填图工作提供技术支撑。

技术领域

本发明属于多元信息融合领域，具体涉及一种适用于花岗岩体岩性划分的综合技术方法。

背景技术

航空放射性数据属于航空物探的一种，其显著的优势是成本低、效益高，并且能在地面难以开展工作的森林、沙漠、沼泽和湖泊等地区开展地质工作。在飞行测量数据采集过程中，由于多种因素的干扰，实际测线会与预期设置有所偏移。航放采样点密度大，一条测线上包含几万至几十万采样点，一个测区往往上百万采样点。由于其测线和测点的不规则性，航放数据是一种不规则的地球物探数据。航空放射性数据包含四种数据：总道(Tc)、钾道(K)、铀道(U)、钍道(Th)。这些数据中蕴含着丰富的信息，既可用于地质填图，又可用于地质找矿，但测量结果受取样岩石(或土壤)类型、土壤湿度、植被覆盖以及测量几何条件等诸多重要因素的影响较大。

1999年4月15日美国航空航天局发射的Landsat－7卫星，用于采集地球图像，进行全球变化的长期连续监测。传感器是增强型专题绘图仪ETM+(Enhanced ThematicMapper)，星上设绝对定标，提高了对地观测分辨率和定位质量，调整了辐射测量精度、范围和灵敏度，通过增益减少了强反射体造成的高亮度饱和效应。该设备包含15m分辨率的全色波段，热红外通道的空间分辨率也达到了60m，每一景覆盖面积为185km＊170km，重叠率为赤道上相邻两景图像旁向重叠率7.3％，轨道方向重叠率为5％，Band6分别具有高、低增益两种图像数据，Band1-5、7增益随季节变化可调整。

遥感图像数据融合是一个对多遥感器的图像数据和其他信息的处理过程，着重于把在空间或时间上冗余或互补的多源数据，按一定的规则或算法进行运算处理，获得比任何单一数据更精确、更丰富的信息，生成一幅具有新的空间、波谱、时间特征的合成图像。它不仅仅是数据间的简单复合，而强调信息的优化，以突出有用的专题信息，消除或抑制无关的信息，改善目标识别的图像环境，从而增加解译的可靠性，减少模糊性(即多义性、不完全性、不确定性和误差)、改善分类、扩大应用范围和效果。

遥感图像解译是根据图像的几何特征和物理性质，进行综合分析，从而揭示出物体或现象的质量和数量特征，以及它们之间的相互关系，进而研究其发生发展过程和分布规律，也就是说根据图像特征来识别它们所代表的物体或现象的性质。

发明内容

本发明要解决的技术问题是提供一种适用于花岗岩体岩性划分的综合技术方法，通过对花岗岩体内不同岩性进行划分，从而为地质填图工作提供技术支撑。

为解决上述技术问题，本发明一种适用于花岗岩体岩性划分的综合技术方法，具体包括以下步骤：

步骤一、对U、Th、K的原始数据进行均一化处理，按照0～255的数值进行线性拉伸，对航空放射性数据进行空间插值，对研究区航空放射性数据的U、Th、K进行插值，并把插值结果进行RGB彩色合成,形成铀钍钾彩色合成图；

步骤二、ETM遥感数据的波段选取，选取ETM遥感数据中的ETM5中红外波段；

步骤三、数据融合，将步骤二中的遥感数据ETM和步骤一中的铀钍钾彩色合成图数据融合，得到融合图；

步骤四、遥感解译，根据步骤三中的融合图主成分波谱融合图像中色彩的差异，进行遥感解译,划分不同岩性的边界。

所述步骤一中，利用Natural neighbor自然邻点插值法对航空放射性数据进行空间插值。