CN105091906B - 高分辨率光学推扫卫星稳态重成像传感器校正方法及系统 - Google Patents
高分辨率光学推扫卫星稳态重成像传感器校正方法及系统 Download PDFInfo
- Publication number
- CN105091906B CN105091906B CN201510386827.5A CN201510386827A CN105091906B CN 105091906 B CN105091906 B CN 105091906B CN 201510386827 A CN201510386827 A CN 201510386827A CN 105091906 B CN105091906 B CN 105091906B
- Authority
- CN
- China
- Prior art keywords
- mrow
- imaging
- msub
- model
- virtual
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 181
- 230000003287 optical effect Effects 0.000 title claims abstract description 33
- 238000012937 correction Methods 0.000 claims abstract description 23
- 238000012892 rational function Methods 0.000 claims abstract description 20
- 238000013507 mapping Methods 0.000 claims abstract description 16
- 239000000523 sample Substances 0.000 claims description 52
- 238000000034 method Methods 0.000 claims description 40
- 239000011159 matrix material Substances 0.000 claims description 20
- 238000001444 catalytic combustion detection Methods 0.000 claims description 18
- 230000010354 integration Effects 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000012952 Resampling Methods 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 2
- 238000013461 design Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 206010044565 Tremor Diseases 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 239000012634 fragment Substances 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 101000694017 Homo sapiens Sodium channel protein type 5 subunit alpha Proteins 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Image Processing (AREA)
Abstract
本发明提供一种高分辨率光学推扫卫星稳态重成像传感器校正方法及系统,包括构建单片TDI CCD非稳态严密几何模型和虚拟CCD稳态成像严密几何模型,根据虚拟CCD稳态成像严密几何模型得到有理函数模型,建立单片TDI CCD 影像与虚拟CCD 重成像影像的一一映射关系,从而根据原始的单片TDI CCD影像生成传感器校正后影像。本发明技术方案结合虚拟CCD成像原理,利用多片TDI CCD非稳态几何模型与虚拟单线阵CCD稳态几何模型定位的一致性,实现多片CCD影像的无缝拼接的同时,校正由平台震颤引起的影像变形,提供用户标准景影像和对应高精度RPC参数,便于后续图像应用。
Description
技术领域
本发明属于遥感卫星数据处理技术领域,更确切的说,本发明是在遥感卫星地面数据处理程中,采用的一种高分辨率光学推扫卫星影像稳态重成像的传感器校正技术方案。
背景技术
随着空间分辨率的逐步提高,高分辨率光学卫星的相机结构变得异常复杂。高分辨率一般指空间分辨率优于5米,为了保证成像质量和成像效率,感光元件由传统的单次积分CCD(电荷藕合器件)更新为多次积分TDI CCD(即时间延迟积分CCD),由单条推扫线阵扩充到多条拼接推扫线阵;同时,空间分辨率的提高也使得光学相机的物理焦距都达到了几米甚至十几米,为了减少成像载荷体积,传统的折射式光学系统设计也由同轴三反和离轴三反系统取而代之。由此可知,高分辨率光学遥感卫星的成像载荷存在分片积分异速、积分时间跳变、CCD不共线、主距与视向量不平行等特点。在卫星数据投入应用前,必须对原始成像数据进行传感器校正处理,得到连续、完整、高精度的影像数据,满足后续国土资源调查、地图测绘的精度需求。
随着分辨率的提高,成像积分时间越来越短,成像频率的提升导致推扫成像相机对平台稳定性愈加敏感,平台震颤也成为影响高分辨率光学卫星影像几何精度的重要因素之一。平台震颤直接导致卫星姿态抖动、影像变形,为了不影响影像产品配准、融合应用,传感器校正除了需要校正相机设计引起的拼接和镜头畸变问题,还需要解决由于卫星平台非稳态运动引起的影像变形问题。因此,高分辨率光学卫星的传感器校正需同时考虑相机设计和线阵推扫非稳态成像问题,将带畸变的分片影像生成稳态匀速成像的完整景影像,并生成与之对应的高精度RFM(有理函数模型),为后续高精度几何校正、影像融合等处理与应用提供标准数据产品。
目前,传感器校正方法大多只考虑相机设计,解决由于多片TDI CCD的拼接问题,并未考虑推扫过程中卫星的稳定性状态;针对分片异速和积分时间跳变的问题,通常采样行频归一化进行图像重采样,容易导致图像过度重采样,降低图像质量。另一方面,在卫星存在震颤时,由于姿态的抖动也难以直接生成高精度的RPCs(RFM多项式系数),只能使用严密几何成像模型,降低了数据应用效率。因此,传统的传感器校正方法已经无法满足高分辨率遥感光学卫星高精度几何处理的要求。
发明内容
针对现有技术中高分辨率光学推扫卫星影像传感器校正方式缺陷,本发明提出一种高分辨率光学卫星稳态重成像的传感器校正技术方案。
本发明技术方案提供一种高分辨率光学推扫卫星稳态重成像的传感器校正方法,包括以下步骤,
步骤1,构建单片TDI CCD非稳态严密几何模型,包括对于像点(s,l),根据在轨几何定标获取的探元指向角模型确定探元s在相机坐标系下的指向,通过轨道模型和姿态内插模型获取每一行扫描行对应时刻的轨道和姿态作为外方位元素,从而建立每一片TDICCD影像每一点(s,l)的严密成像几何模型;其中,s为某单片TDI CCD影像上任意一个探元的探元序号,l为某单片TDI CCD影像上任意一行的行号;
步骤2,构建虚拟CCD稳态成像严密几何模型,包括设定与多个单片TDI CCD拼接长度相等的虚拟单线阵CCD,在实际成像基础上,确定虚拟单线阵CCD的内方位元素,并设定卫星匀速推扫成像每一扫描时间以及稳定的轨道模型和姿态模型,从而建立虚拟CCD影像上任一点(s',l')的严格成像几何模型;其中,s'为虚拟单线阵CCD上任意一个探元的探元序号,l'为虚拟CCD影像上任意一行的行号;
步骤3,采用地形独立法,根据虚拟CCD稳态成像严密几何模型得到有理函数模型,代替虚拟CCD稳态成像严密成像几何模型;
步骤4,虚拟CCD稳态重成像,包括以物方作为参考,通过步骤1和3所得结果,建立单片TDI CCD影像与虚拟CCD重成像影像的一一映射关系,从而根据原始的单片TDICCD影像生成传感器校正后影像。
而且,步骤1中,利用探元指向角模型确定探元s在相机坐标系下的指向计算实现如下,
其中,axi,ayi为多项式系数,i=0,1,2,3,ψx(s)和ψy(s)为与探元序号s相应的相机像方矢量在相机坐标系下的指向角。
而且,步骤1中,所采用的轨道模型是以时间为变量的三次多项式模型,姿态模型采用以时间为变量的拉格朗日模型,时间通过行号l确定。
而且,步骤1中所得每一片TDI CCD影像每一点的严密成像几何模型,建立了像点图像坐标(s,l)与对应物点坐标(X,Y,Z)的关系如下,
其中,λ为比例系数,为成像过程中像方矢量与物方矢量之间的旋转矩阵,与成像时间有关,计算如下,
其中,均为3×3的方阵,分别代表卫星本体坐标系到传感器坐标系、J2000坐标系到卫星本体坐标系以及WGS84坐标系到J2000坐标系的旋转矩阵。
而且,步骤2中,虚拟单线阵CCD沿轨方向放置在所有单片TDI CCD的中间位置,垂轨方向视场覆盖了多个单片TDI CCD共同覆盖的视场范围。
而且,步骤2中,稳定的轨道模型和姿态模型采用以成像时间为变量的三次多项式模型,对于虚拟影像上任意一行l',通过行号l'计算成像时间t',从而内插出成像时刻的轨道和姿态。
而且,步骤3中,所述地形独立法实现方式为,首先将虚拟CCD稳态成像的影像范围划分为m×n个格网,则产生了(m+1)×(n+1)个均匀分布的格网点;然后在三维空间确定若干高程面,通过虚拟CCD稳态成像严密几何模型计算格网点对应高程面上的空间三维坐标,得到多个虚拟控制点;最后采用最小二乘法通过虚拟控制点解算有理函数模型的多项式系数。
而且,步骤4中,建立单片TDI CCD影像与虚拟CCD重成像影像的一一映射关系实现方式为,以物方为基准,首先利用步骤1所得原始的单片TDI CCD影像的严密成像几何模型,计算某一扫描行上某一点(s,l)与地面的交点(B,L,H),H为高程;然后利用步骤3所得有理函数模型将(B,L,H)反算到虚拟扫描景坐标(s',l'),然后通过重采样获取虚拟扫描景灰度值。
本发明还相应提供一种高分辨率光学推扫卫星稳态重成像的传感器校正系统,包括以下模块,
物理单片模型构建模块,用于构建单片TDI CCD非稳态严密几何模型,包括对于像点(s,l),根据在轨几何定标获取的探元指向角模型确定探元s在相机坐标系下的指向,通过轨道模型和姿态内插模型获取每一行扫描行对应时刻的轨道和姿态作为外方位元素,从而建立每一片TDI CCD影像每一点(s,l)的严密成像几何模型;其中,s为某单片TDI CCD影像上任意一个探元的探元序号,l为某单片TDI CCD影像上任意一行的行号;
虚拟线阵模型构建模块,用于构建虚拟CCD稳态成像严密几何模型,包括设定与多个单片TDI CCD拼接长度相等的虚拟单线阵CCD,在实际成像基础上,确定虚拟单线阵CCD的内方位元素,并设定卫星匀速推扫成像每一扫描时间以及稳定的轨道模型和姿态模型,从而建立虚拟CCD影像上任一点(s',l')的严格成像几何模型;其中,s'为虚拟单线阵CCD上任意一个探元的探元序号,l'为虚拟CCD影像上任意一行的行号;
有理函数模型生成模块,用于采用地形独立法,根据虚拟CCD稳态成像严密几何模型得到有理函数模型,代替虚拟CCD稳态成像严密成像几何模型;
虚拟CCD稳态重成像模块,用于以物方作为参考,通过物理单片模型构建模块和虚拟线阵模型构建模块所得结果,建立单片TDI CCD影像与虚拟CCD重成像影像的一一映射关系,从而根据原始的单片TDI CCD影像生成传感器校正后影像。
而且,物理单片模型构建模块中,利用探元指向角模型确定探元s在相机坐标系下的指向计算实现如下,
其中,axi,ayi为多项式系数,i=0,1,2,3,ψx(s)和ψy(s)为与探元序号s相应的相机像方矢量在相机坐标系下的指向角。
而且,物理单片模型构建模块中,所采用的轨道模型是以时间为变量的三次多项式模型,姿态模型采用以时间为变量的拉格朗日模型,时间通过行号l确定。
而且,物理单片模型构建模块中所得每一片TDI CCD影像每一点的严密成像几何模型,建立了像点图像坐标(s,l)与对应物点坐标(X,Y,Z)的关系如下,
其中,λ为比例系数,为成像过程中像方矢量与物方矢量之间的旋转矩阵,与成像时间有关,计算如下,
其中,均为3×3的方阵,分别代表卫星本体坐标系到传感器坐标系、J2000坐标系到卫星本体坐标系以及WGS84坐标系到J2000坐标系的旋转矩阵。
而且,虚拟线阵模型构建模块中,虚拟单线阵CCD沿轨方向放置在所有单片TDICCD的中间位置,垂轨方向视场覆盖了多个单片TDI CCD共同覆盖的视场范围。
而且,虚拟线阵模型构建模块中,稳定的轨道模型和姿态模型采用以成像时间为变量的三次多项式模型,对于虚拟影像上任意一行l',通过行号l'计算成像时间t',从而内插出成像时刻的轨道和姿态。
而且,有理函数模型生成模块中,所述地形独立法实现方式为,首先将虚拟CCD稳态成像的影像范围划分为m×n个格网,则产生了(m+1)×(n+1)个均匀分布的格网点;然后在三维空间确定若干高程面,通过虚拟CCD稳态成像严密几何模型计算格网点对应高程面上的空间三维坐标,得到多个虚拟控制点;最后采用最小二乘法通过虚拟控制点解算有理函数模型的多项式系数。
而且,虚拟CCD稳态重成像模块中,建立单片TDI CCD影像与虚拟CCD重成像影像的一一映射关系实现方式为,以物方为基准,首先利用原始的单片TDI CCD影像的严密成像几何模型,计算某一扫描行上某一点(s,l)与地面的交点(B,L,H),H为高程;然后利用有理函数模型将(B,L,H)反算到虚拟扫描景坐标(s',l'),然后通过重采样获取虚拟扫描景灰度值。
本发明提供一种高分辨率光学推扫卫星影像稳态重成像的传感器校正技术方案,解决高分辨率光学推扫卫星非稳态成像的传感器校正问题。本发明结合虚拟CCD成像原理,利用多片TDI CCD非稳态几何模型与虚拟单线阵CCD稳态几何模型定位的一致性,实现多片CCD影像的无缝拼接的同时,校正由平台震颤引起的影像变形,提供用户标准景影像和对应高精度RPC参数,便于后续图像应用。
附图说明
图1为本发明实施例的内定向模型示意图。
图2为本发明实施例的虚拟单线阵CCD与物理多片CCD的关系示意图。
图3为本发明实施例的虚拟影像与原始影像映射过程。
图4为本发明实施例的高分辨率光学推扫卫星影像稳态重成像的传感器校正原理图。
图5为本发明实施例的高分辨率光学推扫卫星影像稳态重成像的传感器校正流程。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述。
本发明提出一种高分辨率光学卫星稳态重成像的传感器校正方法,利用无畸变的虚拟单线阵CCD代替原始有畸变的分片CCD影像,将原来非稳态(积分跳变、分片异速、姿态抖动)的推扫成像条件转变为稳定、匀速(积分不变、分片同速、姿态平滑)的推扫成像条件,在稳态条件下进行重成像,从而得到无畸变的完整景影像,同时得到稳态成像的高精度RPCs。
具体实施时,本发明技术方案可由本领域技术人员采用计算机软件技术实现自动运行流程。参见图5,实施例的具体实现如下:
步骤1,单片异速严密几何模型构建:即构建单片TDI CCD非稳态严密几何模型,对于像点(s,l),根据在轨几何定标获取的探元指向角模型确定探元s在相机坐标系下的指向,通过高精度轨道和姿态内插模型获取每一行扫描行对应时刻的轨道和姿态(外方位元素),从而建立每一片TDI CCD影像每一点的严密成像几何模型。
实施例首先进行原始分片非稳态成像建模,即单片TDI CCD非稳态严密几何模型构建,是根据轨道数据、姿态数据、相机成像时间数据,通过高精度姿轨内插模型获取每一行扫描行对应时刻的轨道和姿态(外方位元素),即在姿态、轨道高精度建模的基础上进行姿态、轨道高精度内插;然后再利用分片相机文件中的相机参数(内方位元素),并根据影像坐标(s,l)分别建立每一片TDI CCD每一成像行的严密成像几何模型。
对每一片CCD的每一扫描行建立严密成像几何模型。根据摄影测量原理,严密成像几何模型依据共线方程进行构建,包括两个部分:内定向模型(像方矢量确定)和外定向模型(物方矢量确定)。
1)内定向模型
根据在轨几何定标确定的探元指向角模型,给定探元序号即可确定相机像方矢量在相机坐标系O1-X1Y1Z1下的指向角ψx(s)、ψy(s)(如图1所示ψx、ψy),计算公式如下:
其中,s为探元序号,axi,ayi为多项式系数,i=0,1,2,3。需要注意的是,每一片CCD都具有独立的多项式系数。
因此,对于单片TDI CCD影像上的任一点(s,l),l为为单片TDI CCD影像上任意一行的行号,s为列号即单片TDI CCD影像上任意一个探元的探元序号,则根据公式(1)可计算该点在相机坐标系下的探元指向角。进一步可计算像方矢量Vimg。
2)外定向模型
外定向模型包括轨道模型和姿态模型。高分辨率卫星的轨道由GPS测得卫星在WGS84坐标系下的位置(Xs、Ys、Zs)和速度(Vx、Vy、Vz),姿态由星敏感器和陀螺组合定姿测得卫星在J2000惯性坐标系下的姿态四元数[q0 q1 q2 q3]T。实施例采用的高精度轨道模型是以时间为变量的三次多项式模型,姿态模型采用以时间为变量的拉格朗日模型,时间通过行号l可以确定。
轨道模型采用三次多项式拟合得到,以保证轨道的平滑特性。传感器校正以景为单位,取该景影像成像弧段的轨道观测点对该弧段的位置和速度作为观测值,利用最小二乘法进行三次多项式拟合,分别求解位置和速度的多项式系数。那么在该弧段内任意时刻t的位置和速度都可以通过时间内插得到。
姿态模型采用拉格朗日内插模型,避免采用多项式损失姿态观测精度。任意时刻t的姿态角利用最邻近的n个历元处(通常取8个历元)的姿态内插计算。
通过行号l即可获得该扫描行的时间t。考虑到TDI CCD积分时间跳变问题,首先需要确定l的积分时间段n,根据公式(3)计算成像时刻。
其中t0为起始行成像时刻,Δti(i=1,2,...,n)为每个积分时间段扫描行的积分时间,li(i=1,2,...,n-1)为每段积分时间内的扫描行数。基于扫描行时间模型和轨道、姿态模型,根据原始影像上像点所在的行号得到成像时刻,进而内插计算卫星的位置矢量、速度矢量和姿态角元素。
3)分片TDI CCD严密成像几何方程
根据内定向模型和外定向模型可以得到分片严密成像几何模型,建立像点图像坐标(s,l)与其对应物点坐标(X,Y,Z)的关系。
其中,λ为比例系数,为成像过程中像方矢量与物方矢量之间的旋转矩阵,与成像时间有关,即
均为3×3的方阵,分别代表卫星本体坐标系到传感器坐标系、J2000坐标系到卫星本体坐标系以及WGS84坐标系到J2000坐标系的旋转矩阵。表示相机的安装矩阵,由在轨几何定标确定。
涉及天球坐标系和地球坐标系的转换,需要根据成像时刻t计算地球旋转参数,如公式(5)所示。
其中,P(t)、N(t)、R(t)、W(t)分别表示相应的成像时刻t岁差矩阵、章动矩阵、地球旋转矩阵和极移矩阵。
步骤2,虚拟CCD影像几何模型构建:构建虚拟CCD稳态成像严密几何模型,在实际成像基础上,确定虚拟单线阵CCD的内方位元素,并设定卫星匀速推扫成像每一扫描时间以及稳定的轨道模型和姿态模型,从而建立虚拟CCD稳态成像的严格几何模型。
类似单片异速严密成像几何模型,虚拟CCD影像的几何模型也是从内定向和外定向模型入手,建立严密几何模型;不同的是,虚拟CCD不再是多个单片独立的物理CCD,而是与多个单片CCD拼接长度相等的完整虚拟线阵CCD(如图2),且在成像时虚拟CCD外方位元素在真实轨道和姿态上更加平滑、稳定。
1)虚拟CCD内定向模型
如图2所示,虚拟CCD的视场覆盖了多个单片CCD共同覆盖的视场,且虚拟CCD是一个完全理想的线阵,因此也无CCD变形,镜头无畸变,焦距无缩放,相机参数可以完全按照设计值进行设计。为了减小虚拟CCD偏视场角与单片CCD的差异,可设定将虚拟CCD放置在所有CCD的中间位置,例如图2中虚拟单线阵CCD置于两行TDI CCD之间。为了工程应用的便利,同样可以根据虚拟CCD设计位置利用三次多项式模型进行拟合,得到虚拟CCD的指向角模型(如公式1),
即
需要注意的是多项式系axi,ayi数值与物理CCD指向角模型系数不一样。对于虚拟CCD上任意一个探元s'均可通过指向角模型确定其像方指向。
2)虚拟CCD外定向
为了保证虚拟CCD成像的稳定性,令每一扫描的积分时间相同,即不存在积分时间跳变,积分时间可以取该景影像积分时间的均值,那么任意行的时间可有起始行时间通过积分时间累计。同时,轨道和姿态均采用以成像时间为变量的三次多项式模型拟合,得到平滑稳定的运动轨迹和指向。通过这两个方面保证虚拟CCD在完全理想、稳定的状态下成像,从而可以获取无畸变的影像。对于虚拟影像上任意一行l',可通过行号计算成像时间t',从而获取成像时刻的轨道和姿态。结合虚拟CCD的内定向模型,从而则可建立虚拟CCD影像上任一点的严密成像几何模型(形式如公式4),
即
步骤3,稳态成像RPCs生成:采用地形独立法,根据虚拟CCD稳态成像严密几何模型得到有理函数模型(RFM),代替虚拟CCD稳态成像严密几何模型。
RFM的生成方法采用地形独立法,该方法是一种依赖严密成像几何模型的方法。首先将虚拟CCD稳态成像的影像范围按照一定的间隔划分为m×n个格网,则产生了(m+1)×(n+1)个均匀分布的格网点;然后在三维空间确定足够数量一定高程间隔的高程面,通过虚拟CCD稳态成像严格几何成像模型计算格网点对应高程面上的空间三维坐标,即可得到足够数量的虚拟控制点;最后采用最小二乘法即可通过虚拟控制点解算RFM多项式系数RPC。
具体实施时,本领域技术人员可自行预设划分的间隔大小或m、n的数值。
由于虚拟CCD在完全稳定的条件下成像,那么严密成像模型可由RFM高精度替代,从而提高正反算的计算效率。RFM是一个与传感器无关的高分辨率卫星严密几何模型的替代模型,直接建立影像坐标与地面地理坐标的关系,可通过严密几何模型与多个高程面交会得到的点作为虚拟控制点,解算80个RPCs。
为了保证计算的稳定性,RFM将像点图像坐标(s',l')和经纬度坐标(B,L)和椭球高H进行正则化处理,使坐标范围在[-1,1]之间。
其中,ln,sn,U,V,W为l',s',B,L,H归一化结果,LineOff,SampleOff,LonOff,LatOff,HeiOff为归一化偏移系数,LineScale,SampleScale,LonScale,LatScale,HeiScale为归一化比例系数,NumL(U,V,W)和DenL(U,V,W)表示列方向有理多项式分子和分母,NumS(U,V,W)和DenS(U,V,W)表示行方向有理多项式分子和分母,具体形式如下:
NumL(U,V,W)=a1+a2V+a3U+a4W+a5VU+a6VW+a7UW+a8V2+a9U2+
a10W2+a11UVW+a12V3+a13VU2+a14VW2+a15V2U+a16U3+
a17UW2+a18V2W+a19U2W+a20W3
DenL(U,V,W)=b1+b2V+b3U+b4W+b5VU+b6VW+b7UW+b8V2+b9U2+
b10W2+b11UVWb12V3+b13VU2+b14VW2+b15V2U+b16U3+
b17UW2+b18V2W+b19U2W+b20W3
NumS(U,V,W)=c1+c2V+c3U+c4W+c5VU+c6VW+c7UW+c8V2+c9U2+
c10W2+c11UVW+c12V3+c13VU2+c14VW2+c15V2U+c16U3+
c17UW2+c18V2W+c19U2W+c20W3
DenS(U,V,W)=d1+d2V+d3U+d4W+d5VU+d6VW+d7UW+d8V2+d9U2+
d10W2+d11UVW+d12V3+d13VU2+d14VW2+d15V2U+d16U3+
d17UW2+d18V2W+d19U2W+d20W3
其中,ai,bi,ci,di(i=1,2…20)为多项式系数。
步骤4,虚拟稳态扫描景生成:虚拟CCD稳态重成像,以物方作为参考,通过步骤1和3中的几何模型,建立单片TDI CCD影像与虚拟CCD重成像影像的一一映射关系,从而可对原始单片CCD影像进行灰度重采样,生成稳态重成像影像,即传感器校正后影像。
通过上述(1)、(2)所得两个成像模型,即建立了虚拟扫描景与原始影像的像点坐标对应关系,实现了基于物方定位一致性的分片与虚拟影像映射。虚拟扫描景的生成过程实际上是在模拟虚拟CCD线阵推扫成像的过程。图4中,(a)部分为虚拟单线阵CCD成像示意图,(b)部分为非稳态成像示意图,(c)部分为稳态重成像示意图。如图4中(a)所示,虚拟CCD与原始分片物理CCD在同一成像面上,沿飞行方向依次对地成像。在t2时刻,虚拟CCD线阵对地面上一条连续的投影线AB成像,并记录为虚拟扫描景上的一条扫描行。同时,AB也是由t1、t3时刻分片物理CCD对地成像构成。因此二者可以通过几何模型进行一一映射g(s,l)的关联。当成像姿态不稳定时,如虚拟CCD和原始CCD影像共享一套外方位参数,则虚拟CCD影像仍存在影像畸变(如图4中(b)所示);若虚拟CCD采用稳态外方位参数,则通过虚拟CCD重成像可获取无畸变影像(如图4中(c)所示)。
考虑到单片TDI CCD采用严格模型,为了避免严格模型反算不收敛的问题,本发明实现虚拟稳态扫描景影像与原始影像的映射过程如图3所示。首先利用原始影像严密几何模型(即式(4))计算某一扫描行上某一点(s,l)与地面的交点(B,L,H),高程H通过STRM DEM获取。然后利用虚拟稳态成像RFM(即式(6))将(B,L,H)反算到虚拟扫描景坐标(s',l'),然后通过重采样方法(如双线性内插、双三次卷积等)获取虚拟扫描景灰度值,即灰度重采样。如此,得到整景虚拟稳态成像影像,即传感器校正影像。
具体实施时,还可以采用模块化方式提供相应系统。本发明实施例相应提供一种高分辨率光学推扫卫星稳态重成像的传感器校正系统,包括以下模块,
物理单片模型构建模块,用于构建单片TDI CCD非稳态严密几何模型,包括对于像点(s,l),根据在轨几何定标获取的探元指向角模型确定探元s在相机坐标系下的指向,通过轨道模型和姿态内插模型获取每一行扫描行对应时刻的轨道和姿态作为外方位元素,从而建立每一片TDI CCD影像每一点(s,l)的严密成像几何模型;其中,s为某单片TDI CCD影像上任意一个探元的探元序号,l为某单片TDI CCD影像上任意一行的行号;
虚拟线阵模型构建模块,用于构建虚拟CCD稳态成像严密几何模型,包括设定与多个单片TDI CCD拼接长度相等的虚拟单线阵CCD,在实际成像基础上,确定虚拟单线阵CCD的内方位元素,并设定卫星匀速推扫成像每一扫描时间以及稳定的轨道模型和姿态模型,从而建立虚拟CCD影像上任一点(s',l')的严格成像几何模型;其中,s'为虚拟单线阵CCD上任意一个探元的探元序号,l'为虚拟CCD影像上任意一行的行号;
有理函数模型生成模块,用于采用地形独立法,根据虚拟CCD稳态成像严密几何模型得到有理函数模型,代替虚拟CCD稳态成像严密成像几何模型;
虚拟CCD稳态重成像模块,用于以物方作为参考,通过物理单片模型构建模块和虚拟线阵模型构建模块所得结果,建立单片TDI CCD影像与虚拟CCD重成像影像的一一映射关系,从而根据原始的单片TDI CCD影像生成传感器校正后影像。
各模块具体实现可参见相应步骤,本发明不予赘述。
本文中所描述的具体实施例仅仅是对本发明精神做举例说明。本发明技术领域的技术人员可以对所描述的具体实施例做各种各样的修改补充或者采用类似的方式替代,但并不会偏离本发明精神或者超越所附权利要求书所定义的范围。
Claims (12)
1.一种高分辨率光学推扫卫星稳态重成像的传感器校正方法,其特征在于:包括以下步骤,
步骤1,构建单片TDI CCD非稳态严密几何模型,包括对于像点(s,l),根据在轨几何定标获取的探元指向角模型确定探元s在相机坐标系下的指向,通过轨道模型和姿态内插模型获取每一行扫描行对应时刻的轨道和姿态作为外方位元素,从而建立每一片TDI CCD影像每一点(s,l)的严密成像几何模型;其中,s为某单片TDI CCD影像上任意一个探元的探元序号,l为某单片TDI CCD影像上任意一行的行号;
步骤2,构建虚拟CCD稳态成像严密几何模型,包括设定与多个单片TDI CCD拼接长度相等的虚拟单线阵CCD,在实际成像基础上,确定虚拟单线阵CCD的内方位元素,并设定卫星匀速推扫成像每一扫描时间以及稳定的轨道模型和姿态模型,从而建立虚拟CCD影像上任一点(s',l')的严格成像几何模型;其中,s'为虚拟单线阵CCD上任意一个探元的探元序号,l'为虚拟CCD影像上任意一行的行号;
步骤3,采用地形独立法,根据虚拟CCD稳态成像严密几何模型得到有理函数模型,代替虚拟CCD稳态成像严密成像几何模型;
所述地形独立法实现方式为,首先将虚拟CCD稳态成像的影像范围划分为m×n个格网,则产生了(m+1)×(n+1)个均匀分布的格网点;然后在三维空间确定若干高程面,通过虚拟CCD稳态成像严密几何模型计算格网点对应高程面上的空间三维坐标,得到多个虚拟控制点;最后采用最小二乘法通过虚拟控制点解算有理函数模型的多项式系数;
步骤4,虚拟CCD稳态重成像,包括以物方作为参考,通过步骤1和3所得结果,建立单片TDI CCD影像与虚拟CCD重成像影像的一一映射关系,从而根据原始的单片TDI CCD影像生成传感器校正后影像;
建立单片TDI CCD影像与虚拟CCD重成像影像的一一映射关系实现方式为,以物方为基准,首先利用步骤1所得原始的单片TDI CCD影像的严密成像几何模型,计算某一扫描行上某一点(s,l)与地面的交点(B,L,H),H为高程,(B,L)为经纬度坐标;然后利用步骤3所得有理函数模型将(B,L,H)反算到虚拟扫描景坐标(s',l'),然后通过重采样获取虚拟扫描景灰度值。
2.根据权利要求1所述高分辨率光学推扫卫星稳态重成像的传感器校正方法,其特征在于:步骤1中,利用探元指向角模型确定探元s在相机坐标系下的指向计算实现如下,
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&psi;</mi>
<mi>x</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>ax</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>ax</mi>
<mn>1</mn>
</msub>
<mo>&times;</mo>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>ax</mi>
<mn>2</mn>
</msub>
<mo>&times;</mo>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>ax</mi>
<mn>3</mn>
</msub>
<mo>&times;</mo>
<msup>
<mi>s</mi>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&psi;</mi>
<mi>y</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>ay</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>ay</mi>
<mn>1</mn>
</msub>
<mo>&times;</mo>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>ay</mi>
<mn>2</mn>
</msub>
<mo>&times;</mo>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>ay</mi>
<mn>3</mn>
</msub>
<mo>&times;</mo>
<msup>
<mi>s</mi>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,axi,ayi为多项式系数,i=0,1,2,3,ψx(s)和ψy(s)为与探元序号s相应的相机像方矢量在相机坐标系下的指向角。
3.根据权利要求1所述高分辨率光学推扫卫星稳态重成像的传感器校正方法,其特征在于:步骤1中,所采用的轨道模型是以时间为变量的三次多项式模型,姿态模型采用以时间为变量的拉格朗日模型,时间通过行号l确定。
4.根据权利要求2所述高分辨率光学推扫卫星稳态重成像的传感器校正方法,其特征在于:步骤1中所得每一片TDI CCD影像每一点的严密成像几何模型,建立了像点图像坐标(s,l)与对应物点坐标(X,Y,Z)的关系如下,
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mi>x</mi>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mi>y</mi>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mi>&lambda;</mi>
<mo>&CenterDot;</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>o</mi>
<mi>b</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>g</mi>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>X</mi>
<mo>-</mo>
<msub>
<mi>X</mi>
<mi>s</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>Y</mi>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mi>s</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>Z</mi>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mi>s</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,λ为比例系数,(Xs,Ys,Zs)为卫星在WGS84坐标系下的位置,为成像过程中像方矢量与物方矢量之间的旋转矩阵,与成像时间有关,计算如下,
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>o</mi>
<mi>b</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>g</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>B</mi>
<mi>o</mi>
<mi>d</mi>
<mi>y</mi>
</mrow>
<mrow>
<mi>s</mi>
<mi>e</mi>
<mi>n</mi>
<mi>s</mi>
<mi>o</mi>
<mi>r</mi>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>J</mi>
<mn>2000</mn>
</mrow>
<mrow>
<mi>B</mi>
<mi>o</mi>
<mi>d</mi>
<mi>y</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>W</mi>
<mi>G</mi>
<mi>S</mi>
<mn>84</mn>
</mrow>
<mrow>
<mi>J</mi>
<mn>2000</mn>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,均为3×3的方阵,分别代表卫星本体坐标系到传感器坐标系、J2000坐标系到卫星本体坐标系以及WGS84坐标系到J2000坐标系的旋转矩阵。
5.根据权利要求1所述高分辨率光学推扫卫星稳态重成像的传感器校正方法,其特征在于:步骤2中,虚拟单线阵CCD沿轨方向放置在所有单片TDI CCD的中间位置,垂轨方向视场覆盖了多个单片TDI CCD共同覆盖的视场范围。
6.根据权利要求1所述高分辨率光学推扫卫星稳态重成像的传感器校正方法,其特征在于:步骤2中,稳定的轨道模型和姿态模型采用以成像时间为变量的三次多项式模型,对于虚拟影像上任意一行l',通过行号l'计算成像时间t',从而内插出成像时刻的轨道和姿态。
7.一种高分辨率光学推扫卫星稳态重成像的传感器校正系统,其特征在于:包括以下模块,
物理单片模型构建模块,用于构建单片TDI CCD非稳态严密几何模型,包括对于像点(s,l),根据在轨几何定标获取的探元指向角模型确定探元s在相机坐标系下的指向,通过轨道模型和姿态内插模型获取每一行扫描行对应时刻的轨道和姿态作为外方位元素,从而建立每一片TDI CCD影像每一点(s,l)的严密成像几何模型;其中,s为某单片TDI CCD影像上任意一个探元的探元序号,l为某单片TDI CCD影像上任意一行的行号;
虚拟线阵模型构建模块,用于构建虚拟CCD稳态成像严密几何模型,包括设定与多个单片TDI CCD拼接长度相等的虚拟单线阵CCD,在实际成像基础上,确定虚拟单线阵CCD的内方位元素,并设定卫星匀速推扫成像每一扫描时间以及稳定的轨道模型和姿态模型,从而建立虚拟CCD影像上任一点(s',l')的严格成像几何模型;其中,s'为虚拟单线阵CCD上任意一个探元的探元序号,l'为虚拟CCD影像上任意一行的行号;
有理函数模型生成模块,用于采用地形独立法,根据虚拟CCD稳态成像严密几何模型得到有理函数模型,代替虚拟CCD稳态成像严密成像几何模型;
所述地形独立法实现方式为,首先将虚拟CCD稳态成像的影像范围划分为m×n个格网,则产生了(m+1)×(n+1)个均匀分布的格网点;然后在三维空间确定若干高程面,通过虚拟CCD稳态成像严密几何模型计算格网点对应高程面上的空间三维坐标,得到多个虚拟控制点;最后采用最小二乘法通过虚拟控制点解算有理函数模型的多项式系数;
虚拟CCD稳态重成像模块,用于以物方作为参考,通过物理单片模型构建模块和虚拟线阵模型构建模块所得结果,建立单片TDI CCD影像与虚拟CCD重成像影像的一一映射关系,从而根据原始的单片TDI CCD影像生成传感器校正后影像;
建立单片TDI CCD影像与虚拟CCD重成像影像的一一映射关系实现方式为,以物方为基准,首先利用原始的单片TDI CCD影像的严密成像几何模型,计算某一扫描行上某一点(s,l)与地面的交点(B,L,H),H为高程,(B,L)为经纬度坐标;然后利用有理函数模型将(B,L,H)反算到虚拟扫描景坐标(s',l'),然后通过重采样获取虚拟扫描景灰度值。
8.根据权利要求7所述高分辨率光学推扫卫星稳态重成像的传感器校正系统,其特征在于:物理单片模型构建模块中,利用探元指向角模型确定探元s在相机坐标系下的指向计算实现如下,
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&psi;</mi>
<mi>x</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>ax</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>ax</mi>
<mn>1</mn>
</msub>
<mo>&times;</mo>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>ax</mi>
<mn>2</mn>
</msub>
<mo>&times;</mo>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>ax</mi>
<mn>3</mn>
</msub>
<mo>&times;</mo>
<msup>
<mi>s</mi>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&psi;</mi>
<mi>y</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>ay</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>ay</mi>
<mn>1</mn>
</msub>
<mo>&times;</mo>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>ay</mi>
<mn>2</mn>
</msub>
<mo>&times;</mo>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>ay</mi>
<mn>3</mn>
</msub>
<mo>&times;</mo>
<msup>
<mi>s</mi>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,axi,ayi为多项式系数,i=0,1,2,3,ψx(s)和ψy(s)为与探元序号s相应的相机像方矢量在相机坐标系下的指向角。
9.根据权利要求7所述高分辨率光学推扫卫星稳态重成像的传感器校正系统,其特征在于:物理单片模型构建模块中,所采用的轨道模型是以时间为变量的三次多项式模型,姿态模型采用以时间为变量的拉格朗日模型,时间通过行号l确定。
10.根据权利要求8所述高分辨率光学推扫卫星稳态重成像的传感器校正系统,其特征在于:物理单片模型构建模块中所得每一片TDI CCD影像每一点的严密成像几何模型,建立了像点图像坐标(s,l)与对应物点坐标(X,Y,Z)的关系如下,
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mi>x</mi>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mi>y</mi>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mi>&lambda;</mi>
<mo>&CenterDot;</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>o</mi>
<mi>b</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>g</mi>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>X</mi>
<mo>-</mo>
<msub>
<mi>X</mi>
<mi>s</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>Y</mi>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mi>s</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>Z</mi>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mi>s</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,λ为比例系数,(Xs,Ys,Zs)为卫星在WGS84坐标系下的位置,为成像过程中像方矢量与物方矢量之间的旋转矩阵,与成像时间有关,计算如下,
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>o</mi>
<mi>b</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>g</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>B</mi>
<mi>o</mi>
<mi>d</mi>
<mi>y</mi>
</mrow>
<mrow>
<mi>s</mi>
<mi>e</mi>
<mi>n</mi>
<mi>s</mi>
<mi>&omega;</mi>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>J</mi>
<mn>2000</mn>
</mrow>
<mrow>
<mi>B</mi>
<mi>o</mi>
<mi>d</mi>
<mi>y</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>W</mi>
<mi>G</mi>
<mi>S</mi>
<mn>84</mn>
</mrow>
<mrow>
<mi>J</mi>
<mn>2000</mn>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,均为3×3的方阵,分别代表卫星本体坐标系到传感器坐标系、J2000坐标系到卫星本体坐标系以及WGS84坐标系到J2000坐标系的旋转矩阵。
11.根据权利要求7所述高分辨率光学推扫卫星稳态重成像的传感器校正系统,其特征在于:虚拟线阵模型构建模块中,虚拟单线阵CCD沿轨方向放置在所有单片TDI CCD的中间位置,垂轨方向视场覆盖了多个单片TDI CCD共同覆盖的视场范围。
12.根据权利要求7所述高分辨率光学推扫卫星稳态重成像的传感器校正系统,其特征在于:虚拟线阵模型构建模块中,稳定的轨道模型和姿态模型采用以成像时间为变量的三次多项式模型,对于虚拟影像上任意一行l',通过行号l'计算成像时间t',从而内插出成像时刻的轨道和姿态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510386827.5A CN105091906B (zh) | 2015-06-30 | 2015-06-30 | 高分辨率光学推扫卫星稳态重成像传感器校正方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510386827.5A CN105091906B (zh) | 2015-06-30 | 2015-06-30 | 高分辨率光学推扫卫星稳态重成像传感器校正方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105091906A CN105091906A (zh) | 2015-11-25 |
CN105091906B true CN105091906B (zh) | 2018-03-02 |
Family
ID=54572855
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510386827.5A Active CN105091906B (zh) | 2015-06-30 | 2015-06-30 | 高分辨率光学推扫卫星稳态重成像传感器校正方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105091906B (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105466400B (zh) * | 2015-12-16 | 2018-05-04 | 贵州天地通科技有限公司 | 一种利用rpc探测多源卫星影像同名像点成像交会角方法 |
CN107945235B (zh) * | 2017-10-17 | 2022-02-01 | 许昌学院 | 一种高轨大面阵静止卫星影像几何定位仿真方法 |
CN107705267B (zh) * | 2017-10-18 | 2020-06-26 | 中国科学院电子学研究所 | 基于控制矢量的光学卫星影像几何校正方法 |
CN108898565B (zh) * | 2018-07-10 | 2021-09-24 | 中国科学院长春光学精密机械与物理研究所 | Tdi ccd相机摆扫成像影像几何畸变还原的反变换方法 |
CN109696182A (zh) * | 2019-01-23 | 2019-04-30 | 张过 | 一种星载推扫式光学传感器内方位元素定标方法 |
CN110030976B (zh) * | 2019-04-08 | 2020-10-30 | 武汉大学 | 保持原始分辨率的遥感虚拟线阵参数提取及影像拼接方法 |
CN110211054B (zh) * | 2019-04-28 | 2021-01-15 | 张过 | 一种星载推扫式光学传感器无畸变影像制作方法 |
CN110956588B (zh) * | 2019-09-23 | 2022-08-05 | 武汉大学 | 基于加密点最短距离的影像高精度几何纠正方法 |
CN110986998B (zh) * | 2019-10-28 | 2021-09-14 | 武汉大学 | 一种基于有理函数模型的卫星视频相机在轨几何定标方法 |
CN111538051B (zh) * | 2020-04-30 | 2022-08-26 | 中国科学院微小卫星创新研究院 | 一种摆扫大幅宽光学卫星精准处理方法 |
CN111521197B (zh) * | 2020-04-30 | 2022-02-15 | 中国科学院微小卫星创新研究院 | 一种摆扫大幅宽光学卫星传感器校正方法 |
CN111612693B (zh) * | 2020-05-19 | 2023-03-14 | 中国科学院微小卫星创新研究院 | 一种旋转大幅宽光学卫星传感器校正方法 |
CN111862332B (zh) * | 2020-07-30 | 2024-06-18 | 武汉多普云科技有限公司 | 一种卫星影像通用成像模型拟合误差的改正方法及系统 |
CN112212833B (zh) * | 2020-08-28 | 2021-07-09 | 中国人民解放军战略支援部队信息工程大学 | 机械拼接型tdi ccd推扫相机整体几何平差方法 |
CN113160312B (zh) * | 2021-02-07 | 2024-01-26 | 中国测绘科学研究院 | 一种卫星线阵影像rpc模型定位方法 |
CN113393499B (zh) * | 2021-07-12 | 2022-02-01 | 自然资源部国土卫星遥感应用中心 | 一种高分七号卫星全色影像和多光谱影像自动配准方法 |
CN115037420B (zh) * | 2022-08-12 | 2022-11-22 | 中国科学院空天信息创新研究院 | 卫星异常数据处理方法、装置、设备及介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103674063A (zh) * | 2013-12-05 | 2014-03-26 | 中国资源卫星应用中心 | 一种光学遥感相机在轨几何定标方法 |
CN103697864A (zh) * | 2013-12-27 | 2014-04-02 | 武汉大学 | 一种基于大虚拟相机的窄视场双相机影像拼接方法 |
CN103914808A (zh) * | 2014-03-14 | 2014-07-09 | 国家测绘地理信息局卫星测绘应用中心 | 一种资源三号卫星三线阵影像和多光谱影像的拼接方法 |
-
2015
- 2015-06-30 CN CN201510386827.5A patent/CN105091906B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103674063A (zh) * | 2013-12-05 | 2014-03-26 | 中国资源卫星应用中心 | 一种光学遥感相机在轨几何定标方法 |
CN103697864A (zh) * | 2013-12-27 | 2014-04-02 | 武汉大学 | 一种基于大虚拟相机的窄视场双相机影像拼接方法 |
CN103914808A (zh) * | 2014-03-14 | 2014-07-09 | 国家测绘地理信息局卫星测绘应用中心 | 一种资源三号卫星三线阵影像和多光谱影像的拼接方法 |
Non-Patent Citations (3)
Title |
---|
基于虚拟线阵的ZY_102C卫星HR相机内视场拼接方法;潘俊等;《武汉大学学报.信息科学版》;20150430;第40卷(第4期);436-443 * |
资源三号测绘卫星三线阵成像几何模型构建与精度初步验证;唐新明等;《测绘学报》;20120430;第41卷(第2期);194 * |
资源三号测绘卫星传感器校正产品生产方法研究;唐新明等;《武汉大学学报.信息科学版》;20140331;第39卷(第3期);290 * |
Also Published As
Publication number | Publication date |
---|---|
CN105091906A (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105091906B (zh) | 高分辨率光学推扫卫星稳态重成像传感器校正方法及系统 | |
CN110500995B (zh) | 利用rpc参数建立高分辨率卫星影像等效几何成像模型的方法 | |
CN106403902B (zh) | 一种星地协同的光学卫星在轨实时几何定位方法及系统 | |
CN102168972B (zh) | 一种基于rpc的三线阵立体卫星区域网平差改进及检校方法 | |
CN104764443B (zh) | 一种光学遥感卫星严密成像几何模型构建方法 | |
CN103697864B (zh) | 一种基于大虚拟相机的窄视场双相机影像拼接方法 | |
CN106885585B (zh) | 一种基于光束法平差的星载摄影测量系统一体化检校方法 | |
CN110211054A (zh) | 一种星载推扫式光学传感器无畸变影像制作方法 | |
CN108088468A (zh) | 一种基于pos设备与数字航测相机的校验方法及装置 | |
CN102279001B (zh) | 一种星载相机的相移补偿方法 | |
CN111561936A (zh) | 旋转大幅宽光学卫星精准处理方法及系统 | |
CN111538051B (zh) | 一种摆扫大幅宽光学卫星精准处理方法 | |
Hattori et al. | Orientation of high-resolution satellite images based on affine projection | |
CN110986888A (zh) | 一种航空摄影一体化方法 | |
CN103778610B (zh) | 一种星载线阵传感器垂轨摆扫影像的几何预处理方法 | |
KR100520275B1 (ko) | Gps 측량에 의한 수치지도 제작시 촬영기구회전모델을 이용한 선형 스캐닝 영상의 기하 보정방법 | |
CN104976991B (zh) | 一种考虑卫星姿态变化的三线阵影像像方偏差的获取方法 | |
CN112816184A (zh) | 一种光学遥感卫星的无控定标方法及装置 | |
Li et al. | Photogrammetric processing of high-resolution airborne and satellite linear array stereo images for mapping applications | |
JP2011128784A (ja) | Rpc計算装置 | |
CN111667533A (zh) | 一种基于有理函数模型获取卫星严格成像模型参数的方法 | |
CN115854994A (zh) | 一种基于长线列大幅宽双向摆扫影像的dsm生成方法及系统 | |
Chen et al. | Rigorous georeferencing for Formosat-2 satellite images by least squares collocation | |
Chen et al. | Orbit adjustment for EROS A1 high resolution satellite images | |
CN115147313B (zh) | 椭圆轨道遥感图像的几何校正方法、装置、设备及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |