一种艇载高光谱推帚成像数据的几何校正方法
【技术领域】
[0001] 本发明涉及图像处理技术领域，涉及高光谱图像处理技术领域特别，为一种艇载 高光谱推帚成像数据的几何校正方法。
【背景技术】
[0002] 定位定姿系统（position and orientation system, POS)利用全球定位系统和惯 性测量装置直接确定传感器空间位置和姿态的集成技术。
[0003] 20世纪80年代，随着遥感技术的飞速发展，新兴起的高光谱成像技术成为空间对 地表观测的有力工具。高光谱成像数据在为人们提供了有关地表地物的丰富空间信息和光 谱信息的同时，也将卫星遥感数据高精度处理与深度处理需求提升到了一个更高层次。高 光谱成像技术能够获取大量窄波段的准连续光谱数据，每一个像元具有一条几乎连续的光 谱曲线。与传统的多光谱遥感技术相比，高光谱遥感能够提供更加丰富的地物观测信息，可 以用来生成更为复杂的观测模型，提高地表地物状况的综合判别与分析能力。
[0004] 然而，星载遥感平台受轨道的限制，经过带观测区域时间固定，无法实现应急观 测；机载遥感审批程序非常复杂，受天气条件影响较大，且获取影像数据成本较高。飞艇对 地观测平台可提供比星载或机载平台更为丰富的多源同步对地观测能力；可提供星载平台 难以实现的动态对地观测能力和精细对地观测能力；可提供机载平台难以实现的定点悬停 对地观测能力；具有安全、可回收、易维护、可巡航、可更换载荷等有点，将成为未来的主要 遥感平台之一。
[0005] 但飞艇受气流侧风影响较大，自身稳定性较差。而高光谱仪一般采用推帚成像模 式，成像过程中每一扫描线的飞行姿态都在变化，使得高光谱数据产生很大的几何形变。而 这种几何形变与常规图像的几何形变恢复有较大区别，如何对这种图像进行恢复是飞艇遥 感平台需要解决的关键问题。
【发明内容】
[0006] 为了克服上述现有技术的缺点，本发明的目的在于提供一种艇载高光谱推帚成像 数据的几何校正方法，解决了现有算法处理艇载推帚成像数据几何校正效果差的问题。
[0007] 为了实现上述目的，本发明采用的技术方案是：
[0008] 一种艇载高光谱推帚成像数据的几何校正方法，包括如下步骤：
[0009] 步骤1，对于一幅待处理图像数据，读取其相对应的艇载POS数据，在飞艇飞行路 线的经炜度范围内，选择其中一点作为坐标原点，将艇载POS数据的经炜度坐标转换为地 面绝对坐标，将艇载POS数据的海拔高度坐标转换为距地面高度坐标；
[0010] 步骤2,利用艇载POS数据计算高光谱相机扫描线中点的空间坐标；优选地，其实 现方法如下：
[0011] 分别计算由飞艇俯仰、翻滚原因引起的高光谱相机扫描线中点偏移，并根据飞艇 航向角计算高光谱相机扫描线中点的空间坐标。
[0012] 步骤3,利用艇载POS数据计算失真图像全部像素的空间坐标；优选地，其实现方 法如下：
[0013] 首先逐扫描线计算每条扫描线上像素的空间坐标；
[0014] 然后计算每条扫描线上像素相对于扫描线中点的距离从而得到失真图像全部像 素的空间坐标。
[0015] 其中在计算每条扫描线上像素相对于扫描线中点的距离时，考虑由于飞艇不是正 视成像引起的扫描线仿射变换。
[0016] 步骤4,建立失真图像到校正图像的1像素最近邻对应关系；优选地，其实现方法 如下：
[0017] 对失真图像每个像素寻找校正图像中的4个最近邻像素，并建立4对对应关系；
[0018] 对于校正图像的每个像素，在包含该像素的对应关系中寻找距离最小的一对，并 舍弃其他包含该像素对应关系。
[0019] 步骤5,采用最近邻插值与8邻域插值相结合的方式来校正图像，优选地，其实现 方法如下：
[0020] 对校正图像中存在1像素最近邻对应关系的像素进行最近邻插值；
[0021] 对校正图像中其他像素迭代进行8邻域插值，即取像素8邻域内存在像素值的点 的平均值。
[0022] 与现有技术相比，本发明通过利用艇载POS数据计算高光谱相机扫描中心点空间 坐标以及计算失真图像全部像素的空间坐标，建立了从失真图像像素到校正图像像素的对 应关系。利用最近邻插值与8邻域插值相结合的方式来校正图像，快速且不失精度地对图 像的几何形变进行了恢复。
【附图说明】
[0023] 图1为本发明实施例1公开的艇载高光谱推帚成像数据的几何校正方法流程图。
[0024] 图2为艇载POS数据的一行示例。
[0025] 图3为本发明实施例2公开的艇载高光谱推帚成像数据的几何校正方法流程图。
[0026] 图4为飞艇POS数据记录的三个姿态角，翻滚角、俯仰角、偏航角的示意图。
[0027] 图5为相机扫描线受翻滚偏移而产生的仿射变换示意图。
[0028] 图6为建立失真图像到校正图像的1像素最近邻对应关系不意图。
[0029] 图7为失真图像扫描线分布情况的示意图。
【具体实施方式】
[0030] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0031] 本发明公开了一种艇载高光谱推帚成像数据的几何校正方法，以解决现有技术的 算法处理艇载推帚成像数据几何校正效果差的问题。其【具体实施方式】如下所述：
[0032] 实施例一
[0033] 本实施例公开的艇载高光谱推帚成像数据的几何校正方法的流程如图1所示，
[0034] 步骤S11、利用艇载POS数据计算高光谱相机扫描线中点的空间坐标；
[0035] 设（xu yj为通过艇载POS数据记录的经炜度转换的飞艇空间坐标，height为由 飞艇POS数据中的海拔度计算而得的距地面高度，（roll, pitch, yaw)为飞艇POS数据记录 的三个姿态角，翻滚角、俯仰角、偏航角，如图4所示。利用公式1. 1，计算得出高光谱相机扫 描线中点的空间坐标（xM，yM)，
[0036] Offsetroll= height*tan (-roll)
[0037] Offsetpitch= height/cos (roll) *tan (pitch) (0.1 )
[0038] xM= X L+offsetr〇n*cos (yaw) +offsetpitch*sin (yaw)
[0039] yM= y L-〇ffsetroll*sin (yaw) +offsetpitch*cos (yaw)
[0040] 式中，offsetMll是由翻滚角引起的扫描线中点偏移，称之为翻滚偏移；offSet pitch 是由俯仰角引起的扫描线中点偏移，称之为俯仰偏移；写成复数形式即，
[0041] xM+i*yM= X L+i^yL+ofTseUe't+i^ofTsetpitd^e'*1 (0. 2)
[0042] 可以看出，扫描线中点的空间坐标（xM，yM)即由飞艇坐标（?以，加上翻滚偏移和 俯仰偏移得到。
[0043] 步骤S12、利用艇载POS数据计算失真图像全部像素的空间坐标；
[0044] 本发明逐扫描线的计算失真图像全部像素的空间坐标，将高光谱相机扫描线 中点空间坐标（xM，yM)代入如下公式，计算出失真图像一条扫描线上像素的空间坐标 (Xline，Yline) '
【主权项】
1. 一种艇载高光谱推帚成像数据的几何校正方法，其特征在于，包括如下步骤： 步骤1，对于一幅待处理图像数据，读取其相对应的艇载POS数据，在飞艇飞行路线的 经炜度范围内，选择其中一点作为坐标原点，将艇载POS数据的经炜度坐标转换为地面绝 对坐标，将艇载POS数据的海拔高度坐标转换为距地面高度坐标； 步骤2,利用艇载POS数据计算高光谱相机扫描线中点的空间坐标； 步骤3,利用艇载POS数据计算失真图像全部像素的空间坐标； 步骤4,建立失真图像到校正图像的1像素最近邻对应关系； 步骤5,采用最近邻插值与8邻域插值相结合的方式来校正图像。
2. 根据权利要求1所述艇载高光谱推帚成像数据的几何校正方法，其特征在于，所述 步骤2中，分别计算由飞艇俯仰、翻滚原因引起的高光谱相机扫描线中点偏移，并根据飞艇 航向角计算高光谱相机扫描线中点的空间坐标。
3. 根据权利要求1所述艇载高光谱推帚成像数据的几何校正方法，其特征在于，所述 步骤2中，高光谱相机扫描线中点的空间坐标（xM，yM)的计算公式如下： offsetr〇n= height*tan(-roll)offsetpitch= height/cos(roll)*tan (pitch) xM=xL+offsetr〇n*cos(yaw)+offsetpitch*sin(yaw) yM=yL_〇ffsetroll*sin(yaw)+offsetpitch*cos(yaw) 其中，offsetMll是由翻滚角引起的扫描线中点偏移，称之为翻滚偏移；offset^^是 由俯仰角引起的扫描线中点偏移，称之为俯仰偏移，（xuyj为通过艇载POS数据记录的经 炜度转换的飞艇空间坐标，height为由飞艇P0S数据中的海拔度计算而得的距地面高度， (roll,pitch,yaw)为飞艇P0S数据记录的三个姿态角，即翻滚角、俯仰角、偏航角。
4. 根据权利要求1所述的几何校正方法，其特征在于，所述步骤3中： 首先逐扫描线计算每条扫描线上像素的空间坐标； 然后计算每条扫描线上像素相对于扫描线中点的距离从而得到失真图像全部像素的 空间坐标，在计算每条扫描线上像素相对于扫描线中点的距离时，考虑由于飞艇不是正视 成像引起的扫描线仿射变换。
5. 根据权利要求3所述的几何校正方法，其特征在于，所述步骤3中：失真图像一条扫 描线上像素的空间坐标（xliM，yliM)的计算公式如下： h=hpip-ht/r.o.9Trol1 ^) /Tm*t.r.h')
xiine=xM+line*cos(yaw)yime=yM-line*sin(yaw) 其中，FOV为高光谱相机的视角，samples为高光谱相机的采样率，h的物理意义为从高 光谱相机到扫描线中点的距离，a为扫描线上的点相对于相机中心线的角度，line为经过 了仿射变换后，一条扫描线上的点相对于扫描线中心点的距离。
6. 根据权利要求1所述的几何校正方法，其特征在于，所述步骤4中： 对失真图像每个像素寻找校正图像中的4个最近邻像素，并建立4对对应关系； 对于校正图像的每个像素，在包含该像素的对应关系中寻找距离最小的一对，并舍弃 其他包含该像素对应关系。
7. 根据权利要求1所述的几何校正方法，其特征在于，所述步骤5中： 对校正图像中存在1像素最近邻对应关系的像素进行最近邻插值； 对校正图像中其他像素迭代进行8邻域插值，即取像素8邻域内存在像素值的点的平 均值。
【专利摘要】一种艇载高光谱推帚成像数据的几何校正方法，利用飞艇POS数据计算高光谱相机扫描线中点的空间坐标；利用飞艇POS数据计算失真图像全部像素的空间坐标；利用失真图像像素的空间坐标，建立失真图像到校正图像的1像素最近邻对应关系；采用最近邻插值与8邻域插值相结合的方式来校正图像；本发明采用计算失真图像像素空间坐标的方式，建立失真图像与校正图像像素间的对应关系，并采用最近邻插值与8邻域插值相结合的方式快速地校正图像。
【IPC分类】G06T5-00
【公开号】CN104700367
【申请号】CN201510097304
【发明人】庞庆宇, 禹晶, 孙卫东
【申请人】清华大学
【公开日】2015年6月10日
【申请日】2015年3月5日