python——经纬度坐标和平面投影坐标的相互转换

美赛建模需要使用到平面坐标系去做基于时间序列的运动范围模型，所以需要对经纬度进行坐标系转换，转换到XY坐标系下。同时，在机器人的路径规划等方面依然会经常使用到坐标系的转换，所以记录一下。

文章目录

经纬度坐标系和平面投影坐标系
- 地理坐标系 / 大地坐标系
- 平面投影坐标系
墨卡托投影法
- 变量和常量定义
- （1）墨卡托投影正解算公式
- （2）墨卡托投影逆解算公式
Python代码实现
- （1）墨卡托投影正解算公式
- （2）墨卡托投影逆解算公式

经纬度坐标系和平面投影坐标系

地理坐标系 / 大地坐标系

用经纬度表示的是地理坐标系，也称大地坐标系。

大地测量中以参考椭球面为基准面的坐标。地面点P的位置用大地经度L、大地纬度B和大地高H表示。 当点在参考椭球面上时，仅用大地经度和大地纬度表示。大地经度是通过该点的大地子午面与起始大地子午面之间的夹角，大地纬度是通过该点的法线与赤道面的夹角，大地高是地面点沿法线到参考椭球面的距离。

平面投影坐标系

人们比较习惯于使用平面坐标系，平面坐标系用xy表示。

把球体表面的坐标转成平面坐标需要一定的手段，这个手段称为投影。投影是XYZ平面内的。有时候用地理坐标系不够方便，投影方法也不是唯一的，还是为了一个目的，务求使当地的坐标最准确。所以目前就存在了好多投影方法，比如高斯投影、墨卡托投影等。

大地坐标、地理坐标均是球面的，投影坐标是平面的。

墨卡托投影法

若点（B,L）经过墨卡托投影（下图）得到的新坐标为点（X,Y），
其中B0为标准纬度，L0为标准经度，e 为第一偏心率，e ’为第二偏心率，根据墨卡托投影方法我们可以实现坐标系的相互转化。

变量和常量定义

椭圆的长半轴a, 椭圆的短半轴b：

（1）墨卡托投影正解算公式

由经纬度坐标（B,L）计算得到平面投影坐标（X,Y）：

（2）墨卡托投影逆解算公式

由平面投影坐标（X,Y）计算得到经纬度坐标（B,L）。可以用牛顿迭代的方法进行解算，即已知q求B，方法如下：

Python代码实现

（1）墨卡托投影正解算公式

代码还有小问题，需要各位同学帮忙指正。

from math import *def get_coordinate(latitude,longitude):B = latitudeL = longitudeq = log( tan(pi/4 + B/2)*(1-exp(1)*sin(B)) / (1+exp(1)*sin(B))**exp(1/2) )x = K*qy = K*(L-L0)return (x, y)if __name__ == '__main__':latitude = 34longitude = -167a = 6378137.0000 # 单位为mb = 6356752.3142B0 = 0; L0 = 0;e1 = sqrt(pow(a,2) - pow(b,2)) / ae2 = sqrt(pow(a,2) - pow(b,2)) / bK = a* cos(B0) / sqrt(1-pow(exp(2), 2)*pow(sin(B0), 2))data = get_coordinate(latitude,longitude)print(data)

运行示例：

（2）墨卡托投影逆解算公式

公式已有，可以自行复现。

参考论文：

[1]胡正,杨青,卜晓楠,张波,江伟健,宋祥瑞.TDOA定位中经纬度与平面坐标转换方案的研究[J].电子世界,2019(03):36-37.

python——经纬度坐标和平面投影坐标的相互转换相关推荐

平面投影坐标转经纬度坐标
平面投影坐标转经纬度坐标本人在进行开发时,使用geotools的GeoTiffReader类读取卫星影像的顶点坐标,得到的是平面坐标,单位是米,打开arcGIS,加载影像,右下角会有坐标,在视图-& ...
地理坐标（经纬度）转换成投影坐标（XY坐标）
前言:限于需求,项目中所有涉及到经纬度的字段都要转换成XY坐标,面向度娘之后发现都没有可用的,琢磨了之后在以为大佬博客中发现了宝藏使用工具:Proj4 Proj4是一个JavaScript类库,其主 ...
学习笔记——【python】GetGeoTransform()使用，gdal截取图像，使用GDAL进行影像投影坐标、地理坐标、图上坐标的转换
1. GetGeoTransform()使用.gdal截取图像 GetGeoTransform() GeoTransform[0],左上角横坐标(应该是投影坐标) GeoTransform[2],行旋 ...
图像坐标、地理坐标和投影坐标的相互转化（Python版）
本博客系转载~ from osgeo import gdal from osgeo import osr import numpy as npdef getSRSPair(dataset):'''获得 ...
如何判断投影坐标是 3 度带还是 6 度带？如何计算中央子午线经度？
在做坐标转换时通常会涉及到:3 度带.6 度带.中央子午线经度的计算,下面使用两种方法对其进行展开计算: 1.什么是3°带? 三度带就是地球的一周(360度)按照经度方向三度一个划分为120份: 2. ...
supermap iobjects 组件设置坐标，投影坐标之间的转换
组件 :supermap iobjects java 功能使用: 1.设置地理坐标 2.投影坐标转投影坐标 3.地理坐标转投影坐标 4.投影坐标转地理坐标注:注释的根据需要使用 public Res ...
python dataset[trans_科学网—Python GDAL 图像坐标，投影坐标，经纬度坐标三者映射及运行错误解决 - 吴妍潼的博文...
题记: 写该博客是因为自己经常遇到这个问题,而我发现网络上关于这方面浏览量高的一些代码竟然都有误,每次照搬都被虐得很惨.有一些同志在某些博客下方留言说代码有问题,而博主没有回应,也没有更改错误.为了自 ...
python运行报错无法映射_科学网-Python GDAL 图像坐标，投影坐标，经纬度坐标三者映射及运行错误解决-吴妍潼的博文...
题记: 写该博客是因为自己经常遇到这个问题,而我发现网络上关于这方面浏览量高的一些代码竟然都有误,每次照搬都被虐得很惨.有一些同志在某些博客下方留言说代码有问题,而博主没有回应,也没有更改错误.为了自 ...
(Python)卫星RPC有理多项式模型读取与正反投影坐标计算原理与实现
(Python)卫星RPC有理多项式模型读取与正反投影坐标计算原理与实现文章目录 (Python)卫星RPC有理多项式模型读取与正反投影坐标计算原理与实现摘要 RPC几何定位模型介绍 RPC模型库 ...

python——经纬度坐标和平面投影坐标的相互转换

文章目录

经纬度坐标系和平面投影坐标系

地理坐标系 / 大地坐标系

平面投影坐标系

墨卡托投影法

变量和常量定义

（1）墨卡托投影正解算公式

（2）墨卡托投影逆解算公式

Python代码实现

（1）墨卡托投影正解算公式

（2）墨卡托投影逆解算公式

python——经纬度坐标和平面投影坐标的相互转换相关推荐

最新文章

热门文章