互联网常见坐标之间的转换（Python）

最新推荐文章于 2024-08-08 07:08:17 发布
置顶 最新推荐文章于 2024-08-08 07:08:17 发布
阅读量1.1k 收藏 6
点赞数
分类专栏： 地理信息
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/mymhj/article/details/117027453
版权

1.场景描述

随着互联网的兴趣，在企业应用系统的开发过程中，几乎每一个APP（CS、BS、移动端）都会搜集用户或设备的位置数据，然后在相关地图（来自互联网电子地图或自建地图服务）上进行标注展示，如果对空间坐标系不是很了解的话，不管你Coding能力有多强都会被各种位置匹配问题（错误、偏离、飞出等）折磨的晕头转向，不知所措。首先针对几种常见的互联网坐标我们对其做一简单场景梳理

常见的互联网坐标：

地球坐标：WGS84，国际标准、从GPS设备获取的坐标的参考坐标系、国际地图提供商使用的坐标系。OSM、谷歌外国、ArcGISonline

火星坐标：GCJ-02，中国标准、也叫国测局坐标，从国行移动设备中定位获取的坐标的参考坐标系。国内出版的各种地图系统，必须至少采用GCJ-02对地理位置进行首次加密。高德地图、搜搜地图、阿里云地图、腾讯地图。

百度坐标：DB09，百度标准。百度地图、百度空间数据的参考坐标系都使用此坐标系。是在GCJ-02坐标系的基础上做了二次加密

 

备注：天地图使用的是CGC2000，参考椭球非常接近。当前阶段（星历变化可能导致差异变大），变率差异引的椭球面上的维度和高程变化量不大0.1mm,两者相容误差在cm级别，在互联网应用中可以将CGC2000就当作WGS84来用。

 

针对上述互联网坐标系，我将转换场景梳理为以下六种：

WGS84转GCJ02：WGS84_To_GCJ02

GCJ02转WGS84：GCJ02_To_WGS84

GCJ02转BD09：GCJ02_To_BD09

BD09转GCJ02：BD09_To_GCJ02

WGS84转BD09：WGS84_To_BD09

BD09转WGS84：BD09_To_WGS84

 

2.功能实现

 

     参考资料

​​​​​​​

    代码实现

开发工具：VScode    使用语言：Python

只实现经纬度坐标之间的转换，支持点、线、面。Web墨卡托的转换目前没设计到没有实现

 

# -*- coding: utf-8 -*-
import arcpy
import math


# define ellipsoid
x_pi = 3.14159265358979324 * 3000.0 / 180.0
pi = 3.1415926535897932384626  # π
a = 6378245.0  # 长半轴
ee = 0.00669342162296594323  # 偏心率平方


def _transformlat(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 * pi) + 20.0 *
            math.sin(2.0 * lng * pi)) * 2.0 / 3.0
    ret += (20.0 * math.sin(lat * pi) + 40.0 *
            math.sin(lat / 3.0 * pi)) * 2.0 / 3.0
    ret += (160.0 * math.sin(lat / 12.0 * pi) + 320 *
            math.sin(lat * pi / 30.0)) * 2.0 / 3.0
    return ret


def _transformlng(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 * pi) + 20.0 *
            math.sin(2.0 * lng * pi)) * 2.0 / 3.0
    ret += (20.0 * math.sin(lng * pi) + 40.0 *
            math.sin(lng / 3.0 * pi)) * 2.0 / 3.0
    ret += (150.0 * math.sin(lng / 12.0 * pi) + 300.0 *
            math.sin(lng / 30.0 * pi)) * 2.0 / 3.0
    return ret


def out_of_china(lng, lat):
    """
    判断是否在国内，不在国内不做偏移
    :param lng:
    :param lat:
    :return:
    """
    return not (lng > 73.66 and lng < 135.05 and lat > 3.86 and lat < 53.55)


def WGS84_To_GCJ02(lng, lat):
    """
    WGS84转GCJ02(火星坐标系)
    :param lng:WGS84坐标系的经度
    :param lat:WGS84坐标系的纬度
    :return:
    """
    if out_of_china(lng, lat):  # 判断是否在国内
        return [lng, lat]
    dlat = _transformlat(lng - 105.0, lat - 35.0)
    dlng = _transformlng(lng - 105.0, lat - 35.0)
    radlat = lat / 180.0 * pi
    magic = math.sin(radlat)
    magic = 1 - ee * magic * magic
    sqrtmagic = math.sqrt(magic)
    dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi)
    dlng = (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * pi)
    mglat = lat + dlat
    mglng = lng + dlng
    return arcpy.Point(mglng, mglat)


def GCJ02_To_WGS84(lng, lat):
    """
    GCJ02(火星坐标系)转GPS84
    :param lng:火星坐标系的经度
    :param lat:火星坐标系纬度
    :return:
    """
    if out_of_china(lng, lat):
        return [lng, lat]
    dlat = _transformlat(lng - 105.0, lat - 35.0)
    dlng = _transformlng(lng - 105.0, lat - 35.0)
    radlat = lat / 180.0 * pi
    magic = math.sin(radlat)
    magic = 1 - ee * magic * magic
    sqrtmagic = math.sqrt(magic)
    dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi)
    dlng = (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * pi)
    mglat = lat + dlat
    mglng = lng + dlng
    return arcpy.Point(lng * 2 - mglng,lat * 2 - mglat)


def GCJ02_to_BD09(lng, lat):
    """
    火星坐标系(GCJ-02)转百度坐标系(BD-09)
    谷歌、高德——>百度
    :param lng:火星坐标经度
    :param lat:火星坐标纬度
    :return:
    """
    z = math.sqrt(lng * lng + lat * lat) + 0.00002 * math.sin(lat * x_pi)
    theta = math.atan2(lat, lng) + 0.000003 * math.cos(lng * x_pi)
    bd_lng = z * math.cos(theta) + 0.0065
    bd_lat = z * math.sin(theta) + 0.006
    return arcpy.Point(bd_lng, bd_lat)


def BD09_to_GCJ02(bd_lon, bd_lat):
    """
    百度坐标系(BD-09)转火星坐标系(GCJ-02)
    百度——>谷歌、高德
    :param bd_lat:百度坐标纬度
    :param bd_lon:百度坐标经度
    :return:转换后的坐标列表形式
    """
    x = bd_lon - 0.0065
    y = bd_lat - 0.006
    z = math.sqrt(x * x + y * y) - 0.00002 * math.sin(y * x_pi)
    theta = math.atan2(y, x) - 0.000003 * math.cos(x * x_pi)
    gg_lng = z * math.cos(theta)
    gg_lat = z * math.sin(theta)
    return arcpy.Point(gg_lng, gg_lat)


def bd09_to_wgs84(bd_lon, bd_lat):
    newPoint = arcpy.Point()
    newPoint = BD09_to_GCJ02(bd_lon, bd_lat)
    return GCJ02_To_WGS84(newPoint.X, newPoint.Y)


def WGS84_To_BD09(lon, lat):
    newPoint = arcpy.Point()
    newPoint = WGS84_To_GCJ02(lon, lat)
    return GCJ02_to_BD09(newPoint.X, newPoint.Y)


def TransPoint(geom, transType):
    newPoint = arcpy.Point()
    oldPoint = geom.getPart(0)
    if transType == 'WGS84_To_GCJ02':
        newPoint = WGS84_To_GCJ02(oldPoint.X, oldPoint.Y)
    elif transType == 'GCJ02_To_WGS84':
        newPoint = GCJ02_To_WGS84(oldPoint.X, oldPoint.Y)
    elif transType == 'GCJ02_to_BD09':
        newPoint = GCJ02_to_BD09(oldPoint.X, oldPoint.Y)
    elif transType == 'BD09_to_GCJ02':
        newPoint = BD09_to_GCJ02(oldPoint.X, oldPoint.Y)
    elif transType == 'WGS84_To_BD09':
        newPoint = WGS84_To_BD09(oldPoint.X, oldPoint.Y)
    elif transType == 'BD09_To_WGS84':
        newPoint = bd09_to_wgs84(oldPoint.X, oldPoint.Y)
    return newPoint


def TransPolyLine(geom, transType):
    newPoint = arcpy.Point()
    newArray = arcpy.Array()
    array = geom.getPart(0)
    for i in range(0, array.count):
        oldPoint = array[i]
        if transType == 'WGS84_To_GCJ02':
            newPoint = WGS84_To_GCJ02(oldPoint.X, oldPoint.Y)
        elif transType == 'GCJ02_To_WGS84':
            newPoint = GCJ02_To_WGS84(oldPoint.X, oldPoint.Y)
        elif transType == 'GCJ02_to_BD09':
            newPoint = GCJ02_to_BD09(oldPoint.X, oldPoint.Y)
        elif transType == 'BD09_to_GCJ02':
            newPoint = BD09_to_GCJ02(oldPoint.X, oldPoint.Y)
        elif transType == 'WGS84_To_BD09':
            newPoint = WGS84_To_BD09(oldPoint.X, oldPoint.Y)
        elif transType == 'BD09_To_WGS84':
            newPoint = bd09_to_wgs84(oldPoint.X, oldPoint.Y)
        newArray.add(newPoint)
    newLine = arcpy.Polyline(newArray, SR)
    newArray.removeAll()
    return newLine


def TransPolygon(geom, transType):
    newPoint = arcpy.Point()
    newArray = arcpy.Array()
    array = geom.getPart(0)
    for i in range(0, array.count):
        oldPoint = array[i]
        if transType == 'WGS84_To_GCJ02':
            newPoint = WGS84_To_GCJ02(oldPoint.X, oldPoint.Y)
        elif transType == 'GCJ02_To_WGS84':
            newPoint = GCJ02_To_WGS84(oldPoint.X, oldPoint.Y)
        elif transType == 'GCJ02_to_BD09':
            newPoint = GCJ02_to_BD09(oldPoint.X, oldPoint.Y)
        elif transType == 'BD09_to_GCJ02':
            newPoint = BD09_to_GCJ02(oldPoint.X, oldPoint.Y)
        elif transType == 'WGS84_To_BD09':
            newPoint = wgs84_to_bd09(oldPoint.X, oldPoint.Y)
        elif transType == 'BD09_To_WGS84':
            newPoint = bd09_to_wgs84(oldPoint.X, oldPoint.Y)
        newArray.add(newPoint)
    newLine = arcpy.Polyline(newArray, SR)
    newArray.removeAll()
    return newLine


# 输入工作空间
# in_Features = arcpy.GetParameter(0)
# out_Feature_FullPath = arcpy.GetParameter(1)
# transType = arcpy.GetParameter(2)
# arcpy.AddMessage(convert_Type)
transType = "WGS84_To_GCJ02"
inShpFile = "E:\\OutputData\\WGS84_R.shp"
outShpFile = "E:\\OutputData\\WGS_GCJ02_R.shp"
if arcpy.Exists(outShpFile):
    arcpy.Delete_management(outShpFile)
arcpy.CopyFeatures_management(inShpFile, outShpFile)
arcpy.RepairGeometry_management(outShpFile)
# arcpy.AddMessage(out_ShpFile)
cur = arcpy.UpdateCursor(outShpFile)
des = arcpy.Describe(outShpFile)
SR = des.spatialReference
if des.shapeType.upper() == 'POINT':
    for r in cur:
        geom = r.getValue("SHAPE")
        r.setValue("SHAPE", TransPoint(geom, transType))
        cur.updateRow(r)
    del r, cur
elif des.shapeType.upper() == 'POLYLINE':
    for r in cur:
        geom = r.getValue("SHAPE")
        newgeom = TransPolyLine(geom, transType)
        r.setValue("SHAPE", newgeom)
        cur.updateRow(r)
    del r, cur
elif des.shapeType.upper() == 'POLYGON':
    for r in cur:
        geom = r.getValue("SHAPE")
        r.setValue("SHAPE", TransPolygon(geom, transType))
        cur.updateRow(r)
    del r, cur

    效果验证

遥哥_GB_Fore
关注
  • 0
    点赞
  • 踩
  • 6
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
【Python数据分析与可视化】期末复习试题集
风口IT猪的成长录
12-03 2万+
1.Jupyter notebook不具备的功能是 （） Jupyter notebook可以直接生成一份交互式文档 Jupyter notebook可以安装Python库 Jupyter notebook可以导出HTML文件 Jupyter notebook可以将文件分享给他人 right_answers: Jupyter notebook可以将文件分享给他人 2.在Jupyter notebook的命令模式下，要查看所有快捷..
Python各坐标系之间的转换（高斯、WGS84、Web墨卡托、瓦片坐标）
u010410697的博客
06-27 1万+
本文整理了一些地理坐标系之间的转换（python代码） pip install pyproj import math from pyproj import CRS from pyproj import Transformer from pyproj import _datadir, datadir # 加上这个linux打包执行不报错 crs_WGS84 = CRS.from_epsg(4326) # WGS84地理坐标系 crs_WebMercator = CRS.from_epsg(3857)
Python -- 地理/投影坐标系转换
zhangdm的博客
05-15 3672
文章目录Geojson格式坐标转换投影坐标系和地理坐标系使用CRS库转换坐标系（以CGCS2000_3_Degree_GK_Zone_40转WGS84为例）大地两千坐标系配置文件 当我们要在ArcGis中绘制一个图形的时候，我们需要确认三点 1、坐标系参考规则 2、要绘制的点集合 3、点的顺序 有了以上三点，我们就可以在指定的坐标系中，按照点的顺序将点连起来，绘制成一个图形。 在实际开发过程中，我们的的数据来源可能是来源于其他网站，我们所获取的数据可能不能直接满足展示需求，例如我们爬取的网站绘图的坐标系与我
coordTransform_py 使用教程
最新发布
gitblog_00851的博客
08-08 328
coordTransform_py 使用教程 coordTransform_py提供百度坐标系(bd-09)、火星坐标系(国测局坐标系、gcj02)、WGS84坐标系直接的坐标互转，也提供了解析高德地址的方法的python版本项目地址:https://gitcode.com/gh_mirrors/co/coordTransform_py 1. 项目目录结构及介绍 协调变换库 coordTrans...
python坐标系转换库_转载：python库Pyproj进行坐标转换
weixin_40001967的博客
12-09 1146
https://www.baidu.com利用Pyproj进行坐标转换作者：郜科科两个坐标系统的参考椭球不同，实地一个点的不同坐标系的值是不同的，不同的部门采用的坐标系统经常是不一致，所以要转换后才能相互利用。例如目前使用的北京市观测站点位置根据GPS的定位而来，GPS使用的地理坐标系为GCS_WGS_1984，所以其坐标的地理坐标系也为GCS_WGS_1984，而假如需要将这些点显示在Web端的...
python三维坐标转二维_2D坐标系与3D坐标系的相互转换
weixin_39761422的博客
12-13 2119
并不是做关于SLAM方向的，但由于某些任务涉及到，故作此笔记～相机内参矩阵：不同的的深度摄像头具有不同的特征参数，在计算机视觉里，将这组参数设置为相机的内参矩阵C:$$begin{bmatrix} f_x& 0 &c_x \ 0 & f_y & c_y \ 0 & 0 & 1 end{bmatrix}$$fx，fy指相机在x轴和y轴上的焦距，cx，...
实习一 WGS84经纬度坐标与Web墨卡托投影坐标相互转换1
08-08
在IT领域，坐标转换是一种常见的任务，特别是在地理信息系统（GIS）和地图绘制中。本实习主要涉及两种坐标系统之间的转换：WGS84经纬度坐标和Web墨卡托投影坐标。WGS84（World Geodetic System 1984）是一种全球通用...
根据地名地址爬取百度坐标，同时转换坐标系为WGS1984.zip
11-23
在IT领域，地理信息系统（GIS）的坐标转换是一项常见的任务，尤其对于数据分析和地图应用来说至关重要。本项目涉及的关键知识点是地名地址的网络爬虫技术、百度地图API的使用以及坐标系统的转换。 首先，我们要了解...
python实现已知两个点坐标求与正北夹角度数.pdf
09-14
标题中的“python实现已知两个点坐标求与正北夹角度数”是一个常见的地理计算问题，通常在导航、地图应用或者游戏开发中会遇到。这个问题的目的是计算从一个点（起点）沿着直线方向到另一个点（终点）时，这条直线与...
python gps与高德地图poi_poi: 百度以及高德地图POI数据和POI边界坐标数据(目前仅限百度)爬取代码，采用Python编写...
weixin_33642697的博客
03-01 705
这是一个爬取互联网开放的POI数据以及POI边界坐标代码集，包括高德和百度地图，可以爬取指定经纬度范围或者指定城市/省份的POI数据以及POI边界坐标(目前仅限百度)。相关网站代码目录介绍baidu百度地图POI数据爬取(WGS84坐标系)app.py 程序入口文件data 存放爬取的POI数据目录boundary 百度POI边界坐标爬取目录POI边界爬取结果存放目录gaode高德地图POI数据爬...
coordTransform_py-master.zip
09-08
基于python语言的坐标转换代码，可以转换百度坐标，谷歌坐标，火星坐标，实测很准确。坐标转换工具。
python 地图经纬度转换、纠偏
qq_40771567的博客
08-06 2731
# -*- coding: utf-8 -*- import json import urllib import math x_pi = 3.14159265358979324 * 3000.0 / 180.0 pi = 3.1415926535897932384626 # π a = 6378245.0 # 长半轴 ee = 0.00669342162296594323 # 偏心率平方...
使用Python和GDAL给图片进行坐标系转换（GCJ-02转为WGS84）
qq_41440099的博客
03-27 1561
背景 GCJ-02坐标系是最常见的地图坐标系，如何使用gdal把GCJ-02坐标系转为WGS84呢？ 原理 通过已知GCJ-02坐标转换WGS84坐标，添加控制点进行几何校正。 步骤 1.坐标点转换 使用网友的gcj02_to_wgs84函数，（经测试该函数会导致坐标值精度丢失，即相同坐标转换后小数部分数值不一样） https://blog.csdn.net/kaikai136412162/article/details/88894892 2.几何校正 转换后的坐标点作为控制点，进行校正：见 https:/
C#调用python实现shp文件的CGCS2000坐标系向WGS84坐标系转换
AIGIS的博客
07-03 1480
前言：我们项目上最近有个需求，要把shp文件导入到我们的webgis系统中来，展示在前端的网页上。如果是任意的shp想展示在webgis中，我们就使用postgresql提供的ogr2ogr工具，直接将shp转为geojson，数据传给前端就行了。但是现在我们的shp是有固定格式的，我们的数据库中也有对应格式的表，因此问题就转变为了，将shp导入我们的表中。 s...
探秘地理坐标转换神器：coordTransform_py
gitblog_00091的博客
04-10 795
探秘地理坐标转换神器：coordTransform_py coordTransform_py提供百度坐标系(bd-09)、火星坐标系(国测局坐标系、gcj02)、WGS84坐标系直接的坐标互转，也提供了解析高德地址的方法的python版本项目地址:https://gitcode.com/gh_mirrors/co/coordTransform_py 在现代GIS（地理信息系统）应用中，准确无误地...
技术文档|利用Python进行地理坐标系统的转换（WGS-84、GCJ-02、BAIDU-09三者之间转换）
chenxu0136的博客
08-04 8278
本文主要讲述地理坐标系统的原理以及怎么利用Python进行地理坐标系统转换，主要内容包含以下几块： 1. 什么是地理坐标系统？ 2. 常用的地理坐标系统有哪些？常用地图产品分别是什么地理坐标系？ 3. 怎么样利用Python实现地理坐标系的转换？ 4. 介绍一个能够不写代码就进行地理坐标系转换的软件。
GIS数据格式坐标转换（地球坐标WGS84、GCJ-02、火星坐标、百度坐标BD-09、国家大地坐标系CGCS2000）
热门推荐
漏刻有时数据可视化大屏（PHP&ECHARTS智能化开源软件系统）
02-12 1万+
地理信息系统 (GIS) 是一个创建、管理、分析和绘制所有类型数据的系统。GIS 将数据连接到地图，将位置数据（事物所在位置）与所有类型的描述性信息（事物在该位置的情况）集成到一起。这可以为适用于自然科学和几乎所有行业的制图和分析提供基础。GIS 帮助用户了解模式、关系和地理环境。其优势包括改善沟通、提高效率以及更好地管理和决策。解决使用百度webAPI坐标转换时的问题：X→GPS根据法律规定，不支持将任何类型的坐标转换为GPS坐标；
Python坐标系转换
diligent_lee的博客
03-21 8432
1. 相关名词 地心地固坐标系（Earth-Centered, Earth-Fixed，ECEF），简称地心坐标系。 地理坐标系统（Geographic Coordinate System，GCS），坐标系是地心坐标系，用经纬度表示球面上的点。 世界大地测量系统（World Geodetic System, WGS），比如WGS84，是一种地理坐标系统，用于全球定位系统（GPS）。 投影坐标系统（Projection Coordinate System，PCS），在二维平面上用米表示位置。 通用横轴墨卡托
python地图坐标系转换(bd09,gcj02,wgs84三种投影坐标系相互转化)
weixin_38468077的博客
04-14 7572
1.介绍 1.1 GIS之坐标系 坐标系是GIS的重中之重，一般来说，工作底图平面坐标系应采用国家大地坐标系CGCS2000(或相当于精度WGS84坐标系)，投影方式采用高斯-克吕格投影，高程基准采用1985国家高程基准。 1.2 地理坐标系（GCS，Geographic Coordinate System） 地理坐标系其实是用了一个规则的球面来代表地球表面。在球面上画一张经纬网，球面上的点就有了它的经纬度，这就是地球上每一点的坐标。因此，地球坐标系中的坐标是以经纬度来表示的。 1.3 我国常见的GCS-地
地理坐标系转换 python
07-29
在 Python 中，可以使用多种库来进行地理坐标系的转换，其中最常用的是 `pyproj` 和 `geopy`。下面是使用这两个库进行地理坐标系转换的示例代码： 使用 `pyproj` 库进行坐标系转换的示例代码： ```python from pyproj import Transformer # 定义输入和输出坐标系 input_crs = 'EPSG:4326' # WGS84 坐标系 output_crs = 'EPSG:3857' # Web Mercator 坐标系 # 创建转换器 transformer = Transformer.from_crs(input_crs, output_crs, always_xy=True) # 输入经纬度坐标 lon, lat = 116.397458, 39.909716 # 坐标系转换 x, y = transformer.transform(lon, lat) # 输出结果 print(f'转换后的坐标：{x}, {y}') ``` 使用 `geopy` 库进行坐标系转换的示例代码： ```python from geopy.point import Point # 定义输入和输出坐标系 input_crs = 'EPSG:4326' # WGS84 坐标系 output_crs = 'EPSG:3857' # Web Mercator 坐标系 # 创建输入点对象 input_point = Point(39.909716, 116.397458, 0) # 输入经纬度坐标（纬度在前，经度在后） # 坐标系转换 output_point = input_point.to(output_crs) # 输出结果 print(f'转换后的坐标：{output_point.latitude}, {output_point.longitude}') ``` 以上代码中，`EPSG:4326` 表示 WGS84 坐标系，也就是经纬度坐标系，`EPSG:3857` 表示 Web Mercator 坐标系，也就是墨卡托投影坐标系。你可以根据你的实际需求选择不同的输入和输出坐标系。
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值