城市独立坐标系下似中央子午线的确定(投稿).doc

城市独立坐标系下似中央子午线的确定王传江1 ，姚文强1 ，余美义1，王解先2（1. 上海市测绘院，上海200063; 2. 同济大学，上海200092）摘要 城市独立坐标系中央子午线等参数各地不一样，一般不公开，在进行GPS网平差时需当地的中央子午线及北、东方向的加常数等参数，本文就这些问题，提出了一种求解得到近似的城市独立坐标系下中央子午线的方法，通过实例计算，GPS网平差根据似中央子午线等参数得到的结果与根据实际已知的中央子午线等参数计算的结果一致。关键词 GPS网平差 中央子午线 独立坐标系Calculation Of Approximate Central Meridian Under Local Independent Coordinate SystemWang Chuanjiang1 Yao Wenqiang1 Yu Meiyi1 Wang Jiexian2(1 Shanghai Municipal Institute of Surveying and Mapping ; 2 Tongji 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：表1 4个点的大地坐标点名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根据本文方法计算得到的中央子午线（）及坐标两种方法结果坐标差(m)（m）(m)（m）(cm)(cm)13449518.420-11917.4273449518.420-11917.4260.0-0.123449515.310-8739.4453449515.309-8739.4460.10.133458752.017-3968.9933458752.019-3968.991-0.2-0.243458752.0173968.9933458752.0163968.9920.10.1以上所得到的参数为许多组数据，似中央子午线为一个区间范围，根据区间范围内的任一似中央子午线计算出的结果是一致的。为了进一步验证本文的方法，采用某城市独立坐标系越江通道工程的实测数据，由随机软件进行基线解算，在WGS84下进行无约束平差，得到各点的在WGS84下大地坐标及大地高，见表4：表4 某控制网的WGS84坐标点名B（DD.MMSSSSSSSS）L（DDD.MMSSSSSSSS）H(m)200331. *45590737121.*2560013117.0308201331. *41292590121.*3146910913.9908GC0131. *4173825215.3938GC0231. *38650476121. *2216804215.7072JY0331. *25091497121. *2060309017.7013JY0231. *41637900121. *4749394814.8904413431. *57129696121. *1093016015.8785428031. *5740127214.8159428331. *20999727121. *0971256315.6725NH0131. *08704294121. *1468746113.7487NH0231. *36722991121. *0083750617.9135NH0331. *29488206121. *5487045518.1314其中2003、2013、GC01、GC02四点有城市独立坐标系下的坐标，按本文的方法计算该城市独立坐标系的似中央子午线等参数，似中央子午线的范围为，取其中一组数值，用TGPPS网平差软件计算，其结果与BALNET网平差（中央子午线等参数已知，可作为真值）计算的结果比较见表5，最大的差别为3mm。表5 本文方法与已知中央子午线的结果比较点名BALNET网平差（中央子午线等参数已知）TGPPS网平差（计算得到的似中央子午线等参数）两种方法结果坐标差(m)（m）(m)（m）(mm)(mm)4134-76*.874304*.206-76*.876304*.2072-14280149*.618 78*.727149*.615 78*.7303-34283 4*.208240*.280 4*.211240*.281-3-1JY03264*.293290*.047264*.291290*.04423JY02177*.837234*.348177*.839234*.349-2-1NH01 37*.850146*.680 37*.851146*.683-1-3NH02101*.688190*.302101*.688190*.3040-2NH03154*.641220*.811154*.642220*.811-10四、结论对城市独立坐标系，在进行GPS作业时，若不