CN113532377B - 一种利用高分七号激光测高数据辅助区域网平差的方法 - Google Patents

一种利用高分七号激光测高数据辅助区域网平差的方法 Download PDF

Info

Publication number
CN113532377B
CN113532377B CN202110795596.9A CN202110795596A CN113532377B CN 113532377 B CN113532377 B CN 113532377B CN 202110795596 A CN202110795596 A CN 202110795596A CN 113532377 B CN113532377 B CN 113532377B
Authority
CN
China
Prior art keywords
image
point
points
coordinates
laser
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
Application number
CN202110795596.9A
Other languages
English (en)
Other versions
CN113532377A (zh
Inventor
唐新明
刘昌儒
李鸿洲
薛白
李国元
王霞
王洋洋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ministry Of Natural Resources Land Satellite Remote Sensing Application Center
Original Assignee
Ministry Of Natural Resources Land Satellite Remote Sensing Application Center
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Ministry Of Natural Resources Land Satellite Remote Sensing Application Center filed Critical Ministry Of Natural Resources Land Satellite Remote Sensing Application Center
Priority to CN202110795596.9A priority Critical patent/CN113532377B/zh
Publication of CN113532377A publication Critical patent/CN113532377A/zh
Application granted granted Critical
Publication of CN113532377B publication Critical patent/CN113532377B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/02—Affine transformations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/80—Geometric correction
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/50—Depth or shape recovery
    • G06T7/521—Depth or shape recovery from laser ranging, e.g. using interferometry; from the projection of structured light

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Optics & Photonics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Image Processing (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种利用高分七号激光测高数据辅助区域网平差的方法，该方法结合高分七号激光测高数据足印影像，实现了基于激光测高数据的高程控制点自动提取，建立了激光测高数据辅助立体影像的联合区域网平差模型，实现了无控、外业控制以及高程控制条件下的模型统一；将激光测高数据作为高程控制应用到平差计算中，实现区域影像几何精度提升。通过本发明中记载的方法能够在覆盖不同地形类型的区域，仅使用激光测高数据作为高程控制，便可实现区域网高程精度的极大提升。

Description

一种利用高分七号激光测高数据辅助区域网平差的方法
技术领域
本发明涉及卫星数据处理技术领域，尤其涉及一种利用高分七号激光测高数据辅助区域网平差的方法。
背景技术
提高稀少甚至无地面控制点的区域网平差精度，是实现境外和外业测控困难区域高精度测图的核心问题之一，也是主要的技术难点。高分七号激光测高数据具有极高的高程精度，并与立体影像同步获取，误差特性相近、相对精度较高，为了充分利用上述特点，发挥足印影像激光光斑地面信息的作用，亟待需要一种高分七号卫星激光测高数据辅助的立体影像联合区域网平差方法。
发明内容
本发明的目的在于提供一种利用高分七号激光测高数据辅助区域网平差的方法，从而解决现有技术中存在的前述问题。
为了实现上述目的，本发明采用的技术方案如下：
一种利用高分七号激光测高数据辅助区域网平差的方法，包括以下步骤：
S1,采用高分七号激光测高系统获取激光测高原始数据，经过一系列精细改正处理过程后形成可用的激光测高标准产品，包括激光测高点的三维坐标、足印影像以及激光测高点足印影像像点坐标；
S2，将获取的立体影像构建区域网，并进行连接点匹配，获取数量足够、分布均匀的连接点，开展基于有理函数模型的自由网平差，实现区域网影像高精度相对定向，得到相对定向后的立体影像；
S3，根据步骤S1中获得的激光测高点的三维坐标、足印影像以及激光测高点足印影像像点坐标和步骤S2中获得的相对定向后的立体影像，开展基于激光测高点的高程控制点提取过程，获取符合提取要求的激光高程控制点；
S4,利用获取到的激光高程控制点作为高程控制，与相对定向后的立体影像开展联合区域网平差，实现立体影像高程精度提升；
S5，更新立体影像有理函数模型，获取精度提升后的立体影像。
步骤S3中在进行高程控制点提取过程包括A方案和B方案，对于每个激光测高点，首先采用A方案进行高程控制点的提取，当A方案高程控制点提取过程失败之后则采用B方案进行高程控制点的提取过程。
优选的，所述A方案包括：
A31，将激光测高点物方三维坐标代入立体影像有理函数模型，计算得到当前激光测高点在立体影像上的初略位置，设其像点坐标为(p,q),则该点与步骤S1中的激光测高点足印影像像点坐标(x0,y0)构成共轭点对；
A32，以像点坐标(p,q)为中心，以激光光斑半径为最大搜索窗口，按公式(1)逐点计算与以激光测高点足印影像像点坐标(x0,y0)为中心的区域影像相关系数，取最大相关系数点位(x1',y1')作为像素级配准点位：
Figure GDA0003400521640000021
式中，g,g'分别为足印影像和立体影像灰度值，
Figure GDA0003400521640000022
为匹配窗口的平均灰度值，w,h分别表示相关系数搜索窗口的宽和高；
A33，以坐标(x0,y0),(x1',y1')为初值，按公式(2)进行最小二乘匹配，获取子像素级配准点位(x1,y1)’即激光测高点在落点影像上的像点坐标:
Figure GDA0003400521640000031
式中，h0,h1为影像辐射畸变参数，m'j,n'j为影像几何变换参数, j＝0,1,2；
A34，重复步骤A31-A33 ，获取当前激光测高点在所有落点影像上的像点坐标；
A35，重复步骤A31-A34，获取所有激光测高点在立体影像中的像点坐标。
优选的，所述B方案包括：
B31,根据激光测高点三维坐标，通过影像几何模型正变换，获取激光测高点在立体影像上的像点坐标概略位置；
B32,以步骤B31中求得的像点坐标为原点，在激光光斑半径范围内对立体影像与足印影像进行同名点匹配，经过粗差剔除后，得到一定数量的候选点；
B33，利用公式(1)计算所有候选点与激光测高点足印影像像点的相关系数，从中选择离激光光斑中心最近且相关系数最大的点对应的像点坐标作为当前激光测高点在落点影像上的像点坐标，；
B34，重复步骤B31-B33，获取当前激光测高点在所有落点影像上的像点坐标；
B35,重复步骤B31-B34，获取所有激光测高点在立体影像中的像点坐标。
优选的，步骤S1中所述一系列精细改正处理过程包括但不限于全波形处理、几何定位、大气以及潮汐误差改正处理过程。
优选的，步骤S2中所述的自由网平差包括以下步骤：
S21，利用卫星影像自带的RPC文件，构建有理函数模型，表达式如公式 (3)所示：
Figure GDA0003400521640000041
其中，pi为一般多项式，表达式为：
pi＝a0+a1·Y+a2·X+a3·Z+a4·Y·X+a5·Y·Z+a6·X·Z+a7·Y2+ a8·X2+a9·Z2+a10·X·Y·Z+a11·Y3+a12·Y·X2+a13·Y·Z2+ a14·Y2·X+a15·X3+a16·X·Z2+a17·Y2·Z+a18·X2·Z+a19·Z3
其中，(X,Y,Z)为正则化后的地面点三维坐标，(x,y)为对应的正则化像点坐标，ai(i＝0,1,...,18,19)为有理函数系数,；
S22，采用仿射变换模型对公式(3)进行几何定位误差补偿，表达式如公式(4)所示：
Figure GDA0003400521640000042
其中，(Δx,Δy)为像方改正数，(x,y)为利用公式(3)求得的像点坐标，mj,nj(j＝0,1,2)为仿射变换模型参数；
S23，影像仿射变换参数和连接点物方坐标是要求解的两类未知量，将公式(3)-(4)线性化，可得误差方程，如公式(5)所示：
V1＝A1X+B1Y1-l1,P1 (5)
式中，X＝[Δm0 Δm1 Δm2 Δn0 Δn1 Δn2]T，是影像的仿射变换参数改正数，Y1＝[Δlon Δlat Δh]T是连接点物方坐标改正值，A1,B1为对应的系数矩阵，l1为初值计算的常量，P1为权值。
优选的，步骤S4中所述的联合区域网平差具体包括：
S41，利用卫星影像自带的RPC文件，构建有理函数模型，表达式如公式(6)所示：
Figure GDA0003400521640000051
其中，pi为一般多项式，表达式为：
pi＝a0+a1·Y+a2·X+a3·Z+a4·Y·X+a5·Y·Z+a6·X·Z+a7·Y2+ a8·X2+a9·Z2+a10·X·Y·Z+a11·Y3+a12·Y·X2+a13·Y·Z2+ a14·Y2·X+a15·X3+a16·X·Z2+a17·Y2·Z+a18·X2·Z+a19·Z3
其中，(X,Y,Z)为正则化后的地面点三维坐标，(x,y)为对应的正则化像点坐标，ai(i＝0,1,...,18,19)为有理函数系数；
S42，采用仿射变换模型对公式(6)进行几何定位误差补偿，表达式如公式(7)所示：
Figure GDA0003400521640000052
其中，(Δx,Δy)为像方改正数，(x,y)为利用公式(6)求得的像点坐标，mj,nj(j＝0,1,2)为仿射变换模型参数；
S43，影像仿射变换参数和连接点物方坐标是要求解的两类未知量，将公式(6)-(7)线性化，可得误差方程，如公式(8)所示：
V1＝A1X+B1Y1-l1,P1 (8)
式中，X＝[Δm0 Δm1 Δm2 Δn0 Δn1 Δn2]T，是影像的仿射变换参数改正数，Y1＝[Δlon Δlat Δh]T是连接点物方坐标改正值，A1,B1为对应的系数矩阵，l1为初值计算的常量，P1为权值；
S44，对于激光测高点，其高程作为控制，在平差中仅改正其平面坐标，此时，误差方程变为：
V2＝B2Y2-l2,P2 (9)
式中，Y2＝[Δlon Δlat]T是激光测高点物方坐标改正值，B2是相应的系数矩阵，l2是初值计算常量，P2是相应权值。
优选的，步骤S4还包括：当存在其他外业控制点时，开展高程控制点、外业控制点与立体影像的联合区域网平差，实现立体影像高程和平面精度同步提升,此时还需要将外业控制点参与联合区域网平差，所涉及到的误差方程为：
V3＝A2X-l3,P3 (10)
式中，X＝[Δm0 Δm1 Δm2 Δn0 Δn1 Δn2]T是影像的仿射变换参数改正数，A2是相应的系数矩阵，l3是初值计算常量，P3是相应权值。
优选的，步骤S2中的立体影像包括但不限于采用高分七号卫星获取的立体影像。
本发明的有益效果是：
本发明公开了一种利用高分七号激光测高数据辅助区域网平差的方法，该方法结合高分七号激光测高数据中的足印影像，实现了基于激光测高数据的高程控制点自动提取，建立了激光测高数据高程控制点辅助立体影像的联合区域网平差模型，实现了无控、外业控制以及高程控制条件下的模型统一；将激光测高数据高程控制点作为高程控制应用到平差计算中，实现区域影像几何精度提升。通过本发明中记载的方法能够在覆盖不同地形类型的区域，仅使用激光测高数据作为高程控制，便可实现区域网高程精度的极大提升。
附图说明
图1是实施例1中提供的利用高分七号激光测高数据辅助区域网平差的方法流程示意图；
图2是实施例1中测区地形与影像、外业控制点和激光测高点分布示意图；
图3是实施例1中自由网平差平面与高程残差图；
图4是实施例1中激光测高点做高程控制的平面和高程残差图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施方式仅仅用以解释本发明，并不用于限定本发明。
实施例1
本实施例提供了一种利用高分七号激光测高数据辅助区域网平差的方法，如图1所示，包括以下步骤：
S1,采用高分七号激光测高系统获取激光测高原始数据，经过一系列精细改正处理过程后形成可用的激光测高标准产品，包括激光测高点的三维坐标、足印影像以及激光测高点足印影像像点坐标；
S2，将获取的立体影像构建区域网，并进行连接点匹配，获取数量足够、分布均匀的连接点，开展基于有理函数模型的自由网平差，实现区域网影像高精度相对定向，得到相对定向后的立体影像；
S3，根据步骤S1中获得的激光测高点的三维坐标、足印影像以及激光测高点足印影像像点坐标和步骤S2中获得的相对定向后的立体影像，开展基于激光测高点的高程控制点提取过程，获取符合提取要求的激光高程控制点；
S4,利用获取到的激光高程控制点作为高程控制，与相对定向后的立体影像开展联合区域网平差，实现立体影像高程精度提升；
S5，更新立体影像有理函数模型，获取精度提升后的立体影像。
本实施例中的步骤S1中，获取激光测高数据具体为：在工作状态下，高分七号卫星激光测高仪2个激光束以3Hz频率向地面发射激光脉冲，在地面形成沿轨向间隔约2.4km、垂轨向约12.5km的离散激光足印光斑，在幅宽20km的光学影像范围内约有2列共16个激光点。
高分七号卫星激光测高仪同时配置了激光足印相机，用于同步获取激光光斑地面影像，经过后期处理，可得出激光光斑中心在足印影像上的位置。
本实施例中通过高分七号卫星激光测高系统获得的原始数据在经过全波形处理、几何定位以及大气、潮汐等各项改正处理后形成激光测高产品，并提供各测高点的三维坐标及其足印影像坐标；高分七号卫星激光测高数据平面采用WGS84椭球基准，高程为WGS84椭球高。
本实施例中，步骤S2中所述的自由网平差包括以下步骤：
S21，利用卫星影像自带的RPC文件，构建有理函数模型，表达式如公式 (3)所示：
Figure GDA0003400521640000081
其中，pi为一般多项式，表达式为：
pi＝a0+a1·Y+a2·X+a3·Z+a4·Y·X+a5·Y·Z+a6·X·Z+a7·Y2+ a8·X2+a9·Z2+a10·X·Y·Z+a11·Y3+a12·Y·X2+a13·Y·Z2+ a14·Y2·X+a15·X3+a16·X·Z2+a17·Y2·Z+a18·X2·Z+a19·Z3
其中，(X,Y,Z)为正则化后的地面点三维坐标，(x,y)为对应的正则化像点坐标，ai(i＝0,1,...,18,19)为有理函数系数；
S22，采用仿射变换模型对公式(3)进行几何定位误差补偿，表达式如公式(4)所示：
Figure GDA0003400521640000091
其中，(Δx,Δy)为像方改正数，(x,y)为利用公式(3)求得的像点坐标，mj,nj(j＝0,1,2)为仿射变换模型参数；
S23，影像仿射变换参数和连接点物方坐标是要求解的两类未知量，将公式(3)-(4)线性化，可得误差方程，如公式(5)所示：
V1＝A1X+B1Y1-l1,P1 (5)
式中，X＝[Δm0 Δm1 Δm2 Δn0 Δn1 Δn2]T，是影像的仿射变换参数改正数，Y1＝[Δlon Δlat Δh]T是连接点物方坐标改正值，A1,B1为对应的系数矩阵，l1为初值计算的常量，P1为权值。
本实施例中的步骤S3中在进行高程控制点提取过程包括A方案和B 方案，对于每个激光测高点，首先采用A方案进行高程控制点的提取，当 A方案高程控制点提取过程失败之后则采用B方案进行高程控制点的提取过程。
所述A方案包括：
A31，将激光测高点物方三维坐标代入立体影像有理函数模型，计算得到当前激光测高点在立体影像上的初略位置，设其像点坐标为(p,q),则该点与步骤S1中的激光测高点像点坐标(x0,y0)构成共轭点对；
A32，以像点坐标(p,q)为中心，以激光光斑半径为最大搜索窗口，按公式(1)逐点计算与以激光测高点像点坐标(x0,y0)为中心的区域影像相关系数，取最大相关系数点位(x1',y1')作为像素级配准点位：
Figure GDA0003400521640000101
式中，g,g'分别为足印影像和立体影像灰度值，
Figure GDA0003400521640000102
为匹配窗口的平均灰度值，w,h分别表示相关系数搜索窗口的宽和高；
A33，以坐标(x0,y0),(x1',y1')为初值，按公式(2)进行最小二乘匹配，获取子像素级配准点位(x1,y1):
Figure GDA0003400521640000103
式中，h0,h1为影像辐射畸变参数，m'j,n'j为影像几何变换参数, j＝0,1,2；
A34，重复步骤A31-S33，获取当前激光测高点在所有落点影像上的像点坐标；
A35，重复步骤A31-A34，获取所有激光测高点在立体影像中的像点坐标。
所述B方案包括：
B31,根据激光测高点三维坐标，通过影像几何模型正变换，获取激光测高点在立体影像上的像点坐标位置；
B32,以步骤B31中求得的像点坐标为原点，在激光光斑半径范围内对立体影像与足印影像进行同名点匹配，经过粗差剔除后，得到一定数量的候选点；
B33，从候选点中选择离激光光斑中心最近且相关系数最大的点作为最终的激光测高点像点坐标，相关系数由公式(1)计算确定；
B34，重复步骤B31-B33，获取当前激光测高点在所有落点影像上的像点坐标；
B35,重复步骤B31-B34，获取所有激光测高点在立体影像中的像点坐标。
在进行高程控制点的提取过程中，由于A、B方案的混合使用，可能使得一个激光测高点在不同的影像上对应不同的物方地物，也就是一个激光测高点对应多个高程控制点，这就需要对每一个高程控制点进行独立平差操作，如果此时仍作为单个高程控制点进行平差计算，则会产生平差误差，这个误差的计算方法可以采用以下公式进行计算：
Figure GDA0003400521640000111
式中，Δh为高程误差，ye为影像地面分辨率，Δr为像点量测误差， B/H为基高比。
对于高分七号卫星来说，前后视沿轨向夹角分别为26度和5度，基高比为0.575，则影像上高程控制点0.5像素的坐标差异便会引起约0.8米的高程误差，这远远大于激光测高点的标称高程精度，因此要特别重视同一激光测高点衍生多个高程控制点的情况。
步骤S4中所述的联合区域网平差具体包括：
S41，利用卫星影像自带的RPC文件，构建有理函数模型，表达式如公式(6)所示：
Figure GDA0003400521640000112
其中，pi为一般多项式，表达式为：
pi＝a0+a1·Y+a2·X+a3·Z+a4·Y·X+a5·Y·Z+a6·X·Z+a7·Y2+ a8·X2+a9·Z2+a10·X·Y·Z+a11·Y3+a12·Y·X2+a13·Y·Z2+ a14·Y2·X+a15·X3+a16·X·Z2+a17·Y2·Z+a18·X2·Z+a19·Z3
其中，(X,Y,Z)为正则化后的地面点三维坐标，(x,y)为对应的正则化像点坐标，ai(i＝0,1,...,18,19)为有理函数系数；
S42，采用仿射变换模型对公式(6)进行几何定位误差补偿，表达式如公式(7)所示：
Figure GDA0003400521640000121
其中，(Δx,Δy)为像方改正数，(x,y)为利用公式(6)求得的像点坐标，mj,nj为仿射变换模型参数，j＝0,1,2；
S43，影像仿射变换参数和连接点物方坐标是要求解的两类未知量，将公式(6)-(7)线性化，可得误差方程，如公式(8)所示：
V1＝A1X+B1Y1-l1,P1 (8)
式中，X＝[Δm0 Δm1 Δm2 Δn0 Δn1 Δn2]T，是影像的仿射变换参数改正数，Y1＝[Δlon Δlat Δh]T是连接点物方坐标改正值，A1,B1为对应的系数矩阵，l1为初值计算的常量，P1为权值；
S44，对于激光测高点，其高程作为控制，在平差中仅改正其平面坐标，此时，误差方程变为：
V2＝B2Y2-l2,P2 (9)
式中，Y2＝[Δlon Δlat]T是激光测高点物方坐标改正值，B2是相应的系数矩阵，l2是初值计算常量，P2是相应权值。
当存在其他外业控制点时，开展高程控制点、外业控制点与立体影像的联合区域网平差，实现立体影像高程和平面精度同步提升,此时还需要将外业控制点参与联合区域网平差，所涉及到的误差方程为：
V3＝A2X-l3,P3 (10)
式中，X＝[Δm0 Δm1 Δm2 Δn0 Δn1 Δn2]T是影像的仿射变换参数改正数，A2是相应的系数矩阵，l3是初值计算常量，P3是相应权值。
采用上述的方案，选择山东省中部某区域作为实验区域，该区域面积约 2.75万平方公里，整体南北长东西窄。区域整体地形起伏从海拔-5米到 1000米，其中中部是山地和丘陵，周边以平地为主。共收集8轨共80对高分七号立体影像，拍摄时间集中在2020年4月至2021年1月。共收集同轨获取激光测高点5轨602个，以及异轨获取激光测高点1轨67个。经过高程控制点提取，共获得328个有效高程控制点。为对平差结果进行验证分析，在试验区测量了109个高精度外业控制点(平面和高程精度均优于 0.15米)，所有控制点像点坐标均在立体环境下量测，量测精度优于0.3像素。整个测区地形以及控制点、激光测高点分布如图2所示。
根据地形坡度以及高差信息，将山东测区划分为平地、丘陵和山地等3 个不同地形。将激光测高点作为高程控制点与立体影像进行联合区域网平差，验证激光测高点在不同地形条件下对立体影像几何定位精度的影响效果，其结果如下表所示：
山东测区激光测高点对不同地形的平差结果统计表
Figure GDA0003400521640000141
从上表可以看出，在没有任何外部控制条件下，高分七号立体影像整体平面定位精度在9.5米左右，其中山地区域的平面精度好于平地和丘陵，并且最大误差也较小，而平地和丘陵区域精度相当。高程方面，自由网平差后的高程精度在7.9米左右，虽然山地区域高程精度好一些，但是系统误差比较明显，全区域高程系统误差达到了5米。从图3也可以看出，测区平面和高程误差呈现明显的系统误差。经过激光测高点辅助区域网平差后，测区整体平面精度没有较大变化，平面精度依然在8.7米左右，平面最大误差也与自由网平差时相近，同时也没有明显的地形区别，但是高程精度迅速提升，整体达到0.76米，特别是山地区域，高程精度达到0.46米，而高程最大误差也降到了1.72米，相比于自由网平差，高程精度提升非常明显。
从图4可以看出，因为没有平面约束，联合平差后平面误差依然显示为系统误差，整体误差趋势与图3相似，但是高程误差规律不明显，数值降低很多。以上实验结果可以看出本发明提供的方法使用激光测高点对高程精度具有较好的提升作用，特别是山地区域，效果更加明显。
通过采用本发明公开的上述技术方案，得到了如下有益的效果：
本发明公开了一种利用高分七号激光测高数据辅助区域网平差的方法，该方法结合高分七号激光测高数据中的足印影像，实现了基于激光测高数据的高程控制点自动提取，建立了激光测高数据高程控制点辅助立体影像的联合区域网平差模型，实现了无控、外业控制以及高程控制条件下的模型统一；将激光测高数据高程控制点作为高程控制应用到平差计算中，实现区域影像几何精度提升。通过本发明中记载的方法能够在覆盖不同地形类型的区域，仅使用激光测高数据作为高程控制，便可实现区域网高程精度的极大提升。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视本发明的保护范围。

Claims (6)

1.一种利用高分七号激光测高数据辅助区域网平差的方法，其特征在于，包括以下步骤：
S1,采用高分七号激光测高系统获取激光测高原始数据，经过一系列精细改正处理过程后形成可用的激光测高标准产品，包括激光测高点的三维坐标、足印影像以及激光测高点足印影像像点坐标；
S2，将获取的立体影像构建区域网，并进行连接点匹配，获取数量足够、分布均匀的连接点，开展基于有理函数模型的自由网平差，实现区域网影像高精度相对定向，得到相对定向后的立体影像；
S3，根据步骤S1中获得的激光测高点的三维坐标、足印影像以及激光测高点足印影像像点坐标和步骤S2中获得的相对定向后的立体影像，开展基于激光测高点的高程控制点提取过程，获取符合提取要求的激光高程控制点；
S4,利用获取到的激光高程控制点作为高程控制，与相对定向后的立体影像开展联合区域网平差，实现立体影像高程精度提升；
S5，更新立体影像有理函数模型，获取精度提升后的立体影像；
步骤S3中在进行高程控制点提取过程包括A方案和B方案，对于每个激光测高点，首先采用A方案进行高程控制点的提取，当A方案高程控制点提取过程失败之后则采用B方案进行高程控制点的提取过程；
所述A方案包括：
A31，将激光测高点物方三维坐标代入立体影像有理函数模型，计算得到当前激光测高点在立体影像上的初略位置，设其像点坐标为(p,q),则该点与步骤S1中的激光测高点足印影像像点坐标(x0,y0)构成共轭点对；
A32，以像点坐标(p,q)为中心，以激光光斑半径为最大搜索窗口，按公式(1)逐点计算与以激光测高点足印影像像点坐标(x0,y0)为中心的区域影像相关系数，取最大相关系数点位(x1',y1')作为像素级配准点位：
Figure FDA0003400521630000021
式中，g,g'分别为足印影像和立体影像灰度值，
Figure FDA0003400521630000022
为匹配窗口的平均灰度值，w,h分别表示相关系数搜索窗口的宽和高；
A33，以坐标(x0,y0),(x1',y1')为初值，按公式(2)进行最小二乘匹配，获取子像素级配准点位(x1,y1)，即激光测高点在落点影像上的像点坐标：
Figure FDA0003400521630000023
式中，h0,h1为影像辐射畸变参数，m'j,n'j(j＝0,1,2)为影像几何变换参数；
A34，重复步骤A31-A33 ，获取当前激光测高点在所有落点影像上的像点坐标；
A35，重复步骤A31-A34，获取所有激光测高点在立体影像中的像点坐标；所述B方案包括：
B31,根据激光测高点三维坐标，通过影像几何模型正变换，获取激光测高点在立体影像上的像点坐标概略位置；
B32,以步骤B31中求得的像点坐标为原点，在激光光斑半径范围内对立体影像与足印影像进行同名点匹配，经过粗差剔除后，得到一定数量的候选点；
B33，利用公式(1)计算所有候选点与激光测高点足印影像像点的相关系数，从中选择离激光光斑中心最近且相关系数最大的点对应的像点坐标作为当前激光测高点在落点影像上的像点坐标；
B34，重复步骤B31-B33，获取当前激光测高点在所有落点影像上的像点坐标；
B35,重复步骤B31-B34，获取所有激光测高点在立体影像中的像点坐标。
2.根据权利要求1所述的利用高分七号激光测高数据辅助区域网平差的方法，其特征在于，步骤S1中所述一系列精细改正处理过程包括全波形处理、几何定位、大气以及潮汐误差改正处理过程。
3.根据权利要求1所述的利用高分七号激光测高数据辅助区域网平差的方法，其特征在于，步骤S2中所述的自由网平差包括以下步骤：
S21，利用卫星影像自带的RPC文件，构建有理函数模型，表达式如公式(3)所示：
Figure FDA0003400521630000031
其中，pi为一般多项式，表达式为：
pi＝a0+a1·Y+a2·X+a3·Z+a4·Y·X+a5·Y·Z+a6·X·Z+a7·Y2+a8·X2+a9·Z2+a10·X·Y·Z+a11·Y3+a12·Y·X2+a13·Y·Z2+a14·Y2·X+a15·X3+a16·X·Z2+a17·Y2·Z+a18·X2·Z+a19·Z3
其中，(X,Y,Z)为正则化后的地面点三维坐标，(x,y)为对应的正则化像点坐标，ai(i＝0,1,...,18,19)为有理函数系数；
S22，采用仿射变换模型对公式(3)进行几何定位误差补偿，表达式如公式(4)所示：
Figure FDA0003400521630000032
其中，(Δx,Δy)为像方改正数，(x,y)为利用公式(3)求得的像点坐标，mj,nj(j＝0,1,2)为仿射变换模型参数；
S23，影像仿射变换参数和连接点物方坐标是要求解的两类未知量，将公式(3)-(4)线性化，可得误差方程，如公式(5)所示：
V1＝A1X+B1Y1-l1,P1 (5)
式中，X＝[Δm0 Δm1 Δm2 Δn0 Δn1 Δn2]T，是影像的仿射变换参数改正数，Y1＝[Δlon Δlat Δh]T是连接点物方坐标改正值，A1,B1为对应的系数矩阵，l1为初值计算的常量，P1为权值。
4.根据权利要求1所述的利用高分七号激光测高数据辅助区域网平差的方法，其特征在于，步骤S4中所述的联合区域网平差具体包括：
S41，利用卫星影像自带的RPC文件，构建有理函数模型，表达式如公式(6)所示：
Figure FDA0003400521630000041
其中，pi为一般多项式，表达式为：
pi＝a0+a1·Y+a2·X+a3·Z+a4·Y·X+a5·Y·Z+a6·X·Z+a7·Y2+a8·X2+a9·Z2+a10·X·Y·Z+a11·Y3+a12·Y·X2+a13·Y·Z2+a14·Y2·X+a15·X3+a16·X·Z2+a17·Y2·Z+a18·X2·Z+a19·Z3
其中，(X,Y,Z)为正则化后的地面点三维坐标，(x,y)为对应的正则化像点坐标，ai(i＝0,1,...,18,19)为有理函数系数；
S42，采用仿射变换模型对公式(6)进行几何定位误差补偿，表达式如公式(7)所示：
Figure FDA0003400521630000042
其中，(Δx,Δy)为像方改正数，(x,y)为利用公式(6)求得的像点坐标，mj,nj(j＝0,1,2)为仿射变换模型参数；
S43，影像仿射变换参数和连接点物方坐标是要求解的两类未知量，将公式(6)-(7)线性化，可得误差方程，如公式(8)所示：
V1＝A1X+B1Y1-l1,P1 (8)
式中，X＝[Δm0 Δm1 Δm2 Δn0 Δn1 Δn2]T，是影像的仿射变换参数改正数，Y1＝[Δlon Δlat Δh]T是连接点物方坐标改正值，A1,B1为对应的系数矩阵，l1为初值计算的常量，P1为权值；
S44，对于激光测高点，其高程作为控制，在平差中仅改正其平面坐标，此时，误差方程变为：
V2＝B2Y2-l2,P2 (9)
式中，Y2＝[Δlon Δlat]T是激光测高点物方坐标改正值，B2是相应的系数矩阵，l2是初值计算常量，P2是相应权值。
5.根据权利要求1所述的利用高分七号激光测高数据辅助区域网平差的方法，其特征在于，步骤S4还包括：当存在其他外业控制点时，开展高程控制点、外业控制点与立体影像的联合区域网平差，实现立体影像高程和平面精度同步提升,此时还需要将外业控制点参与联合区域网平差，所涉及到的误差方程为：
V3＝A2X-l3,P3 (10)
式中，X＝[Δm0 Δm1 Δm2 Δn0 Δn1 Δn2]T是影像的仿射变换参数改正数，A2是相应的系数矩阵，l3是初值计算常量，P3是相应权值。
6.根据权利要求1所述的利用高分七号激光测高数据辅助区域网平差的方法，其特征在于，步骤S2中的立体影像包括采用高分七号卫星获取的立体影像。
CN202110795596.9A 2021-07-14 2021-07-14 一种利用高分七号激光测高数据辅助区域网平差的方法 Active CN113532377B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110795596.9A CN113532377B (zh) 2021-07-14 2021-07-14 一种利用高分七号激光测高数据辅助区域网平差的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110795596.9A CN113532377B (zh) 2021-07-14 2021-07-14 一种利用高分七号激光测高数据辅助区域网平差的方法

Publications (2)

Publication Number Publication Date
CN113532377A CN113532377A (zh) 2021-10-22
CN113532377B true CN113532377B (zh) 2022-02-22

Family

ID=78127907

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110795596.9A Active CN113532377B (zh) 2021-07-14 2021-07-14 一种利用高分七号激光测高数据辅助区域网平差的方法

Country Status (1)

Country Link
CN (1) CN113532377B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114137567B (zh) * 2021-11-26 2024-07-09 自然资源部国土卫星遥感应用中心 一种基于卫星遥感数据的冰川高程变化监测方法
CN117647232B (zh) * 2024-01-29 2024-04-16 航天宏图信息技术股份有限公司 激光高程点转刺到卫星立体影像的方法、装置及设备

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104931022B (zh) * 2015-04-21 2018-03-16 国家测绘地理信息局卫星测绘应用中心 基于星载激光测高数据的卫星影像立体区域网平差方法
US10053218B2 (en) * 2016-03-31 2018-08-21 General Electric Company System and method for positioning an unmanned aerial vehicle
CN106960174B (zh) * 2017-02-06 2021-06-08 中国测绘科学研究院 高分影像激光雷达高程控制点提取及其辅助定位方法
CN107504981B (zh) * 2017-07-25 2020-01-03 国家测绘地理信息局卫星测绘应用中心 一种基于激光测高数据的卫星姿态误差修正方法及设备
CN111174753B (zh) * 2019-12-28 2021-05-11 同济大学 基于有理函数模型的光学影像与激光测高数据平差方法

Also Published As

Publication number Publication date
CN113532377A (zh) 2021-10-22

Similar Documents

Publication Publication Date Title
CN113532377B (zh) 一种利用高分七号激光测高数据辅助区域网平差的方法
CN104931022B (zh) 基于星载激光测高数据的卫星影像立体区域网平差方法
KR101965965B1 (ko) 위성영상과 제공 rpc로부터 제작된 수치표고모델의 자동 오차보정 방법
CN102506824B (zh) 一种城市低空无人机系统生成数字正射影像图的方法
Li et al. Rigorous photogrammetric processing of HiRISE stereo imagery for Mars topographic mapping
Pan et al. Basic Products of the ZiYuan-3 Satellite and Accuracy Evaluation.
CN102735216B (zh) Ccd立体相机三线阵影像数据平差处理方法
CN113899387B (zh) 一种基于验后补偿的光学卫星遥感影像区域网平差方法及系统
CN108226982B (zh) 单线阵卫星激光联合高精度定位处理方法
CN113538595B (zh) 利用激光测高数据辅助提升遥感立体影像几何精度的方法
CN111473802A (zh) 一种基于线阵推扫的光学传感器内方位元素定标方法
CN102855628A (zh) 多源多时相高分辨率卫星遥感影像自动匹配方法
Chrysoulakis et al. Validation of ASTER GDEM for the Area of Greece
Ganas et al. An investigation into the spatial accuracy of the IKONOS 2 orthoimagery within an urban environment
CN108489468A (zh) 三线阵影像外方位元素平滑方程自适应光束法平差方法
CN114998399B (zh) 一种异源光学遥感卫星影像立体像对预处理方法
CN108447100A (zh) 一种机载三线阵ccd相机的偏心矢量和视轴偏心角标定方法
CN113238228B (zh) 基于水准约束的三维地表形变获取方法、系统及装置
CN110030968A (zh) 一种基于星载立体光学影像的地面遮挡物仰角测量方法
CN105571598A (zh) 一种卫星激光高度计足印相机姿态的测定方法
CN116045920A (zh) 一种基于gf7号立体像对提取dem的方法
CN115346128A (zh) 一种光学立体卫星dem高程改正和融合方法
Zhou et al. DEM generation using Ziyuan-3 mapping satellite imagery without ground control points
CN113624207A (zh) 一种基于轨道约束的卫星影像立体区域网平差方法
CN113066181A (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