电子表格进行高斯投影换算GPS坐标的方法
【GPS】Excel在GPS坐标转换计算中的应用
4
转换公式(4)
平面直角坐标系(X,Y,Z)
图 1:坐标转换流程图
1
各转换步骤对应的公式分别如下所示: (1) GPS定位得到的大地经纬度(B84,L84,H84)坐标通过公式[3](1)转换成空间直角坐标
(X84,Y84,Z84):
⎧X = (N + h)cosBcosL
则 GPS 坐标转换成 54 坐标只需将第一步算的直角坐标代入公式(2)。
(3) 根据 54 坐标系的椭球参数,将第二步所得的(X54,Y54,Z54)换算为大地坐标(B54,L54,H54), 详见公式[3](3):
⎧L = arctan(Y / X)
⎪⎪⎨B = arctan[(Z + Ne2 sin B) / X 2 + Y 2 ]
随着我国各地gps差分台站的不断建立以及美国sa政策的取消使得单机定位的精度大大提高有的已经达到了分米级精度能够满足国土资源调查土地利用更新遥感监测海域使用权清查等工作的应用由于我国测绘成果使用的是1954年北京坐标系或1980年西安坐标系以下分别简称54系和80系而gps测定的坐标是wgs84坐标系坐标需要进行坐标系转换
Excel 在 GPS 坐标转换计算中的应用
冯锦霞,陈立 (中南大学 国土与测绘工程系,湖南 长沙 410083)
摘要:对利用 Excel 电子表格进行 WGS84 坐标系与 54 坐标系转换并进行高斯投影换算的方 法进行了介绍,给出了在 Excel 中完成整个过程的简单方法,实现了用 Excel 完成坐标转换 并进行高斯投影。 关键字: Excel VBA 坐标变换 高斯投影 GPS 坐标
对利用EXCEL电子表格进行高斯投影换算的方法进行了较详细的介绍资料
对利用EXCEL电子表格进行高斯投影换算的方法进行了较详细的介绍,对如何进行GPS坐标系转换进行了分析,提出了一种简单实用的坐标改正转换方法,介绍了用EXCEL完成转换的思路。
[关键字] 电子表格;GPS;坐标转换作为尖端技术GPS,能方便快捷性地测定出点位坐标,无论是操作上还是精度上,比全站仪等其他常规测量设备有明显的优越性。
随着我国各地GPS差分台站的不断建立以及美国SA政策的取消,使得单机定位的精度大大提高,有的已经达到了亚米级精度,能够满足国土资源调查、土地利用更新、遥感监测、海域使用权清查等工作的应用。
在一般情况下,我们使用的是1954年北京坐标系或1980年西安坐标系(以下分别简称54系和80系),而GPS测定的坐标是WGS-84坐标系坐标,需要进行坐标系转换。
对于非测量专业的工作人员来说,虽然GPS定位操作非常容易,但坐标转换则难以掌握,EXCEL是比较普及的电子表格软件,能够处理较复杂的数学运算,用它来进行GPS坐标转换、面积计算会非常轻松自如。
要进行坐标系转换,离不开高斯投影换算,下面分别介绍用EXCEL进行换算的方法和GPS坐标转换方法。
一、用EXCEL进行高斯投影换算从经纬度BL换算到高斯平面直角坐标XY(高斯投影正算),或从XY换算成BL(高斯投影反算),一般需要专用计算机软件完成,在目前流行的换算软件中,存在一个共同的不足之处,就是灵活性较差,大都需要一个点一个点地进行,不能成批量地完成,给实际工作带来许多不便。
笔者发现,用EXCEL可以很直观、方便地完成坐标换算工作,不需要编制任何软件,只需要在EXCEL的相应单元格中输入相应的公式即可。
下面以54系为例,介绍具体的计算方法。
完成经纬度BL到平面直角坐标XY的换算,在EXCEL中大约需要占用21列,当然读者可以通过简化计算公式或考虑直观性,适当增加或减少所占列数。
在EXCEL中,输入公式的起始单元格不同,则反映出来的公式不同,以公式从第2行第1列(A2格)为起始单元格为例,各单元格的公式如下:单元格单元格内容说明A2输入中央子午线,以度.分秒形式输入,如115度30分则输入115.30 起算数据L0B2=INT(A2)+(INT(A2*100)-INT(A2)*100)/60+(A2*10000-INT(A2*100)* 100)/3600把L0化成度C2以度小数形式输入纬度值,如38°14′20″则输入38.1420起算数据BD2以度小数形式输入经度值起算数据LE2=INT(C2)+(INT(C2*100)-INT(C2)*100)/60+(C2*10000-INT(C2*100)* 100)/3600把B化成度F2=INT(D2)+(INT(D2*100)-INT(D2)*100)/60+(D2*10000-INT(D2*100)* 100)/3600把L化成度G2=F2-B2L-L0H2=G2/57.2957795130823 化作弧度I2=TAN(RADIANS(E2)) Tan(B)J2=COS(RADIANS(E2)) COS(B)K2=0.006738525415*J2*J2 L2=I2*I2M2=1+K2N2=6399698.9018/SQRT(M2)O2=H2*H2*J2*J2P2=I2*J2Q2=P2*P2R2=(32005.78006+Q2*(133.92133+Q2*0.7031))S2=6367558.49686*E2/57.29577951308-P2*J2*R2+((((L2-58)*L2+61)*O2/30+(4*K2+5)*M2-L2)*O2/12+1)*N2*I2*O2/2计算结果XT2=((((L2-18)*L2-(58*L2-14)*K2+5)*O2/20+M2-L2)*O2/6+1)*N2*(H2*J 2)计算结果Y表中公式的来源及EXCEL软件的操作方法,请参阅有关资料,这里不再赘述。
excel在高斯投影坐标正算公式中的应用
excel在高斯投影坐标正算公式中的应用
Excel在改变世界:应用高斯投影正算公式
Excel在高斯投影坐标正算公式中的应用
Excel在当今的商务统计、财务会计、商业报表处理等工作中已经成为一种必不可少的办公软件。
在空间信息科学领域,它也正在发挥其重要作用。
例如,它在高斯投影正算公式中的应用。
一、高斯投影正算的概念
高斯投影正算是根据高斯投影原理,利用投影面内的算法和参数,将地球大地坐标转换为投影坐标,或者将投影坐标转换为大地坐标。
二、Excel应用详解
1. 填写Excel工作表
使用Excel进行高斯投影正算时,要先将原始数据(如:地球大地坐标经度、纬度)填入Excel工作表,以便更方便快捷地进行计算。
2. 高斯投影系数求取
接下来,使用高斯投影坐标转换功能,根据起始的大地坐标经度、纬度,提取出系数A、B。
3. 运用公式实现高斯投影正算
有了其他要素的求取,运用Excel的高斯投影正算函数就可以直接完成投影面原点坐标转换到地球大地坐标或者投影面大地坐标转换到地球大地坐标的任务。
可以使用函数=PRJG2BL(A,B,X,Y),其中 A 、B 代表高斯投影的参数,X 、Y 代表地球大地坐标的经纬度,该函数的返回值为(经度,维度),同样可以使用函数=PRJB2LL(A,B,X,Y),输入该函数的参数 X 、Y 就是投影坐标。
三、总结
由以上,可以看出,Excel在高斯投影正算中的应用十分简明:充分利用Excel提供的工具,可以快捷有效地完成地理坐标之间的转换工作,
极大地提升了地理数据处理流程的效率和合理性,极大地为地理数据处理任务提供了可靠的基础支持。
在Excel中实现统一椭球的高斯投影正反算
在Excel中实现统一椭球的高斯投影正反算
覃军;骆奇峰
【期刊名称】《测绘与空间地理信息》
【年(卷),期】2009(32)2
【摘要】<控制测量学>介绍了对于不同的椭球都必须输入不同系数进行高斯投影正反算的方法,由于正反算公式的系数多、系数数字位长,这种方法显得很繁琐.目前,计算机应用非常广泛,可以推导出不同椭球的统一计算公式.在Excel中输入不同的椭球参数和统一的计算公式即可进行高斯正反算和换带计算.另外,笔者在公式的推导过程中和算例重算中发现<控制测量学>(第三版)中存在一些印刷排版错误,在此将一一列出,供同行参考.
【总页数】4页(P205-208)
【作者】覃军;骆奇峰
【作者单位】湖南省张家界市国土资源局永定分局,湖南,张家界,427000;广东省国土资源厅测绘院,广东,广州,510500
【正文语种】中文
【中图分类】P226
【相关文献】
1.Excel在高斯投影坐标正算公式中的应用 [J], 崔学敏;赵雪松;姚继华
2.椭球变换后的高斯投影正反算算法 [J], 邓兴升;汤仲安;花向红;舒玉芝
3.高斯投影正、反算及换带计算程序 [J], 孔达;孟庆武;赵力娟;;;
4.利用高斯投影地图反算图斑椭球面积的误差 [J], WANG Huiran;GAO Ziqiang
5.一种不同参考椭球下高斯投影地形图坐标转换及实现 [J], 王会然
因版权原因,仅展示原文概要,查看原文内容请购买。
用EXCLE完成坐标转化的简易方法
用EXCEL完成GPS坐标转换的简易方法EXCEL, GPS, 坐标用EXCEL完成GPS坐标转换的简易方法[摘要] 对利用EXCEL电子表格进行高斯投影换算的方法进行了较详细的介绍,对如何进行GPS坐标系转换进行了分析,提出了一种简单实用的坐标改正转换方法,介绍了用EXCEL完成转换的思路。
[关键字] 电子表格;GPS;坐标转换作为尖端技术GPS,能方便快捷性地测定出点位坐标,无论是操作上还是精度上,比全站仪等其他常规测量设备有明显的优越性。
随着我国各地GPS差分台站的不断建立以及美国SA政策的取消,使得单机定位的精度大大提高,有的已经达到了亚米级精度,能够满足国土资源调查、土地利用更新、遥感监测、海域使用权清查等工作的应用。
在一般情况下,我们使用的是1954年北京坐标系或1980年西安坐标系(以下分别简称54系和80系),而GPS测定的坐标是WGS-84坐标系坐标,需要进行坐标系转换。
对于非测量专业的工作人员来说,虽然GPS定位操作非常容易,但坐标转换则难以掌握,EXCEL是比较普及的电子表格软件,能够处理较复杂的数学运算,用它来进行GPS坐标转换、面积计算会非常轻松自如。
要进行坐标系转换,离不开高斯投影换算,下面分别介绍用EXCEL进行换算的方法和GPS坐标转换方法。
一、用EXCEL进行高斯投影换算从经纬度BL换算到高斯平面直角坐标XY(高斯投影正算),或从XY换算成BL(高斯投影反算),一般需要专用计算机软件完成,在目前流行的换算软件中,存在一个共同的不足之处,就是灵活性较差,大都需要一个点一个点地进行,不能成批量地完成,给实际工作带来许多不便。
笔者发现,用EXCEL可以很直观、方便地完成坐标换算工作,不需要编制任何软件,只需要在EXCEL的相应单元格中输入相应的公式即可。
下面以54系为例,介绍具体的计算方法。
完成经纬度BL到平面直角坐标XY的换算,在EXCEL中大约需要占用21列,当然读者可以通过简化计算公式或考虑直观性,适当增加或减少所占列数。
工程测量中高斯-克吕格投影换带计算的应用
工程测量中高斯-克吕格投影换带计算的应用在测量工程中,有些测区刚好处于投影带边缘,甚至有些工程横跨两个或两个以上投影带,如交通、水利、电力等较长的线路,为了坐标统一的需要,可以进行坐标换带,将相邻带的坐标换成同一系统的数据。
坐标换带有直接换带计算法和间接换带计算法两种,间接换带计算法就是根据第一带的平面坐标某1,y1和中央子午线的经度L1,按高斯-克吕格投影坐标反算公式求得大地坐标B、L,然后根据B,L和第二带的中央子午线经度L2,按高斯-克吕格投影坐标正算公式求得在第二带中的平面坐标某2、y2。
由于在换带计算中,把椭球面上的大地坐标作为过渡坐标,因而称为间接换带法。
这种方法理论上严密,精度高,而且通用性强,虽然计算量较大,但可用电子计算机计算来克服,已成为坐标换带中最根本的方法。
2、换带计算公式用a表示椭球长半轴,b表示椭球短半轴,f=为扁率,e=为第一偏心率,eˊ=为第二偏心率,N=为卯酉圈曲率半径,R=为子午圈曲率半径,B表示经度,L表示纬度。
2.1高斯-克吕格投影反算公式:B=Bf-[-(5+3Tf+Cf-9TfCf)+(61+90Tf+45Tf2)]L=L0+[D-(1+2Tf+Cf)+(5+28Tf+6Cf+8TfCf+24Tf2)]式中:Nf==,Rf=Bf=+(-)in2+(-)in4+in6,e1==,Tf=tg2Bf,Cf=eˊ2co2Bf,D=2.2高斯-克吕格投影正算公式:某N=k0{M+NtgB[+(5-T+9C+4C2)]+(61-58T+T2+270C-330TC)YE=FE+k0N[A+(1-T+C)+(5-18T+T2+14C-58TC)]式中:k0=1,T=tg2B,C=eˊ2co2B,A=(L-L0)coB,N==M=a[(1---)B-(--)in2B+(-)in4B-in6B东西偏移量FE=500000米+带号某10000003.1进行换带计算的步骤分析通常建立独立坐标系的方法是以一个国家控制点和方位角作为起算数据,观测边投影到平均高程面。
EXCEL公式进行经纬度与XY坐标的相互转换
EXCEL公式进行经纬度与XY坐标的相互转换•一、用EXCEL进行高斯投影换算从经纬度B、L换算到高斯平面直角坐标X、Y(高斯投影正算),或从X、Y换算成B、L(高斯投影反算),一般需要专用计算机软件完成。
在目前流行的换算软件中不足之处,就是灵活性较差,大都需要一个点一个点地进行,不能成批量地完成,给实际工作带来许多不便。
而用EXCEL可以很直观、方便地完成坐标换算工作,不需要编制任何软件,只需要在EXCEL的相应单元格中输入相应的公式即可。
下面以1954年北京坐标系为例,介绍具体的计算方法。
上图为编辑好的EXCEL表(红色为输入数据项)完成经纬度B、L到平面直角坐标X、Y的换算,在EXCEL中大约需要占用21列,当然读者可以通过简化计算公式或考虑直观性,适当增加或减少所占列数。
在EXCEL中以公式从第3行第1列(A3格)为起始单元格为例,各单元格的公式如下:(1)单元格A3输入中央子午线,以度、分、秒形式输入,如107度0分则输入107.00 (2)单元格B3公式如上图,把L0化成度形式。
(3)单元格C3以度小数形式输入纬度值,如23°44′01″则输入23.4401。
(4)单元格D3以度小数形式输入经度值,如107°42′48″则输入107.4248。
(5)单元格E3公式如上图,把纬度B化成度形式。
(6)单元格F3公式如上图,把经度L化成度形式。
(7)各个单元格输入公式如下:表中计算公式见由孔祥元等主编、武汉大学2002年出版的《控制测量学》,EXCEL软件的操作方法请参阅有关资料。
按上面表格中的公式输入到相应单元格后,就可方便地由经纬度求得平面直角坐标。
当输入完所有的经纬度后,用鼠标下拉即可得到所有的计算结果。
表中的许多单元格公式为中间过程,可以用EXCEL的列隐藏功能把这些没有必要显示的列隐藏起来,表面上形成标准的计算报表,使整个计算表简单明了,可计算的数据量是无限的,当第一次输入公式后,相当于自己完成了一软件的编制,可存储起来供今后重复使用。
利用EXCEL行高斯投影正反算
利用EXCEL行高斯投影正反算————————————————————————————————作者:————————————————————————————————日期:利用EXCEL 进行高斯投影正反算在工作中常需要将大量经纬度转换成高斯平面坐标、将高斯平面坐标转换成经纬度、将6度带坐标转换成3度带坐标等坐标转换问题。
面对这些问题,我们希望能找个坐标转换软件进行批量转换从而降低劳动强度、提高工作效率。
然而我们通常用的软件对批量转换往往有一定的限制而且对转换数据的格式要求比较严格不容易掌握和使用。
实际上我们通常用的办公软件EXCEL 就可以完成这项工作。
EXCEL 办公软件操作简单方便、易于掌握。
想要用EXCEL 实现正反算,我们必须知道高斯投影正反算数学公式。
(高斯投影正算实际就是把大地坐标通过高斯投影数学模型转换为平面坐标,反之则为高斯投影反算)下面就是高斯投影公式:正算公式:"2322"4""4sin cos sin cos (59)22N N x X B Bl B B t l ηρρ=++-+ "322"3524"5""3"5cos cos (1)cos (118)6120NN N y X Bl B t l B t t l ηρρρ=++-++-+反算公式:23222424635(539)(619045)224720f f f f f f f f f f ff f f ft t N B B y t t y t t y M N M N M N ηη=-+-+---+ 22324535111(12)(52824)cos 6cos 120cos f f f f f f f f f fl y t y t t y N B N B N B η=-+++++ 我们只需要把上述两个公式用EXCEL 函数写于EXCEL 表格中就可以进行计算了。
利用EXCEL进行高斯投影正反算
利用EXCEL 进行高斯投影正反算在工作中常需要将大量经纬度转换成高斯平面坐标、将高斯平面坐标转换成经纬度、将6度带坐标转换成3度带坐标等坐标转换问题。
面对这些问题,我们希望能找个坐标转换软件进行批量转换从而降低劳动强度、提高工作效率。
然而我们通常用的软件对批量转换往往有一定的限制而且对转换数据的格式要求比较严格不容易掌握和使用。
实际上我们通常用的办公软件EXCEL 就可以完成这项工作。
EXCEL 办公软件操作简单方便、易于掌握。
想要用EXCEL 实现正反算,我们必须知道高斯投影正反算数学公式。
(高斯投影正算实际就是把大地坐标通过高斯投影数学模型转换为平面坐标,反之则为高斯投影反算)下面就是高斯投影公式:正算公式:"2322"4""4sin cos sin cos (59)22N N x X B Bl B B t l ηρρ=++-+ "322"3524"5""3"5cos cos (1)cos (118)6120NN N y X Bl B t l B t t l ηρρρ=++-++-+ 反算公式:23222424635(539)(619045)224720f f f f f f f f f f ff f f ft t N B B y t t y t t y M N M N M N ηη=-+-+---+ 22324535111(12)(52824)cos 6cos 120cos f f f f f f f f f fl y t y t t y N B N B N B η=-+++++ 我们只需要把上述两个公式用EXCEL 函数写于EXCEL 表格中就可以进行计算了。
其实,我们用的其他软件中高斯正反算基本全都是用这两个公式实现的,然而由于各种程序语言以及软件计算数据时精度不同所以最后结果也会有很小的误差。
GPS经纬坐标到直角坐标系的转换
GPS经纬坐标到直角坐标系的转换GPS经纬坐标到直角坐标系的转换用EXCEL完成GPS坐标转换的简易方法[摘要] 对利用EXCEL电子表格进行高斯投影换算的方法进行了较详细的介绍,对如何进行GPS坐标系转换进行了分析,提出了一种简单实用的坐标改正转换方法,介绍了用EXCEL完成转换的思路。
[关键字] 电子表格;GPS;坐标转换作为尖端技术GPS,能方便快捷性地测定出点位坐标,无论是操作上还是精度上,比全站仪等其他常规测量设备有明显的优越性。
随着我国各地GPS差分台站的不断建立以及美国SA政策的取消,使得单机定位的精度大大提高,有的已经达到了亚米级精度,能够满足国土资源调查、土地利用更新、遥感监测、海域使用权清查等工作的应用。
在一般情况下,我们使用的是1954年北京坐标系或1980年西安坐标系(以下分别简称54系和80系),而GPS测定的坐标是WGS-84坐标系坐标,需要进行坐标系转换。
对于非测量专业的工作人员来说,虽然GPS定位操作非常容易,但坐标转换则难以掌握,EXCEL是比较普及的电子表格软件,能够处理较复杂的数学运算,用它来进行GPS坐标转换、面积计算会非常轻松自如。
要进行坐标系转换,离不开高斯投影换算,下面分别介绍用EXCEL进行换算的方法和GPS坐标转换方法。
一、用EXCEL进行高斯投影换算从经纬度BL换算到高斯平面直角坐标XY(高斯投影正算),或从XY换算成BL(高斯投影反算),一般需要专用计算机软件完成,在目前流行的换算软件中,存在一个共同的不足之处,就是灵活性较差,大都需要一个点一个点地进行,不能成批量地完成,给实际工作带来许多不便。
笔者发现,用EXCEL可以很直观、方便地完成坐标换算工作,不需要编制任何软件,只需要在EXCEL的相应单元格中输入相应的公式即可。
下面以54系为例,介绍具体的计算方法。
完成经纬度BL到平面直角坐标XY的换算,在EXCEL中大约需要占用21列,当然读者可以通过简化计算公式或考虑直观性,适当增加或减少所占列数。
Excel 2003在高斯投影坐标换算中的应用
Excel 2003在高斯投影坐标换算中的应用
谷超杰
【期刊名称】《测绘与空间地理信息》
【年(卷),期】2008(31)6
【摘要】对Excel 2003函数功能进行充分的应用,使其成为批量坐标换算的一个真实有效的工具.针对坐标换算处理的整个过程,给出坐标换算的新方法.重点论述运用Excel 2003函数功能实现高斯投影坐标变换,可正算、反算、换带、求面积.同时结合实例体现Excel 2003处理数据的特点,给出一个批量坐标换算的高效率、便捷和直观的方法.
【总页数】3页(P190-192)
【作者】谷超杰
【作者单位】黑龙江省第六地质勘察院,黑龙江,佳木斯,154002
【正文语种】中文
【中图分类】P282.1
【相关文献】
1.Excel在高斯投影坐标正算公式中的应用 [J], 崔学敏;赵雪松;姚继华
2.矿山测量中高斯投影坐标计算的Excel实现 [J], 付永波;仲树强;陈刚;王洪涛
3.基于matlab的高斯投影正反算与相邻带坐标换算程序设计 [J], 徐翰;周强波
4.Excel在平面坐标换算中的应用 [J], 洪建军;洪盛远
5.用CASIO fx—4500P进行高斯投影带坐标换算 [J], 张崇安
因版权原因,仅展示原文概要,查看原文内容请购买。
EXCEL公式进行经纬度与XY坐标的相互转换
EXCEL公式进行经纬度与XY坐标的相互转换•一、用EXCEL进行高斯投影换算从经纬度B、L换算到高斯平面直角坐标X、Y(高斯投影正算),或从X、Y换算成B、L(高斯投影反算),一般需要专用计算机软件完成。
在目前流行的换算软件中不足之处,就是灵活性较差,大都需要一个点一个点地进行,不能成批量地完成,给实际工作带来许多不便。
而用EXCEL可以很直观、方便地完成坐标换算工作,不需要编制任何软件,只需要在EXCEL的相应单元格中输入相应的公式即可。
下面以1954年北京坐标系为例,介绍具体的计算方法。
上图为编辑好的EXCEL表(红色为输入数据项)完成经纬度B、L到平面直角坐标X、Y的换算,在EXCEL中大约需要占用21列,当然读者可以通过简化计算公式或考虑直观性,适当增加或减少所占列数。
在EXCEL中以公式从第3行第1列(A3格)为起始单元格为例,各单元格的公式如下:(1)单元格A3输入中央子午线,以度、分、秒形式输入,如107度0分则输入107.00 (2)单元格B3公式如上图,把L0化成度形式。
(3)单元格C3以度小数形式输入纬度值,如23°44′01″则输入23.4401。
(4)单元格D3以度小数形式输入经度值,如107°42′48″则输入107.4248。
(5)单元格E3公式如上图,把纬度B化成度形式。
(6)单元格F3公式如上图,把经度L化成度形式。
(7)各个单元格输入公式如下:表中计算公式见由孔祥元等主编、武汉大学2002年出版的《控制测量学》,EXCEL软件的操作方法请参阅有关资料。
按上面表格中的公式输入到相应单元格后,就可方便地由经纬度求得平面直角坐标。
当输入完所有的经纬度后,用鼠标下拉即可得到所有的计算结果。
表中的许多单元格公式为中间过程,可以用EXCEL的列隐藏功能把这些没有必要显示的列隐藏起来,表面上形成标准的计算报表,使整个计算表简单明了,可计算的数据量是无限的,当第一次输入公式后,相当于自己完成了一软件的编制,可存储起来供今后重复使用。
用EXCEL批量计算坐标转换
用EXCEL批量计算坐标转换在Excel中进行坐标转换,可以通过使用VBA(Visual Basic for Applications)编写宏来实现批量计算。
下面是一个示例,可以将经纬度转换为高斯投影坐标:1. 在Excel中创建一个新的工作表。
2.在第一列中输入经度值,例如A1单元格输入:经度。
3.在第二列中输入纬度值,例如B1单元格输入:纬度。
4.在第三列中输入转换后的X坐标的标题,例如C1单元格输入:高斯X坐标。
5.在第四列中输入转换后的Y坐标的标题,例如D1单元格输入:高斯Y坐标。
8.在模块中输入以下VBA代码:```vbaOption ExplicitSub ConvertCoordinatesDim ws As WorksheetDim lng As Double, lat As Double '经纬度Dim x As Double, y As Double '高斯坐标Dim i As Long'指定要处理的工作表Set ws = ThisWorkbook.Sheets("Sheet1")'从第二行开始遍历每一行For i = 2 To ws.Cells(ws.Rows.Count, 1).End(xlUp).Row'读取经纬度值lng = ws.Cells(i, 1).Valuelat = ws.Cells(i, 2).Value'进行坐标转换,可以使用适当的转换算法'以下是一个示例,将经度值转换为高斯X坐标,纬度值转换为高斯Y坐标x = lng * 200y = lat * 100'将转换后的坐标值写入对应的单元格ws.Cells(i, 3).Value = xws.Cells(i, 4).Value = yNext iEnd Sub```10. 返回Excel界面,并在工作表上运行宏。
选择开发选项卡中的宏,找到ConvertCoordinates并点击运行。
Excel 2003在高斯投影坐标换算中的应用
c se tera z t no oria o v r o a s po c o s gte u ci s f x e 2 0 .t el e efn t no oiv u ss h l a o f odn t c n es n i G u s r e t nu i n t n cl 0 3 Irai dt c o f s i ei i c e i n j i n h f o oE z h u i p te a dn gt ec l l i f or ia ,o ria a s r t n b te nt opo ci n s a dae lu t n I em a t e i n eai a ua o o o n t c odnt t nf ma o e e r e t nz e ,n rac c l i .nt e ni , v c tn c d e er o i w w j o o a ao h m t
d s u s st e c a a tr t s a d a v n a e o rc s ig d t sn x e 0 3 c mb n d wi n e a l a d p t f r a d a meh d ic s e h h cei i n d a tg f o e s aa u ig E c l 0 o ie t a x mpe, n us o w r to r sc p n 2 h o ac o dn t o v rin wi h g frb t h c r i ae c n eso t ih—e iin y a d c n e i n e h f ce c n o v n e c .
形投影 的一种 。
测量工作 中 , 量 数 据 的 处理 是 测 量 工 作 中 的核 心 测 工作 , 量坐标 的换 算 是 数 据处 理 工 作 中最 耗 时 的工 作 测 内容之一 , 高 坐标 换 算 的 速度 是 我 们测 绘 工 作 者 的 目 提
Excel VBA在GPS坐标转换计算中的应用
Excel VBA在GPS坐标转换计算中的应用摘 要为了使GPS的观测成果在实际中得到应用,必须把GPS观测得到的WGS-84坐标转换成实际需要的国家大地坐标或地方独立坐标。
本文利用Excel软件的宏程序VBA进行编程,实现不同坐标系之间坐标的相互转换。
并且给出一些关于坐标转换的算例,结果表明:所编程序可以方便、快捷、可靠、准确地解决不同坐标间的转换问题,能够满足设计计算工作的要求,是可行的。
关键词Excel VBA;GPS;坐标转换中图分类号P2 文献标识码 A 文章编号1673-9671-(2012)011-0165-01GPS测量直接提供的坐标(B,L,H)是1984年世界大地坐标系(Word Geodetic System 1984即WGS-84)的坐标。
而在实际应用中,我国采用的是1954北京坐标系或者1980西安坐标系下的高斯投影坐标(x,y)。
如何实现WGS-84坐标系与北京54坐标系或者1980西安坐标系的转换,一直是GPS应用中人们关心的热点。
1 坐标系统1.1 地形面、参考椭球面和大地水准面地形面是地球的自然表面,包含了地球外部的陆地地形和海底地形,是一个极度不规则的面。
大地水准面是一个物理参考面,是地球的一个重力等位面。
由于陆地质量的影响,存在着局部的重力异常,它们将使大地水准面的形状发生扭曲。
参考椭球面是一个几何或数学参考面,是一个与大地水准面非常接近的旋转椭球面。
1.2 坐标、坐标系、基准和坐标参照系坐标是用于在一个给定维数的空间中相对一个参照系来确定点的位置的一组数。
坐标系指的是描述空间位置的表达形式。
基准指的是为描述空间位置而定义的一系列点、线、面。
在大地测量中的基准一般是指为确定点在空间中的位置,而采用的地球椭球或参考椭球的几何参数和物理参数,及其在空间的定位、定向方式,以及在描述空间位置时所采用的单位长度的定义。
1.3 常用坐标系的表现形式1)空间直角坐标系。
用EXCEL转换经纬度坐标到北京54坐标的方法
用EXCEL转换经纬度坐标到北京54坐标的方法[复制链接]lihongbing lihongbing 当前离线注册时间2007-10-27在线时间72 小时帖子133精华0威望6 点土地券1618 亩活跃度231 度看IP卡狗仔卡版主主题7华夏券0 张串个门加好友打招呼发消息电梯直达1楼发表于2008-4-1 15:51:23 |只看该作者|倒序浏览分享到:用EXCEL进行高斯投影换算从经纬度BL换算到高斯平面直角坐标XY(高斯投影正算),或从XY换算成BL(高斯投影反算),一般需要专用计算机软件完成,在目前流行的换算软件中,存在一个共同的不足之处,就是灵活性较差,大都需要一个点一个点地进行,不能成批量地完成,给实际工作带来许多不便。
笔者发现,用EXCEL可以很直观、方便地完成坐标换算工作,不需要编制任何软件,只需要在EXCEL的相应单元格中输入相应的公式即可。
下面以54系为例,介绍具体的计算方法。
完成经纬度BL到平面直角坐标XY的换算,在EXCEL中大约需要占用21列,当然读者可以通过简化计算公式或考虑直观性,适当增加或减少所占列数。
在EXCEL中,输入公式的起始单元格不同,则反映出来的公式不同,以公式从第2行第1列(A2格)为起始单元格为例,各单元格的公式如下:单元格单元格内容说明A2输入中央子午线,以度.分秒形式输入,如115度30分则输入115.30起算数据L0B2=INT(A2)+(INT(A2*100)-INT(A2)*100)/60+(A2*10000-INT(A2*100)*100)/3600把L0化成度C2以度小数形式输入纬度值,如38°14′20″则输入38.1420起算数据BD2以度小数形式输入经度值起算数据LE2=INT(C2)+(INT(C2*100)-INT(C2)*100)/60+(C2*10000-INT(C2*100)*100)/3600 把B化成度F2=INT(D2)+(INT(D2*100)-INT(D2)*100)/60+(D2*10000-INT(D2*100)*100)/3600 把L化成度G2=F2-B2L-L0H2=G2/57.2957795130823化作弧度I2=TAN(RADIANS(E2))Tan(B)J2=COS(RADIANS(E2))COS(B)K2=0.006738525415*J2*J2L2=I2*I2M2=1+K2N2=6399698.9018/SQRT(M2)O2=H2*H2*J2*J2P2=I2*J2Q2=P2*P2R2=(32005.78006+Q2*(133.92133+Q2*0.7031))S2=6367558.49686*E2/57.29577951308-P2*J2*R2+((((L2-58)*L2+61)*O2/30+(4*K2+5)*M2-L2)*O2/12+1)*N2*I2*O2/2计算结果XT2=((((L2-18)*L2-(58*L2-14)*K2+5)*O2/20+M2-L2)*O2/6+1)*N2*(H2*J2)计算结果Y表中公式的来源及EXCEL软件的操作方法,请参阅有关资料,这里不再赘述。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
利用EXCEL电子表格进行高斯投影换算GPS坐标的方法
2009-06-13 10:05
[摘要] 对利用EXCEL电子表格进行高斯投影换算的方法进行了较详细的介绍,对如何进行GPS坐标系转换进行了分析,提出了一种简单实用的坐标改正转换方法,介绍了用EXCEL完成转换的思路。
[关键字] 电子表格;GPS;坐标转换
作为尖端技术GPS,能方便快捷性地测定出点位坐标,无论是操作上还是精度上,比全站仪等其他常规测量设备有明显的优越性。
随着我国各地GPS差分台站的不断建立以及美国SA政策的取消,使得单机定位的精度大大提高,有的已经达到了亚米级精度,能够满足国土资源调查、土地利用更新、遥感监测、海域使用权清查等工作的应用。
在一般情况下,我们使用的是1954年北京坐标系或1980年西安坐标系(以下分别简称54系和80系),而GPS测定的坐标是WGS-84坐标系坐标,需要进行坐标系转换。
对于非测量专业的工作人员来说,虽然GPS定位操作非常容易,但坐标转换则难以掌握,EXCEL是比较普及的电子表格软件,能够处理较复杂的数学运算,用它来进行GPS坐标转换、面积计算会非常轻松自如。
要进行坐标系转换,离不开高斯投影换算,下面分别介绍用EXCEL进行换算的方法和GPS坐标转换方法。
一、用EXCEL进行高斯投影换算
从经纬度BL换算到高斯平面直角坐标XY(高斯投影正算),或从XY换算成BL(高斯投影反算),一般需要专用计算机软件完成,在目前流行的换算软件中,存在一个共同的不足之处,就是灵活性较差,大都需要一个点一个点地进行,不能成批量地完成,给实际工作带来许多不便。
笔者发现,用EXCEL可以很直观、方便地完成坐标换算工作,不需要编制任何软件,只需要在EXCEL的相应单元格中输入相应的公式即可。
下面以54系为例,介绍具体的计算方法。
完成经纬度BL到平面直角坐标XY的换算,在EXCEL中大约需要占用21列,当然读者可以通过简化计算公式或考虑直观性,适当增加或减少所占列数。
在EXCEL中,输入公式的起始单元格不同,则反映出来的公式不同,以公式从第2行第1列(A2格)为起始单元格为例,各单元格的公式如下:
单元格
单元格内容
说明
A2
输入中央子午线,以度.分秒形式输入,如115度30分则输入115.30
起算数据L0
B2
=INT(A2)+(INT(A2*100)-INT(A2)*100)/60+(A2*10000-INT(A2*100)*100)/3600
把L0化成度
C2
以度小数形式输入纬度值,如38°14′20″则输入38.1420
起算数据B
D2
以度小数形式输入经度值
起算数据L
E2
=INT(C2)+(INT(C2*100)-INT(C2)*100)/60+(C2*10000-INT(C2*100)*100)/3600 把B化成度
F2
=INT(D2)+(INT(D2*100)-INT(D2)*100)/60+(D2*10000-INT(D2*100)*100)/3600 把L化成度
G2
=F2-B2
L-L0
H2
=G2/57.2957795130823
化作弧度
I2
=TAN(RADIANS(E2))
Tan(B)
J2
=COS(RADIANS(E2))
COS(B)
K2
=0.006738525415*J2*J2
L2
=I2*I2
M2
=1+K2
N2
=6399698.9018/SQRT(M2)
O2
=H2*H2*J2*J2
P2
=I2*J2
Q2
=P2*P2
R2
=(32005.78006+Q2*(133.92133+Q2*0.7031))
S2
=6367558.49686*E2/57.29577951308-P2*J2*R2+((((L2-58)*L2+61)*
O2/30+(4*K2+5)*M2-L2)*O2/12+1)*N2*I2*O2/2
计算结果X
T2
=((((L2-18)*L2-(58*L2-14)*K2+5)*O2/20+M2-L2)*O2/6+1)*N2*(H2*J2)
计算结果Y
表中公式的来源及EXCEL软件的操作方法,请参阅有关资料,这里不再赘述。
按上面表格中的公式输入到相应单元格后,就可方便地由经纬度求得平面直角坐标。
当输入完所有的经纬度后,用鼠标下拉即可得到所有的计算结果。
表中的许多单元格公式为中间过程,可以用EXCEL的列隐藏功能把这些没有必要显示的列隐藏起来,表面上形成标准的计算报表,使整个计算表简单明了。
从理论上讲,可计算的数据量是无限的,当第一次输入公式后,相当于自己完成了一软件的编制,可另存起来供今后重复使用,一劳永逸。
二、GPS坐标转换方法与面积计算
GPS所采用的坐标系是美国国防部1984世界坐标系,简称WGS-84,它是一个协议地球参考系,坐标系原点在地球质心。
GPS的测量结果与我国的54系或80系坐标相差几十米至一百多米,随区域不同,差别也不同,经粗落统计,我国西部相差70米左右,东北部140米左右,南部75米左右,中部45米左右。
由此可见,必须将WGS-84坐标进行坐标系转换才能供标图使用。
坐标系之间的转换一般采用七参数法或三参数法,其中七参数为X平移、Y平移、Z平移、X旋转、Y旋转、Z旋转以及尺度比参数,若忽略旋转参数和尺度比参数则为三参数方法,三参数法为七参数法的特例。
这里的Z、Y、Z是空间大地直角坐标系坐标,为转换过程的中间值。
在实际工作中我们常用的是平面直角坐标,是否可以跳过空间直角坐标系,省略复杂的运算,进行简单转换呢?为此,笔者进行了长期的实践,证明是可行的。
其在原理是:不把GPS所测定的WGS-84坐标当作WGS-84坐标,而是当作具有一定系统性误差的54系坐标值,然后通过国家已知点纠正,消除该系统误差。
我们暂把该方法称作坐标改正法,下面以WGS-84坐标转换成54系坐标为例,介绍数据处理方法:
首先,在测区附近选择一国家已知点,在该已知点上用GPS测定WGPS-84坐标系经纬度B和L,把此坐标视为有误差的54系坐标,利用54系EXCEL将经纬度BL转换成平面直角坐标X’Y’,然后与已知坐标比较则可计算出偏移量:
△X=X-X’
△Y=Y-Y’
式中的X、Y为国家控制点的已知坐标,X’、Y’为测定坐标,△X和△Y为偏移量。
求得偏移量后,就可以用此偏移量纠正测区内的其他测量点了。
把其他GPS测量点的经纬度测量值,转换成平面坐标X’Y’,在此XY坐标值上直接加上偏移值就得到了转换后的54系坐标:
X=X’+△X
Y=Y’+△Y
在上述EXCEL计算表的最后两列,附加上求得的改正数并分别与计算出来的XY相加后,即得到转换结果。
若测量路线是一闭合区域的话,可把计算结果按路线顺序排列起来,再输入相应的计算公式,即可计算出该区域的面积。
有关用坐标计算面积的原理与公式,这里不再叙述,读者可参阅有关资料。
需要说明的是,面积的计算精度基本上不受坐标转换精度的影响,若只需要求算面积的话,可不进行坐标系转换这一步,只需要把BL化成XY就行了。
就1:1万比例尺成图而言,在一般的县行政区范围内(如40Km×40Km),用此简单的坐标改正法进行转换与较复杂的七参数法没有多大差别。
能否满足1:1万比例尺变更调查的要求,主要取决于GPS接收机本身的精度,与转换方法的选择关系不大。
当面积较大时,使用该方法可能会使误差增大,这时可考虑分区域转换。