热门标签 | HotTags
当前位置:  开发笔记 > 前端 > 正文

手持GPS坐标系转换的心得体会

手持GPS坐标系转换的心得体会(福建省水利水电勘测设计研究院;福建省闽江学院)李桂炎陈海金李青对于坐标系的转换,GPS的使用者造成一知半解
手持GPS坐标系转换的心得体会
(福建省水利水电勘测设计研究院;福建省闽江学院)

李桂炎 陈海金 李青

对于坐标系的转换,GPS的使用者造成一知半解,尤其是对于接触不久的人,搞不明白竟然是怎么一回事。我经过多年来使用手持GPS来抛砖引玉,希望能引出更多的高手来为我院风电等项目导航。

归纳常见的坐标转换问题,多数为WGS84转换成西安1980坐标系或1954年北京坐标系。其中WGS84坐标系属于大地坐标,就是我们常说经纬度的坐标,而西安1980坐标系和1954年北京坐标系属于平面直角坐标。对于什么是大地坐标,什么是平面直角坐标,以及它们如何建立,本人将分段来介绍.

一、坐标怎样转换?

GPS卫星星历是以WGS84坐标系为根据而建立的,我国目前应用的地形图却属于1954年北京坐标系或1980年国家大地坐标系;因为不同坐标系之间存在着平移和旋转关系(WGS84坐标系与我国应用的坐标系之间的误差约为75~80m),所以在我国应用GPS进行绝对定位必须进行坐标转换,转换后的绝对定位精度可由75~80提高到5~10m。精确一句话,“减少误差,提高精度”。
二、如何在WGS84坐标系和1954年北京坐标系之间进行转换呢?
学习测量的人都知道,地球是一个近似的椭球体。因此为了研究方便,科学家们根据各自的理论建立了不同的椭球模型来模拟地球的形状。而且我们刚才讨论了半天的各种坐标系也是建立在这些椭球基准之上的。比如北京54坐标系采用的就是克拉索夫斯基椭球模型。而对应于 WGS84坐标系有一个WGS84椭球,其常数采用 IUGG第17届大会大地测量常数的推荐值。WGS84椭球两个最常用的几何常数:长半轴:6378137±2(m);扁率:1:298.257223563

之所以说到半长轴和扁率倒数是因为要在不同的坐标系之间转换,就需要转换不同的椭球基准。这就需要两个很重要的转换参数DA、DF。其中DA的含义是两个椭球基准之间半长轴的差;DF的含义是两个椭球基准之间扁率倒数的差。

在进行坐标转换时,这两个转换参数是固定的,这里,我们给出在WGS84转换成1954年北京坐标系进行参数如下:
例如:BJ54坐标转换WGS84坐标参数:DA:-108; DF:0.0000005

DX:22; DY:-118; DZ:-30.5,同理,也可以把西安1980坐标转换WGS84坐标参数。

福建省地区GPS坐标转换参数(福建省二类调查GPS坐标转换参数)

WGS84转BJ54坐标参数 (表1)

序号

DX DY DZ DA DF 适宜地区

1 -21 -125 -44 -108 0.0000005 福州、宁德

2 -19 -116 -57 -108 0.0000005 漳州、泉州、厦门

3 -19 -120 -47 -108 0.0000005 南平、三明

4 -20 -123 -42 -108 0.0000005 龙岩

5 -21 -120 -54 -108 0.0000005 莆田

三、椭球的基准转换过来了,还要那些计算?
椭球的基准转换过来了,那么由于建立椭球的原点还是不一致的,还需要在DX、DY、DZ这三个空间平移参量,来将两个不同的椭球原点重合,这样一来才能使两个坐标系的椭球完全转换过来。而由于各地的地理位置不同,所以在各个地方的这三个坐标轴平移参量也是不同的,因此需要用当地的已知点来计算这三个参数。具体的计算方法是:
第一步:搜集应用区域内GPS“B”级网三个以上网点WGS84坐标系B、L、H值及我国坐标系(BJ54或西安80)B、L、h、x值。(注:B、L、H分别为大地坐标系中的大地纬度、大地经度及大地高,h、x分别为大地坐标系中的高程及高程异常。各参数可以通过各我省测绘局或测绘院具有“A”级、“B”级网的单位获得。)
第二步:计算不同坐标系三维直角坐标值。计算公式如下:
X=(N+H)cosBcosL
Y=(N+H)cosBsinL
Z=[N(1-e2)+H]sinB

不同坐标系对应椭球的有关常数详见下表:

(表2)

项 目 WGS84坐标系 BJ54坐标系 西安80坐标系

A 6378137m 6378245m 6378140m

E2 0.00669437999013 0.006693427 0.006694385

F 1/298.257223563 1/298.3 1/298.257

(注:X、Y、Z为大地坐标系中的三维直角坐标;A为大地坐标系对应椭球之长半轴;e2为大地坐标系对应椭球第一偏心率; N为该点的卯酉圈曲率半径,N=A/(1-e2sin2B)1/2;H=h+x,该处H为BJ54或西安80坐标系中的大地高)。
第三步:求出DX,DY,DZ。即利用WGS84坐标系的X、Y、Z值,减去我国坐标系的对应值,得出实现坐标系统转换的三个参数。(应算出WGS84与北京

和西安坐标系两套参数。)

第四步:参数验证。参数计算之后必须对其进行验证。验证的方法是在应用区域内选择5个以上水准点进行实测,实测值与测绘部门提供的理论值对比,

如果最大误差不大于15米,平均误差不大于10米,则计算出的参数可以使用,否则要重新计算或查找出现问题的原因。

对了,还有一个很重要的事情,要在位置格式的地方,选择用户自定义方式,输入如下参数:

中央经线: 福建省各地区根据省测绘局提供的中央经线来输入计算

投影比例:1;

东西偏差:500000;

四、福建省风海上钻探K点坐标实例计算:

(根据测绘局海测图提供经度、纬度如下)
(表3)
点 经度 纬度
度 分 秒 度 分 秒

K1 119 26 25.2 25 42 23.0

K2 119 22 44.6

25 42 43.6

K3 119 27 36.3

25 43 15.2

高斯正算计算表(BJ54)

X前缀=0米 Y前缀=500000米 椭球:长半轴=6378245米 1/扁率=298.3 第 1页

序号 点名 纬度B 经度L 中央子午线L0 纵坐标X 横坐标Y
1

K1 25.42230 119.26252 120 2844476.292 443827.392

2 K2 25.42436 119.22446 120 2845137.741 437679.903
3 K3 25.43152 119.27363 120 2846074.476 445816.265


计算表(西安80)

X前缀=0米 Y前缀=500000米 椭球:长半轴= 6378140米 1/扁率=298.257 第 2页

序号 点名 纬度B 经度L 中央子午线L0 纵坐标X 横坐标Y
1

K1 25.42230 119.26252 120 2844426.972 443828.312

2 K2 25.42436 119.22446 120 2845088.410 437680.923
3 K3 25.43152 119.27363 120 2846025.129 445817.152


计算表(WGS84)

X前缀=0米 Y前缀=500000米 椭球:长半轴= 6378137米 1/扁率=298.257223563 第 3页

序号 点名 纬度B 经度L 中央子午线L0 纵坐标X 横坐标Y
1

K1 25.42230 119.26252 120 2844425.647 443828.338

2 K2 25.42436 119.22446 120 2845087.085 437680.952
3 K3 25.43152 119.27363 120 2846023.803 445817.178
从高斯正算计算表中可以看出西安80与WGS84差值比较小,以上在进行坐标转换的五个参数都已经得到了,那么可以将BJ54坐标参数输入到GPS手持机中,到野外进行点位放样与测量了。

转:https://www.cnblogs.com/bluespot/archive/2009/08/03/1537477.html



推荐阅读
  • 本文详细介绍了如何使用ActionScript 3.0 (AS3) 连接并操作MySQL数据库。通过具体的代码示例和步骤说明,帮助开发者理解并实现这一过程。 ... [详细]
  • 在即将迎来26岁生日之际,作者的人生陷入了低谷。经过近三年的硕士学习后,最终决定退学,并且面临没有工作经验的困境。尽管如此,作者依然坚定地选择为自己的人生负责。 ... [详细]
  • 使用Python在SAE上开发新浪微博应用的初步探索
    最近重新审视了新浪云平台(SAE)提供的服务,发现其已支持Python开发。本文将详细介绍如何利用Django框架构建一个简单的新浪微博应用,并分享开发过程中的关键步骤。 ... [详细]
  • 本文详细介绍了美国最具影响力的十大财团,包括洛克菲勒、摩根、花旗银行等。这些财团在历史发展过程中逐渐形成,并对美国的经济、政治和社会产生深远影响。 ... [详细]
  • 基于KVM的SRIOV直通配置及性能测试
    SRIOV介绍、VF直通配置,以及包转发率性能测试小慢哥的原创文章,欢迎转载目录?1.SRIOV介绍?2.环境说明?3.开启SRIOV?4.生成VF?5.VF ... [详细]
  • 深入探讨CPU虚拟化与KVM内存管理
    本文详细介绍了现代服务器架构中的CPU虚拟化技术,包括SMP、NUMA和MPP三种多处理器结构,并深入探讨了KVM的内存虚拟化机制。通过对比不同架构的特点和应用场景,帮助读者理解如何选择最适合的架构以优化性能。 ... [详细]
  • 本题通过将每个矩形视为一个节点,根据其相对位置构建拓扑图,并利用深度优先搜索(DFS)或状态压缩动态规划(DP)求解最小涂色次数。本文详细解析了该问题的建模思路与算法实现。 ... [详细]
  • 在网页开发中,页面加载速度是一个关键的用户体验因素。为了提升加载效率,避免在PageLoad事件中进行大量数据绑定操作,可以采用异步加载和特定控件来优化页面加载过程。 ... [详细]
  • 本文介绍了一种根据用户选择动态切换屏幕界面的方法,通过定义不同的选择块(Selection Block),实现灵活的用户交互体验。 ... [详细]
  • 本题探讨了在一个有向图中,如何根据特定规则将城市划分为若干个区域,使得每个区域内的城市之间能够相互到达,并且划分的区域数量最少。题目提供了时间限制和内存限制,要求在给定的城市和道路信息下,计算出最少需要划分的区域数量。 ... [详细]
  • 本文详细探讨了HTML表单中GET和POST请求的区别,包括它们的工作原理、数据传输方式、安全性及适用场景。同时,通过实例展示了如何在Servlet中处理这两种请求。 ... [详细]
  • 在现代Web应用中,当用户滚动到页面底部时,自动加载更多内容的功能变得越来越普遍。这种无刷新加载技术不仅提升了用户体验,还优化了页面性能。本文将探讨如何实现这一功能,并介绍一些实际应用案例。 ... [详细]
  • 解决SVN图标显示异常问题的综合指南
    本文详细探讨了SVN图标无法正常显示的问题,并提供了多种有效的解决方案,涵盖不同环境下的具体操作步骤。通过本文,您将了解如何排查和修复这些常见的SVN图标显示故障。 ... [详细]
  • 磁盘健康检查与维护
    在计算机系统运行过程中,硬件或电源故障可能会导致文件系统出现异常。为确保数据完整性和系统稳定性,定期进行磁盘健康检查至关重要。本文将详细介绍如何使用fsck和badblocks工具来检测和修复文件系统及硬盘扇区的潜在问题。 ... [详细]
  • 本文将探讨Java编程语言中对象和类的核心概念,帮助读者更好地理解和应用面向对象编程的思想。通过实际例子和代码演示,我们将揭示如何在Java中定义、创建和使用对象。 ... [详细]
author-avatar
sherklock
这个家伙很懒,什么也没留下!
PHP1.CN | 中国最专业的PHP中文社区 | DevBox开发工具箱 | json解析格式化 |PHP资讯 | PHP教程 | 数据库技术 | 服务器技术 | 前端开发技术 | PHP框架 | 开发工具 | 在线工具
Copyright © 1998 - 2020 PHP1.CN. All Rights Reserved | 京公网安备 11010802041100号 | 京ICP备19059560号-4 | PHP1.CN 第一PHP社区 版权所有