一种航空高光谱岩矿探测多元数据处理方法技术

本发明专利技术属于遥感技术领域，具体公开一种航空高光谱岩矿探测多元数据处理方法，该方法具体包括以下步骤：步骤(1)高光谱图像几何校正；步骤(2)大气校正，得到多元光谱数据；步骤(3)对上述步骤(2)中得到的多元光谱数据进行辐射定标与校正；步骤(4)提取对上述步骤(3)中辐射定标与校正后的多元光谱数据。本发明专利技术的方法使高光谱数据的空间位置准确，在实现数据精度的相对最优化。

A Multivariate Data Processing Method for Aerial Hyperspectral Rock and Mineral Exploration

【技术实现步骤摘要】
一种航空高光谱岩矿探测多元数据处理方法
本专利技术属于遥感
，具体涉及一种航空高光谱岩矿探测多元数据处理方法。
技术介绍
测谱学是地质领域的重要技术手段之一，能够无损、快速、大面积对地质构造、地层界线和典型矿物进行探测和识别。航空高光谱遥感是光谱学的重要数据来源，其所包含的空间、辐射和光谱三重信息，在算法的支持下，实现了岩矿相关信息定量遥感目标。传统上，各个单位在实施航空高光谱技术进行岩矿探测时，研究重点多集中在后端的高光谱信息提取工作，处理多采用传统生态环境评价的处理方法。在几何校正方面，采用遴选控制点的方法，造成条带间偏差；在辐射校正方面，缺乏仪器的标定方法；大气校正的参数常常缺失；没有充分发挥辅助数据等问题。容易忽视岩矿微弱光谱特征信息，无法实现多元数据处理的精细要求，造成所获取的航空高光谱数据常含有一些病态数据，影响到后续岩矿提取的精度。因此，亟待需要对岩矿航空高光谱数据的处理技术进行规范化，设计一套适合岩矿探测的多元数据处理技术。
技术实现思路
本专利技术的目的在于提供一种航空高光谱岩矿探测多元数据处理方法，该方法使高光谱数据的空间位置准确，在实现数据精度的相对最优化。实现本专利技术目的的技术方案：一种航空高光谱岩矿探测多元数据处理方法，该方法具体包括以下步骤：步骤(1)高光谱图像几何校正；步骤(2)大气校正，得到多元光谱数据；步骤(3)对上述步骤(2)中得到的多元光谱数据进行辐射定标与校正；步骤(4)提取对上述步骤(3)中辐射定标与校正后的多元光谱数据。所述的步骤(1)的具体步骤如下：步骤(1.1)定姿定位系统POS数据处理；步骤(1.2)对高光谱图像进行Bundle校正；步骤(1.3)对高光谱图像进行几何变换。所述的步骤(1.1)中的定姿定位系统POS数据包括机载GPS数据、BundleGCP数据、地面GPS基站数据和IMU数据管理。所述的步骤(2)的具体步骤如下：步骤(2.1)基于FLAASH模型的大气校正；步骤(2.2)获取大气校正参数；步骤(2.3)应用地面同步数据，得到多元光谱数据。所述的步骤(2.2)中采用太阳光度计CE318测得大气校正参数。所述的步骤(2.2)中的大气校正参数包括大气透过率、消光光学厚度、气溶胶光学厚度、大气水汽柱总量和臭氧总量。所述的步骤(2.3)中的多元数据信息包括日光强度、云量、气温、能见度、风向、风速、测量高度、天顶角、方位角和时间。所述的步骤(3)的具体步骤如下：步骤(3.1)光谱定标；步骤(3.2)辐射校正。本专利技术的有益技术效果在于：本专利技术的方法通过对所采集的航空高光谱数据进行一系列处理，得到高精度岩矿光谱数据，按照几何校正、大气校正、辐射定标与校正和辅助信息提取的步骤，逐渐去除航空高光谱数据中的干扰信息，使高光谱数据的空间位置准确，在辅助信息提取的基础上，实现数据精度的相对最优化。本专利技术的方法为实现航空高光谱岩矿探测多元数据处理的系统化、自动化和规范化提供了技术支持。结合岩矿探测的实际需求，设计一套涵盖几何校正、大气校正、辐射定标与校正和辅助信息提取的技术体系，为现有相关工作的提升提供了决策依据。本专利技术的方法不限定航空高光谱的波段范围，适用于可见光、近红外、中红外和热红外等多个传感器。通过一系列的数据处理技术，有利于从大噪声背景中提取岩矿微弱光谱信息，增强岩矿的特征波段，克服现有高光谱数据被云、植被和山体阴影遮挡等多个不利因素。为相关岩矿探测多元数据处理工作提供了新的技术支持。附图说明图1为本专利技术所提供的一种航空高光谱岩矿探测多元数据处理方法的流程图。具体实施方式下面结合附图和实施例对本专利技术作进一步详细说明。如图1所示，本专利技术所提供的一种航空高光谱岩矿探测多元数据处理方法，该方法具体包括以下步骤：步骤(1)高光谱图像几何校正，具体步骤如下：对航空高光谱图像进行几何校正，从而消除机载高光谱影像附上真实的地理坐标几何畸变。步骤(1.1)定姿定位系统POS数据处理定姿定位系统POS数据包括机载GPS数据、BundleGCP数据、地面GPS基站数据和IMU数据管理。按照项目分割原则，对定姿定位系统每个项目的POS数据进行存储。再利用定姿定位系统POSPacMMS5.3软件处理上述POS数据，导出相应项目的定姿定位系统POS数据，快速处理和分析校正POS数据效果。处理结束后，由于POS数据中记录了每一个GPS时间的姿态数据和位置数据，通过与基站数据差分后处理最终得到结果文件sbet.out，该文件中包含了每一个GPS时间的六方位元素(XYZ和RPH)。步骤(1.2)对高光谱图像进行Bundle校正由POS系统计算出的结果中的六方位元素是POS系统坐标下的，而不是高光谱图像坐标系统下的，因此与高光谱图像的外方位元素有一定的偏移，为了精确获取该偏移量，需要对高光谱图像进行Bundle校正。选取一块很小的飞行区域，在该区域内测量地面控制点，并在高光谱图像上找到对应的点，计算得出高光谱图像的外方位元素的六个偏移量。步骤(1.3)对高光谱图像进行几何变换高光谱原始图像中记录了每一帧图像的GPS时间，机载高光谱数据生成了一个时间标志记录文件*.att文件，该文件中记录了每一帧图像的GPS时间，该GPS时间都包含在POS系统结果文件的时间段内，通过时间比对和坐标投影变换就可以获得每一帧图像的姿态数据和位置数据，从而完成高光谱图像的几何变换。步骤(2)大气校正，得到多元光谱数据，具体步骤如下：由于电磁波在“太阳-大气-地物-大气-传感器”的传播过程中，受大气的散射、反射和吸收作用，机载高光谱传感器所接收的地物目标电磁波谱与地物标准光谱存在失真。因此，需要进行大气校正。步骤(2.1)基于FLAASH模型的大气校正基于太阳波谱范围内(不含热辐射)和平面郎伯体，在传感器处接收的像元光谱辐射亮度公式为：式中，L为传感器像元接收到的总辐射亮度；ρ为像素表面反射率；ρe为像素周围平均表面反射率；S为大气球面反照率；Lα为大气后向散射辐射率(大气程辐射)；A、B为大气条件和几何条件的两个系数。参数A、B、S和Lα的值是通过辐射传输模型MODTRAN的计算获取的，需要用到观测视场角、太阳角度、平均海拔高度，以及假设的大气模型、气溶胶类型、能见度范围。空间平均反射率用于计算大气点扩散函数，它是描述地表未处于视线路径上的部分点对目标像素的辐射贡献关系。气溶胶厚度的反演应用暗目标法。步骤(2.2)获取大气校正参数在地面同步测量中，采用太阳光度计CE318，太阳光度计CE318能自动跟踪太阳作直射辐射测量，还能进行太阳高度角、太阳主平面和极化测量。通过太阳光度计CE318测得的直射太阳辐射数据可反演出的大气校正参数数据有：大气透过率、消光光学厚度、气溶胶光学厚度、大气水汽柱总量和臭氧总量。大气校正参数数据可反演大气气溶胶粒子尺度谱分布和气溶胶相函数。步骤(2.3)应用地面同步数据，得到多元光谱数据将测区的日光强度、云量、气温、能见度、风向、风速、测量高度、天顶角、方位角、时间等多元数据信息进行处理，为大气校正提供数据源。将上述多元数据信息反演得到在不同波段上的气溶胶光学厚度、水汽含量数据、大气浑浊度和光谱卷积算法反演后的气溶胶光学厚度。步骤(3)对上述步骤(2)中得到的多元光谱数据进行辐射本文档来自技高网...

【技术保护点】
1.一种航空高光谱岩矿探测多元数据处理方法，其特征在于，该方法具体包括以下步骤：步骤(1)高光谱图像几何校正；步骤(2)大气校正，得到多元光谱数据；步骤(3)对上述步骤(2)中得到的多元光谱数据进行辐射定标与校正；步骤(4)提取对上述步骤(3)中辐射定标与校正后的多元光谱数据。

【技术特征摘要】
1.一种航空高光谱岩矿探测多元数据处理方法，其特征在于，该方法具体包括以下步骤：步骤(1)高光谱图像几何校正；步骤(2)大气校正，得到多元光谱数据；步骤(3)对上述步骤(2)中得到的多元光谱数据进行辐射定标与校正；步骤(4)提取对上述步骤(3)中辐射定标与校正后的多元光谱数据。2.根据权利要求1所述的一种航空高光谱岩矿探测多元数据处理方法，其特征在于：所述的步骤(1)的具体步骤如下：步骤(1.1)定姿定位系统POS数据处理；步骤(1.2)对高光谱图像进行Bundle校正；步骤(1.3)对高光谱图像进行几何变换。3.根据权利要求2所述的一种航空高光谱岩矿探测多元数据处理方法，其特征在于：所述的步骤(1.1)中的定姿定位系统POS数据包括机载GPS数据、BundleGCP数据、地面GPS基站数据和IMU数据管理。4.根据权利要求3所述的一种航空高光谱岩矿探测多元数据处理方法，其特征在于：所述的步骤(2)的...

【专利技术属性】
技术研发人员：张东辉，赵英俊，李璇琼，秦凯，赵宁博，张玉燕，喻翔，
申请(专利权)人：核工业北京地质研究院，
类型：发明
国别省市：北京,11