卫星影像的正射校正方法和系统转让专利
申请号 : CN201510703289.8
文献号 : CN105387846B
文献日 : 2017-11-10
发明人 : 朱德海 , 叶思菁 , 岳彦利 , 刘帝佑 , 杜振博 , 李林 , 熊全 , 姚晓闯 , 李鹏山
申请人 : 中国农业大学
摘要 :
权利要求 :
1.一种卫星影像的正射校正方法,其特征在于,包括:
S1:在地理坐标系下划定正射校正变换范围;获取所述变换范围内的数字高程模型数据,将所述数字高程模型数据转换到所述地理坐标系下;
S2:获取待校正卫星影像的有理函数系数,根据所述有理函数系数建立有理函数模型,通过所述有理函数模型正变换估算所述待校正卫星影像的顶点像元物方坐标;
S3:根据所述待校正卫星影像的顶点像元物方坐标,提取与所述待校正卫星影像对应的所述转换到地理坐标系下的数字高程模型数据;
S4:将与所述待校正卫星影像对应的所述数字高程模型数据通过有理函数模型正变换,重新计算所述待校正卫星影像的顶点像元物方坐标,从中提取物方坐标的最大值与最小值;
S5:根据投影变换关系,计算所述物方坐标的最大值与最小值对应的目标平面坐标;
S6:根据所述物方坐标最大值与最小值对应的目标平面坐标与所述待校正卫星影像的像元空间分辨率,建立所述目标平面坐标系与所述待校正卫星影像的像方坐标系之间的对应关系。
2.根据权利要求1所述的卫星影像的正射校正方法,其特征在于,所述步骤S1还包括:在所述变换范围内建立地理格网,基于所述地理格网对所述数字高程模型数据执行重采样、裁剪、镶嵌,将所述数字高程模型数据分成若干分块栅格文件,使所述地理格网的每个网格对应唯一的一个分块栅格文件。
3.根据权利要求1或2所述的卫星影像的正射校正方法,其特征在于,所述步骤S3包括:将所述顶点像元物方坐标的最大值增加0.01°,将所述顶点像元物方坐标的最小值减小0.01°,作为提取所述数字高程模型数据的边界范围,并根据所述边界范围建立高程二维数组,建立所述高程二维数组的数组元素与所述数字高程模型数据的对应关系,将所述数字高程模型数据赋值到数组。
4.根据权利要求3所述的卫星影像的正射校正方法,其特征在于,所述高程二维数组为[xls,yls]:xls=[((lonmax+0.01)-(lonmin-0.01))/0.00025]+1yls=[((latmax+0.01)-(latmin-0.01))/0.00025]+1其中,“[]”为向下取整符号,(lonmax,latmax)为顶点像元物方坐标的最大值,(lonmin,latmin)为顶点像元物方坐标的最小值。
5.根据权利要求1或2所述的卫星影像的正射校正方法,其特征在于,步骤S5还包括:建立地理坐标系与目标平面坐标系之间的投影变换关系,具体包括如下步骤:在所述目标平面坐标系下针对数据处理区域建立格网;
取所述格网的至少三个顶点的坐标;
计算所述三个顶点的坐标在地理坐标系下的同名点坐标;
根据所述三个顶点的坐标以及相应的同名点坐标,计算变换系数;
根据所述变换系数建立投影逆变换线性方程组。
6.根据权利要求5所述的卫星影像的正射校正方法,其特征在于,所述投影逆变换线性方程组为:B=Ba+θy(y-ya)+θx(x-xa)
其中,(B,L),(x,y)为待投影要素分别在地理坐标系与目标平面坐标系下的坐标;
θx,θy分别为不同的变换系数,所述变换系数的计算公式为:θx=(Bb-Ba)/(xb-xa);θy=(Bd-Ba)/(yd-ya);
其中,取所述至少三个顶点的坐标A(xa,ya)、B(xb,yb)、D(xd,yd),且xa=xd,ya=yb,所述顶点坐标A(xa,ya)、B(xb,yb)、D(xd,yd)在地理坐标系下的同名点坐标为A’(Ba,La),B’(Bb,Lb),D’(Bd,Ld)。
7.根据权利要求1或2所述的卫星影像的正射校正方法,其特征在于,所述步骤S6具体包括:根据所述物方坐标最大值与最小值对应的目标平面坐标与所述待校正卫星影像的像元空间分辨率建立校正结果图像二维数组,每个所述校正结果图像二维数组对应所述待校正卫星影像的一个光谱波段;
针对所述校正结果图像二维数组逐元素提取元素中心点的目标平面坐标;
根据所述投影变换关系,计算所述中心点的目标平面坐标对应的地理坐标系下的物方坐标,并从所述数字高程模型数据中提取所述物方坐标位置的高程值,通过所述有理函数模型计算所述物方坐标对应的待校正卫星影像的像方坐标。
8.根据权利要求7所述的卫星影像的正射校正方法,其特征在于,在步骤S6之后还包括:提取与所述像方坐标相邻的若干像元值;
通过双线性内插计算所述像方坐标对应的物方坐标的像元值;
将所述像元值作为对应的所述校正 结果图像二维数组的数组元素;
在目标平面坐标系下,根据所述校正 结果图像二维数组创建多光谱栅格文件。
9.一种卫星影像的正射校正系统,其特征在于,包括:
重组单元,用于在地理坐标系下划定正射校正变换范围;获取所述变换范围内的数字高程模型数据,将所述数字高程模型数据转换到所述地理坐标系下;
估算单元,用于获取待校正卫星影像的有理函数系数,根据所述有理函数系数建立有理函数模型,通过所述有理函数模型正变换估算所述待校正卫星影像的顶点像元物方坐标;
提取单元,用于根据所述待校正卫星影像的顶点像元物方坐标,提取与所述待校正卫星影像对应的所述数字高程模型数据;
重新计算单元,用于将与所述待校正卫星影像对应的所述数字高程模型数据通过有理函数模型正变换,重新计算所述待校正卫星影像的顶点像元物方坐标,从中提取物方坐标最大值与最小值;
变换单元,用于建立所述地理坐标系与目标平面坐标系之间的投影变换关系,根据所述投影变换关系,计算所述物方坐标最大值与最小值对应的目标平面坐标;
校正单元,用于根据所述物方坐标最大值与最小值对应的目标平面坐标与所述待校正卫星影像的像元空间分辨率,建立所述目标平面坐标系与所述待校正卫星影像的像方坐标系之间的对应关系。
说明书 :
卫星影像的正射校正方法和系统
技术领域
背景技术
机(简称PMS,组合幅宽60km),四台16米分辨率多光谱相机(简称WFV,组合幅宽800km),其中WFV影像实现四天的覆盖周期。其中,GF-1 LV1A影像产品具备详细的辐射定标系数、光谱响
应参数、有理多项式系数,以其良好的时空分辨率特性在农情资源与灾害监测等领域具有
广泛的应用前景。GF-1 LV1A卫星影像产品不具备空间参考,且幅宽较大,其几何畸变随着
像元位置与星下点距离增大而逐渐严重,同时考虑到影像更新快速、数据量大等特征,如何
快速执行正射校正,改正地形起伏和传感器误差引起的卫星影像几何畸变,成为GF-1 LV1A
数据处理的基础性工作之一。
量中的共线条件方程为代表,综合考虑传感器成像时引起几何畸变的地形起伏、大气折射、
相机透镜畸变及传感器平台的位置、姿态变化等因素,建立像方和物方的数学关系,成像理
论严密,定位精度较高,但解算复杂,且因为某些商业原因传感器的核心信息和卫星轨道参
数并未公开而难以广泛适用。通用传感器模型是一种依赖于传感器物理模型或地面控制点
建立的独立于不同传感器的通用数学模型,直接采用数学函数形式描述地面点物方坐标和
相应像方坐标间的几何关系,其中发展最快、应用最广和最具代表性的便是有理函数模型
(rational polynomial coefficient,RPC)。国内外许多学者采用RPC模型进行了SPOT-5卫
星、QuickBird卫星、IKONOS卫星等高分辨率卫星影像的正射校正实验,得到了较高精度的
实验结果。在此基础上,有学者对比了多种形式下的RPC模型参数求解精度,证明了RPC模型
对SAR卫星传感器成像几何拟合的一致性;还有学者提出了“地形相关”和“地形无关”的RPC模型参数求解方式,研究了用最小二乘方法求解RPC模型参数的算法;相关的研究还包括
RPC模型精度与误差补偿,控制点动态提取,基于GPU(通用计算编程模型)提高校正效率等
等。
间,这无疑降低了数据处理效率,所以提供一种快速的卫星影像的正射校正方法具有深远
的意义。
发明内容
最小值;
的对应关系。
的每个网格对应唯一的一个分块栅格文件。
二维数组,建立所述高程二维数组的数组元素与所述数字高程模型数据的对应关系,将所
述数字高程模型数据赋值到数组。
待校正卫星影像的一个光谱波段;
函数模型计算所述物方坐标对应的待校正卫星影像的像方坐标。
坐标;
坐标最大值与最小值;
坐标系之间的对应关系。
校正影像范围快速读取DEM数据进入内存数组以辅助正射校正计算;另一方面将基于数值
法的快速投影变换过程集成进入正射校正计算过程,通过减少中间文件的读写提高计算效
率,并可以适应计算区域范围的变化,满足工程化应用要求,输出平面坐标参考下卫星影像
正射校正产品;此外还可以减少硬盘存储空间的压力,减少硬盘损耗。
附图说明
具体实施方式
下步骤:S1:在地理坐标系下划定正射校正变换范围;获取所述变换范围内的数字高程模型
数据,将所述数字高程模型数据转换到所述地理坐标系下;
最小值;
的对应关系。
剪、镶嵌,将所述数字高程模型数据分成若干分块栅格文件,使所述地理格网的每个网格对
应唯一的一个分块栅格文件。
“度”单位,其经度跨度为△lon,纬度跨度为△lat;最小经度为lon0,最大经度为lon0+△lon,最小纬度为lat0,最大纬度为lat0+△lat,均可以被1°整除。在变换范围内建立0.5°*
0.5°的地理格网,将地理格网视为相对独立的DEM数据分块存储单元;针对每个地理格网边
界坐标,其经度最大值、最小值必有一个可以被1°(度)整除,纬度最大值、最小值必有一个
可以被1°(度)整除;依据地理格网边界坐标最小值对地理格网编码,编码方式采用“(纬度
最小值*2)-(经度最小值*2)”。将DEM数据集转换到地理坐标系下,将DEM数据的栅格分辨率
统一重采样到0.00025°,并应用所述地理格网划分DEM数据集,通过裁剪、镶嵌处理使每个
地理格网对应唯一的、边界与地理格网相同的DEM分块栅格文件。所述DEM分块栅格文件相
互独立,行、列数均为2000并且文件名称与格网编码一致。在此基础上,在执行正射校正过
程中可以根据GF-1 LV1A多光谱影像数据边界坐标迅速提取相交的DEM分块栅格文件。
具体地,针对单景GF-1 LV1A多光谱影像数据提取RPC系数,建立RPC模型。通过RPC模型正变
换估算影像顶点像元物方坐标E(lone,late),F(lonf,latf),G(long,latg),H(lonh,lath),从中可以提取物方坐标最大值(lonmax,latmax)与最小值(lonmin,latmin)。
0.01°,将所述顶点像元物方坐标的最小值减小0.01°,作为提取所述数字高程模型数据的
边界范围,计算与该边界范围相交的地理格网在经度方向的数量Gxn以及在纬度方向的数
量Gyn,其中,Gxn与Gyn的计算公式如下:
方向的像元数量为yls,所述高程二维数组为[xls,yls]:
述数字高程模型数据的分块栅格文件像元的对应关系,将所述数字高程模型数据赋值到数
组。如图3所示,所述DEM分块栅格文件在空间分布上无缝邻接,可以根据其在空间上的排列
顺序,以西北顶点为像元起点读取DEM分块栅格文件,建立初始高程二维数组OrigDemArray
[Gox,Goy]。其中,Gox=Gxn*2500;Goy=Gyn*2500。
DemArray[u,w]=OrigDemArray[iniX+u,iniY+w];
小值;具体的,将所述高程二维数组DemArray应用于RPC模型正变换,重新计算待校正单景
GF-1 LV1A多光谱影像顶点像元物方坐标K(lonk,latk),M(lonm,latm),N(lonn,latn),P(lonp,latp),从中提取地理坐标最大值(lonmax’,latmax’)与最小值(lonmin’,latmin’)。
与目标平面坐标系之间的投影变换关系,具体包括如下步骤:在所述目标平面坐标系下针
对数据处理区域建立等公里跨度格网;每个格网经度跨度△x,纬度跨度△y,如图4所示,将每个格网视为相对独立的处理单元取所述格网的至少三个顶点的坐标;计算所述三个顶点
的坐标在地理坐标系下的同名点坐标;根据所述三个顶点的坐标以及相应的同名点坐标,
计算变换系数;根据所述变换系数建立投影逆变换线性方程组;较优的,可以取顶点坐标A
(xa,ya)、B(xb,yb)、C(xc,yc)D(xd,yd),且xa=xd,xb=xc,ya=yb,yc=yd,所述顶点A、B、D在地理坐标系下的同名点坐标为A’(Ba,La),B’(Bb,Lb),D’(Bd,Ld),则所述变换系数的计算公式为:
的对应关系。其中较优的,步骤S6具体包括:根据所述物方坐标最大值与最小值对应的目标
平面坐标与所述待校正卫星影像的像元空间分辨率建立n个校正结果图像二维数组Sarray
[zxn,zyn],每个所述校正结果图像二维数组对应所述待校正卫星影像的一个光谱波段;
zxn=[xmax-xmin/Sxy]+1;zyn=[ymax-ymin/Sxy]+1;其中“[]”为向下取整符号。针对所
述校正结果图像二维数组逐元素提取元素中心点的目标平面坐标;根据所述投影变换关
系,计算所述中心点的目标平面坐标对应的地理坐标系下的物方坐标,并从所述数字高程
模型数据中提取所述物方坐标位置的高程值,通过所述有理函数模型计算所述物方坐标对
应的待校正卫星影像的像方坐标。以一个校正结果图像二维数组Sarray0为例,设Sarray0与
所述待校正GF-1 LV1A影像的第1波段对应,遍历全部数组元素:对于某一数组元素Sarray0
[x0,y0],x0={0,1,...zxn-1},y0={0,1,…zyn-1},提取元素中心点平面坐标(Xg,Yg),应用所述基于格网的线性投影逆变换体系计算对应的地理坐标(Bg,Lg),在此基础上从所述
高程二维数组DemArray提取所述地理坐标位置的高程Hg。将(Bg,Lg,Hg)应用于RPC模型计
算物方坐标对应的待校正单景GF-1 LV1A影像的像方坐标(s,l)。具体的,s=Tx*Scale+
Soff;l=Ty*Lcale+Loff;
述校准结果图像二维数组的数组元素;在目标平面坐标系下,根据所述校准结果图像二维
数组创建多光谱栅格文件。具体的,从所述待校正GF-1 LV1A影像的第1波段数据中,提取与
所述像方坐标(s,l)相邻的若干像元值,应用双线性内插计算所述地理坐标(Bg,Lg)对应的
像元值,并赋值给所述数组元素Sarray0[x0,y0]。对n个校正结果图像二维数组Sarray执行
上述赋值过程。在目标平面坐标参考系下,根据所述n个校正结果图像二维数组Sarray创建
多光谱栅格文件,完成单景GF-1 LV1A影像正射校正过程。
换范围内的数字高程模型数据,将所述数字高程模型数据转换到所述地理坐标系下;估算
单元,用于获取待校正卫星影像的有理函数系数,根据所述有理函数系数建立有理函数模
型,通过所述有理函数模型正变换估算所述待校正卫星影像的顶点像元物方坐标;提取单
元,用于根据所述待校正卫星影像的顶点像元物方坐标,提取与所述待校正卫星影像对应
的所述数字高程模型数据;重新计算单元,用于将与所述待校正卫星影像对应的所述数字
高程模型数据通过有理函数模型正变换,重新计算所述待校正卫星影像的顶点像元物方坐
标,从中提取物方坐标最大值与最小值;变换单元,用于建立所述地理坐标系与目标平面坐
标系之间的投影变换关系,根据所述投影变换关系,计算所述物方坐标最大值与最小值对
应的目标平面坐标;校正单元,用于根据所述物方坐标最大值与最小值对应的目标平面坐
标与所述待校正卫星影像的像元空间分辨率,建立所述目标平面坐标系与所述待校正卫星
影像的像方坐标系之间的对应关系。
校正影像范围快速读取DEM数据进入内存数组以辅助正射校正计算;另一方面将基于数值
法的快速投影变换过程集成进入正射校正计算过程,通过减少中间文件的读写提高计算效
率,并可以适应计算区域范围的变化,满足工程化应用要求,输出平面坐标参考下GF-1
LV1A影像正射校正产品;此外还可以减少硬盘存储空间的压力,减少硬盘损耗。
所限定的范围之内。