一种适用于花岗岩体岩性划分的综合技术方法技术

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

An integrated technical method for lithologic classification of granite bodies

The invention belongs to the field of multi-element information fusion, and specifically relates to an integrated technical method applicable to the lithology division of granite bodies. The following steps are as follows: Step 1, spatial interpolation of airborne radioactive data to form a uranium and thorium color synthesis map; step two, ETM remote sensing data band selection, and select the remote sensing data of the middle infrared ETM5 band ETM; step three, data fusion, step 2 remote sensing data and step one uranium and thorium color color Composite map data fusion, get the fusion diagram; step four, remote sensing interpretation. The invention provides a technical support for geological mapping work by dividing different lithologies in the granite body.

【技术实现步骤摘要】
一种适用于花岗岩体岩性划分的综合技术方法
本专利技术属于多元信息融合领域，具体涉及一种适用于花岗岩体岩性划分的综合技术方法。
技术介绍
航空放射性数据属于航空物探的一种，其显著的优势是成本低、效益高，并且能在地面难以开展工作的森林、沙漠、沼泽和湖泊等地区开展地质工作。在飞行测量数据采集过程中，由于多种因素的干扰，实际测线会与预期设置有所偏移。航放采样点密度大，一条测线上包含几万至几十万采样点，一个测区往往上百万采样点。由于其测线和测点的不规则性，航放数据是一种不规则的地球物探数据。航空放射性数据包含四种数据：总道(Tc)、钾道(K)、铀道(U)、钍道(Th)。这些数据中蕴含着丰富的信息，既可用于地质填图，又可用于地质找矿，但测量结果受取样岩石(或土壤)类型、土壤湿度、植被覆盖以及测量几何条件等诸多重要因素的影响较大。1999年4月15日美国航空航天局发射的Landsat－7卫星，用于采集地球图像，进行全球变化的长期连续监测。传感器是增强型专题绘图仪ETM+(EnhancedThematicMapper)，星上设绝对定标，提高了对地观测分辨率和定位质量，调整了辐射测量精度、范围和灵敏度，通过增益减少了强反射体造成的高亮度饱和效应。该设备包含15m分辨率的全色波段，热红外通道的空间分辨率也达到了60m，每一景覆盖面积为185km＊170km，重叠率为赤道上相邻两景图像旁向重叠率7.3％，轨道方向重叠率为5％，Band6分别具有高、低增益两种图像数据，Band1-5、7增益随季节变化可调整。遥感图像数据融合是一个对多遥感器的图像数据和其他信息的处理过程，着重于把在空间或时间上冗余或互补的多源数据，按一定的规则或算法进行运算处理，获得比任何单一数据更精确、更丰富的信息，生成一幅具有新的空间、波谱、时间特征的合成图像。它不仅仅是数据间的简单复合，而强调信息的优化，以突出有用的专题信息，消除或抑制无关的信息，改善目标识别的图像环境，从而增加解译的可靠性，减少模糊性(即多义性、不完全性、不确定性和误差)、改善分类、扩大应用范围和效果。遥感图像解译是根据图像的几何特征和物理性质，进行综合分析，从而揭示出物体或现象的质量和数量特征，以及它们之间的相互关系，进而研究其发生发展过程和分布规律，也就是说根据图像特征来识别它们所代表的物体或现象的性质。
技术实现思路
本专利技术要解决的技术问题是提供一种适用于花岗岩体岩性划分的综合技术方法，通过对花岗岩体内不同岩性进行划分，从而为地质填图工作提供技术支撑。为解决上述技术问题，本专利技术一种适用于花岗岩体岩性划分的综合技术方法，具体包括以下步骤：步骤一、对U、Th、K的原始数据进行均一化处理，按照0～255的数值进行线性拉伸，对航空放射性数据进行空间插值，对研究区航空放射性数据的U、Th、K进行插值，并把插值结果进行RGB彩色合成,形成铀钍钾彩色合成图；步骤二、ETM遥感数据的波段选取，选取ETM遥感数据中的ETM5中红外波段；步骤三、数据融合，将步骤二中的遥感数据ETM和步骤一中的铀钍钾彩色合成图数据融合，得到融合图；步骤四、遥感解译，根据步骤三中的融合图主成分波谱融合图像中色彩的差异，进行遥感解译,划分不同岩性的边界。所述步骤一中，利用Naturalneighbor自然邻点插值法对航空放射性数据进行空间插值。本专利技术的有益技术效果在于：本专利技术建立了从空间插值到遥感解译的一系列流程，提出了新的技术思路和方法，为地质填图工作提供了依据。本专利技术可用于我国花岗岩体地质填图工作。本专利技术基于明确的地质、遥感、物探方法组合，可操作性强。相关工作成果，已由野外实地勘查验证，对地质填图工作的效率提升较高。具体实施方式下面结合实施例对本专利技术作进一步详细说明。本专利技术一种适用于花岗岩体岩性划分的综合技术方法，具体包括以下步骤：步骤一、对U、Th、K的原始数据进行均一化处理，按照0～255的数值进行线性拉伸，然后，利用Naturalneighbor自然邻点插值法对航空放射性数据U、Th、K进行空间插值，并把插值结果进行RGB彩色合成为铀钍钾彩色合成图，其中，K、K/U、K/Th分别对应RGB三基色，在航空放射性数据空间插值结果中，以K、K/U、K/Th与ETM5融合后的影像效果最好，颜色存在明显差异，这种色彩差异反映了不同岩性的差异；步骤二、ETM遥感数据的波段选取，在ETM+数据各个波段中，中红外ETM5波段(1.55—1.75μm)被认为是所有波段中最佳的一个，用于分辨道路、裸露土壤、水，它还能在不同植被之间有好的对比度，并且有较好的穿透大气、云雾的能力，因此，选择中红外ETM5波段的遥感数据ETM；步骤三、数据融合方法，数据融合属于一种属性融合，它是同一地区多种遥感平台、多波段、多时相遥感数据之间以及遥感与非遥感数据之间的复合，数据融合将多源数据加以智能化合成，产生比单一信息源更精确、更完全、更可靠的估计和判断。无论是遥感信息，还是航空放射性信息，单独解决地质问题的能力都是有限的，且具有一定的偶然性和多解性，把两种信息进行融合处理，得到既具有航放的异常信息，又具有遥感丰富的地物纹理信息的综合图像，增强岩性的识别区分能力；将步骤二中的遥感数据ETM和步骤一中的铀钍钾彩色合成图数据融合，得到融合图，数据融合的方法为主成分波谱融合，主成分波谱融合是对图像依次进行主成分变换、直方图匹配和主成分逆变换的一系列过程，从几何意义来看，主成分变换后的主分量空间坐标系与变换前的多光谱空间坐标系相比旋转了一个角度，主成分变换的第一分量PCI包含信息最多，而后面的分量信息含量逐渐减少，它主要实现了在几乎不丢失信息量的同时对数据进行了压缩，同时主成分变换后的前几个主分量信噪比大，噪声相对小，因此突出了主要信息，也达到了增强图像的目的；具体为主成分波谱融合主要步骤包括三步：第一步，对步骤一中的K、K/U、K/Th彩色合成图内插后进行主成分变换；第二步，对步骤二中的ETM5数据进行直方图匹配，然后用匹配后的数据替换第一主成分波段；第三步，进行主成分逆变换。函数使用最近邻域法运算，将K、K/U、K/Th彩色合成图重采样到ETM5遥感数据像元尺寸；步骤四、遥感解译，根据步骤三中的主成分波谱融合图像中色彩的差异，进行遥感解译,划分不同岩性的边界，遥感图像解译是根据图像的几何特征和物理性质，进行综合分析，从而揭示出物体或现象的质量和数量特征，以及它们之间的相互关系，进而研究其发生发展过程和分布规律，也就是说根据图像特征来识别它们所代表的物体或现象的性质。本文档来自技高网...

【技术保护点】
一种适用于花岗岩体岩性划分的综合技术方法，其特征在于：具体包括以下步骤：步骤一、对U、Th、K的原始数据进行均一化处理，按照0～255的数值进行线性拉伸，对航空放射性数据进行空间插值，对研究区航空放射性数据的U、Th、K进行插值，并把插值结果进行RGB彩色合成,形成铀钍钾彩色合成图，然后对K、U、Th进行比值运算，得到K/U、K/Th，将K、K/U、K/Th对应RGB进行彩色合成，得到K、K/U、K/Th彩色合成图；步骤二、ETM遥感数据的波段选取，选取ETM遥感数据中的ETM5中红外波段；步骤三、数据融合，将步骤二中的遥感数据ETM和步骤一中的铀钍钾彩色合成图数据融合，得到融合图；步骤四、遥感解译，根据步骤三中的融合图主成分波谱融合图像中色彩的差异，进行遥感解译,划分不同岩性的边界。

【技术特征摘要】
1.一种适用于花岗岩体岩性划分的综合技术方法，其特征在于：具体包括以下步骤：步骤一、对U、Th、K的原始数据进行均一化处理，按照0～255的数值进行线性拉伸，对航空放射性数据进行空间插值，对研究区航空放射性数据的U、Th、K进行插值，并把插值结果进行RGB彩色合成,形成铀钍钾彩色合成图，然后对K、U、Th进行比值运算，得到K/U、K/Th，将K、K/U、K/Th对应RGB进行彩色合成，得到K、K/U、K/Th彩色合成图；步骤二、E...

【专利技术属性】
技术研发人员：孔维豪，何建国，张松，高洪雷，田明明，李平，
申请(专利权)人：核工业北京地质研究院，
类型：发明
国别省市：北京,11