投影于抵偿高程面上的坐标计算方法及其公式推导

《投影于抵偿高程面上的坐标计算方法及其公式推导》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、投影于抵偿高程面上的坐标计算方法及其公式推导摘要：讨论、分析投影于（任意）抵偿高程面上的平面坐标计算方法及其计算公式的推导。供同行们讨论与参考。关键词：交通工程；公路控制测量；投影于（任意）抵偿高程面上的平面坐标计算方法。0前言国家有关规范规定，在大、中型工程测量中，其控制网必须与国家控制点联测，或采用国家坐标系统，以达到测量资源共享、成果共用的目的。国家坐标系统是采用高斯－克吕格正形投影（简称“高斯投影”），即先由大地面投影到参考椭球面，再由参考椭球面投影到高斯平面；而高程面则是投影到大地水准面上。公路测量常用的处理方法是，采用分带形式，以减小高斯投影产生的长度变形；而高程面的投影

2、，因为测区平均高程面与大地水准面的差值和地球曲率半径相比微不足道，故忽略不计。然而，随着公路建设的不断扩大与发展，公路（特别是高速公路）从平原微丘区向山岭重丘区（乃至高原地区）延伸，测区高程面由数10m增加到数百米乃至数千米；由于高程面的不同所产生的长度变形对工程建设的影响是我们必须考虑的问题。如，据有关计算表明，当大地高程面H＝700m时，其长度变形为11cm/km，远大于规范允许值，这对于重要工程的测量是一个不可忽略的小数。本文通过分析讨论，提出在（任意）选定的抵偿高程面上的平面坐标的计算方法来解决长度变形问题。1独立坐标系中投影于抵偿高程面上的坐标换算在独立坐标系中，原有坐标X

3、、Y投影高程面为H0，测区平均高程面为H，为使实测边长与成图平面上的边长相一致，不致产生过大的长度投影变形，需将测区平均高程面H作为抵偿高程面（简称投影面），从而建立新的地方独立坐标系统。利用原有坐标X、Y换算成新的投影面（抵偿高程面）上的独立坐标Xˊ、Yˊ，一般取测区中心或附近点为投影原点（X0、Y0），换算过程中不考虑椭球面正形投影到高斯平面上长度改化变形因素对坐标换算的影响，公式推导如下。如图1所示：图1R-投影区地球平均曲率半径H0-原坐标投影面高程H-新坐标投影面（抵偿高程面）高程X0、Y0-投影原点坐标X、Y-原坐标Xˊ、Yˊ-投影于抵偿高程面上的新坐标因为：(Xˊ-X0

4、)/（R+H）=（X-X0）/（R+H0）Xˊ=X0+（X-X0）（R+H）/（R+H0）所以：Xˊ=X0+(X-X0)〔1+（H-H0）/(R+H0)（1）同理：Yˊ=Y0+(Y-Y0)〔1+（H-H0）/(R+H0)〕以上（1）式即为新老坐标投影换算公式2投影于抵偿高程面上的高斯平面坐标换算将1954年北京坐标换算为投影于地方独立抵偿高程面上的高斯平面坐标，按以下两种方法考虑：2.1方法一该方法是以抵偿高程面作为地方独立参考椭球面，通过计算任一实测边长D0投影于地方独立参考椭球面(即抵偿高程面)并经高斯正形投影改化后，得到新的地方独立坐标系抵偿高程Hp上的高斯平面边长度D2ˊ，与

5、该实测边长在1954年北京坐标系统高斯平面上长度D2（实测边长D0投影于克氏参考椭球上经高斯正形投影改化后在高斯平面上的长度）的比值D2ˊ/D2关系式，推导出新旧投影面上的高斯平面坐标换算公式。如图2所示：图2D0-实测边长长度R-投影区地球平均曲率半径D1-克氏椭球面上长度D1ˊ-地方独立椭球面（抵偿高程面）上长度D2-1954年北京坐标系高斯平面上长度D2ˊ---地方独立坐标系高斯平面上长度（D1ˊ经高斯正形投影改化后得到的长度）Hm-实测边长两端点平均高程Hp-抵偿高程面高程（即投影面高程）Ym-实测边长两端点距1954年北京坐标系中央子午线之平均距离Ymˊ-实测边长两端点距地

6、方独立坐标系中央子午线之平均距离X、Y-1954年北京坐标Xˊ、Yˊ-投影于抵偿高程面上的高斯平面坐标因为：D1=D0[1-Hm/（R+Hm）]D1ˊ=D1(1+Hp/R)D2=D1(1+Ym2/2R2)-------高斯正形投影改化公式D2ˊ=D1ˊ[1+Ymˊ2/2(R+Hp)2]=D1(1+Hp/R)[1+Ymˊ2/2(R+Hp)2]---------高斯正形投影改化公式所以：任一边长投影于抵偿高程面上的高斯平面长度与其投影于1954年北京坐标系高斯平面上的长度之比=D2ˊ/D2=(1+Hp/R)[1+Ymˊ2/2(R+Hp)2]/[1+Ym2/2R2]由于长度投影变形与方向无

7、关，且坐标与长度成正比例线性关系。因此：Yˊ=Y(1+Hp/R)[1+Yˊ2/2(R+Hp)2]/[1+Y2/2R2]Xˊ=X(1+Hp/R)[1+Yˊ2/2(R+Hp)2]/[1+Y2/2R2]以上（2）式即为投影于抵偿高程面上的高斯平面坐标换算公式，其中Yˊ坐标计算可以通过解一元二次方程求得，或者取Yˊ=Y(1+Hp/R)式计算，误差较小。2.2方法二该方法是首先将原坐标X、Y（1954年北京坐标）进行高斯投影反算，求得大地坐标经纬度L、B。然后在考虑