CN102221690B - 高分辨率卫星遥感定位的组合岭估计光束法平差方法 - Google Patents
高分辨率卫星遥感定位的组合岭估计光束法平差方法 Download PDFInfo
- Publication number
- CN102221690B CN102221690B CN 201010148266 CN201010148266A CN102221690B CN 102221690 B CN102221690 B CN 102221690B CN 201010148266 CN201010148266 CN 201010148266 CN 201010148266 A CN201010148266 A CN 201010148266A CN 102221690 B CN102221690 B CN 102221690B
- Authority
- CN
- China
- Prior art keywords
- ridge
- estimated
- parameter
- mountain range
- prime
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Image Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
Description
技术领域
本方法涉及一种高分辨率卫星遥感定位的光束法平差法,尤其是涉及一种高分辨率卫星遥感定位的组合岭估计光束法平差方法。
背景技术
在空中三角测量中,光束法平差在理论上是最严密的方法,因为这种方法直接使用了像点坐标观测值,但是实验结果表明光束法平差并不是精度最高的,主要是由于还存在残余的系统误差,解决的办法主要是自检校法。李德仁等学者将光束法平差应用于高分辨率卫星影像定位中,但是由于法方程的病态性导致计算的结果不甚理想,因为未知外方位线角元素强相关,现有的解决方法主要有四种:增加虚拟观测方程、合并强相关项、线角分求法、岭估计。在处理高分辨率遥感影像定位法方程病态问题时岭估计特别是广义岭估计近来受到越来越多的重视。但是广义岭估计首先需要获取最小二乘解来解算岭参数,而当法方程系数阵严重病态时,可能无法获得最小二乘解,因为法方程系数阵无法求逆。本方法针对该问题结合文献的方法提出一种组合岭估计方法进行单线阵高分辨率影像的光束法平差,取得了良好的效果,可达0.5m的平面定位精度和0.7m的高程定位精度。
方法内容
本方法的目的就是为了克服上述现有技术存在的缺陷而提供的一种能提高定位精度的高分辨率卫星遥感定位的组合岭估计光束法平差方法
本方法的目的可以通过以下技术方案来实现:
一种高分辨率卫星遥感定位的组合岭估计光束法平差方法,其特征在于,包括以下步骤:
1)读取像对数据、控制数据、传感数据,并设定限差与迭代次数;
2)数据预处理;
3)通过L曲线法计算岭参数;
7)打印计算结果并进行精度分析。
所述的步骤2)数据预处理包括空间后方交会、空间前方交会。
所述的步骤3)中的岭参数计算过程如下:
通过L曲线法计算岭参数K值,即选择L曲线上曲率最大的点相对应的K值作为所求的参数,根据岭估计的估计准则可得和都是岭参数K的函数,选择不同的K值并以为横坐标、为纵坐标画图,得到L曲线, 则L曲线上的点的曲率q的计算公式为:
对上式求最大值,就可以得到最大曲率所对应的K值即为岭参数。
与现有技术相比,本方法具有以下优点:
1、解决了由于法方程系数严重病态而无法求得最小二乘解,从而不能进行广义岭估计的问题;
2、提高了定位精度,用本方法进行单线阵高分辨率影像的光束法平差,取得了良好的效果,可达0.5m的平面定位精度和0.7m的高程定位精度。
附图说明
图1为本发明的流程图;
图2为前方交会误差分布图(放大20倍);
图3为基于岭估计的光束法平差误差分布图(放大500倍);
图4为基于组合岭估计的光束法平差误差分布图(放大600倍)。
具体实施方式
下面结合附图和具体实施例对本方法进行详细说明。
如图1所示,一种高分辨率卫星遥感定位的组合岭估计光束法平差方法,包括以下步骤:
1)读取像对数据、控制数据、传感数据,并设定限差与迭代次数;
2)对数据进行空间后方交会、空间前方交会处理;
3)通过L曲线法计算岭参数;
7)打印计算结果并进行精度分析。
单线阵光束法平差模型
现有大部分商用高分辨率卫星影像(QuickBird、IKONOS、SPOT等)都是单线阵沿飞行方向推扫成像。与航片不同的是,卫星影像每一行的外方位元素都不一样,不过每一扫描行的影像与被摄物体之间仍然满足共线方程。
假设飞行方向为像坐标x轴方向,扫描行方向为像坐标y轴方向,若第i扫描行的外方位元素为Xsi,Ysi,Zsi,ωsi,κsi,(X,Y,Z)为地面点的物方空间坐标;(Xsi,Ysi,Zsi)为第i行投影中心的物方空间坐标。其共线方程的形式:
设外方位角元素分别为ωi,κi,由于卫星飞行相对稳定,可以认为外方位元素为时间t的函数(一般认为取2次足够保证精度),故并可用下述数学模型进行描述:
将(1)式按泰勒公式展开至一次项,得误差方程:
令:
则(2)式可表示为:
Vi=AiXi-Li (3)
式中a11,a12,…,a26的具体表达式为:
其中:
X=(ATPA)-1ATPL
基于组合岭估计的光束法平差
由于线角元素高度相关,法方程系数ATPA严重病态,故采用有偏估计,有偏估计中比较实用的是广义岭估计。
广义岭估计
解算广义岭估计值的关键是确定岭参数K,其基本原则是寻找使均方差达到最小的对角阵K。而岭估计中K是个正常数。
组合岭估计
由于广义岭估计求解K需要首先算出最小二乘解,但是当法方程严重病态时往往最小二乘解无解,本文提出应用L曲线法先计算出岭估计参数K值,即选择L曲线上曲率最大的点相对应的K值作为所求的参数。根据岭估计的估计准则知道和都是岭参数K的函数,选择不同的K值并以为横坐标、为纵坐标画图,得到L曲线。令
则L曲线上的点的曲率q的计算公式为:
对上式求最大值,就可以得到最大曲率所对应的K值即为所求值。
实施例
本方法采用QuickBird卫星的某地区异轨立体像对,为Basic产品。影像分辨率为0.7m,标称地面定位精度为23m,覆盖范围大约为18km×18km。控制点坐标采用基于某市虚拟参考网络VRS(Virtual reference station)的GPSRTK((real time kinematic)测得,精度在5cm以内,共有69个测点。所有测点的高度变化在100m以内。从中选择任选23个已知点对立体定位结果进行检验。采用VC+编制程序进行光束法平差解算得到的结果如表1所示:
表1各种计算方法的误差比较
为了更清楚的展示误差的分布,利用matlab绘制了误差图,前方交会误差分布图如图2,基于岭估计的光束法平差误差分布图如图3,而基于组合岭估计的光束法平差误差分布图如图4,图2清楚的表明前方交会解算的结果具有很强的系统误差,而基于岭估计和组合岭估计的光束法平差计算出的结果则表现出杂乱无章的误差分布,亦即没有明显的系统误差,并且本文提出的基于组合岭估计方法能获取更高精度的加密点坐标。另外控制点的分布对加密点的坐标也有很大的影响,加密点周围控制点多的比周围控制点少的精度高。
Claims (5)
2.根据权利要求1所述的一种高分辨率卫星遥感定位的组合岭估计光束法平差方法,其特征在于,所述的步骤2)数据预处理包括空间后方交会和空间前方交会。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010148266 CN102221690B (zh) | 2010-04-15 | 2010-04-15 | 高分辨率卫星遥感定位的组合岭估计光束法平差方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010148266 CN102221690B (zh) | 2010-04-15 | 2010-04-15 | 高分辨率卫星遥感定位的组合岭估计光束法平差方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102221690A CN102221690A (zh) | 2011-10-19 |
CN102221690B true CN102221690B (zh) | 2013-02-13 |
Family
ID=44778292
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010148266 Expired - Fee Related CN102221690B (zh) | 2010-04-15 | 2010-04-15 | 高分辨率卫星遥感定位的组合岭估计光束法平差方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102221690B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104123457B (zh) * | 2014-07-21 | 2017-03-08 | 湖北工业大学 | 一种稳健的卫星遥感影像有理函数模型参数估计方法 |
CN104978761B (zh) * | 2015-06-01 | 2018-08-24 | 北京航空航天大学 | 一种rpc模型校正系数获取方法 |
CN107490800B (zh) * | 2017-08-07 | 2019-09-03 | 桂林电子科技大学 | 一种卫星导航快速定位方法、装置和卫星导航接收机 |
CN110389365B (zh) * | 2018-04-23 | 2021-09-10 | 中移物联网有限公司 | 一种卫星导航定位方法及装置、终端、存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101226057A (zh) * | 2008-02-01 | 2008-07-23 | 武汉朗视软件有限公司 | 一种数字近景摄影测量方法 |
CN101604018A (zh) * | 2009-07-24 | 2009-12-16 | 中国测绘科学研究院 | 高分辨率遥感影像数据处理方法及其系统 |
CN101644595A (zh) * | 2009-09-01 | 2010-02-10 | 南京大学 | 一种复杂水位过程的拟合方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3196683B2 (ja) * | 1997-03-18 | 2001-08-06 | トヨタ自動車株式会社 | はんだ付部の検査方法及びその装置 |
-
2010
- 2010-04-15 CN CN 201010148266 patent/CN102221690B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101226057A (zh) * | 2008-02-01 | 2008-07-23 | 武汉朗视软件有限公司 | 一种数字近景摄影测量方法 |
CN101604018A (zh) * | 2009-07-24 | 2009-12-16 | 中国测绘科学研究院 | 高分辨率遥感影像数据处理方法及其系统 |
CN101644595A (zh) * | 2009-09-01 | 2010-02-10 | 南京大学 | 一种复杂水位过程的拟合方法 |
Non-Patent Citations (4)
Title |
---|
JP特开平10-261084A 1998.09.29 |
刘世杰.高分辨率卫星遥感影像成像模型与定位技术研究.《中国优秀硕士学位论文全文数据库》.2008,全文. * |
岳庆兴等.SPOT影像光束法平差复共线性消除.《海洋测绘》.2008,第28卷(第02期),全文. * |
黄幼才.岭估计及其应用.《武汉大学学报(信息科学版)》.1987,第12卷(第04期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN102221690A (zh) | 2011-10-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11988781B2 (en) | Extrinsic calibration method of multiple 3D LiDAR sensors for autonomous navigation system | |
Nuth et al. | Co-registration and bias corrections of satellite elevation data sets for quantifying glacier thickness change | |
CN109917356B (zh) | 一种机载激光扫描系统误差标定方法 | |
CN107504981B (zh) | 一种基于激光测高数据的卫星姿态误差修正方法及设备 | |
CN106885576B (zh) | 一种基于多点地形匹配定位的auv航迹偏差估计方法 | |
CN102353378B (zh) | 一种矢量形式信息分配系数的组合导航系统自适应联邦滤波方法 | |
CN103913181A (zh) | 一种基于参数辨识的机载分布式pos传递对准方法 | |
CN102636159B (zh) | 多镜头航天线阵相机系统在轨几何自检校方法 | |
CN102221690B (zh) | 高分辨率卫星遥感定位的组合岭估计光束法平差方法 | |
CN104931022A (zh) | 基于星载激光测高数据的卫星影像立体区域网平差方法 | |
CN105224737A (zh) | 一种空间目标轨道改进初值修正方法 | |
CN102208108B (zh) | 摄像机大视场高精度快速现场全局标定方法 | |
CN113466890B (zh) | 基于关键特征提取的轻量化激光雷达惯性组合定位方法和系统 | |
CN103927744B (zh) | 一种基于指向姿态的敏捷卫星观测目标条带分割方法 | |
CN105225240A (zh) | 一种基于视觉特征匹配与拍摄角度估计的室内定位方法 | |
CN102706352A (zh) | 一种航空用线状目标的矢量地图匹配导航方法 | |
CN107991676A (zh) | 星载单航过InSAR系统对流层误差校正方法 | |
Giles et al. | A method for sub-pixel scale feature-tracking using Radarsat images applied to the Mertz Glacier Tongue, East Antarctica | |
CN104978761A (zh) | 一种rpc模型校正系数获取方法 | |
CN102147249B (zh) | 基于直线特征的星载光学线阵影像精确纠正处理方法 | |
CN103575298A (zh) | 基于自调节的ukf失准角初始对准方法 | |
CN105004321A (zh) | 顾及曝光不同步的无人机gps辅助光束法平差方法 | |
CN105389819A (zh) | 一种鲁棒的半标定下视图像极线校正方法及系统 | |
CN105353359A (zh) | 异步组网雷达系统误差校正方法 | |
Dachauer et al. | Aerodynamic roughness length of crevassed tidewater glaciers from UAV mapping |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130213 Termination date: 20160415 |