CN104102844A - 基于11范数约束的有理函数模型参数求解和全参数优化方法 - Google Patents
基于11范数约束的有理函数模型参数求解和全参数优化方法 Download PDFInfo
- Publication number
- CN104102844A CN104102844A CN201410352962.3A CN201410352962A CN104102844A CN 104102844 A CN104102844 A CN 104102844A CN 201410352962 A CN201410352962 A CN 201410352962A CN 104102844 A CN104102844 A CN 104102844A
- Authority
- CN
- China
- Prior art keywords
- rational function
- function model
- delta
- model
- reference mark
- 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.)
- Granted
Links
Landscapes
- Image Processing (AREA)
Abstract
基于l1范数约束的有理函数模型参数求解和全参数优化方法:在地形无关方案中,当卫星影像内部畸变较大时,仍可以得到可靠的有理函数模型参数;在地形相关方案中,可利用少于39个控制点求解78个有理函数模型系数;在已知有理函数模型系数的条件下,利用少量地面控制点对有理函数模型进行全参数优化。本发明提高了有理函数模型参数求解的稳定性,减少了所需控制点的个数,可以在少于39个控制点的条件下求出有理函数模型参数的可靠解,并提高了影像精校正的精度。
Description
技术领域
本发明涉及有理函数模型参数的求解和优化方法,能够自适应地精简有理函数模型,减少所需控制点的数量,提高图像校正精度。可应用于遥感、摄影测量、测绘、图像处理等领域。
背景技术
有理函数模型(Rational Function Model,RFM)是一种通用的几何成像模型,它具有良好的内插特性和连续性并独立于传感器和平台,此外还具有拟合精度高、成像参数及模型保密、计算速度快等优势,因而被广泛的应用于卫星遥感影像的几何校正中。
有理函数模型将像点坐标表示为以相应地面点空间坐标为自变量的有理多项式的比值,其标准方程为如式(1)所示:
其中l和s是像点在像平面上的行列值标准化后的结果[OGC(OpenGIS Consortium),1999.The OpenGISTM Abstract Specification-Topic7:The Earth Imagery Case.],(X,Y,Z)是物方点空间坐标标准化后的结果,
Nl(X,Y,Z)=a0+a1Z+a2Y+a3X+a4ZY+a5ZX+a6YX+a7Z2
+a8Y2+a9X2+a10ZYX+a11Z2Y+a12Z2X+a13ZY2
+a14Y2X+a15ZX2+a16YX2+a17Z3+a18Y3+a19X3
Dl(X,Y,Z)=b0+b1Z+b2Y+b3X+b4ZY+b5ZY+b6YX+b7Z2
+b8Y2+b9X2+b10ZYX+b11Z2Y+b12Z2X+b13ZY2
+b14Y2X+b15ZX2+b16YX2+b17Z3+b18Y3+b19X3
Ns(X,Y,Z)=c0+c1Z+c2Y+c3X+c4ZY+c5ZX+c6YX+c7Z2
+c8Y2+c9X2+c10ZYX+c11Z2Y+c12Z2X+c13ZY2
+c14Y2X+c15ZX2+c16YX2+c17Z3+c18Y3+c19X3
Ds(X,Y,Z)=d0+d1Z+d2Y+d3X+d4ZY+d5ZX+d6YX+d7Z2
+d8Y2+d9X2+d10ZYX+d11Z2Y+d12Z2X+d13ZY2
+d14Y2X+d15ZX2+d16YX2+d17Z3+d18Y3+d19X3
式中,ai,bi,ci,di(i=0,1,…,19)即为有理多项式系数(Rational Polynomial Coefficients,RPCs),其中b0和d0的值为1。
对影像坐标行列值的标准化公式为:
其中LINE_OFF、LINE_SCALE、SAMP_OFF和SAMP_SCALE为影像坐标的标准化因子。Line表示影像行坐标,Sample表示影像的列坐标。
对地面点的标准化公式为:
其中LONG_OFF、LONG_SCALE、LAT_OFF、LAT_SCALE、HEIGHT_OFF、HEIGHT_SCALE为地面坐标标准化因子。Longitude表示地面点的经度、Latitude表示地面点的纬度、Height表示地面点的高程。
有理函数模型的使用主要包括两大类:地形无关的方案和地形相关的方案:
1)地形无关的方案是通过卫星遥感影像的严格成像模型生成多个高程面的虚拟控制点格网,再利用这些虚拟控制点计算地形无关的有理函数模型参数(RPCs)。最后利用真实地面控制点对地形无关的有理函数模型进行像方或物方补偿,实现影像的精校正。例如像方多项式模型进行误差补偿的具体形式如(4)所示:
Δl=a0+al·l′+as·s′+als·l′·s′
+al2·l′2+as2·s′2+…
(4)
Δs=b0+bl·l′+bc·s′+bls·l′·s′
+bl2·l′2+bs2·s′2+…
其中,Δl和Δs是像平面坐标计算值与实际值之间的偏差,l′和s′是平面坐标的计算值,a0,al,as,…和b0,bl,bs,…是多项式模型的系数。较为常用的是像方仿射变换模型,最少需要3个控制点就能取得比较理想的改正效果。
2)地形相关方案则是直接利用大量的真实控制点解算有理函数模型参数,得到地形相关的有理函数模型,然后利用地形相关的有理函数模型进行影像的精校正。
其中地形无关的方案应用较为广泛,而地形相关方案的应用相对较少,主要原因是需要大量的控制点且模型的解算很不稳定。
近年来发射的和将来要发射的高分辨率卫星(如ALOS、GeoEye、Pleiades、SPOT-6/7、ZY-3、GF-1/2等)均提供地形无关的RPCs来描述其影像的成像模型,用户可以十分方便地利用这些参数对影像进行几何处理,在精确控制点和DEM的支持下,还能进行高精度的正射校正。
尽管有理函数模型具有很多的优势,但同时还存在一些缺点:1)模型包含78个RPCs,至少需要39个控制点才能求解;2)78个RPCs之间存在很强的相关性,这使得RPCs的求解成为一个病态性问题;3)有理函数模型是纯数学模型,缺乏几何意义;4)有理函数模型在窄视场角的高分辨率影像中应用很成功,但是对于宽视场角和定向参数不稳定的情况下适用性较差。
法方程的病态性是求解RPCs的主要困难,常用通过正则化(例如Tikhonov正则化、岭估计、Levenberg-Marquardt方法等)或奇异值截断的方法将法方程修正为良态,然后进行求解。这些方法可以在一定程度上改善方程的病态性、防止“过拟合”现象,但是仍然需要大量的控制点来解算RPCs。当控制点的数量不足够多时,得到的解往往很不稳定。此外,对于HJ-1A/B之类的国产卫星,由于传感器视场角大、平台姿态不稳定,图像的内部畸变比较大,在解算RPCs时更容易出现模型不稳定的现象,且影像的精确纠正很难通过简单的像方改正来实现,这是的有理函数模型无法应用于这类卫星遥感影像。
鉴于有理函数模型的通用性、简单性、保密性、计算量小等优势,克服现有方法中存在的问题,从而提高有理函数模型计算的稳定性及其适用范围具有重要的实用价值。
发明内容
本发明的目的在于解决现有技术的不足,提出了一种基于l1范数约束的有理函数模型求解方法和一种对有理函数模型全参数进行优化的精纠正方法。在地形无关方案中,当卫星影像内部畸变较大时,仍可以得到可靠的有理函数模型参数;在地形相关方案中,可利用少于39个控制点求解78个有理函数模型系数;在已知有理函数模型系数的条件下,利用少量地面控制点对有理函数模型进行全参数优化。本发明提高了有理函数模型参数求解的稳定性,减少了所需控制点的个数,可以在少于39个控制点的条件下求出有理函数模型参数的可靠解,并提高了影像精校正的精度。
本发明的技术方案如下:
一种基于l1范数约束的有理函数模型求解方法,所述方法利用控制点进行求解,控制点可以是真实地面控制点,也可以是通过严格成像模型生成的虚拟控制点;包括以下步骤,
步骤1.1,对所有控制点的像平面坐标和地面坐标进行标准化;
步骤1.2,通过简单的变形将有理函数模型转化为关于各RPCs的线性模型,并将控制点所列的误差方程用矩阵形式来表示;
步骤1.3,对线性模型模型加入l1范数约束,将问题转化为带l1范数约束的最优化问题;
步骤1.4,利用LARS(Least Angle Regression)算法[Efron,Bradley,Trevor Hastie,lain Johnstone,and Robert Tibshirani.″Least angle regression.″The Annals of statistics32,no.2(2004):407499.]求解带l1范数约束的最优化问题;
步骤1.5,得到有理函数模型参数RPCs,生成有理函数模型。
一种对有理函数模型全参数进行优化的精纠正方法,所述方法中已有的有理函数模型参数从卫星影像自带的RPC文件中获取或通过地形无关方案利用严格成像模型生成;需要利用真实地面控制点对有理函数模型参数进行优化;包括以下步骤,
步骤2.1,对所有真实地面控制点的像平面坐标和地面坐标进行标准化;
步骤2.2,通过简单的变形将有理函数模型转化为关于各RPCs的线性模型,然后控制点所列的误差方程可以用矩阵形式来表示;
步骤2.3,将有理函数模型参数的初始值带入误差方程,得到关于有理函数模型参数改正量的误差方程;
步骤2.4,对模型加入l1范数约束,将问题转化为带l1范数约束的最优化问题;
步骤2.5,利用LARS算法求解有理函数模型参数改正量的估计值;
步骤2.5,修正有理函数模型参数,生成改正后有理函数模型。
附图说明
图1是本发明实施例一流程图;
图2是本发明实施例二流程图;
图3是本发明实施例三流程图。
具体实施方式
本发明提出两种方法:一种是基于l1范数约束的有理函数模型求解方法,另一种是对有理函数模型全参数进行优化的精纠正方法。第一种方法可以分别应用于地形无关方案的有理函数参数求解和地形相关方案的有理函数模型参数求解,第二种方法可以用于已知初始有理函数模型参数前提下的有理函数模型精化。以下结合三个实施例和附图1~3详细说明本发明 技术方案。
(1)实施例一提供基于地形无关方案的有理函数模型参数求解方法,参见附图1:
步骤1.1,根据卫星遥感影像的严格成像模型生成虚拟控制点格网;
卫星遥感影像的严格成像模型是根据严格的物理成像过程建立的,用来描述像素坐标(行列号)与大地经纬度坐标(经度、纬度和高程)之间的映射关系,根据像素坐标和大地经纬度坐标的映射方向可分为正算模型(式(5))和反算模型(式(6))。
(X,Y)=F(l,s,Z) (5)
(l,s)=F-1(X,Y,Z) (6)
式(5)和式(6)中l表示像平面点的行坐标,s表示像平面点的列坐标,X表示经度,Y表示纬度,Z表示高程,F(·)表示正算模型,F-1(·)表示反算模型。
首先利用严格成像模型计算影像的覆盖范围,再根据全球DEM数据和影像覆盖范围统计影像覆盖范围内的最大和最小高程值,并在最大和最小高程之间建立7个高程面,在每个高程面按一定的间隔利用严格成像模型的正算模型(X,Y)=F(l,s,Z)生成21*21的均匀分布的控制点格网,这样就建立了7个不同高程面的三维控制点格网(n个控制点)。
步骤1.2,按式(2)和式(3)对n个控制点的像平面坐标和地面坐标进行标准化;
步骤1.3,通过变形将有理函数模型(1)转化为关于各RPCs的线性模型,如式(7)所示
然后n个控制点所列的误差方程可以矩阵形式来表示,如式(8)所示,
L=Ax+V (8)
其中
L=[l1,l2,…,ln,s1,s2,…,sn]T,
x=[a0,…,a19,b1,…b19,c0,…,c19,d1,…,d19]T,
V为随机误差,
A=[A1,A2,…,An]T,
(ll,sl)为第i个控制点标准化后的像平面行列值,
(Xi,Yi,Zi)为第i个控制点标准化后的地面坐标(经度、纬度和高程),
i=1,2,…,n;
步骤1.4,对模型(8)加入l1范数约束,将问题转化为带l1范数约束的最优化问题:
其中
α为l1范数约束因子,通常可取值为10-5;
||·||1和||·||2分别表示l1范数和l2范数,例如对于任意向量x:
||x||1=∑i|xl|表示向量x的l1范数,
表示向量x的l2范数;
步骤1.5,利用LARS(Least Angle Regression)算法[Efron,Bradley,Trevor Hastie,lain Johnstone,and Robert Tibshirani.″Least angle regression.″The Annals of statistics32,no.2(2004):407-499.]求解得到未知向量x的估计值
步骤1.6,根据估计值得到有理函数模型参数RPCs,生成有理函数模型。
(2)实施例二提供基于地形相关方案的有理函数模型参数求解方法,参见附图2:
步骤2.1,按式(2)和式(3)对n个真实控制点的像平面坐标和地面坐标进行标准化;
步骤2.2,通过变形将有理函数模型(1)转化为关于各RPCs的线性模型,如式(7)所示,然后将n个控制点所列的误差方程以矩阵形式来表示,如式(8)所示,该步骤与13相同;
步骤2.3,对模型(8)加入l1范数约束,将问题转化为带l1范数约束的最优化问题;
步骤2.4,利用LARS算法求解得到未知向量x的估计值
步骤2.5,根据估计值得到有理函数模型参数RPCs,生成有理函数模型。
(3)实施例三提供已知初始有理函数模型参数条件下的有理函数模型精化,参见附图3:
步骤3.1,按式(2)和式(3)对n个真实控制点的像平面坐标和地面坐标进行标准化;
步骤3.2,通过简单的变形将有理函数模型(1)转化为关于各RPCs的线性模型,如式(7)所示,然后n个控制点所列的误差方程可以矩阵形式来表示,如式(8)所示,该步骤与1.3相同;
步骤3.3,将向量x分解为两个分量的和
x=x0+Δx (10)
其中
然后模型(8)可转换为
L′=AΔx+V (11)
其中L′=L-Ax0,V为随机误差。
步骤3.4,对模型(11)加入l1范数约束,将问题转化为带l1范数约束的最优化问题:
其中α为l1范数约束因子,通常可取值为10-5;
步骤3.5,利用LARS(Least Angle Regression)算法求解得到改正向量Δx的估计值
步骤3.6,根据估计值修正有理函数模型参数,得到有理函数模型参数的估计值
然后生成改正后的有理函数模型。
Claims (2)
1.一种基于l1范数约束的有理函数模型求解方法,所述方法利用n个控制点进行求解,控制点可以是真实地面控制点,也可以是通过严格成像模型生成的虚拟控制点;其特征是:包括以下步骤,
步骤1.1,对n个控制点的像平面坐标和地面坐标进行标准化;
步骤1.2,通过变形将有理函数模型(1)转化为关于各RPCs的线性模型(2),
然后n个控制点所列的误差方程可以矩阵形式来表示,如式(3)所示,
L=Ax+V (3)
其中
L=[l1,l2,…,ln,s1,s2,…,sn]T,
x=[a0,…,a19,b1,…b19,c0,…,c19,d1,…,d19]T,
V为随机误差,
A=[A1,A2,…,An]T,
(li,si)为第i个控制点标准化后的像平面行坐标和列坐标,
(Xi,Yi,Zi)为第i个控制点标准化后的地面经度、纬度和高程,
i=1,2,…,n;
步骤1.3,对模型(3)加入l1范数约束,将问题转化为带l1范数约束的最优化问题:
其中α为l1范数约束因子;
步骤1.4,利用LARS算法求解得到未知向量x的估计值
步骤1.5,根据估计值得到有理函数模型参数RPCs,生成有理函数模型。
2.一种对有理函数模型全参数进行优化的精纠正方法,所述方法中已有的有理函数模型参数从卫星影像自带的RPC文件中获取或通过地形无关方案利用严格成像模型生成;需要利用n个真实地面控制点对有理函数模型参数进行优化;其特征为:包括以下步骤,
步骤2.1,对n个真实控制点的像平面坐标和地面坐标进行标准化;
步骤2.2,通过变形将有理函数模型(1)转化为关于各RPCs的线性模型(2),然后将n个控制点所列的误差方程以矩阵形式来表示,如式(3)所示,该步骤与1.2相同;
步骤2.3,将向量x分解为两个分量的和
x=x0+Δx (5)
其中
然后模型(3)可转换为
L′=AΔx+V (6)
其中L′=L-Ax0,V为随机误差。
步骤2.4,对模型(6)加入l1范数约束,将问题转化为带l1范数约束的最优化问题:
其中α为l1范数约束因子;
步骤2.5,利用LARS算法求解得到改正向量Δx的估计值
步骤2.6,根据估计值修正有理函数模型参数,得到有理函数模型参数的估计值
然后生成改正后的有理函数模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410352962.3A CN104102844B (zh) | 2014-07-24 | 2014-07-24 | 基于11范数约束的有理函数模型参数求解和全参数优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410352962.3A CN104102844B (zh) | 2014-07-24 | 2014-07-24 | 基于11范数约束的有理函数模型参数求解和全参数优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104102844A true CN104102844A (zh) | 2014-10-15 |
CN104102844B CN104102844B (zh) | 2017-11-17 |
Family
ID=51670989
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410352962.3A Active CN104102844B (zh) | 2014-07-24 | 2014-07-24 | 基于11范数约束的有理函数模型参数求解和全参数优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104102844B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104537614A (zh) * | 2014-12-03 | 2015-04-22 | 中国资源卫星应用中心 | 一种环境一号卫星ccd影像正射校正方法 |
CN104978761A (zh) * | 2015-06-01 | 2015-10-14 | 北京航空航天大学 | 一种rpc模型校正系数获取方法 |
CN105761228A (zh) * | 2016-03-09 | 2016-07-13 | 中国测绘科学研究院 | 通过微修正rpc参数实现卫星遥感影像高精度几何纠正方法 |
CN106777728A (zh) * | 2016-12-23 | 2017-05-31 | 南方电网科学研究院有限责任公司 | 一种复共轭有理函数对的实现电路及实现方法 |
CN109541593A (zh) * | 2018-10-30 | 2019-03-29 | 北京航空航天大学 | 一种改进的最小费用流InSAR相位解缠方法 |
CN115100079A (zh) * | 2022-08-24 | 2022-09-23 | 中国科学院空天信息创新研究院 | 遥感图像的几何校正方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101608914A (zh) * | 2009-07-23 | 2009-12-23 | 武汉大学 | 基于复共线性分析的rpc参数优选方法 |
CN101719118A (zh) * | 2009-11-24 | 2010-06-02 | 中国科学院对地观测与数字地球科学中心 | 一种用于克服有理函数模型病态性的改进的奇异值修正方法 |
CN101894270A (zh) * | 2010-07-26 | 2010-11-24 | 中国科学院遥感应用研究所 | 面向遥感影像分类的样本自动选取方法 |
CN103383773A (zh) * | 2013-03-26 | 2013-11-06 | 中国科学院遥感与数字地球研究所 | 一种动态提取图像控制点的遥感卫星图像自动正射纠正的框架和方法 |
CN103823981A (zh) * | 2014-02-28 | 2014-05-28 | 武汉大学 | 一种数字高程模型辅助的卫星影像区域网平差方法 |
-
2014
- 2014-07-24 CN CN201410352962.3A patent/CN104102844B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101608914A (zh) * | 2009-07-23 | 2009-12-23 | 武汉大学 | 基于复共线性分析的rpc参数优选方法 |
CN101719118A (zh) * | 2009-11-24 | 2010-06-02 | 中国科学院对地观测与数字地球科学中心 | 一种用于克服有理函数模型病态性的改进的奇异值修正方法 |
CN101894270A (zh) * | 2010-07-26 | 2010-11-24 | 中国科学院遥感应用研究所 | 面向遥感影像分类的样本自动选取方法 |
CN103383773A (zh) * | 2013-03-26 | 2013-11-06 | 中国科学院遥感与数字地球研究所 | 一种动态提取图像控制点的遥感卫星图像自动正射纠正的框架和方法 |
CN103823981A (zh) * | 2014-02-28 | 2014-05-28 | 武汉大学 | 一种数字高程模型辅助的卫星影像区域网平差方法 |
Non-Patent Citations (1)
Title |
---|
LU ZHANG等: "Rational function modeling for spaceborne SAR datasets", 《ISPRS JOURNAL OF PHOTOGRAMMETRY AND REMOTE SENSING》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104537614A (zh) * | 2014-12-03 | 2015-04-22 | 中国资源卫星应用中心 | 一种环境一号卫星ccd影像正射校正方法 |
CN104978761A (zh) * | 2015-06-01 | 2015-10-14 | 北京航空航天大学 | 一种rpc模型校正系数获取方法 |
CN104978761B (zh) * | 2015-06-01 | 2018-08-24 | 北京航空航天大学 | 一种rpc模型校正系数获取方法 |
CN105761228A (zh) * | 2016-03-09 | 2016-07-13 | 中国测绘科学研究院 | 通过微修正rpc参数实现卫星遥感影像高精度几何纠正方法 |
CN106777728A (zh) * | 2016-12-23 | 2017-05-31 | 南方电网科学研究院有限责任公司 | 一种复共轭有理函数对的实现电路及实现方法 |
CN109541593A (zh) * | 2018-10-30 | 2019-03-29 | 北京航空航天大学 | 一种改进的最小费用流InSAR相位解缠方法 |
CN115100079A (zh) * | 2022-08-24 | 2022-09-23 | 中国科学院空天信息创新研究院 | 遥感图像的几何校正方法 |
CN115100079B (zh) * | 2022-08-24 | 2022-11-22 | 中国科学院空天信息创新研究院 | 遥感图像的几何校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104102844B (zh) | 2017-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104102844A (zh) | 基于11范数约束的有理函数模型参数求解和全参数优化方法 | |
CN107689064B (zh) | 顾及光行差改正的卫星光学严格几何成像模型构建方法 | |
CN102168972B (zh) | 一种基于rpc的三线阵立体卫星区域网平差改进及检校方法 | |
CN105510913A (zh) | 基于类光学像方改正的异源光学和sar遥感影像联合定位方法 | |
Shen et al. | Correcting bias in the rational polynomial coefficients of satellite imagery using thin-plate smoothing splines | |
CN103759714A (zh) | 一种三线阵卫星影像区域网平差方法 | |
CN108007476B (zh) | 一种天基干涉成像雷达高度计的干涉定标方法及系统 | |
CN107316280B (zh) | 离岛卫星影像rpc模型高精度几何定位方法 | |
CN103632349A (zh) | 一种降低tdi-ccd相机图像模糊度的方法 | |
CN101598785B (zh) | 卫星影像各级产品的有理函数成像模型生成方法 | |
CN104123457A (zh) | 一种稳健的卫星遥感影像有理函数模型参数估计方法 | |
CN111650579A (zh) | 一种岩移参数自适应获取的InSAR矿区三维形变估计方法、装置及介质 | |
CN111724465A (zh) | 基于平面约束优选虚拟控制点的卫星影像平差方法及装置 | |
CN103778612A (zh) | 一种基于全色影像的卫星颤振探测与补偿方法 | |
CN109188483B (zh) | 一种时序化高精度外方位元素自动定标方法 | |
Qiu et al. | A high-precision calibration approach for Camera-IMU pose parameters with adaptive constraints of multiple error equations | |
CN101598786B (zh) | 卫星影像二/三级产品的有理函数成像模型生成方法 | |
CN106157258B (zh) | 一种星载sar图像几何校正方法 | |
Li et al. | Improve the ZY-3 height accuracy using ICESat/GLAS laser altimeter data | |
CN110046430B (zh) | 基于椭球分层大气折射的光学卫星影像精确对地定位方法 | |
CN105678716A (zh) | 一种地基sar大气干扰相位校正方法及装置 | |
CN115097450A (zh) | 跨轨道高分三号sar偏移量大梯度滑坡形变估计方法 | |
CN111738135B (zh) | 一种顾及斜距和地理坐标映射偏差的sar影像特征提取方法 | |
CN108957500B (zh) | 一种传感器观测视线与地球表面交点的计算方法 | |
Shen et al. | An improved method for transforming GPS/INS attitude to national map projection frame |
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 |