Python实现常见的4种坐标互相转换
作者:费弗里 发布时间:2023-11-21 06:22:38
标签:Python,坐标,转换
一、简介
主流被使用的地理坐标系并不统一,常用的有WGS84、GCJ02(火星坐标系)、BD09(百度坐标系)以及百度地图中保存矢量信息的web墨卡托,本文利用Python编写相关类以实现4种坐标系统之间的互相转换。
二、代码及说明
import math
class LngLatTransfer():
def __init__(self):
self.x_pi = 3.14159265358979324 * 3000.0 / 180.0
self.pi = math.pi # π
self.a = 6378245.0 # 长半轴
self.es = 0.00669342162296594323 # 偏心率平方
pass
def GCJ02_to_BD09(self, gcj_lng, gcj_lat):
"""
实现GCJ02向BD09坐标系的转换
:param lng: GCJ02坐标系下的经度
:param lat: GCJ02坐标系下的纬度
:return: 转换后的BD09下经纬度
"""
z = math.sqrt(gcj_lng * gcj_lng + gcj_lat * gcj_lat) + 0.00002 * math.sin(gcj_lat * self.x_pi)
theta = math.atan2(gcj_lat, gcj_lng) + 0.000003 * math.cos(gcj_lng * self.x_pi)
bd_lng = z * math.cos(theta) + 0.0065
bd_lat = z * math.sin(theta) + 0.006
return bd_lng, bd_lat
def BD09_to_GCJ02(self, bd_lng, bd_lat):
'''
实现BD09坐标系向GCJ02坐标系的转换
:param bd_lng: BD09坐标系下的经度
:param bd_lat: BD09坐标系下的纬度
:return: 转换后的GCJ02下经纬度
'''
x = bd_lng - 0.0065
y = bd_lat - 0.006
z = math.sqrt(x * x + y * y) - 0.00002 * math.sin(y * self.x_pi)
theta = math.atan2(y, x) - 0.000003 * math.cos(x * self.x_pi)
gcj_lng = z * math.cos(theta)
gcj_lat = z * math.sin(theta)
return gcj_lng, gcj_lat
def WGS84_to_GCJ02(self, lng, lat):
'''
实现WGS84坐标系向GCJ02坐标系的转换
:param lng: WGS84坐标系下的经度
:param lat: WGS84坐标系下的纬度
:return: 转换后的GCJ02下经纬度
'''
dlat = self._transformlat(lng - 105.0, lat - 35.0)
dlng = self._transformlng(lng - 105.0, lat - 35.0)
radlat = lat / 180.0 * self.pi
magic = math.sin(radlat)
magic = 1 - self.es * magic * magic
sqrtmagic = math.sqrt(magic)
dlat = (dlat * 180.0) / ((self.a * (1 - self.es)) / (magic * sqrtmagic) * self.pi)
dlng = (dlng * 180.0) / (self.a / sqrtmagic * math.cos(radlat) * self.pi)
gcj_lng = lat + dlat
gcj_lat = lng + dlng
return gcj_lng, gcj_lat
def GCJ02_to_WGS84(self, gcj_lng, gcj_lat):
'''
实现GCJ02坐标系向WGS84坐标系的转换
:param gcj_lng: GCJ02坐标系下的经度
:param gcj_lat: GCJ02坐标系下的纬度
:return: 转换后的WGS84下经纬度
'''
dlat = self._transformlat(gcj_lng - 105.0, gcj_lat - 35.0)
dlng = self._transformlng(gcj_lng - 105.0, gcj_lat - 35.0)
radlat = gcj_lat / 180.0 * self.pi
magic = math.sin(radlat)
magic = 1 - self.es * magic * magic
sqrtmagic = math.sqrt(magic)
dlat = (dlat * 180.0) / ((self.a * (1 - self.es)) / (magic * sqrtmagic) * self.pi)
dlng = (dlng * 180.0) / (self.a / sqrtmagic * math.cos(radlat) * self.pi)
mglat = gcj_lat + dlat
mglng = gcj_lng + dlng
lng = gcj_lng * 2 - mglng
lat = gcj_lat * 2 - mglat
return lng, lat
def BD09_to_WGS84(self, bd_lng, bd_lat):
'''
实现BD09坐标系向WGS84坐标系的转换
:param bd_lng: BD09坐标系下的经度
:param bd_lat: BD09坐标系下的纬度
:return: 转换后的WGS84下经纬度
'''
lng, lat = self.BD09_to_GCJ02(bd_lng, bd_lat)
return self.GCJ02_to_WGS84(lng, lat)
def WGS84_to_BD09(self, lng, lat):
'''
实现WGS84坐标系向BD09坐标系的转换
:param lng: WGS84坐标系下的经度
:param lat: WGS84坐标系下的纬度
:return: 转换后的BD09下经纬度
'''
lng, lat = self.WGS84_to_GCJ02(lng, lat)
return self.GCJ02_to_BD09(lng, lat)
def _transformlat(self, lng, lat):
ret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + \
0.1 * lng * lat + 0.2 * math.sqrt(math.fabs(lng))
ret += (20.0 * math.sin(6.0 * lng * self.pi) + 20.0 *
math.sin(2.0 * lng * self.pi)) * 2.0 / 3.0
ret += (20.0 * math.sin(lat * self.pi) + 40.0 *
math.sin(lat / 3.0 * self.pi)) * 2.0 / 3.0
ret += (160.0 * math.sin(lat / 12.0 * self.pi) + 320 *
math.sin(lat * self.pi / 30.0)) * 2.0 / 3.0
return ret
def _transformlng(self, lng, lat):
ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + \
0.1 * lng * lat + 0.1 * math.sqrt(math.fabs(lng))
ret += (20.0 * math.sin(6.0 * lng * self.pi) + 20.0 *
math.sin(2.0 * lng * self.pi)) * 2.0 / 3.0
ret += (20.0 * math.sin(lng * self.pi) + 40.0 *
math.sin(lng / 3.0 * self.pi)) * 2.0 / 3.0
ret += (150.0 * math.sin(lng / 12.0 * self.pi) + 300.0 *
math.sin(lng / 30.0 * self.pi)) * 2.0 / 3.0
return ret
def WGS84_to_WebMercator(self, lng, lat):
'''
实现WGS84向web墨卡托的转换
:param lng: WGS84经度
:param lat: WGS84纬度
:return: 转换后的web墨卡托坐标
'''
x = lng * 20037508.342789 / 180
y = math.log(math.tan((90 + lat) * self.pi / 360)) / (self.pi / 180)
y = y * 20037508.34789 / 180
return x, y
def WebMercator_to_WGS84(self, x, y):
'''
实现web墨卡托向WGS84的转换
:param x: web墨卡托x坐标
:param y: web墨卡托y坐标
:return: 转换后的WGS84经纬度
'''
lng = x / 20037508.34 * 180
lat = y / 20037508.34 * 180
lat = 180 / self.pi * (2 * math.atan(math.exp(lat * self.pi / 180)) - self.pi / 2)
return lng, lat
整个模块的使用方式可用下面的导图概括,其中每个函数都只需要传入经纬度坐标信息:
来源:https://www.cnblogs.com/feffery/p/11023673.html
0
投稿
猜你喜欢
- image.jsp------------------------------生成随机验证码图片的Jsp页面 代码如下: <
- 首先我们介绍使用自己的数据库查询多个手机号码,那还是建议你拥有一个自己的的手机号码数据库。正常情况下,只是满足一般查询的话,你不需要去购买专
- 在 玉伯 的文章 《一道大题目,嘿嘿》 中有这样一段代码:[] == ![]也许很多同学迷惑:咦,这个如何转换呢?首先,我们了解下逻辑 NO
- 第一个示例: 简单的jsp自定标签获取内容: 首先创建一个jsp实例类然后继承SimpleTagSupport类 然后实现父类的doTag(
- 下拉框连动JavaScript代码,市区二级联动多级联动下拉选择框,动态获取下一级琥珀无限级联动菜单-JavaScript版 <htm
- 支持CSS属性Safari和WebKit实施大子的CSS 2.1规格所界定的万维网联盟( W3C ) ,以及部分的CSS 3规格。 。这个C
- 注意:安装时要保证Oracle安装目录不能带有中文字符(如果第一次安装出现“加载数据库错误areasQueries”的错误,一般是因为Ora
- 一、两层结构的ASP应用有何缺点 在Browser/Server 应用程序开发领域,微软公司的IIS/ASP以其强大的功能,良好的扩展能力,
- 平常我们只听说过ADO等一类ASP对象,但在ASP中还有个鲜为人知的专门SQL Server的ASP访问对象,它就是SQLOLE.SQLSe
- 首先声明,在这组里我是个绝对的菜鸟。再次声明,小爝这个菜鸟在“网页设计”这个圈里混了快1年了。 摘要:我知道我有多少底,所以我在总结我的成长
- 本文实例为大家分享了js实现黑白div块画空心图形的具体代码,供大家参考,具体内容如下<!DOCTYPE html><ht
- 内容摘要:近年来,视频点播VOD(Video on Demand)的名字在媒体上出现得越来越多。VOD技术使人们可以根据
- 为什么要页面静态化?1.动态文件执行过程:语法分析-编译-运行2.静态文件,不需要编译,减少了服务器脚本运行的时间,降低了服务器的响应时间,
- 无头模式添加,可以让selenium模拟登录,进入到后台运行这里以登录打开公司内网下载数据为例,因为涉及私密问题,所以有些地方我们进行覆盖,
- 1,file_put_contents()函数 2,使用PHP内置缓存机制实现页面静态化:output_bufferingphp中outpu
- 在后台添加中一般我们都设置了一个“关键字”,而且一般输入可能如下: 娃娃,毛绒娃娃,日本充气娃娃 我们将根据这个来做,具体实现如下: 先看一
- 此文译自Fred Wilson 2010年2月在迈阿密举行的Web未来应用的年会上的演讲谢谢青云推荐了这篇这么好的演说谢谢卓和百忙中抽空帮我
- 从我们论坛中收集了这段HTML制作页面需要最大化、最小化时可以借鉴参考。最大化效果:<OBJECT id="max
- PHP原型模式Prototype Pattern是什么原型模式是一种创建型模式,它可以通过复制现有对象来创建新的对象,而无需知道具体的创建过
- 我一直很难理解Javascript语言的继承机制。它没有"子类"和"父类"的概念,也没有"