[发明专利]一种工测数据与三维数字地球高程拟合的方法

摘要

本发明提供了一种工测数据与三维数字地球高程拟合的方法，包括步骤101，工测数据相对坐标与工测控制点投影坐标相互修正；步骤102，根据工测投影坐标和三维数字地球坐标计算七参数；步骤103，使用七参数将工测投影坐标中每个工测点的坐标转换到目标数字地球的坐标系统下获得新的工测投影坐标，以新的工测投影坐标进行地形拟合，将工测数据高程调整到三维数字地球高程附近；步骤104，以新的工测投影坐标修正三维数字地球高程坐标。本发明可以在确保工测数据工测点相对位置关系不变的前提下将工测数据和由此产生的输电线路设计成果提交到三维数字地球进行显示，实现工测数据与三维数字地球的融合。

主权项

一种工测数据与三维数字地球高程拟合的方法，其特征在于，包括：步骤101，工测数据相对坐标与工测控制点投影坐标相互修正；步骤102，根据工测投影坐标和三维数字地球坐标计算七参数；步骤103，使用七参数将工测投影坐标中每个工测点的坐标转换到目标数字地球的坐标系统下获得新的工测投影坐标，以新的工测投影坐标进行地形拟合，将工测数据高程调整到三维数字地球高程附近；步骤104，以新的工测投影坐标修正三维数字地球高程坐标；所述步骤101，进一步包括：获得工测点相对坐标；判断工测点相对坐标与工测控制点投影坐标是否吻合；如果不吻合，则根据工测点相对坐标推算新的投影坐标；利用平差法修正新的投影坐标；所述获得工测相对坐标，包括：已知有n个桩分别为Z1、Z2、……、Zn；档距分别为l1、l2、……、ln‑1；转角分别为α1、α2、……、αn；高程分别为h1、h2、……、hn；测得桩的投影坐标分别为(x1,y1,z1)、(x2,y2,z2)、……、(xn,yn,zn)；所述根据工测点相对坐标推算新的投影坐标，利用平差法修正新的投影坐标，包括：步骤1：生成中间坐标：中间坐标记为(x′1,y′1,z′1)、(x′2,y′2,z′2)、……、(x′n,y′n,z′n)，则生成方法如下：记ΔH＝z1‑h1，中间坐标递推公式如下：(x1′,y1′)=(x1,y1)z1′=h1+ΔH=z1(u1,v1)=(x2-x1,y2-y1)l]]>(x2′,y2′)=(x1′,y1′)+(u1,v1)·l1z2′=h2+ΔH(u2,v2)=(u1,v1)·cos(α2)+(-v1,u1)·sin(α2)......(xn′,yn′)=(xn-1′,yn-1′)+(un-1,vn-1)·ln-1zn′=hn+ΔH(un,vn)=(un-1,vn-1)·cos(αn)+(-vn-1,un-1)·sin(αn)]]>其中，ΔH为第一个桩投影坐标高程与工测坐标高程之差，ui,vi为中间坐标下第i档的单位方向向量即第i个桩到第i+1个桩的单位方向向量，u1,v1由投影坐标确定，第i档的单位方向向量由第i‑1档的单位方向向量及第i档起始桩转角确定；步骤2、进行平差处理生成新投影坐标：新投影坐标记为(x″1,y″1,z″1)、(x″2,y″2,z″2)、……、(x″n,y″n,z″n)，则新投影坐标与中间坐标的转换公式如下：xi′′yi′′=cos(θ)-sin(θ)sin(θ)cos(θ)xi′yi′+ΔxΔyzi′′=zi′+Δh]]>其中i∈{1,2,...,n}，θ、Δx、Δy、Δh为待求参数；为了求解θ、Δx、Δy、Δh，需使的以下两个函数值最小：f(θ,Δx,Δy)=Σi=1n[(xi′′-xi)2+(yi′′-yi)2]g(Δh)=Σi=1n(zi′′-zi)2]]>使用最小二乘法求解θ、Δx、Δy、Δh，其中θ是坐标系旋转角度，Δx、Δy、Δh为坐标的平移值；所述步骤102根据工测投影坐标和三维数字地球坐标计算七参数，进一步包括：工测投影坐标中工测点的WGS84经纬度高程为(B1,L1,H1)、(B2,L2,H2)、……、(Bn,Ln,Hn)；把测得的WGS84经纬度高程转换为WGS84椭球下的空间直角坐标(a1,b1,c1)、(a2,b2,c2)、……、(an,bn,cn)；把新投影坐标转换为其所对应椭球下的空间直角坐标(a′1,b′1,c′1)、(a′2,b′2,c′2)、……、(a′n,b′n,c′n)；七参数转换公式为：其中i∈{1,2,...,n}，Dx、Dy、Dz、Rx、Ry、Rz、m为待求参数；为了求解Dx、Dy、Dz、Rx、Ry、Rz、m，需使的以下函数值最小：使用最小二乘法求解七参数，包括：Dx、Dy、Dz、Rx、Ry、Rz、m。