城市独立坐标系下似中央子午线的确定王传江1，姚文强1，余美义1，王解先2（1.上海市测绘院，上海200063;2.同济大学，上海200092）摘要城市独立坐标系中央子午线等参数各地不一样，一般不公开，在进行GPS网平差时需当地的中央子午线及北、东方向的加常数等参数，本文就这些问题，提出了一种求解得到近似的城市独立坐标系下中央子午线的方法，通过实例计算，GPS网平差根据似中央子午线等参数得到的结果与根据实际已知的中央子午线等参数计算的结果一致。关键词GPS网平差中央子午线独立坐标系CalculationOfApproximateCentralMeridianUnderLocalIndependentCoordinateSystemWangChuan激ang1YaoWenqiang1YuMeiyi1Wang激exian2(1ShanghaiMunicipalInstituteofSurveyingandMapping;2Tong激University)一、前言随着GPS定位精度的不断提高，GPS技术在测量中的应用也越来越广泛。由于GPS观测解算得到的定位结果表示在WGS84全球坐标系中。我国测绘成果普遍表示在北京54坐标系中，它以克拉索夫斯基椭球为参考椭球，投影方式为Gauss投影，以3度或6度带划分整个中国所在区域。我国大部分城市为了避免Gauss投影变形带来的不便，而采用地方独立坐标系。地方独立坐标系的建立仍采用克拉索夫斯基椭球，中央子午线定在城市中央，投影面定为城市平均高度。这些原因使得我国的平面坐标较复杂，而城市中央子午线等参数各地不一样，一般不公开，在进行GPS网平差前时，需要知道当地的中央子午线及北、东向的加常数等信息。本文提出一种求解得到城市独立坐标系下近似中央子午线（似中央子午线）的方法，从而完成WGS84成果到城市独立坐标成果的转换，对工程施工应用有参考意义。二、数学模型对GPS网，在WGS84椭球下进行平差计算，得到各点WGS84椭球下的大地经纬度和大地高，通过以下过程转换成平面坐标：1)根据WGS84的椭球参数，由（1）式将换算至空间直角坐标：（1）其中：2)根据54椭球的椭球参数，由（2）式通过迭代将换算至大地坐标形式：(2)3)给定一个中央子午线、投影面高程，由（3）式进行高斯投影，将投影为Gauss坐标：（3）式中的符号见[3]。若GPS测定的点中，已知点的平面坐标为，则这些点的平面坐标与已知坐标之间的关系：（4）其中：为北、东向加常数；为缩放尺度；为旋转矩阵，为旋转角。（4）式是平面转换的公式，它是线形模型，由于Gauss投影变形是非线性的，它的一次项与成正比，因此平面转换模型只适合50km范围内较小的工程使用。在给定的中央子午线L0时，对2个以上公共点，由最小二乘法计算，即可以求出四参数、及。因此对城市独立坐标系，在确定中央子午线后，北、东向加常数可按上述（4）式解算。一般独立坐标系的中央子午线位于城市的中部，在完成GPS测量后，以所在测区区域中部的平均子午线L’为中心，在区间[L’-,L’+]内，以1秒步长变化子午线L’值，对每一L’值，经过以上（1）～（4）步骤的计算，可得到一组、及四参数，故对应与多个似中央子午线L’，可以得到多组四参数、、，再根据得到的其中任一组中央子午线及加常数，由GPS网平差软件完成GPS网平差计算。三、数据验证与应用为了验证本文提出的方法，给定四点的WGS84经纬度坐标及大地高，如下表1所示，按给定不同中央子午线及加常数等参数计算得到相应的独立坐标系下的坐标，然后用本文给定的方法计算不同似中央子午线及相应的转换后坐标，具体见表2、表3：表14个点的大地坐标点名B（DD.MMSS）L（DD.MMSS）H（m）131.1000122.200050231.1000122.220050331.1500122.250050431.1500122.300050表2独立坐标和转换坐标比较1点名独立坐标系中央子午线()，北向加常数0m,东向加常数500Km根据本文方法计算得到的中央子午线（～）及坐标两种方法结果坐标差(m)（m）(m)（m）(cm)(cm)13449511.722499205.5053449511.722499205.5060.0-0.123449511.961502383.4853449511.960502383.4840.10.133458753.694507144.1883458753.697507144.189-0.3-0.143458762.080515082.1813458762.079515082.1800.10.1表3独立坐标和转换坐标比较2点名独立坐标系中央子午线()，北向加常数0m,东向加常数0m根据本文方法计算得到的中央子午线（～）及坐标两种方法结...