探讨RTK测量中中央经线的选择.docVIP

探讨RTK测量中中央经线的选择 　　摘要：RTK测量的坐标系统设置和参数计算过程中，椭球基准的选择和中央经线的设置对参数计算的的影响和作用。本文选择WGS-84椭球进行投影计算和地方坐标系的转换参数计算来说明影响和在应用中作用。 　　关键字：高斯投影中央经线参数计算 　　中图分类号：O221.8 文献标识码：A 文章编号： 　　1引言 　　RTK技术以其高精度、操作简单等优势在测绘工程领域得到广泛应用。但RTK测量的是WGS-84坐标系下的坐标，需要将其转换到当地坐标系下的坐标。WGS-84坐标系和地方坐标系转换有两种思路，空间转换模型和平面转换模型,由于地方坐标系椭球基准和投影高程面的不确定性，在局部区域测量时一般采用平面转换模型。在测量过程中，椭球基准和投影中央子午线的设置对坐标转换参数产生一定的影响。本文通过实际数据的解算来说明这一点，希望能对RTK的应用有所帮助。 　　2 高斯投影 　　高斯投影是一种横轴椭圆柱面正形投影。在RTK坐标转换过程中，是将WGS-84经纬度坐标，在设置的投影基准上进行投影计算，然后利用四参数进行平面的转换.而投影基准的选择上，不同的设置有不同的投影。例如Trimble在坐标系统设置时，一般情况下选择（1）无投影/无基准（2）键入参数两种设置，在选择无投影/无基准时椭球基准为WGS-84椭球；在键入参数时，椭球基准根据输入设置而定，可以选择北京54，也可以选择WGS-84椭球。下面以WGS-84椭球基准来进行高斯投影计算，公式如下： 　　x=X+NsinBcosB+NsinBcosB(5-t+9+4)+NsinBcosB(61- 58t+t)(1) 　　y=NcosB+NcosB(1-t+)+NcosB(5-18t+t+14-58t)(2) 　　X=c[+( +++)]=1-+-+ 　　=-1 　　=-+ 　　=-+ 　　= 　　 　　利用（1）、（2）计算高斯投影平面坐标x,y的计算精度可达0.001m.在公式中=L-L0,其中L表示WGS-84坐标系中某一点的经度，L0表示该点所在区域的中央经线的经度,B表示大地纬度，N表示卯酉圈曲率半径，t=tanB,表示第二偏心率，=cosB。 　　3 平面转换 　　在进行平面转换时，在工作选择不少于2对高精度的控制点坐标对，两个点最好能贯穿整个测区。转换模型，通常涉及2个平移参数，1个旋转参数和1个尺度缩放参数共计4个转换参数。现已知5个坐标对，选择3个求解4参数。 　　 （3） 　　（3）式中，（，）是坐标平移量，是缩放尺度比参数，是旋转参数，这就是待求的四参数。利用3个公共点（X,Y）可列出6个方程，误差方程式为： 　　 （4） 　　将3个公共点的坐标代入（4）式，应用最小二乘原理求出坐标转换参数。求得转换参数后，再利用上述模型进行坐标转换，转换后的坐标与已知值的残差大小反映转换参数的精度。 　　4 程序代码 　　利用Visual Basic 语言实现WGS-84大地坐标到平面坐标的转换，并根据已知对应点求解四参数，部分程序原代码如下： 　　4.1坐标转换 　　Public sub zh() 　　Dim a#, b1#, c#, n#, l#, b#, x#, y#, e1#, e2#, w#, dd#, mm#, ss# 　　a = 6378137: b1 = a * 297.257223563 / 298.257223563 　　dd = Fix(Val(Trim(Text1.Text))):mm = Fix((Val(Trim(Text1.Text)) - dd) * 100) 　　ss = ((Val(Trim(Text1.Text)) - dd) * 100 - mm) * 100:b = (dd + mm / 60 + ss / 3600) / 57.295779513 　　dd = Fix(Val(Trim(Text3.Text))):mm = Fix((Val(Trim(Text3.Text)) - dd) * 100) 　　ss = ((Val(Trim(Text3.Text)) - dd) * 100 - mm) * 100:zyjx = (dd + mm / 60 + ss / 3600) 　　dd = Fix(Val(Trim(Text2.Text))):mm = Fix((Val(Trim(Text2.Text)) - dd) * 100) 　　ss = ((Val(Trim(Text2.Text)) - dd) * 100 - mm) * 100:l = ((dd + mm / 60 + ss / 3600) - zyjx) / 57.2957795