GPS高精度工程平面控制网的建立.docVIP

GPS高精度工程平面控制网的建立 【摘 要】 本文介绍了一种利用GPS技术建立高精度工程平面控制网的方法，既可以保留相对坐标的精度，又可以满足工程建设的要求。着重论述了如何将WGS-84中的大地坐标换算到高精度工程平面控制网的坐标。 【关键词】 GPS 高精度平面控制网 高斯-克吕格投影 GPS是一门先进的空间科学技术，利用GPS技术建立高精度工程平面控制网，采集的数据是WGS-84椭球下的大地坐标，在平面控制网中，要将其转换到平面坐标。对于建立高精度平面控制网，考虑到投影变形，一般采用独立坐标系，所有的工程几何量（角度、距离等），都在这个独立坐标系中计算得到，达到施工所需的精度要求。 1 步骤 按照如下步骤来建立该独立坐标系控制网：（1）在处理数据时，应先将网中心点的WGS-84系大地坐标作为固定点，进行无约束平差，检验GPS网本身的符合精度、以及有没有明显的系统差，并剔除有粗差的观测数据。（2）将WGS-84系大地坐标进行高斯-克吕格投影，选择区域中心经度为中央子午线，计算平面坐标。（3）为了减少投影变形，满足高精度工程建设的需要，可以将投影面设为区域的平均大地高程面上。 以上步骤可以减小投影变形，提高控制网的实际精度。 2 理论基础与数学模型 GPS解算出来的无约束成果为WGS-84坐标系统内的大地坐标，我们需要的是平面坐标，需要对大地坐标进行投影转换。 （1）GPS的大地坐标高斯投影面，应用以下公式进行坐标转换求得高斯平面坐标（X，Y）： 式中为椭球面上P点与中央子午线的经差，P点在中央子午线之东为正，在西为负。B为P点的大地纬度，X为由赤道到纬度为B的子午线弧长。当P点的（L，B）为已知时（由于中央子午线之经度L0是已知的，则即可算出），即可按式1计算P点投影后的高斯平面坐标（x，y）。 式1所表示的（x，y）和（L，B）的函数关系，即确定了式1中F1和F2的具体形式。 当 3.5°时，式1换算的精度为±0.1m，欲要换算精确至0.001m的坐标式，可将1式继续扩充，现直接写出如下： 由式2中看出：当B 0时，x 0，而y随而变，故知赤道投影为直线，且为y轴（横坐标轴）。又 0时，y 0，而x X，故知中央子午线投影为直线，且为x轴（纵坐标轴）。而中央子午线与赤道之交点投影后为平面坐标之原点。 在式2中，令 常数，则得子午线投影曲线以纬度B为参数的参数方程。由此方程可知，当B值增大时，y值减小，而x值增加。又因Cos（-B） Cos（B），故无论B为正或为负，y的值都不变。所以子午线的投影曲线弯向中央子午线并向两极收敛，同时还对称于中央子午线和赤道。同样，令B 常数，则得平行圈投影曲线为向上和向下弯曲的曲线，它也对称于中央子午线和赤道。此外，我们还可以知道，离开中央子午线愈远的子午线，投影后弯曲愈厉害，长度变形也愈大，而任一曲线投影后均有向中央子午线和两极弯曲的趋势。 （2）如果需要将投影面设在该区域平均大地高程面上，则需进行下列变换。 如图1所示，地面P地的WGS-84大地坐标为（B，L），P地点对参考椭球的法线为椭球面交于P点，P点的WGS-84大地坐标为（B，L），现假设有一为WGS-84参考椭球相似的椭球，其相应投影平面在测区的平均高程面处，我们称这个椭球为相似参考椭球，相似参考椭球面为参考椭球面平行，即过P地点对应两椭球的法线重合，P地在相似参考椭球面上P的大地坐标为（B，L），两椭球的，，，相等。 如图2所示Nm为工程独立坐标系原点处对应的椭球卯酉圈曲率半径， 所以Hm a即da Hm 然后可以有由GPS工程平面控制网中每个点在平均椭球即相似参考椭球的大地坐标（Bm，Lm）以所选择的中央子午线（n为网中点的总数）并应用上面所提供的高斯投影计算公式来求每个点在大地高高程面上的平面直角坐标（x’，y’）。 3 程序设计基本思路和流程 （1）运行流程如图3： （2）基本思路：窗体Form1用于计算主界面，Cmd1主要是保存数据，Cmd2主要是用来计算平面坐标，首先调用子程序hhm，输入经纬度，然后作出选择是否将投影面设在工程区的平均WGS-84大地高程面上，如果是，则计算出在工程区的平均WGS-84大地高程面上的坐标，如果否，就直接跳出这一部分程序直接计算平面坐标（x，y）。 模块Module来定义bs，cs是布尔型数据，文件名是字符型数据 最后调用子程序ABCDE计算弧长XX，再计算x，y坐标，最后输出x，y坐标。 （3）部分代码 Subhhm（） DimlAsDouble Diml01AsDouble l01 0 hm 0 Fori 1ToForm1.MSFl.Rows-1 将L转化为秒 l Val（Left（Form1.MSFl.TextMatrix（i，2），Len（Form1.MSFl.TextMat