解手持GPS参数实例.docVIP

解手持GPS参数程序和实例 作者：陈平 手持GPS一般内部设置五个转换参数，只要计算出五个参数（DX，DY，DZ，DA，DF），并在仪器中输入这五个参数即可完成了坐标转换工作。 1、文件名称：“84-54-80-XYZ” “Q=1=84,Q=2=54,Q=3=80”?Q If Q=1:Then 6378137-A:0.00669437999013-E: IfEnd If Q=2:Then 6378245-A:0.006693427-E: IfEnd If Q=3:Then 6378140-A:0.006694385-E: IfEnd “H”?H:”B”?B:”L”?L:A÷√(1-E(sin(B)2)-N:(N+H)cos(B)cos(L)-X: (N+H)cos(B)sin(L)-Y:(N(1-E)+H)sin(B)-Z:”X=”:X◢ ”Y=”:Y◢ ”Z=”:Z◢ 编写此文件的数学公式依据为下： 不同坐标系对应的地球椭球参数，表1： 坐标系名称 长半轴a（m） 第一偏心率平方e2 扁率f WGS-84 6378137 0. 00669437999013 1/ 298. 257223563 北京54 6378245 0. 006693427 1/ 298. 3 西安80 6378140 0. 006694385 1/ 298. 257 对于对于同一空间点,大地坐标系与空间直角坐标系有下列转换关系： N = A/ (1 – E2sinB2) 1/ 2 在“84-54-80-XYZ”程序中，A为相应坐标系里的a，E为相应坐标系里的e2 2、实例：三个高等级控制点坐标，表2：（下列坐标为北京54坐标系，y717里的38为带号，以3度带划分的。中心子午线经度为：38°*3=114°） 点名 x y h 001 2710395.717717 245.046 002 2709699.810500 62.428 003 2710797.460600 57.355 用GPS实测这三个控制点的WGS-84坐标系统下大地坐标（B，L）和大地高H，表3： 点名 B L H 001 24°29′50.38″ 113°44′18.94″ 248.249 002 24°29′26.98″ 113°37′59.10″ 65.648 003 24°30′02.59″ 113°37′33.64″ 60.520 用“QH1-4”程序（在覃辉编著的《CASIO fx-5800P公路与铁路施工测量程序》一书中可以查到）将表2里的1、2、3号点通过高斯反算（反算时去带号，比717就输入473387.717）成大地坐标（B，L），反算的结果如下，表4： 点名 B L H 001 24°29′48.7084″ 113°44′14.8053″ 245.046 002 24°29′25.2982″ 113°37′54.9763″ 62.428 003 24°30′0.9093″ 113°37′29.5165″ 57.355 将表3和表4里的B、L、H值分别输入“84-54-80-XYZ”程序里，可得出2组数据： 由表3里的B、L、H值算出的数据：（WGS-84坐标系） 点名 X2 Y2 Z2 001 -2337913.315 5316182.665 2628611.08 002 -2328172.572 5320600.196 2627880.177 003 -2327331.879 5320467.041 2628875.098 由表4里的B、L、H值算出的数据：（北京54坐标系） 点名 X1 Y1 Z1 001 -2337853.553 5316335.964 2628609.761 002 -2328112.854 5320753.365 2627878.553 003 -2327272.171 5320620.230 2628873.548 将上一步计算中所得到的成果X1、Y1、Z1和X2、Y2、Z2代入到以下函数模型中计算， 可解出3组DX、DY、DZ，取平均数。五参数中，DA、DF是常值，而DX、DY、DZ在不同的地区参数值都不一样。其中，相对于1954年北京坐标系而言，DA=-108m，DF=0.0000005；而对于1980年国家大地坐标系而言，DA=-3m，DF=-0.0000000025，DA为需要转换的两个坐标系之间椭球半长轴之差。DF为需要转换的两个坐标系之间椭球扁率之差。（所以，84相对于北京54：DA=6378137-6378245=-108，DF=1/ 298. 257223563-1/ 29