一种全球数字高程模型约束的影像区域网平差方法技术

本发明专利技术公开了一种全球数字高程模型约束的影像区域网平差方法，首先构建立体影像平差区域网，布设密集连接点。随后获取测区内全球数字高程模型数据，将测区划分为平地、丘陵和山地三种地形区域，设置利用全球数字高程模型数据针对平差计算获得的不同地形区域内连接点对应地面点高程值改正数的约束规则；最后开展区域网平差计算，在每次迭代计算时均利用该约束规则对连接点对应地面点的高程值改正数进行约束，直至平差结束，并根据平差结果为每张影像生成新的成像几何模型参数文件。本发明专利技术通过在不同地形区域采用不同规则使用全球数字高程模型对区域网平差进行高程约束，实现无地面控制条件下卫星遥感影像高程精度的大幅度提升。度提升。度提升。

【技术实现步骤摘要】
一种全球数字高程模型约束的影像区域网平差方法


[0001]本专利技术涉及摄影测量
，尤其涉及一种全球数字高程模型约束的影像区域网平差方法。

技术介绍

[0002]数字高程模型是用一组有序数值阵列形式表示地面高程的一种实体地面模型。目前，通过互联网可以公开免费获取的全球范围数字高程模型产品主要包括：
①
SRTM
‑
DEM(全称Space Shuttle Radar Topography Mission)：由美国太空总署等制作并免费提供，其数据覆盖范围为北纬60
°
到南纬56
°
之间的所有陆地区域，占地球陆地表面的80％，全球最高分辨率为30米，官方标称高程中误差10米。
②
ASTER GDEM(全称Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model)：由美国NASA和日本联合制作并免费提供，其数据覆盖范围为北纬83...

【技术保护点】

【技术特征摘要】
1.一种全球数字高程模型约束的影像区域网平差方法，其特征在于：包括如下步骤，S1、将测区内连续覆盖的遥感卫星立体影像构建一个平差区域网，并采用像方仿射变换补偿的基于有理函数模型的区域网平差模型；有理函数模型即RFM；S2、在平差区域网内立体影像上布设密集的连接点；S3、选取并从互联网上下载测区内一种或多种免费全球数字高程模型产品，并根据地表高程起伏情况将测区概略地划分为平地、丘陵和山地三种地形区域，结合全球数字高程模型产品在不同地形区域的高程精度经验值，分别为平差计算时求解的位于三种地形区域的连接点对应地面点高程值改正数设置不同的高程约束规则；S4、确定连接点对应地面点的平面坐标初值、高程值初值以及影像RFM仿射变换模型参数初值；S5、针对每个连接点，将影像RFM的仿射变换参数和连接点对应地面点的三维坐标作为未知数，逐点构建连接点的误差方程；S6、对误差方程进行法化，获取法方程；采用最小二乘平差原理对法方程进行整体平差计算，计算每张影像RFM的仿射变换参数以及连接点对应地面点的三维坐标的改正数；S7、每完成一次平差计算，判断是否满足平差收敛条件，若是，则平差计算结束；若否，则利用平差计算获取的改正数更新影像RFM的仿射变换参数以及连接点对应地面点的三维坐标，同时利用高程约束规则对更新后的不同地形区域连接点对应地面点的高程值进行修改，重复步骤S5到S6，再次进行平差迭代计算，直到满足收敛条件为止；S8、利用最后一次平差迭代计算得到的影像RFM的仿射变换参数，对影像RFM进行像方仿射变换补偿，为每张影像生成新的RFM。2.根据权利要求1所述的全球数字高程模型约束的影像区域网平差方法，其特征在于：步骤S1中的RFM是利用有理多项式建立影像的像方与其对应的物方之间的数学映射关系，RFM方程定义为，RFM方程定义为，其中，(P,L,H)为归一化后的地面坐标；(x,y)为归一化后的影像坐标；Num
S
(P,L,H)、Den
S
(P,L,H)、Num
L
(P,L,H)和Den
L
(P,L,H)分别表示一般多项式，多项式中每一个变量P,L,H的幂均不超过3次，所有变量的幂之和也不超过3次；影像RFM的误差补偿采用像方仿射变换模型，则RFM描述的像点坐标(x,y)和地面点坐标(P,L,H)之间的关系被修正为，标(P,L,H)之间的关系被修正为，其中，(
△
x,
△
y)为影像坐标(x,y)的系统误差像方补偿值，其值为，
△
y＝a0+a1·
y+a2·
x
△
x＝b0+b1·
y+b2·
x其中，(a0,a1,a2,b0,b1,b2)表示仿射变换模型的参数。3.根据权利要求1所述的全球数字高程模型约束的影像区域网平差方法，其特征在于：步骤S2中先采用影像自动匹配算法匹配连接点，当影像自动匹配算法匹配的连接点分布不均匀、数量不足或维度不够时，采用人工判读方法补测连接点。4.根据权利要求1所述的全球数字高程模型约束的影像区域网平差方法，其特征在于：步骤S3中将测区划分为平地、丘陵和山地三种地形区域时，利用全球数字高程模型计算各立体影像覆盖范围的地面坡度角和高程差，当立体影像覆盖范围内地面坡度角小于2度或高程差小于80米，则该立体影像覆盖区域为平地地形区域；当立体影像覆盖范围内地面坡度角小于6度或高程差小于300米，则该立体影像覆盖区域为丘陵地形区域；否则则为山地地形区域。5.根据权利要求1所述的全球数字高程模型约束的影像区域网平差方法，其特征在于：步骤S3中，根据所采用的全球数字高程模型产品在不同地形区域的高程精度经验值，分别设置每次平差迭代时计算平地、丘陵和山地地形区域连接点对应地面点高程改正值的高程约束规则，具体如下，h
plan_corr
∈[(h
init
‑
σ
plan
)，(h
init
+σ
plan
)]h
hill_corr
∈[(h
init
‑
σ
hill
)，(h
init
+σ
hill
)]h
mount_corr
∈[(h
init
‑
σ
mount
)，(h
init
+σ
mount
)]其中，h
plan...

【专利技术属性】
技术研发人员：周平，唐新明，王懿哲，王霞，刘昌儒，王艺颖，岳明宇，郭莉，李丹丹，
申请(专利权)人：自然资源部国土卫星遥感应用中心，
类型：发明
国别省市：