CN107314763A - 一种基于制约函数非线性估计的卫星影像区域网平差方法 - Google Patents
一种基于制约函数非线性估计的卫星影像区域网平差方法 Download PDFInfo
- Publication number
- CN107314763A CN107314763A CN201710583487.4A CN201710583487A CN107314763A CN 107314763 A CN107314763 A CN 107314763A CN 201710583487 A CN201710583487 A CN 201710583487A CN 107314763 A CN107314763 A CN 107314763A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- msup
- mtr
- mtd
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
- G01C11/04—Interpretation of pictures
- G01C11/30—Interpretation of pictures by triangulation
- G01C11/34—Aerial triangulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
- G06T7/73—Determining position or orientation of objects or cameras using feature-based methods
- G06T7/74—Determining position or orientation of objects or cameras using feature-based methods involving reference images or patches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10016—Video; Image sequence
- G06T2207/10021—Stereoscopic video; Stereoscopic image sequence
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
Abstract
本发明涉及一种基于制约函数非线性估计的卫星影像区域网平差方法,包括以下步骤,S01:构建基于RFM的像方空间影像定位的函数模型;S02:引入约束条件;S03:对待估参数进行求解。其优点在于,在常规区域网平差数学模型基础上引入约束条件,避免无控制点情况下由于必要起算数据缺失引起的解不唯一问题,同时采用制约函数法对未知参数进行最优估计,保留函数模型的原始非线性形式,将带约束条件的非线性问题转化为求解无约束函数极值的问题,避免采用线性最优估计算法带来的模型误差问题,从而保证无控制点情况下卫星影像区域网平差解的高精度。
Description
技术领域
本发明涉及卫星摄影测量技术领域,具体地说,是一种基于制约函数非线性估计的卫星影像区域网平差方法。
背景技术
航天技术、计算机技术、网络技术、信息处理技术的进步推动遥感对地观测技术不断向前发展,高分辨率测绘卫星系统不断涌现,利用高分辨率遥感卫星进行航天摄影测量具有效率高、覆盖范围大、不受空间管制政策限制的优势,因此高分辨率卫星影像成为获取中小比例尺地理空间信息的重要手段之一。影像的高精度几何定位是高分辨率卫星影像广泛应用的重要前提,是利用卫星影像生产4D数字产品(DTM、DOM、DLG、DRG)的基础。目前“立体影像--地面控制点--空三”是摄影测量非常成熟的几何定位处理过程,利用具有一定视差的立体影像对,匹配出足够数量同名连接点,以构像模型为基础,补偿系统误差,对未知数进行最优估计实现高精度对地定位目的,然而这种方法依赖地面控制点,在境外,岛屿、沙漠等控制点难以获取区域受到限制,而且环保、海洋、矿产等领域对影像内部拼接精度要求高于对绝对精度的要求,因此,影像无控制点区域网平差技术的研究具有重要意义,不仅可以有效解决无控区域测绘问题,还能提高影像间内部符合精度,提高影像间拼接精度。目前卫星应用商在元数据中提供有理多项式系数(Rational Polynomial Coefficient,RPC)文件,可以用来直接计算地面点三维坐标,然而定轨测姿误差会延续到RPC文件中,导致直接解算精度降低,而且通常采用高斯-牛顿法将非线性问题转换为线性问题进行参数最优估计,RFM模型分子分母皆为三次项,强制线性化会带来模型误差,影响解的精度。
中国发明专利CN201010588994.5,公开日为2011.08.31,公开了一种基于RPC的三线阵立体卫星区域网平差改进及验校方法,包括以下步骤:(1)由卫星下传辅助数据提取卫星的轨道和姿态,并根据提取的卫星的轨道和姿态确定地面点与像点之间的投影关系,即严格成像模型;(2)根据严格成像模型,求解出三线阵各个影像的RPC参数,通过影像匹配找出三线阵之间的连接点;(3)列出三线阵影像的前视、正视、后视三个仿射变换公式;(4)根据步骤(3)列出的三个仿射变换公式建立误差方程,将步骤(2)找出的连接点坐标以及地面控制点库中的控制点坐标代入误差方程进行解算,获得仿射变换参数的改正数,利用该改正数修正步骤(3)的仿射变换公式;(5)利用相机参数和修正后的正视仿射变换公式,再次列写前视和后视仿射变换公式;(6)将步骤(4)修正后的前视、后视仿射变换公式与步骤(5)列写前视和后视仿射变换公式进行联立,获得相机参数的改正数,修正相机参数,完成三线阵立体卫星区域网平差改进及检校。但是这种方法依赖控制点,无法适用于无控区域测绘问题。
中国发明专利CN201110091756.8,公开日为2011.10.12,公开了一种基于RFM模型的多源星载SAR影像自动匹配方法,首先计算影像各自的RPC参数,然后,在每个金字塔影像层,利用影像RPC参数进行待匹配点初始点位的预测、近似核线几何建立约束匹配、匹配窗口影像几何粗纠正,并对每层金字塔影像匹配结果采用基于RFM模型区域网平差删除错误的匹配点,同时实现影像RPC参数的精化和匹配点物方坐标的计算。逐层精化匹配结果直到原始影像层,最后利用最小二乘影像匹配方法精化匹配结果,实现多源星载SAR影像同名点的自动可靠匹配。本发明将RFM模型引入到多源星载SAR影像自动匹配,并将RFM模型区域网平差融入到各层金字塔影像匹配过程中,可以有效删除匹配中错误的匹配点,有效降低同名点人工量测的工作量。这种方法采用高斯-牛顿法将非线性问题转换为线性问题进行参数最优估计,RFM模型分子分母皆为三次项,强制线性化会带来模型误差,影响解的精度。
中国发明专利CN201410025159.9,公开日为2014.04.30,公开了一种三线阵卫星影像区域网平差方法,包括以下步骤:1)读取卫星影像,并将各类点的坐标信息存成相应的点文件;2)读取各文件,得到卫星影像的RPC模型以及控制点和检查点的地面坐标、影像坐标和连接点的像点坐标;3)确定地面点坐标和像方变换模型的初值;4)逐点构建误差方程并进行法化;5)对法方程利用谱修正法迭代估计求解改正参数和地面点坐标;6)计算此次平差迭代后所能达到的物方精度;7)计算此次平差迭代后所能达到的像方精度,通过像方残差并利用信息扩散模型确定下一次迭代像点坐标观测值的权值P;8)当满足定向参数中的平移参数均小于阈值时时,平差迭代结束;当不满足预设条件时继续迭代,直至迭代结束;9)当整个平差迭代结束时,输出最平差物方精度和像方精度的精度报告。该方法利用谱修正法迭代估计求解改正参数和地面点坐标,设定迭代次数以及阈值,但是这种方法无法控制精度,只能趋近于相应点,需要与地面点进行匹配,当被测区域无控制点时,该方法就失去作用。
因此,亟需一种能够有效降低误差、提高精度、适用于无控区域,对地面控制点依赖较低的卫星影像区域网平差方法,而目前关于这种方法还未见报道。
发明内容
本发明的目的是针对现有技术中的不足,提供一种基于制约函数非线性估计的卫星影像区域网平差方法。
为实现上述目的,本发明采取的技术方案是:
一种基于制约函数非线性估计的卫星影像区域网平差方法,所述基于制约函数非线性估计的卫星影像区域网平差方法包括以下步骤:
S01:构建基于RFM的像方空间影像定位的函数模型;
S011:获取立体影像及RPC文件;
S012:匹配获取足够数量的同名连接点;
S013:构建RFM模型:
利用RPC文件构建描述物方点与对应像点之间定量关系的RFM模型;
S014:构建系统误差补偿模型:
对每一个像点进行系统误差补偿,补偿模型如(1)所示,
其中,Δp、Δr表示像方系统误差,(1,s)表示像点坐标,a0,as,a1,b0,bs,b1表示每景影像对应的6个仿射变换系数,作为待求定向参数,其中a0,as,a1修正影像行方向误差,b0,bs,b1修正影像列方向误差;
S015:建立区域网平差数学模型:
将Δp、Δr代入RFM模型,即可得到像方空间影像区域网平差的函数模型,如(2)所示:
其中,(P,L,H)表示地面点大地坐标经平移和缩放后的标准化坐标,取值位于(-1.0~+1.0)之间,NumL、DenL、NumS、DenS表示有理多项式,LINE_OFF、SAMPLE_OFF为像点标准化的平移参数,LINE_SCALE、SAMPLE_SCALE为像点标准化的比例参数;
S016:构建目标函数:
根据观测值与未知数之间的函数关系(式(2))构建目标函数,将像点理论值与实测值间的残差平方和最小作为目标函数,如(3)所示:
f(X)=min(VTPV) (3)
其中,f(X)为目标函数,X为所有待估参数,包括连接点三维坐标及每景影像的6个定向参数,V表示像点理论值与实测值间的残差向量,P表示观测值的权;
S02:引入约束条件:
约束条件为平差后网中重心点的位置保持不变,平差后网中各连接点至重心点的矢径角变化的带权总和为0、平差后网中各点至重心点的矢径距变化的带权总和为0,如(4)所示;
STx=O (4)
其中,x表示连接点坐标的改正数;
S03:对待估参数进行求解;
S031:构建罚函数:
联合目标函数(式(3))与约束条件(式(4)),构建罚函数,如(5)所示:
其中,P(X,M)为罚函数,M为罚因子;
S032:选取初始罚因子;
S033:计算对应的待估值:
利用罚函数对未知参数求一阶导数,并令一阶导数为0,得到带有M的极值的罚函数的解,将初始罚因子代入,得到解;
S034:判断待估值与约束集的距离是否小于限差:
若得到的解不满足约束条件,变更M的值,重复步骤S033;若满足约束条件,则进入步骤S035;
S035:获得影像定向参数及加密点坐标:
满足约束条件的罚函数如(6)所示:
minP(X,Mk)=P(X(k),Mk) (6)
其中,Mk为罚函数最小时的罚因子,X(k)为对应的最优待估参数值;
获得的解,即X(k),为区域网平差解算得到最优连接点三维坐标及影像定向参数。
作为一种优选的技术方案,所述的步骤S012中,利用最小二乘影像匹配算法对立体影像序列进行匹配,获取足够数量的同名连接点,并建立连接点文件,所述的连接点文件记录的信息包括每个连接点所在影像的编号、像点坐标等。
作为一种优选的技术方案,所述的步骤S013中,在RFM模型中,像点坐标(1,s)表示为以相应地面点空间坐标(φ,λ,h)为自变量的多项式的比值,如(7)所示:
其中,(ln,sn)和(P,L,H)分别表示像点坐标(1,s)和地面点坐标(φ,λ,h)经平移和缩放后的标准化坐标,取值位于(-1.0~+1.0)之间,其变换关系如(8)~(16)所示:
其中,LAT_OFF、LONG_OFF、HEIGHT_OFF、LINE_OFF、SAMP_OFF为标准化的平移参数,LAT_SCALE、LONG_SCALE、HEIGHT_SCALE、LINE_SCALE、SAMP_SCALE为标准化的比例参数,ai,bi,ci,di(i=0,1,2,…,19)为RPC文件。
作为一种优选的技术方案,所述的步骤S033中,利用罚函数对未知参数求一阶导数,并令一阶导数为0,如(17)所示:
其中,X1,X2,X3,...为未知参数。
作为一种优选的技术方案,所述的步骤S034中,若得到的解不满足约束条件,则增大M的值,重复步骤S033。
本发明优点在于:
1、将平差后网中重心点的位置保持不变、网中各连接点至重心点的矢径角变化的带权总和为0、网中各点至重心点的矢径距变化的带权总和为0作为约束条件,引入到常规基于RFM模型的区域网平差数学模型中,避免由于必要起算数据缺失引起的解不唯一问题,使得无控制点条件下卫星影像区域网平差解算成为可能,流程易于设计与实现;
2、采用制约函数法对未知参数进行最优估计,保留函数模型的原始非线性形式,将带约束条件的非线性问题转化为求解无约束函数极值的问题,降低了高斯-牛顿法强制线性化引入的模型误差,从而保证无控制点情况下卫星影像区域网平差解的高精度;
3、适用于境外、岛屿、沙漠等控制点难以获取的区域,满足环保、海洋、矿产等领域对影像内部拼接精度高于绝对精度的要求。
附图说明
附图1是本发明的一种基于制约函数非线性估计的卫星影像区域网平差方法的总架构图。
附图2是本发明的一种基于制约函数非线性估计的卫星影像区域网平差方法的详细流程图。
附图3是本发明的构建基于RFM的像方空间影像定位的函数模型的流程图。
附图4是本发明的对待估参数进行求解的流程图。
具体实施方式
下面结合附图对本发明提供的具体实施方式作详细说明。
实施例1
参照图1,图1是本发明的一种基于制约函数非线性估计的卫星影像区域网平差方法的总架构图。
本发明的一种基于制约函数非线性估计的卫星影像区域网平差方法,包括以下步骤:
S01:构建基于RFM的像方空间影像定位的函数模型;
获取影像元数据中提供有理多项式系数文件(Rational PolynomialCoefficient),即RPC文件,用于确定物方点与对应像点之间的定量关系;对影像做形如仿射变换的几何纠正补偿系统误差,将系统误差表示为像点坐标的仿射变换,附加在有理函数模型中进行区域网平差,实现像方空间影像定位。
像点坐标系统误差的公式如(1)所示:
其中,Δp、Δr表示像方系统误差,(l,s)表示像点坐标,a0,as,a1,b0,bs,b1表示每景影像对应的6个仿射变换系数,作为待求定向参数,其中a0,as,a1修正影像行方向误差,b0,bs,b1修正影像列方向误差。
将Δp、Δr代入有理函数模型(Rational Function Model,RFM)即可得到像方空间影像区域网平差的函数模型,如(2)所示:
其中,(P,L,H)表示地面点大地坐标经平移和缩放后的标准化坐标,取值位于(-1.0~+1.0)之间,NumL、DenL、NumS、DenS表示有理多项式,LINE_OFF、SAMPLE_OFF为像点标准化的平移参数,LINE_SCALE、SAMPLE_SCALE为像点标准化的比例参数;地面点、像点坐标标准化的平移、缩放参数及有理多项式系数由卫星应用商以RPC文件的形式附加在影像元数据中提供。
根据观测值与未知数之间的函数关系(式(2))构建目标函数,将像点理论值与实测值间的残差平方和最小作为目标函数,如(3)所示:
f(X)=min(VTPV) (3)
其中,f(X)为目标函数,X为所有待估参数,包括连接点三维坐标及每景影像的6个定向参数,V表示像点理论值与实测值间的残差向量,P表示观测值的权。
S02:引入约束条件;
将平差后网中重心点的位置保持不变、网中各连接点至重心点的矢径角变化的带权总和为0、网中各点至重心点的矢径距变化的带权总和为0作为约束条件,如(4)所示;
STx=0 (4)
其中,x表示连接点坐标的改正数。
S03:对待估参数进行求解;
联合目标函数(式(3))与约束条件(式(4)),构建罚函数,如(5)所示:
其中,P(X,M)为罚函数,M为罚因子。
将约束非线性问题转化为求解无约束的罚函数极值问题,其计算公式如(6)所示:
minP(X,Mk)=P(X(k),Mk) (6)
其中,Mk为罚函数最小时的罚因子,X(k)为对应的最优待估参数值。
利用罚函数对未知参数求一阶导数,并令一阶导数为0,得到带有M的极值罚函数的解;取不同的M值,得到不同的解,并判断解是否满足约束条件,当满足约束条件,停止迭代,区域网平差结束,获得的解即为区域网平差解算得到的最优连接点及影像定向参数。
本发明的优点在于,将平差后网中重心点的位置保持不变、网中各连接点至重心点的矢径角变化的带权总和为0、网中各点至重心点的矢径距变化的带权总和为0作为约束条件,引入到常规基于RFM模型的区域网平差数学模型中,避免由于必要起算数据缺失引起的解不唯一问题,使得无控制点条件下卫星影像区域网平差解算成为可能,流程易于设计与实现;采用制约函数法对未知参数进行最优估计,保留函数模型的原始非线性形式,将带约束条件的非线性问题转化为求解无约束函数极值的问题,降低了高斯-牛顿法强制线性化引入的模型误差,从而保证无控制点情况下卫星影像区域网平差解的高精度;适用于境外、岛屿、沙漠等控制点难以获取的区域,满足环保、海洋、矿产等领域对影像内部拼接精度高于绝对精度的要求。
实施例2
参照图2,图2是一种基于制约函数非线性估计的卫星影像区域网平差方法的详细流程图。本发明的本发明的一种基于制约函数非线性估计的卫星影像区域网平差方法的具体流程如下。
S01:构建基于RFM的像方空间影像定位的函数模型
参照图3,图3是本发明的构建基于RFM的像方空间影像定位的函数模型的流程图。
S011:获取立体影像及RPC文件
获取立体影像以及影像元数据中的RPC文件。
S012:匹配获取足够数量的同名连接点
利用最小二乘影像匹配算法对立体影像序列进行匹配,获取足够数量的同名连接点,并建立连接点文件,连接点文件记录的信息包括每个连接点所在影像的编号、像点坐标等。
S013:构建RFM模型
利用RPC文件,构建描述物方点与对应像点之间定量关系的RFM模型,在RFM模型中,像点坐标(1,s)表示为以相应地面点空间坐标(φ,λ,h)为自变量的多项式的比值,如(7)所示:
其中,(ln,sn)和(P,L,H)分别表示像点坐标(1,s)和地面点坐标(φ,λ,h)经平移和缩放后的标准化坐标,取值位于(-1.0~+1.0)之间,其变换关系如(8)~(16)所示:
其中,LAT_OFF、LONG_OFF、HEIGHT_OFF、LINE_OFF、SAMP_OFF为标准化的平移参数,LAT_SCALE、LONG_SCALE、HEIGHT_SCALE、LINE_SCALE、SAMP_SCALE为标准化的比例参数,ai,bi,ci,di(i=0,1,2,…,19)为RPC文件。
S013:构建系统误差补偿模型
对每一个像点进行系统误差补偿,如(17)所示:
其中,Δp、Δr表示像方系统误差,(l,s)表示像点坐标,a0,as,a1,b0,bs,b1表示每景影像对应的6个仿射变换系数,作为待求定向参数,其中a0,as,a1修正影像行方向误差,b0,bs,b1修正影像列方向误差。
S014:建立区域网平差数学模型
将Δp、Δr代入RFM模型,即可得到像方空间影像区域网平差的函数模型,如(18)所示:
S015:构建目标函数
根据观测值与未知数之间的函数关系(式(18))构建目标函数,将像点理论值与实测值间的残差平方和最小作为目标函数,如(19)所示:
f(X)=min(VTPV) (19)
其中,f(X)为目标函数,X为所有待估参数,包括连接点三维坐标及每景影像的6个定向参数,V表示像点理论值与实测值间的残差向量,P表示观测值的权。
S02:引入约束条件
对每个连接点附加约束条件,约束条件为表示平差后网中重心点的位置保持不变、网中各连接点至重心点的矢径角变化的带权总和为0、网中各点至重心点的矢径距变化的带权总和为0,如(20)所示:
STx=0 (20)
其中,x表示连接点坐标的改正数。
S03:对待估参数进行求解
参照图4,图4是本发明的对待估参数进行求解的流程图
S031:构建罚函数
联合目标函数(式(19))与约束条件(式(20)),构建罚函数,如(21)所示:
其中,P(X,M)为罚函数,M为罚因子。
S032:选取初始罚因子
S033:计算对应的待估值
利用罚函数对未知参数求一阶导数,并令一阶导数为0,如(22)所示:
其中,X1,X2,X3,...为未知参数。
得到带有M的极值罚函数的解,将初始罚因子代入,得到对应的待估值。
S034:判断待估值与约束集的距离是否小于限差
将待估值与约束集进行对比,若得到的解不满足约束条件,变更M的值,重复步骤S033;
若满足约束条件,区域网平差结束,进入步骤S035。
S035:获得影像定向参数及加密点坐标
满足约束条件的罚函数如(23)所示:
minP(X,Mk)=P(X(k),Mk) (23)
其中,Mk为罚函数最小时的罚因子,X(k)为对应的最优待估参数值。
获得的解,即X(k),为区域网平差解算得到最优连接点三维坐标及影像定向参数。
本发明的优点在于,将平差后网中重心点的位置保持不变、网中各连接点至重心点的矢径角变化的带权总和为0、网中各点至重心点的矢径距变化的带权总和为0作为约束条件,可有效避免由于必要起算数据缺失引起的解不唯一问题,使得无控制点条件下卫星影像区域网平差解算成为可能,流程易于设计与实现。
实施例3
本发明的本发明的一种基于制约函数非线性估计的卫星影像区域网平差方法的具体实施方式如下。
S01:构建基于RFM的像方空间影像定位的函数模型
S011:获取立体影像及RPC文件
获取立体影像以及影像元数据中的RPC文件。
S012:匹配获取足够数量的同名连接点
利用最小二乘影像匹配算法对立体影像序列进行匹配,获取足够数量的同名连接点,并建立连接点文件,连接点文件记录的信息包括每个连接点所在影像的编号、像点坐标等。
S013:构建RFM模型
利用RPC文件,构建描述物方点与对应像点之间定量关系的RFM模型,在RFM模型中,像点坐标(l,s)表示为以相应地面点空间坐标(φ,λ,h)为白变量的多项式的比值,如(7)所示:
其中,(ln,sn)和(P,L,H)分别表示像点坐标(1,s)和地面点坐标(φ,λ,h)经平移和缩放后的标准化坐标,取值位于(-1.0~+1.0)之间,其变换关系如(8)~(16)所示:
其中,LAT_OFF、LONG_OFF、HEIGHT_OFF、LINE_OFF、SAMP_OFF为标准化的平移参数,LAT_SCALE、LONG_SCALE、HEIGHT_SCALE、LINE_SCALE、SAMP_SCALE为标准化的比例参数,ai,bi,ci,di(i=0,1,2,…,19)为RPC文件。
S013:构建系统误差补偿模型
对每一个像点进行系统误差补偿,如(17)所示:
其中,Δp、Δr表示像方系统误差,(1,s)表示像点坐标,a0,as,a1,b0,bs,b1表示每景影像对应的6个仿射变换系数,作为待求定向参数,其中a0,as,a1修正影像行方向误差,b0,bs,b1修正影像列方向误差。
S014:建立区域网平差数学模型
将Δp、Δr代入RFM模型,即可得到像方空间影像区域网平差的函数模型,如(18)所示:
S015:构建目标函数
根据观测值与未知数之间的函数关系(式(18))构建目标函数,将像点理论值与实测值间的残差平方和最小作为目标函数,如(19)所示:
f(X)=min(VTPV) (19)
其中,f(X)为目标函数,X为所有待估参数,包括连接点三维坐标及每景影像的6个定向参数,V表示像点理论值与实测值间的残差向量,P表示观测值的权。
S02:引入约束条件
对每个连接点附加约束条件,约束条件为表示平差后网中重心点的位置保持不变、网中各连接点至重心点的矢径角变化的带权总和为0、网中各点至重心点的矢径距变化的带权总和为0,如(20)所示:
STx=0 (20)
其中,x表示连接点坐标的改正数。
S03:对待估参数进行求解
S031:构建罚函数
联合目标函数(式(19))与约束条件(式(20)),构建罚函数,如(21)所示:
其中,P(X,M)为罚函数,M为罚因子。
S032:选取初始罚因子
S033:计算对应的待估值
利用罚函数对未知参数求一阶导数,并令一阶导数为0,如(22)所示:
其中,X1,X2,X3,...为未知参数。
得到带有M的极值罚函数的解,将初始罚因子代入,得到对应的待估值。
S034:判断待估值与约束集的距离是否小于限差
将待估值与约束集进行对比,若得到的解不满足约束条件,增大M的值,重复步骤S033;
若满足约束条件,区域网平差结束,进入步骤S035。
S035:获得影像定向参数及加密点坐标
满足约束条件的罚函数如(23)所示:
minP(X,Mk)=P(X(k),Mk) (23)
其中,Mk为罚函数最小时的罚因子,X(k)为对应的最优待估参数值。
获得的解,即X(k),为区域网平差解算得到最优连接点三维坐标及影像定向参数。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员,在不脱离本发明方法的前提下,还可以做出若干改进和补充,这些改进和补充也应视为本发明的保护范围。
Claims (5)
1.一种基于制约函数非线性估计的卫星影像区域网平差方法,其特征在于,所述基于制约函数非线性估计的卫星影像区域网平差方法包括以下步骤:
S01:构建基于RFM的像方空间影像定位的函数模型;
S011:获取立体影像及RPC文件;
S012:匹配获取足够数量的同名连接点;
S013:构建RFM模型:
利用RPC文件构建描述物方点与对应像点之间定量关系的RFM模型;
S014:构建系统误差补偿模型:
对每一个像点进行系统误差补偿,补偿模型如(1)所示:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>p</mi>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mi>s</mi>
</msub>
<mo>&CenterDot;</mo>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mi>l</mi>
</msub>
<mo>&CenterDot;</mo>
<mi>l</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>r</mi>
<mo>=</mo>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mi>s</mi>
</msub>
<mo>&CenterDot;</mo>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>b</mi>
<mi>l</mi>
</msub>
<mo>&CenterDot;</mo>
<mi>l</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Δp、Δr表示像方系统误差,(l,s)表示像点坐标,a0,as,a1,b0,bs,b1表示每景影像对应的6个仿射变换系数,作为待求定向参数,其中a0,as,a1修正影像行方向误差,b0,bs,b1修正影像列方向误差;
S015:建立区域网平差数学模型:
将Δp、Δr代入RFM模型,即可得到像方空间影像区域网平差的函数模型,如(2)所示:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mi>&Delta;</mi>
<mi>p</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mi>N</mi>
<mi>u</mi>
<mi>m</mi>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>,</mo>
<mi>L</mi>
<mo>,</mo>
<mi>H</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>D</mi>
<mi>e</mi>
<mi>n</mi>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>,</mo>
<mi>L</mi>
<mo>,</mo>
<mi>H</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mi>L</mi>
<mi>I</mi>
<mi>N</mi>
<mi>E</mi>
<mo>_</mo>
<mi>S</mi>
<mi>C</mi>
<mi>A</mi>
<mi>L</mi>
<mi>E</mi>
<mo>+</mo>
<mi>L</mi>
<mi>I</mi>
<mi>N</mi>
<mi>E</mi>
<mo>_</mo>
<mi>O</mi>
<mi>F</mi>
<mi>F</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>=</mo>
<mi>&Delta;</mi>
<mi>r</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mi>N</mi>
<mi>u</mi>
<mi>m</mi>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>,</mo>
<mi>L</mi>
<mo>,</mo>
<mi>H</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>D</mi>
<mi>e</mi>
<mi>n</mi>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>,</mo>
<mi>L</mi>
<mo>,</mo>
<mi>H</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mi>S</mi>
<mi>A</mi>
<mi>M</mi>
<mi>P</mi>
<mi>L</mi>
<mi>E</mi>
<mo>_</mo>
<mi>S</mi>
<mi>C</mi>
<mi>A</mi>
<mi>L</mi>
<mi>E</mi>
<mo>+</mo>
<mi>S</mi>
<mi>A</mi>
<mi>M</mi>
<mi>P</mi>
<mi>L</mi>
<mi>E</mi>
<mo>_</mo>
<mi>O</mi>
<mi>F</mi>
<mi>F</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,(P,L,H)表示地面点大地坐标经平移和缩放后的标准化坐标,取值位于(-1.0~+1.0)之间,NumL、DenL、NumS、DenS表示有理多项式,LINE_OFF、SAMPLE_OFF为像点标准化的平移参数,LINE_SCALE、SAMPLE_SCALE为像点标准化的比例参数;
S016:构建目标函数:
根据观测值与未知数之间的函数关系(式(2))构建目标函数,将像点理论值与实测值间的残差平方和最小作为目标函数,如(3)所示:
f(X)=min(VTPV) (3)
其中,f(X)为目标函数,X为所有待估参数,包括连接点三维坐标及每景影像的6个定向参数,V表示像点理论值与实测值间的残差向量,P表示观测值的权;
S02:引入约束条件:
约束条件为平差后网中重心点的位置保持不变,平差后网中各连接点至重心点的矢径角变化的带权总和为0、平差后网中各点至重心点的矢径距变化的带权总和为0,如(4)所示:
STx=0 (4)
其中,x表示连接点坐标的改正数;
S03:对待估参数进行求解;
S031:构建罚函数:
联合目标函数(式(3))与约束条件(式(4)),构建罚函数,如(5)所示:
<mrow>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>,</mo>
<mi>M</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>M</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>l</mi>
</munderover>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msup>
<msub>
<mi>S</mi>
<mi>j</mi>
</msub>
<mi>T</mi>
</msup>
<mi>x</mi>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,P(X,M)为罚函数,M为罚因子;
S032:选取初始罚因子;
S033:计算对应的待估值:
利用罚函数对未知参数求一阶导数,并令一阶导数为0,得到带有M的极值罚函数的解,将初始罚因子代入,得到解;
S034:判断待估值与约束集的距离是否小于限差:
若得到的解不满足约束条件,变更M的值,重复步骤S033;若满足约束条件,则进入步骤S035;
S035:获得影像定向参数及加密点坐标:
满足约束条件的罚函数如(6)所示:
min P(X,Mk)=P(X(k),Mk) (6)
其中,Mk为罚函数最小时的罚因子,X(k)为对应的最优待估参数值;
获得的解,即X(k),为区域网平差解算得到最优连接点三维坐标及影像定向参数。
2.根据权利要求1所述的基于制约函数非线性估计的卫星影像区域网平差方法,其特征在于,所述的步骤S012中,利用最小二乘影像匹配算法对立体影像序列进行匹配,获取足够数量的同名连接点,并建立连接点文件,所述的连接点文件记录的信息包括每个连接点所在影像的编号、像点坐标等。
3.根据权利要求2所述的基于制约函数非线性估计的卫星影像区域网平差方法,其特征在于,所述的步骤S013中,在RFM模型中,像点坐标(l,s)表示为以相应地面点空间坐标(φ,λ,h)为自变量的多项式的比值,如(7)所示:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>l</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>N</mi>
<mi>u</mi>
<mi>m</mi>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>,</mo>
<mi>L</mi>
<mo>,</mo>
<mi>H</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>D</mi>
<mi>e</mi>
<mi>n</mi>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>,</mo>
<mi>L</mi>
<mo>,</mo>
<mi>H</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>s</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>N</mi>
<mi>u</mi>
<mi>m</mi>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>,</mo>
<mi>L</mi>
<mo>,</mo>
<mi>H</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>D</mi>
<mi>e</mi>
<mi>n</mi>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>,</mo>
<mi>L</mi>
<mo>,</mo>
<mi>H</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,(ln,sn)和(P,L,H)分别表示像点坐标(1,s)和地面点坐标(φ,λ,h)经平移和缩放后的标准化坐标,取值位于(-1.0~+1.0)之间,其变换关系如(8)~(16)所示:
<mrow>
<mi>P</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&phi;</mi>
<mo>-</mo>
<mi>L</mi>
<mi>A</mi>
<mi>T</mi>
<mo>_</mo>
<mi>O</mi>
<mi>F</mi>
<mi>F</mi>
</mrow>
<mrow>
<mi>L</mi>
<mi>A</mi>
<mi>T</mi>
<mo>_</mo>
<mi>S</mi>
<mi>C</mi>
<mi>A</mi>
<mi>L</mi>
<mi>E</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>L</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&lambda;</mi>
<mo>-</mo>
<mi>L</mi>
<mi>O</mi>
<mi>N</mi>
<mi>G</mi>
<mo>_</mo>
<mi>O</mi>
<mi>F</mi>
<mi>F</mi>
</mrow>
<mrow>
<mi>L</mi>
<mi>O</mi>
<mi>N</mi>
<mi>G</mi>
<mo>_</mo>
<mi>S</mi>
<mi>C</mi>
<mi>A</mi>
<mi>L</mi>
<mi>E</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>H</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>h</mi>
<mo>-</mo>
<mi>H</mi>
<mi>E</mi>
<mi>I</mi>
<mi>G</mi>
<mi>H</mi>
<mi>T</mi>
<mo>_</mo>
<mi>O</mi>
<mi>F</mi>
<mi>F</mi>
</mrow>
<mrow>
<mi>H</mi>
<mi>E</mi>
<mi>I</mi>
<mi>G</mi>
<mi>H</mi>
<mi>T</mi>
<mo>_</mo>
<mi>S</mi>
<mi>C</mi>
<mi>A</mi>
<mi>L</mi>
<mi>E</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>l</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>l</mi>
<mo>-</mo>
<mi>L</mi>
<mi>I</mi>
<mi>N</mi>
<mi>E</mi>
<mo>_</mo>
<mi>O</mi>
<mi>F</mi>
<mi>F</mi>
</mrow>
<mrow>
<mi>L</mi>
<mi>I</mi>
<mi>N</mi>
<mi>E</mi>
<mo>_</mo>
<mi>S</mi>
<mi>C</mi>
<mi>A</mi>
<mi>L</mi>
<mi>E</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>s</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>s</mi>
<mo>-</mo>
<mi>S</mi>
<mi>A</mi>
<mi>M</mi>
<mi>P</mi>
<mo>_</mo>
<mi>O</mi>
<mi>F</mi>
<mi>F</mi>
</mrow>
<mrow>
<mi>S</mi>
<mi>A</mi>
<mi>M</mi>
<mi>P</mi>
<mo>_</mo>
<mi>S</mi>
<mi>C</mi>
<mi>A</mi>
<mi>L</mi>
<mi>E</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>N</mi>
<mi>u</mi>
<mi>m</mi>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>,</mo>
<mi>L</mi>
<mo>,</mo>
<mi>H</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mi>L</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mi>P</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>4</mn>
</msub>
<mi>L</mi>
<mi>P</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>5</mn>
</msub>
<mi>L</mi>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>6</mn>
</msub>
<mi>P</mi>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>7</mn>
</msub>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>8</mn>
</msub>
<msup>
<mi>P</mi>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>9</mn>
</msub>
<msup>
<mi>H</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>10</mn>
</msub>
<mi>P</mi>
<mi>L</mi>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>11</mn>
</msub>
<msup>
<mi>L</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>12</mn>
</msub>
<msup>
<mi>LP</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>13</mn>
</msub>
<msup>
<mi>LH</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>14</mn>
</msub>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
<mi>P</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>15</mn>
</msub>
<msup>
<mi>P</mi>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>16</mn>
</msub>
<msup>
<mi>PH</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>17</mn>
</msub>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>18</mn>
</msub>
<msup>
<mi>P</mi>
<mn>2</mn>
</msup>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>19</mn>
</msub>
<msup>
<mi>H</mi>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>D</mi>
<mi>e</mi>
<mi>n</mi>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>,</mo>
<mi>L</mi>
<mo>,</mo>
<mi>H</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mi>L</mi>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<mi>P</mi>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>3</mn>
</msub>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>4</mn>
</msub>
<mi>L</mi>
<mi>P</mi>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>5</mn>
</msub>
<mi>L</mi>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>6</mn>
</msub>
<mi>P</mi>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>7</mn>
</msub>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>8</mn>
</msub>
<msup>
<mi>P</mi>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>9</mn>
</msub>
<msup>
<mi>H</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>10</mn>
</msub>
<mi>P</mi>
<mi>L</mi>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>11</mn>
</msub>
<msup>
<mi>L</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>12</mn>
</msub>
<msup>
<mi>LP</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>B</mi>
<mn>13</mn>
</msub>
<msup>
<mi>LH</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>14</mn>
</msub>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
<mi>P</mi>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>15</mn>
</msub>
<msup>
<mi>P</mi>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>16</mn>
</msub>
<msup>
<mi>PH</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>17</mn>
</msub>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>18</mn>
</msub>
<msup>
<mi>P</mi>
<mn>2</mn>
</msup>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>19</mn>
</msub>
<msup>
<mi>H</mi>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>N</mi>
<mi>u</mi>
<mi>m</mi>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>,</mo>
<mi>L</mi>
<mo>,</mo>
<mi>H</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mi>L</mi>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<mi>P</mi>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>4</mn>
</msub>
<mi>L</mi>
<mi>P</mi>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>5</mn>
</msub>
<mi>L</mi>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>6</mn>
</msub>
<mi>P</mi>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>7</mn>
</msub>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>8</mn>
</msub>
<msup>
<mi>P</mi>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>9</mn>
</msub>
<msup>
<mi>H</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>10</mn>
</msub>
<mi>P</mi>
<mi>L</mi>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>11</mn>
</msub>
<msup>
<mi>L</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>12</mn>
</msub>
<msup>
<mi>LP</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>13</mn>
</msub>
<msup>
<mi>LH</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>14</mn>
</msub>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
<mi>P</mi>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>15</mn>
</msub>
<msup>
<mi>P</mi>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>16</mn>
</msub>
<msup>
<mi>PH</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>17</mn>
</msub>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>18</mn>
</msub>
<msup>
<mi>P</mi>
<mn>2</mn>
</msup>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>19</mn>
</msub>
<msup>
<mi>H</mi>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>N</mi>
<mi>u</mi>
<mi>m</mi>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>,</mo>
<mi>L</mi>
<mo>,</mo>
<mi>H</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mi>L</mi>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
<mi>P</mi>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>3</mn>
</msub>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>4</mn>
</msub>
<mi>L</mi>
<mi>P</mi>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>5</mn>
</msub>
<mi>L</mi>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>6</mn>
</msub>
<mi>P</mi>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>7</mn>
</msub>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>8</mn>
</msub>
<msup>
<mi>P</mi>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>9</mn>
</msub>
<msup>
<mi>H</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>10</mn>
</msub>
<mi>P</mi>
<mi>L</mi>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>11</mn>
</msub>
<msup>
<mi>L</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>12</mn>
</msub>
<msup>
<mi>LP</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>13</mn>
</msub>
<msup>
<mi>LH</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>14</mn>
</msub>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
<mi>P</mi>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>15</mn>
</msub>
<msup>
<mi>P</mi>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>16</mn>
</msub>
<msup>
<mi>PH</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>17</mn>
</msub>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>18</mn>
</msub>
<msup>
<mi>P</mi>
<mn>2</mn>
</msup>
<mi>H</mi>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>19</mn>
</msub>
<msup>
<mi>H</mi>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,LAT_OFF、LONG_OFF、HEIGHT_OFF、LINE_OFF、SAMP_OFF为标准化的平移参数,LAT_SCALE、LONG_SCALE、HEIGHT_SCALE、LINE_SCALE、SAMP_SCALE为标准化的比例参数,ai,bi,ci,di(i=0,1,2,…,19)为RPC文件。
4.根据权利要求3所述的基于制约函数非线性估计的卫星影像区域网平差方法,其特征在于,所述的步骤S033中,利用罚函数对未知参数求一阶导数,并令一阶导数为0,如(17)所示:
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>P</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>X</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>P</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>X</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>P</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>X</mi>
<mn>3</mn>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
<mo>...</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,X1,X2,X3,...为未知参数。
5.根据权利要求4所述的基于制约函数非线性估计的卫星影像区域网平差方法,其特征在于,所述的步骤S034中,若得到的解不满足约束条件,则增大M的值,对M进行迭代,重复步骤S033。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710583487.4A CN107314763B (zh) | 2017-07-18 | 2017-07-18 | 一种基于制约函数非线性估计的卫星影像区域网平差方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710583487.4A CN107314763B (zh) | 2017-07-18 | 2017-07-18 | 一种基于制约函数非线性估计的卫星影像区域网平差方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107314763A true CN107314763A (zh) | 2017-11-03 |
CN107314763B CN107314763B (zh) | 2019-05-07 |
Family
ID=60178489
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710583487.4A Active CN107314763B (zh) | 2017-07-18 | 2017-07-18 | 一种基于制约函数非线性估计的卫星影像区域网平差方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107314763B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108919319A (zh) * | 2018-05-15 | 2018-11-30 | 中国人民解放军战略支援部队信息工程大学 | 海岛礁卫星影像无地面控制点定位方法及系统 |
CN109752958A (zh) * | 2019-02-20 | 2019-05-14 | 成都市勘察测绘研究院 | 一种地面控制网更新维护的自适应网型约束平差方法 |
CN111369453A (zh) * | 2020-02-26 | 2020-07-03 | 武汉大学 | 一种基于平均高程面的影像快速几何预处理方法 |
CN112906198A (zh) * | 2021-01-21 | 2021-06-04 | 西南交通大学 | 任意平面测量控制网非线性最小二乘平差数据处理方法 |
CN112966612A (zh) * | 2021-03-10 | 2021-06-15 | 广东海洋大学 | 基于牛顿积分神经动力学的北极海冰遥感图像提取方法 |
CN113759402A (zh) * | 2020-06-03 | 2021-12-07 | 千寻位置网络有限公司 | Gnss形变监测方法、系统 |
CN113899386A (zh) * | 2021-09-27 | 2022-01-07 | 武汉大学 | 基于立体基准网的多源光学卫星遥感影像协同区域网平差方法及系统 |
CN114562982A (zh) * | 2022-03-09 | 2022-05-31 | 北京市遥感信息研究所 | 一种光学和sar异源卫星影像联合平差的定权方法和装置 |
CN117151991A (zh) * | 2023-07-06 | 2023-12-01 | 北京数慧时空信息技术有限公司 | 一种遥感影像rpc修复方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102168972A (zh) * | 2010-12-15 | 2011-08-31 | 中国资源卫星应用中心 | 一种基于rpc的三线阵立体卫星区域网平差改进及检校方法 |
CN102213762A (zh) * | 2011-04-12 | 2011-10-12 | 中交第二公路勘察设计研究院有限公司 | 基于rfm模型的多源星载sar影像自动匹配方法 |
CN103759714A (zh) * | 2014-01-20 | 2014-04-30 | 国家测绘地理信息局卫星测绘应用中心 | 一种三线阵卫星影像区域网平差方法 |
CN104112078A (zh) * | 2014-07-28 | 2014-10-22 | 中国科学院遥感与数字地球研究所 | 一种自适应几何关系强弱的有理函数模型区域网平差方法 |
-
2017
- 2017-07-18 CN CN201710583487.4A patent/CN107314763B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102168972A (zh) * | 2010-12-15 | 2011-08-31 | 中国资源卫星应用中心 | 一种基于rpc的三线阵立体卫星区域网平差改进及检校方法 |
CN102213762A (zh) * | 2011-04-12 | 2011-10-12 | 中交第二公路勘察设计研究院有限公司 | 基于rfm模型的多源星载sar影像自动匹配方法 |
CN103759714A (zh) * | 2014-01-20 | 2014-04-30 | 国家测绘地理信息局卫星测绘应用中心 | 一种三线阵卫星影像区域网平差方法 |
CN104112078A (zh) * | 2014-07-28 | 2014-10-22 | 中国科学院遥感与数字地球研究所 | 一种自适应几何关系强弱的有理函数模型区域网平差方法 |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108919319B (zh) * | 2018-05-15 | 2021-01-08 | 中国人民解放军战略支援部队信息工程大学 | 海岛礁卫星影像无地面控制点定位方法及系统 |
CN108919319A (zh) * | 2018-05-15 | 2018-11-30 | 中国人民解放军战略支援部队信息工程大学 | 海岛礁卫星影像无地面控制点定位方法及系统 |
CN109752958B (zh) * | 2019-02-20 | 2021-12-03 | 成都市勘察测绘研究院 | 一种地面控制网更新维护的自适应网型约束平差方法 |
CN109752958A (zh) * | 2019-02-20 | 2019-05-14 | 成都市勘察测绘研究院 | 一种地面控制网更新维护的自适应网型约束平差方法 |
CN111369453A (zh) * | 2020-02-26 | 2020-07-03 | 武汉大学 | 一种基于平均高程面的影像快速几何预处理方法 |
CN113759402A (zh) * | 2020-06-03 | 2021-12-07 | 千寻位置网络有限公司 | Gnss形变监测方法、系统 |
CN112906198A (zh) * | 2021-01-21 | 2021-06-04 | 西南交通大学 | 任意平面测量控制网非线性最小二乘平差数据处理方法 |
CN112966612A (zh) * | 2021-03-10 | 2021-06-15 | 广东海洋大学 | 基于牛顿积分神经动力学的北极海冰遥感图像提取方法 |
CN113899386A (zh) * | 2021-09-27 | 2022-01-07 | 武汉大学 | 基于立体基准网的多源光学卫星遥感影像协同区域网平差方法及系统 |
CN113899386B (zh) * | 2021-09-27 | 2023-11-21 | 武汉大学 | 基于立体基准网的多源光学卫星遥感影像协同区域网平差方法及系统 |
CN114562982A (zh) * | 2022-03-09 | 2022-05-31 | 北京市遥感信息研究所 | 一种光学和sar异源卫星影像联合平差的定权方法和装置 |
CN114562982B (zh) * | 2022-03-09 | 2023-09-26 | 北京市遥感信息研究所 | 一种光学和sar异源卫星影像联合平差的定权方法和装置 |
CN117151991A (zh) * | 2023-07-06 | 2023-12-01 | 北京数慧时空信息技术有限公司 | 一种遥感影像rpc修复方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107314763B (zh) | 2019-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107314763A (zh) | 一种基于制约函数非线性估计的卫星影像区域网平差方法 | |
Hu et al. | Understanding the rational function model: methods and applications | |
CN102506824B (zh) | 一种城市低空无人机系统生成数字正射影像图的方法 | |
CN112017224B (zh) | Sar数据区域网平差处理方法和系统 | |
CN105046251B (zh) | 一种基于环境一号卫星遥感影像的自动正射校正方法 | |
CN107192376B (zh) | 基于帧间连续性的无人机多帧图像目标定位校正方法 | |
CN112393714B (zh) | 一种基于无人机航拍与卫星遥感融合的影像校正方法 | |
CN109727278B (zh) | 一种机载LiDAR点云数据与航空影像的自动配准方法 | |
CN106403902A (zh) | 一种星地协同的光学卫星在轨实时几何定位方法及系统 | |
CN109709551B (zh) | 一种星载合成孔径雷达影像的区域网平面平差方法 | |
CN102410831B (zh) | 多条带扫描成像模型的设计及定位方法 | |
CN106127697A (zh) | 无人机机载成像高光谱几何校正方法 | |
CN110111274B (zh) | 一种星载推扫式光学传感器外方位元素定标方法 | |
CN102168972A (zh) | 一种基于rpc的三线阵立体卫星区域网平差改进及检校方法 | |
CN107330927B (zh) | 机载可见光图像定位方法 | |
CN110030968B (zh) | 一种基于星载立体光学影像的地面遮挡物仰角测量方法 | |
CN111003214B (zh) | 基于云控制的国产陆地观测卫星姿轨精化方法 | |
CN105004321B (zh) | 顾及曝光不同步的无人机gps辅助光束法平差方法 | |
CN104864853B (zh) | 一种高分辨率三线阵卫星沿轨方向姿态颤振的探测方法 | |
CN109188483B (zh) | 一种时序化高精度外方位元素自动定标方法 | |
CN113096047B (zh) | 广义云驱动与辐射协同遥感影像几何精校正方法及系统 | |
CN110986888A (zh) | 一种航空摄影一体化方法 | |
CN113255162A (zh) | 基于非刚性概率模型的车载激光点云自动误差纠正方法 | |
CN116758234A (zh) | 一种基于多点云数据融合的山地地形建模方法 | |
CN110443837B (zh) | 一种直线特征约束下的城区机载激光点云与航空影像配准方法和系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |