CN101907705B - 通用的多源遥感影像几何校正模型联合平差方法 - Google Patents
通用的多源遥感影像几何校正模型联合平差方法 Download PDFInfo
- Publication number
- CN101907705B CN101907705B CN 201010242711 CN201010242711A CN101907705B CN 101907705 B CN101907705 B CN 101907705B CN 201010242711 CN201010242711 CN 201010242711 CN 201010242711 A CN201010242711 A CN 201010242711A CN 101907705 B CN101907705 B CN 101907705B
- Authority
- CN
- China
- Prior art keywords
- centerdot
- partiald
- model
- reference mark
- unknown
- 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
Landscapes
- Image Processing (AREA)
Abstract
本发明从遥感影像的正算几何模型出发,建立了一种通用的多源遥感影像成像模型的区域网平差,并采用改进的Levenberg-Marquardt算法解算所建立的平差方程,克服了模型病态性对解的影响。该方法与遥感影像的具体成像模型无关,可以对不同的成像模型进行联合平差。在多源遥感数据的区域网平差过程中,可以针对各数据源的特点选择最合适的成像模型,充分发挥各种模型在严密性、定位准确性、保密性等方面的特点。
Description
技术领域
本发明建立了一种通用的多源遥感影像成像模型的区域网平差,该方法可以对不同的成像模型进行联合平差。在多源遥感数据的区域网平差过程中,可以针对各数据源的特点选择最合适的成像模型,充分发挥各种模型在严密性、定位准确性、保密性等方面的特点。
技术背景
卫星遥感影像的成像几何模型是描述地面点三维空间坐标与相应像点在像平面坐标系的二维坐标之间的映射关系的数学模型,它是进行影像几何校正和目标定位的基础。随着航天技术、传感器技术的发展,对于卫星遥感影像,已经出现了各种各样的成像模型,它们在严密性、复杂性及定位准确性方面都有着各自的特点。目前最常用的几何模型有严格成像模型、多项式模型、仿射变换模型和有理函数模型等,对于不同的卫星遥感数据源,严格成像模型又具有各自不同的表达形式。
对于各种成像模型,现在已有相应的区域网平差方法,但这些方法或者是针对具体的成像几何模型进行的,或者是基于有理函数模型进行的。针对具体成像模型的区域网平差方法缺乏通用性,不能用于多源遥感影像的联合平差;基于有理函数模型进行的区域网平差方法利用了有理函数模型的通用性,可以用于多源遥感影像的联合平差,但是对于一些严格成像模型已知的数据源(如Landsat系列、Spot系列等),有理函数模型并不是最佳的成像模型,这种联合平差的方法就不能发挥不同成像几何模型的优势。本文将针对遥感影像成像几何模型建立一种统一的区域网平差方法,实现不同成像模型的联合平差。
发明内容
本发明将不同的遥感影像成像几何模型描述为通用的正算模型,建立了一种统一的区域网平差方法,实现了不同成像模型的联合平差。
按照地面点三维空间坐标与相应像点二维平面坐标之间的映射顺序,卫星遥感影像的成像几何模型可以分为两类:反算模型和正算模型。正算模型表示为:
其中,(X,Y,Z)表示控制点的地面坐标,(x,y)表示控制点在影像上的量测坐标,t=(t1,t2,…tn)T表示传感器几何校正模型的参数。
不同卫星的严格成像模型、仿射变换模型、多项式模型和有理函数模型等均可表示成通用的正算模型,因此可以将多种成像模型联合起来进行区域网平差。模型线性化后,每个控制点的误差方程为:
其中,
将(1.2)表示成矩阵形式:
V=Adt+Bds-L (1-3)
其中,
dt=(dt1,dt2,…,dtn)T,ds=(dX,dY,dZ)T,V=(vx,vy)T,L=(lx,ly)T。
这里没有求解函数偏导的具体解析形式,而是采用数值计算的方法来逼近各函数偏导,如:
现有m景影像(可以为多源数据),n个控制点,其中nc个已知控制点,nu个未知控制点(即连接点),所有控制点所对应的像平面点个数为n′,设第l景影像的模型参数个数为pl。
单个控制点的线性化误差方程可以表示为:
Vi=Aidti+Bidsi-Li (1-5)
其中,
dsi=(dXi,dYi,dZi)T
1)若第i个像平面点在第j景影像上,对应于第k个未知控制点,则:
2)若第i个像平面点在第j景影像上,对应于一个已知控制点,则:
联立所有控制点的误差方程,多模型的区域网平差方程可以记为:
V=Adt+Bds-L (1-6)
其中,
多模型区域网平差方程(1-6)的解算:首先采用改化法方程的方法,消去多模型的区域网平差方程中所有连接点的地面坐标未知数,仅保留各传感器几何模型的未知参数,然后利用改进的LM(Levenberg-Marquardt)算法求解出各未知参数。
在传统的LM算法中,第k次迭代计算的未知数改变量ΔX为:
式中I是单位矩阵,J(k)是系数矩阵的Jacobian矩阵,μ(k)>0为阻尼系数。
对LM算法的改进如下:
其中
则(1-7)式改为:
附图说明
图1是使用本发明方法的具体实例处理过程示意图。
具体实施过程
下文结合说明书附图1,以Landsat-5和CBERS卫星遥感影像数据为例,对本发明的具体实施方式作详细说明。其中Landsat-5数据采用严格成像模型,CBERS数据采用有理函数模型。本发明所述方法包含但不限于所举实例。
步骤1:读取已经标准化后的卫星遥感影像数据;
步骤2:输入已知控制点和未知控制点(及连接点);
步骤3:利用卫星影像参数获取模型参数的初始值,并计算未知控制点地面坐标的初始值,当前迭代次数置为0;
步骤4:若当前迭代次数小于最大迭代次数,转至步骤5;若当前迭代次数达到最大迭代次数,转至步骤9;
步骤5:依次对输入的卫星遥感影像数据建立误差方程,模型的线性化过程中采用数值计算的方法逼近模型对各参数的偏导;
步骤6:将所有卫星遥感影像数据的误差方程联立为整体的区域网平差方程,采用改化法方程的方法,消去其中所有连接点的地面坐标未知数,仅保留各传感器几何模型的未知参数;
步骤7:使用改进的LM算法求出未知参数的改正值向量,并更新参数的值;
步骤8:用求得的改正值向量的模与规定的限差比较,若小于限差,转至步骤9;否则,当前迭代次数增加1,转至步骤4;
步骤9:解算完成,输出各影像模型参数。
Claims (1)
1.一种通用的多源遥感影像几何校正模型联合平差的方法,其主要步骤是:
步骤(1):读取已经标准化后的卫星遥感影像数据;
步骤(2):输入已知控制点、未知控制点及连接点;
步骤(3):利用卫星影像参数获取模型参数的初始值,并计算未知控制点地面坐标的初始值,当前迭代次数置为0;
步骤(4):若当前迭代次数小于最大迭代次数,转至步骤(5);若当前迭代次数达到最大迭代次数,转至步骤(9);
步骤(5):将不同的遥感影像成像几何模型描述为通用的正算模型:
其中,(X,Y,Z)表示控制点的地面坐标,(x,y)表示控制点在影像上的量测坐标,t=(t1,t2,…tn)T表示传感器几何校正模型的参数;
依次对输入的卫星遥感影像数据建立误差方程,模型的线性化过程中采用数值计算的方法逼近模型对各参数的偏导,模型线性化后,每个控制点的误差方程为:
其中,
步骤(6):将各误差方程联立成整体的区域网平差方程:
V=Adt+Bds-L (1-3)
式中,
m为进行区域网平差的影像个数;
n为控制点的个数;
nu为未知控制点的个数;
n′为所有控制点所对应的像平面点个数;
设第i个像平面点在第j景影像上,
其中pi为第i景影像的模型参数个数;
①若第i个像平面点对应于第k个未知控制点,则:
②若第i个像平面点对应于一个已知控制点,则:
采用改化法方程的方法,消去其中所有连接点的地面坐标未知数,仅保留各传感器几何模型的未知参数;
步骤(7):使用改进的LM算法求出未知参数的改正值向量,并更新参数的值;
其中对LM算法的改进如下:
在传统的LM算法中,第k次迭代计算的未知数改正量ΔX为:
式中I是单位矩阵,J(k)是系数矩阵的Jacobian矩阵,μ(k)>0为阻尼系数;
设 的对角线元素组成向量
其中
将(1-4)式改为:
则(1-5)式是改进的LM算法中第k次迭代计算的未知数改正量;
步骤(8):用求得的改正值向量的模与规定的限差比较,若小于限差,转至步骤(9);否则,当前迭代次数增加1,转至步骤(4);
步骤(9):解算完成,输出各影像模型参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010242711 CN101907705B (zh) | 2010-08-03 | 2010-08-03 | 通用的多源遥感影像几何校正模型联合平差方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010242711 CN101907705B (zh) | 2010-08-03 | 2010-08-03 | 通用的多源遥感影像几何校正模型联合平差方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101907705A CN101907705A (zh) | 2010-12-08 |
CN101907705B true CN101907705B (zh) | 2013-07-31 |
Family
ID=43263219
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010242711 Active CN101907705B (zh) | 2010-08-03 | 2010-08-03 | 通用的多源遥感影像几何校正模型联合平差方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101907705B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102322863B (zh) * | 2011-07-26 | 2013-07-24 | 武汉大学 | 一种遥感卫星多星联合逆向定轨定姿方法 |
CN103413272B (zh) * | 2013-07-22 | 2016-12-28 | 中国科学院遥感与数字地球研究所 | 低空间分辨率多源遥感图像空间一致性校正方法 |
CN103759714B (zh) * | 2014-01-20 | 2015-12-09 | 国家测绘地理信息局卫星测绘应用中心 | 一种三线阵卫星影像区域网平差方法 |
CN104123457B (zh) * | 2014-07-21 | 2017-03-08 | 湖北工业大学 | 一种稳健的卫星遥感影像有理函数模型参数估计方法 |
CN104764443B (zh) * | 2015-04-24 | 2017-07-18 | 国家测绘地理信息局卫星测绘应用中心 | 一种光学遥感卫星严密成像几何模型构建方法 |
CN112987024B (zh) * | 2019-08-23 | 2023-10-13 | 长沙天仪空间科技研究院有限公司 | 一种基于合成孔径激光雷达的成像装置及方法 |
CN114562982B (zh) * | 2022-03-09 | 2023-09-26 | 北京市遥感信息研究所 | 一种光学和sar异源卫星影像联合平差的定权方法和装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1834689A (zh) * | 2006-04-05 | 2006-09-20 | 中国科学院上海技术物理研究所 | 基于统计约束模型的线阵扫描图像非均匀性校正方法 |
EP1382980B1 (en) * | 2002-07-18 | 2010-02-24 | Itt Manufacturing Enterprises, Inc. | Active remote sensing using lock-in amplifiers and beam steering |
-
2010
- 2010-08-03 CN CN 201010242711 patent/CN101907705B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1382980B1 (en) * | 2002-07-18 | 2010-02-24 | Itt Manufacturing Enterprises, Inc. | Active remote sensing using lock-in amplifiers and beam steering |
CN1834689A (zh) * | 2006-04-05 | 2006-09-20 | 中国科学院上海技术物理研究所 | 基于统计约束模型的线阵扫描图像非均匀性校正方法 |
Non-Patent Citations (1)
Title |
---|
陈朋山等.抗差LM算法求解遥感影像严格物理模型.《科学技术与工程》.2009,第9卷(第16期),第4614-4618页. * |
Also Published As
Publication number | Publication date |
---|---|
CN101907705A (zh) | 2010-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101907705B (zh) | 通用的多源遥感影像几何校正模型联合平差方法 | |
CN102663767B (zh) | 视觉测量系统的相机参数标定优化方法 | |
CN103759714B (zh) | 一种三线阵卫星影像区域网平差方法 | |
CN101149836B (zh) | 一种三维重构的双摄像机标定方法 | |
CN101216555B (zh) | Rpc模型参数提取方法和几何纠正方法 | |
CN105866735B (zh) | 基于mds模型的修正代价函数的到达时间差迭代定位方法 | |
CN110969665B (zh) | 一种外参标定方法、装置、系统及机器人 | |
CN102778224B (zh) | 一种基于极坐标参数化的航空摄影测量光束法平差的方法 | |
CN107871327A (zh) | 基于特征点线的单目相机位姿估计和优化方法及系统 | |
CN106772524A (zh) | 一种基于秩滤波的农业机器人组合导航信息融合方法 | |
CN102609940A (zh) | 利用地面激光扫描技术进行测量对象表面重建时点云配准误差处理方法 | |
CN108919304B (zh) | 一种基于参考平面的移动测量系统中pos误差补偿方法 | |
CN108917772A (zh) | 基于序列图像的非合作目标相对导航运动估计方法 | |
Zheng et al. | Minimal solvers for 3d geometry from satellite imagery | |
CN104729481A (zh) | 一种基于pnp透视模型的合作目标位姿精度测量方法 | |
CN105046046A (zh) | 一种集合卡尔曼滤波局地化方法 | |
CN105353351A (zh) | 一种基于多信标到达时间差改进型定位方法 | |
CN107991676A (zh) | 星载单航过InSAR系统对流层误差校正方法 | |
CN106340059A (zh) | 一种基于多体感采集设备三维建模的自动拼接方法 | |
Ren et al. | High precision calibration algorithm for binocular stereo vision camera using deep reinforcement learning | |
CN104318566A (zh) | 可返回多个高程值的新型多视影像铅垂线轨迹匹配方法 | |
CN103513247B (zh) | 合成孔径雷达图像和光学图像同名点匹配的方法 | |
Wanchun et al. | Closed form algorithm of double-satellite TDOA+ AOA localization based on WGS-84 model | |
CN101718548A (zh) | 基于平面标志物的位姿处理方法 | |
CN106023146A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |