一种机载摆扫式高光谱数据快速几何校正方法技术

本发明专利技术公开了一种机载摆扫式高光谱数据快速几何校正方法，涉及高光谱数据几何校正及应用技术领域。该方法利用DEM数据在无地面控制点的情况下，实现了对机载高光谱数据的系统级几何校正，同时，采用局部匹配算法有效优化了系统级几何校正的处理效率，进一步提升了航空高光谱数据几何校正的效率和精度，而且采用本发明专利技术提供的方法，可以实现对多景航空高光谱数据进行工程化批量化处理，为工程化应用提供了解决方案。

A Fast Geometric Correction Method for Airborne Sweeping Hyperspectral Data

The invention discloses an airborne swing hyperspectral data fast geometric correction method, which relates to the technical field of hyperspectral data geometric correction and application. The method realizes system-level geometric correction of Airborne Hyperspectral Data by using DEM data without ground control points. At the same time, the processing efficiency of system-level geometric correction is effectively optimized by using local matching algorithm, which further improves the efficiency and accuracy of geometric correction of Airborne Hyperspectral data. Moreover, the method provided by the invention can achieve high aeronautical accuracy for multi-scene. The engineering batch processing of spectral data provides a solution for engineering applications.

【技术实现步骤摘要】
一种机载摆扫式高光谱数据快速几何校正方法
本专利技术涉及高光谱数据几何校正及应用
，尤其涉及一种机载摆扫式高光谱数据快速几何校正方法。
技术介绍
航空遥感具有自主性强、精度高、效率高、灵活方便等优点，成为快速获取高精度遥感数据的有效手段。随着人们对于遥感数据的需求逐渐向着高空间、高时间、高光谱的“三高”方向发展，航空高光谱方面的应用变得越来越引人关注。近年来，航空高光谱技术已逐渐深入到各领域中，并且在定量化应用方面取越来越多的成果。几何校正作为遥感数据处理的重要过程，其精度对于定量化应用效果具有决定性作用。一般来说，航空遥感的工作高度在1000-3000m左右，剧烈变化的大气气流使得航空高光谱数据获取条件十分复杂，成像过程涉及参数无法被准确估计，此外，高光谱数据量大的特点也逐渐成为几何校正效率的难点，如何高效进行几何校正并获得高精度数据产品成为高光谱数据工程化应用所面临的挑战。
技术实现思路
本专利技术的目的在于提供一种机载摆扫式高光谱数据快速几何校正方法，从而解决现有技术中存在的前述问题。为了实现上述目的，本专利技术采用的技术方案如下：一种机载摆扫式高光谱数据快速几何校正方法，包括如下步骤：S1，从原始高光谱影像中提取高信噪比的单波段影像；S2，利用DEM数据和GPS数据，生成系统级几何校正控制点：S3，利用系统级几何校正控制点对所述单波段影像进行校正，得到单波段系统级几何校正影像；S4，将高分辨率影像作为参考影像，将单波段系统级几何校正影像作为待校正影像，利用待校正影像与参考影像同名点地理信息以及S2中得到的系统级几何校正控制点，采用逆变换方式生成几何精校正控制点；S5，利用几何精校正控制点对高光谱原始数据进行处理，得到几何校正后的高光谱影像。优选地，S2包括如下步骤：S201，按照如下公式计算传感器在高度H条件下的扫描向量其中，为初始坐标系下扫描向量，为机体坐标系下扫描向量，R、P、H分别代表翻滚角ω、俯仰角偏航角κ的变换矩阵，上式中各项展开后变形为：其中，s为扫描方向，θ为扫描向量与星下点方向夹角；S202，按照如下公式计算传感器扫描行中各个像元的真实扫描向量其中，h为传感器相对地面的高度值；S203，根据DEM数据和传感器GPS数据，计算得到任一DEM像元点对于传感器的扫描向量S204，按照如下公式寻找与最接近的扫描向量得到传感器扫描像元对应DEM像元的真实坐标信息：S205，按照S204的方法遍历所有的像元，得到原始影像的地理查找表，获得系统级几何校正控制点。优选地，S204采用如下原理的基于滑框的局部匹配优化的方法进行处理：在完成首次最优解求取并获得初始像元坐标之后，以初始像元为中心生成边长L＝2R+1的正方形滑框，以步长为1逐像元遍历DEM，每次求解范围将缩减到L×L范围内。优选地，S4包括如下步骤：S401，利用待校正影像与参考影像同名点地理信息，按照如下公式计算多项式校正系数：其中，u、v和x、y分别为变换前后的影像坐标，aij、bij分别为多项式校正系数，n表示多项式次数；S402，利用S2中得到的系统级几何校正控制点以及S401中得到的多项式校正系数，按照S401中的公式，逐像元计算生成几何精校正控制点。优选地，S5之后还包括步骤：重复S1至S5对多景航空高光谱数据进行工程化批量化处理。本专利技术的有益效果是：本专利技术提供的机载摆扫式高光谱数据快速几何校正方法，利用DEM数据在无地面控制点的情况下，实现了对机载高光谱数据的系统级几何校正，同时，采用局部匹配算法有效优化了系统级几何校正的处理效率，进一步提升了航空高光谱数据几何校正的效率和精度，而且采用本专利技术提供的方法，可以实现对多景航空高光谱数据进行工程化批量化处理，为工程化应用提供了解决方案。附图说明图1为本专利技术提供的机载摆扫式高光谱数据快速几何校正方法流程示意图；图2为传感器姿态参数示意图；图3为基于DEM的系统级几何校正原理示意图；图4为基于滑框的局部匹配算法原理示意图；图5是新疆罗布泊地区HyMap机载高光谱原始影像示意图；图6是对新疆罗布泊地区HyMap机载高光谱原始影像进行校正后的影像示意图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施方式仅仅用以解释本专利技术，并不用于限定本专利技术。如图1所示，本专利技术实施例提供了一种机载摆扫式高光谱数据快速几何校正方法，包括如下步骤：S1，从原始高光谱影像中提取高信噪比的单波段影像；S2，利用DEM数据和GPS数据，生成系统级几何校正控制点：S3，利用系统级几何校正控制点对所述单波段影像进行校正，得到单波段系统级几何校正影像；S4，将高分辨率影像作为参考影像，将单波段系统级几何校正影像作为待校正影像，利用待校正影像与参考影像同名点地理信息以及S2中得到的系统级几何校正控制点，采用逆变换方式生成几何精校正控制点；S5，利用几何精校正控制点对高光谱原始数据进行处理，得到几何校正后的高光谱影像。现有技术中，一般是通过分步骤的方法对高光谱数据进行几何校正，需要对所有波段进行两次校正。而由于高光谱影像具有上百个波段，所以，采用现有技术的方法对高光谱数据进行几何校正，速度会非常慢。本专利技术中，依据高光谱影像的上百个波段中，各个波段间几何信息是近似一致的特点，采用的上述方法中，首先从上百个波段中选择一个信噪比高的波段进行后续控制点提取，然后再利用最终获得的控制点对高光谱影像的所有的上百个波段进行一次性几何校正，所以，节省了很多时间，有效的提高了高光谱影像几何校正过程的效率。本专利技术的一个优选实施例中，S2可以包括如下步骤：S201，按照如下公式计算传感器在高度H条件下的扫描向量其中，为初始坐标系下扫描向量，为机体坐标系下扫描向量，R、P、H分别代表翻滚角ω、俯仰角偏航角κ的变换矩阵，上式中各项展开后变形为：其中，s为扫描方向，θ为扫描向量与星下点方向夹角；上述方法在实施过程中，具体可以为，首先假定在某一高度H时传感器各个飞行姿态角度都为0的状态为初始状态，此时扫描向量可记作对于摆扫型传感器来说，其横轴Xi的变化将取决于扫描的方向和角度s·tanθ，扫描方式决定其纵轴方向无变化，所以Yi为0，Zi为-1。如图2所示，当传感器在获取数据时，各个飞行姿态角度是变化的，此时扫描向量记作通过获取翻滚角(ω)、俯仰角及偏航角(κ)三个姿态参数可以构建机体坐标系与初始坐标系下扫描向量的变换关系。在高度h的条件下计算的扫描向量还需将其向量的模长转换为实际值，从而得到传感器扫描行中各个像元的真实扫描向量具体的，S202，可以按照如下公式计算传感器扫描行中各个像元的真实扫描向量其中，h为传感器相对地面的高度值；S203，如图3所示，根据DEM数据和传感器GPS数据，计算得到任一DEM像元点对于传感器的扫描向量S204，可以按照如下公式寻找与最接近的扫描向量得到传感器扫描像元对应DEM像元的真实坐标信息：该步骤在具体实施过程中，由于数据量较大，采用现有的匹配过程效率极其低下，因此，本专利技术中，采用了基于滑框的局部匹配算法对上述过程进行优化。如图4所示，为一个[m×n]的矩阵A，m和n分别代表DEM的长和宽，为一个[ns×nl]的矩阵B，ns和n本文档来自技高网...

【技术保护点】
1.一种机载摆扫式高光谱数据快速几何校正方法，其特征在于，包括如下步骤：S1，从原始高光谱影像中提取高信噪比的单波段影像；S2，利用DEM数据和GPS数据，生成系统级几何校正控制点：S3，利用系统级几何校正控制点对所述单波段影像进行校正，得到单波段系统级几何校正影像；S4，将高分辨率影像作为参考影像，将单波段系统级几何校正影像作为待校正影像，利用待校正影像与参考影像同名点地理信息以及S2中得到的系统级几何校正控制点，采用逆变换方式生成几何精校正控制点；S5，利用几何精校正控制点对高光谱原始数据进行处理，得到几何校正后的高光谱影像。

【技术特征摘要】
1.一种机载摆扫式高光谱数据快速几何校正方法，其特征在于，包括如下步骤：S1，从原始高光谱影像中提取高信噪比的单波段影像；S2，利用DEM数据和GPS数据，生成系统级几何校正控制点：S3，利用系统级几何校正控制点对所述单波段影像进行校正，得到单波段系统级几何校正影像；S4，将高分辨率影像作为参考影像，将单波段系统级几何校正影像作为待校正影像，利用待校正影像与参考影像同名点地理信息以及S2中得到的系统级几何校正控制点，采用逆变换方式生成几何精校正控制点；S5，利用几何精校正控制点对高光谱原始数据进行处理，得到几何校正后的高光谱影像。2.根据权利要求1所述的机载摆扫式高光谱数据快速几何校正方法，其特征在于，S2包括如下步骤：S201，按照如下公式计算传感器在高度H条件下的扫描向量其中，为初始坐标系下扫描向量，为机体坐标系下扫描向量，R、P、H分别代表翻滚角ω、俯仰角偏航角κ的变换矩阵，上式中各项展开后变形为：其中，s为扫描方向，θ为扫描向量与星下点方向夹角；S202，按照如下公式计算传感器扫描行中各个像元的真实扫描向量其中，h为传感器相对地面的高度值；S203，根据DEM数据和传感器GPS数据，计算得到任一DEM像元点对于传感器的扫...

【专利技术属性】
技术研发人员：于峻川，李逸川，闫柏琨，刘文良，马燕妮，
申请(专利权)人：中国国土资源航空物探遥感中心，
类型：发明
国别省市：北京,11