地心地固坐标系WGS84经纬度与笛卡尔直角坐标系相互转换的推导、理解与代码实现(C++和matlab)

本文介绍了中国、美国和日本的大地坐标系发展历史,并详细解析了经纬度坐标与直角坐标系之间的转换方法,包括数学推导及C++、Matlab代码实现。
该文章已生成可运行项目,



一、大地坐标系现状简析

1.1 我国

  解放后在我国大陆地区曾先后使用了三种坐标系。其一,于1954年引入了苏联1942年普尔科沃坐标系(PULK0V01942),经局部平差建立了“1954年北京坐标系(BJ54)”,使用克拉索夫斯基(Krassowsky)椭球,原点位于前苏联普尔科沃;其二,在20世纪70年代末期,利用天文大地测量资料通过统一整体平差,建立了“1980年国家大地坐标系”,使用UGG75椭球,大地原点位于陕西省西安市北泾阳县永乐镇,又称“1980年西安大地坐标系”;其三,在1980年国家大地坐标系基础上建立“新1954年北京坐标系”,恢复使用了克拉索夫斯基椭球,有利于新坐标成果的使用和新、旧图的拼接。
  未来几年,我国将逐渐使用新一代大地坐标系“2000国家大地坐标系”(China Geodetic Coordinate System2000,CGCS2000)。其定义与协议地球参考系的定义一致,即原点是包括海洋和大气的整个地球质心,初始定向由1984.0时国际时间局定向给定,定向的时间演化使得地壳无整体旋转。

1.2 美国

  美国国防部以军事目的获取全球地面定位信息,相继建立了多个全球地心坐标系。美国国防部于1984年继WGS60、WGS66、WGS72之后建立起来第四个世界大地坐标系世界大地坐标系(World Geodetic System1984,WGS84)。WGS84规定其原点在地球质心、Z轴指向BH1984.0定义的协议地球极方向、X轴指向BH1984.0的零子午面和协议地球极赤道的交点,Y轴与X、Z轴构成右手地心地固ECEF直角坐标系。

1.3 日本

  日本在1918年之前,使用的是“日本东京1892年大地坐标系(JTD1892)”,又称“旧东京坐标系”,坐标原点在东京麻布旧东京天文台(旧海军观象台)子午仪中心点,采用1841年贝塞尔椭球,a=6377397.155米,1/f=299.1528。由于东京坐标系与WGS84坐标系的差异产生的问题日益突出。于是,2001年6月日本新测量法明确规定了其新的大地基准应该与WGS84保持一致,同时新建立的“日本大地基准2000(JGD2000)”取代了“东京坐标系”。东京坐标系通过中间坐标系Tokyo97与JGD2000大地坐标系(水平)相连,Tokyo97采用东京坐标系原点和贝塞尔参考椭球,轴向与1TRF94一致。东京坐标系坐标首先通过平差或网格内插变换至Tokyo97后,然后再变换至TRF94。这样便形成了JGD2000大地坐标系(水平),它采用GRS80椭球。

二、经纬度坐标(L,B,H)与直角坐标系坐标(X,Y,Z)相互转换

  以目前使用最为广泛的WGS84的经纬度坐标进行转换,是基于同一基准的大地坐标系转换。

2.1 正解(L,B,H)——>(X,Y,Z)

2.1.1 数学推导

在这里插入图片描述
  如上示意图,求大地坐标系中点P(L,B,H)在大地空间直角坐标系的坐标(X,Y,Z)分为三步:

  1. 首先,求得P点在其子午线平面坐标系O-xy内坐标(x,y)与纬度B的关系:

在这里插入图片描述

  1. 其次,求得子午线平面坐标(x,y)在大地空间直角坐标系坐标(X,Y,

本文章已经生成可运行项目
评论 18
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

可见一班

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值