常见GIS坐标系的坐标转换

最新推荐文章于 2024-05-21 22:44:12 发布
最新推荐文章于 2024-05-21 22:44:12 发布
阅读量305 收藏 8
点赞数 3
分类专栏： Java 文章标签： python java
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xhp312098226/article/details/138225964
版权

几种常见的 GIS 坐标系的简单介绍及坐标数据的相互转换方法。

一、坐标系介绍

1. WGS-84：大地坐标系

  • 国际通用坐标系，也叫地球坐标系。
  • 各 GPS 设备、境外地图（包括但不限于谷歌地图、高德地图、百度地图的境外版）等使用的坐标系。

2. GCJ-02：国测局坐标系

  • 国内通用坐标系，也叫火星坐标系。
  • 各境内地图（高德地图、腾讯地图等）等使用的坐标系。

3. BD-09：百度坐标系

  • GCJ-02 基础上的二次加密。
  • 百度标准，百度 SDK，百度地图，百度 GeoCoding 等使用的坐标系。

二、坐标数据转换方法-Java 版

1. 准备工作

1.1 定义依赖的一些常量和转换方法

    private static final double PI = Math.PI;
    private static final double X_PI = PI * 3000.0 / 180.0;
    private static final double AA = 6378245.0;
    private static final double EE = 0.00669342162296594323;

    /**
     * 经纬度偏移转换-经度
     *
     * @param lng
     * @param lat
     * @return
     */
    private static double transformLng(double lng, double lat) {
        double ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + 0.1 * lng * lat + 0.1 * Math.sqrt(Math.abs(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;
    }

    /**
     * 经纬度偏移转换-维度
     *
     * @param lng
     * @param lat
     * @return
     */
    private static double transformLat(double lng, double lat) {
        double ret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + 0.1 * lng * lat + 0.2 * Math.sqrt(Math.abs(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;
    }

    /**
     * 判断坐标是否超出境内
     *
     * @param lng
     * @param lat
     * @return
     */
    private static boolean outOfChina(double lng, double lat) {
        // 纬度3.86~53.55,经度73.66~135.05；非精确判断，可能误判
        return !(73.66 < lng && lng < 135.05 && 3.86 < lat && lat < 53.55);
    }

2. WGS-84 和 GCJ-02 的相互转换

2.1 WGS-84 转 GCJ-02

    /**
     * WGS-84转GCJ-02
     *
     * @param lng
     * @param lat
     * @return
     */
    public static double[] wgs84ToGcj02(double lng, double lat) {
        if (outOfChina(lng, lat)) {
            return new double[]{lng, lat};
        } else {
            double dlat = transformLat(lng - 105.0, lat - 35.0);
            double dlng = transformLng(lng - 105.0, lat - 35.0);
            double radlat = lat / 180.0 * PI;
            double magic = Math.sin(radlat);
            magic = 1 - EE * magic * magic;
            double sqrtMagic = Math.sqrt(magic);
            dlat = (dlat * 180.0) / ((AA * (1 - EE)) / (magic * sqrtMagic) * PI);
            dlng = (dlng * 180.0) / (AA / sqrtMagic * Math.cos(radlat) * PI);
            double mglat = lat + dlat;
            double mglng = lng + dlng;
            return new double[]{mglng, mglat};
        }
    }

2.2 GCJ-02 转 WGS-84

    /**
     * GCJ-02转WGS-84
     *
     * @param lng
     * @param lat
     * @return
     */
    public static double[] gcj02ToWgs84(double lng, double lat) {
        if (outOfChina(lng, lat)) {
            return new double[]{lng, lat};
        } else {
            double dlat = transformLat(lng - 105.0, lat - 35.0);
            double dlng = transformLng(lng - 105.0, lat - 35.0);
            double radlat = lat / 180.0 * PI;
            double magic = Math.sin(radlat);
            magic = 1 - EE * magic * magic;
            double sqrtMagic = Math.sqrt(magic);
            dlat = (dlat * 180.0) / ((AA * (1 - EE)) / (magic * sqrtMagic) * PI);
            dlng = (dlng * 180.0) / (AA / sqrtMagic * Math.cos(radlat) * PI);
            double mglat = lat + dlat;
            double mglng = lng + dlng;
            return new double[]{lng * 2 - mglng, lat * 2 - mglat};
        }
    }

3. GCJ-02 和 BD-09 的相互转换

3.1 GCJ-02 转 BD-09

    /**
     * GCJ-02转BD-09
     *
     * @param lng
     * @param lat
     * @return
     */
    public static double[] gcj02ToBd09(double lng, double lat) {
        double z = Math.sqrt(lng * lng + lat * lat) + 0.00002 * Math.sin(lat * X_PI);
        double theta = Math.atan2(lat, lng) + 0.000003 * Math.cos(lng * X_PI);
        double bdLng = z * Math.cos(theta) + 0.0065;
        double bdLat = z * Math.sin(theta) + 0.006;
        return new double[]{bdLng, bdLat};
    }

3.2 BD-09 转 GCJ-02

    /**
     * BD-09转GCJ-02
     *
     * @param lng
     * @param lat
     * @return
     */
    public static double[] bd09ToGcj02(double lng, double lat) {
        double x = lng - 0.0065;
        double y = lat - 0.006;
        double z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * X_PI);
        double theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * X_PI);
        double gcjLng = z * Math.cos(theta);
        double gcjLat = z * Math.sin(theta);
        return new double[]{gcjLng, gcjLat};
    }

4. WGS-84 和 BD-09 的相互转换

4.1 WGS-84 转 BD-09

    /**
     * WGS-84转BD-09
     *
     * @param lng
     * @param lat
     * @return
     */
    public static double[] wgs84ToBd09(double lng, double lat) {
        double[] gcjPoint = wgs84ToGcj02(lng, lat);
        return gcj02ToBd09(gcjPoint[0], gcjPoint[1]);
    }

4.2 BD-09 转 WGS-84

    /**
     * BD-09转WGS-84
     *
     * @param lng
     * @param lat
     * @return
     */
    public static double[] bd09ToWgs84(double lng, double lat) {
        double[] gcjPoint = bd09ToGcj02(lng, lat);
        return gcj02ToWgs84(gcjPoint[0], gcjPoint[1]);
    }

三、坐标数据转换方法-Python 版

1. 准备工作

1.1 依赖的一些常量和转换方法

import math

# 常量定义
PI = math.pi
x_PI = float(PI * float(3000.0) / float(180.0))
aa = float(6378245.0)
ee = 0.00669342162296594323

# 经纬度偏移转换
def transform_lat(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 transform_lng(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):
    # 纬度3.86~53.55,经度73.66~135.05；非精确判断，可能误判
    if 73.66 < lng < 135.05 and 3.86 < lat < 53.55:
        return False

2. WGS-84 和 GCJ-02 的相互转换

2.1 WGS-84 转 GCJ-02

# WGS-84转GCJ-02
def wgs84_to_gcj02(lng, lat):
    if out_of_china(lng, lat):
        return [lng, lat]
    else:
        dlat = transform_lat(lng - 105.0, lat - 35.0)
        dlng = transform_lng(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) / ((aa * (1 - ee)) / (magic * sqrtmagic) * PI)
        dlng = (dlng * 180.0) / (aa / sqrtmagic * math.cos(radlat) * PI)
        mglat = lat + dlat
        mglng = lng + dlng
        return [mglng, mglat]

2.2 GCJ-02 转 WGS-84

# GCJ-02转WGS-84
def gcj02_to_wgs84(lng, lat):
    if out_of_china(lng, lat):
        return [lng, lat]
    else:
        dlat = transform_lat(lng - 105.0, lat - 35.0)
        dlng = transform_lng(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) / ((aa * (1 - ee)) / (magic * sqrtmagic) * PI)
        dlng = (dlng * 180.0) / (aa / sqrtmagic * math.cos(radlat) * PI)
        mglat = lat + dlat
        mglng = lng + dlng
        return [lng * 2 - mglng, lat * 2 - mglat]

3. GCJ-02 和 BD-09 的相互转换

3.1 GCJ-02 转 BD-09

# GCJ-02转BD-09
def gcj02_to_bd09(lng, lat):
    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 [bd_lng, bd_lat]

3.2 BD-09 转 GCJ-02

# BD-09转GCJ-02
def bd09_to_gcj02(bd_lon, bd_lat):
    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 [gg_lng, gg_lat]

4. WGS-84 和 BD-09 的相互转换

4.1 WGS-84 转 BD-09

# WGS-84转BD-09
def wgs84_to_bd09(lng, lat):
    gcj_point = wgs84_to_gcj02(lng, lat)
    bd_point = gcj02_to_bd09(gcj_point[0], gcj_point[1])
    return [bd_point[0], bd_point[1]]

4.2 BD-09 转 WGS-84

# BD-09转WGS-84
def bd09_to_wgs84(lng, lat):
    gcj_point = bd09_to_gcj02(lng, lat)
    wgs_point = gcj02_to_wgs84(gcj_point[0], gcj_point[1])
    return [wgs_point[0], wgs_point[1]]
-XHP-
关注
  • 3
    点赞
  • 踩
  • 8
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
coordzbzh-GIS坐标转换.rar
04-22
支持单点坐标不同坐标系（空间、平面经纬度、投影坐标系间的转换） 支持中英子午线的设置 可以定义文件格式后批量转换输出 gis数据处理必须工具
坐标系转换专题.ppt
03-28
坐标系转换专题
GIS坐标系
bobozai86的博客
11-27 1509
坐标系简介
GIS中坐标系的基本概念
若海软件科技
12-02 2271
GIS坐标系的分类 1、地理坐标系（Geographic Coordinate System, GCS） 2、投影坐标系（Projected Coordinate System, PCS） 平面极坐标系，这个在初中的时候都学过。GIS中的坐标系可以从这个概念入门学习。可以理解为“为确定两个值之间相对变化而产生的确定值“。这两个值在平面上。 在理解一下球坐标系，“为确定三个值之间相对变化而产生的确定值“。这三个值球面上。 现在再来看下GIS中的地理坐标系，就是在GIS学科中的球坐标系。以地球的质心为圆心，
经纬度及其与坐标系的转换
kissgoodbye2012的博客
07-13 1万+
经纬度和地理坐标系的关系
常见地图坐标系以及转换方法、转换工具
泛、不系之舟
06-08 4760
因为国内对地理坐标有特殊的政策，GPS获取的原始坐标不能在互联网上对外提供，所有电子地图必须对位置做偏移处理，这就导致从底层接口得到的坐标展示在地图上会有偏移。不同的坐标系之间可能有几十到几百米的偏移，所以在地图上绘制矢量图形、做地理数据可视化时，需要修正不同坐标系之间的偏差。 包含常见地图坐标系以及转换方法、转换工具......
GIS基础（2）坐标系的转换
ysqzhatian的博客
07-24 3667
1 引言 上篇记录了GIS坐标系的一些知识，主要是地理坐标系和投影坐标系相关内容。坐标系的转换其实我一直以来都是比较懵的，什么四参数，七参数什么的都不太清楚。正好就此机会全部盘点一下。 2 三个坐标系 前面说过：数据的空间分析必须在统一空间参考下进行，否则结果并不可靠。所以，我们需要转换坐标。因为很多空间数据的操作都与ArcGIS有关，在讨论坐标系转换前我们需要了解三个坐标系： 2.1 真实坐标系 数据的真实坐标系是数据记录本身所对应的坐标系，加载数据后可以从Extent里面看出数据所记录的真实数据是怎么样
GIS：深圳独立坐标系与国家2000坐标系互转教程
热门推荐
Oruizn的博客
08-02 2万+
经过多年的发展，在深圳独立坐标系下的地理数据已经积累很多了，如今国家要求统一转换到国家2000坐标系。对于北京54、西安80这类坐标系下的数据，借助ArcGIS软件能够很轻松地转换为国家2000，但是由于缺少深圳独立坐标系的具体参数，从而存在转换难度。详细的转换参数掌握在深圳市测绘相关部门，对于机关单位能够无偿地进行转换，而对于社会企业则有一定的收费标准。因此产生了一定的需求，如何用较少的成本实现深圳独立坐标系转国家2000坐标系？.........
GIS地图坐标系转换方法，快快收藏
layliangbo的博客
09-08 2451
开发地图的时最烦的就是不同的地图工具以及数据来源坐标系不同造成偏移甚至报错。以下整理了部分的地图坐标转换，一个文件直接搞定所有情况。
meethigher-Gis坐标系4326与3857及高德百度坐标系转换
meethigher的博客
05-22 8931
参考 GIS基础知识 - 坐标系、投影、EPSG:4326、EPSG:3857 GIS基础教程之坐标系 - 知乎 Axis Order — GeoTools 28-SNAPSHOT User Guide EPSG.io: Coordinate Systems Worldwide 区域面积-距离/面积计算-示例中心-JS API 示例 | 高德地图API 一、坐标系 分类 坐标系分为两种 地理坐标系（Geographic Coordinate System, GCS） 投影坐标系（Projected C
坐标系转换工具【支持批量】
12-31
【GIS工具箱】集成了坐标拾取器、坐标系转换工具和随机坐标生成工具；已被300人的公司广泛使用。良心资源，童叟无欺！ 打包出售，更实惠！ https://download.csdn.net/download/jingqia/75970880
geojson文件坐标转换工具
05-24
geojson 格式文件坐标系转换工具，c#编写。 支持wgs、gcj、百度坐标之间互相转换。 内附git地址。
python中用“*”运算符复制列表时子列表实际上引用的是同一个列表，列表和嵌套列表的坑
ThreeS_tones的博客
05-21 116
当你修改了一个子列表的元素时，所有的子列表都会随之改变。因为他们实际上都是同一个子列表的引用。输出如下，可以看到每个子列表中的第5个元素相同，并没有实现从0.8-1.2分布。先看一个例子，我想改变load的子列表中第五个元素，从0.8-1.2分布。的时候，Python实际上。
JVM学习-执行引擎
丢爸
05-21 421
当然是否需要启动JIT编译器将字节码直接编译为对应平台的本地机器指令，则需要根据代码被调用。大部分的程序代码转换为物理机的目标代码或虚拟机能执行的指令集之前，都需要经过下图各个步骤。，将其直接编译为对应平台的本地机器指令，以此提升Java程序的执行性能。”，JIT编译器在运行时会针对那些频繁被调用的热点代码做出。Java代码编译是由Java源码编译器来完成，流程图如下。而定，关于那些需要被编译为本地代码的字节码，称为“Java字节码的执行是由JVM执行引擎来完成。
使用Python操作excel单元格——在单元格中插入公式
xll_007的博客
05-16 358
使用Python操作excel单元格——在单元格中插入公式。 通过使用Python的openpyxl库，来操作excel单元格，在单元格中插入公式的操作。把学习的过程分享给大家。大佬勿喷！
MaxKB，基于 LLM 大语言模型的知识库问答系统，开箱即用，支持快速嵌入到第三方业务系统。
最新发布
AIGC Studio：分享AIGC前沿知识和好玩应用，公众号同名。
05-21 70
MaxKB 是一款基于 LLM 大语言模型的知识库问答系统。MaxKB = Max Knowledge Base，旨在成为企业的最强大脑。开箱即用：支持直接上传文档、自动爬取在线文档，支持文本自动拆分、向量化、RAG（检索增强生成），智能问答交互体验好；无缝嵌入：支持零编码快速嵌入到第三方业务系统；
python—使用排序算法按照字符串长度排序
2401_83954530的博客
05-18 145
del a[0]”，排完序之后，字符串的长度（数字）小的排在前面，列表最前面的元素即为长度最小的字符串，将它打印出来，然后删除这个最短的字符串。【问题描述】首先输入字符串的个数n,再逐行输入每个字符串。c=b ”，逐个计算列表中字符串的长度。然后将字符串的长度代表字符串本身长度，方便排序。a.append(b)”，逐行输入每个字符串，然后将输入的字符串以元素的形式存入列表。“n=eval(input())”，导入需要计算的字符串个数。“a.sort()”，对列表中进行排序。
jQuery知识点（二）
shsjssnn的博客
05-21 466
1.获取节点： 通过各种选择器或选择器的组合来获取想要的网页元素。 2.创建节点： 元素节点通过工厂函数$()直接创建，文本节点和属性节点可以随着元素节点创建时一起创建，也可以通过先创建元素节点，再使用操作内容和属性的方法产生。 3.插入节点： append()、prepend(): 是父元素的方法 appendTo()、prependTo(): 是要插入的节点的方法 before()、after(): 是要插入
gis坐标转换为经纬度在线小程序
07-03
### 回答1： GIS坐标转换为经纬度的在线小程序是一种能够帮助用户实现坐标转换的工具。通过该小程序，用户可以输入不同的GIS坐标（如高斯投影坐标、UTM坐标等），并将其转换为经纬度坐标。具体实现方法如下所述： 首先，用户需要在小程序的输入框中输入待转换的GIS坐标。小程序可以支持多种常用的GIS坐标格式，例如输入格式可以是（X,Y）或者X,Y形式的坐标。 接下来，小程序将根据用户输入的GIS坐标类型，使用相应的算法进行坐标转换。例如，如果用户输入的是高斯投影坐标，小程序将使用高斯投影反算公式进行计算，将其转换为经纬度坐标。 在计算完成后，小程序将会在界面上显示转换后的经纬度坐标。用户可以随时调整输入的GIS坐标，实时查看对应的经纬度坐标。 此外，小程序还可以提供一些额外的功能，例如批量转换、坐标定位等。用户可以通过这些功能进行一次性转换多个坐标，或者在地图上标注坐标点。 综上所述，GIS坐标转换为经纬度的在线小程序可以方便用户在移动设备上进行坐标转换操作。帮助用户快速获取需要的经纬度坐标信息，适用于需要频繁进行GIS坐标转换的场景。同时，小程序的便携性也使得用户可以随时随地进行坐标转换，提高工作效率。 ### 回答2： GIS坐标转换为经纬度在线小程序通常可以通过以下步骤完成： 1. 首先，需要获取用户输入的GIS坐标，包括坐标系类型以及坐标数值。可以通过输入框或地图选择的方式来实现用户输入。 2. 然后，根据用户选择的GIS坐标系类型，进行相应的转换算法。常见的GIS坐标系包括UTM坐标系、高斯-克吕格投影等。针对不同的坐标系，可以编写不同的转换代码。 3. 在转换算法中，需要考虑坐标系的原点、投影参数、椭球体参数等。这些参数通常根据用户选择的坐标系类型进行设置。 4. 转换完成后，将经纬度数值显示给用户。可以以文本形式或地图标注的形式展示。 5. 同时，还可以添加一些附加功能，比如点击地图获取GIS坐标、批量转换多个坐标等。 6. 最后，为了提升用户体验，可以设计一些可选项，如选择常用的GIS坐标系、保存历史转换记录等。 需要注意的是，GIS坐标转换为经纬度涉及到复杂的数学计算和坐标系知识，因此在编写小程序时需要进行充分的测试和验证，确保转换结果准确可靠。同时，还需要确保小程序的界面友好、功能简洁，以便用户能够方便地进行坐标转换操作。 ### 回答3： GIS坐标转换为经纬度的在线小程序是一个方便用户将不同坐标系下的坐标转换为经纬度的工具。用户可以通过输入GIS坐标，选择当前坐标系，然后点击转换按钮进行转换操作。 这个在线小程序的实现需要用到一些基本的数学公式和算法，例如投影算法、反三角函数等。通过这些算法，可以将用户输入的GIS坐标转换为经纬度值。 在实现过程中，先要确定用户当前使用的坐标系，根据不同的坐标系选择不同的转换算法。然后，根据用户输入的GIS坐标，应用相应的算法进行计算，得出经纬度值。 在小程序界面上，可以设置一个输入框用于用户输入GIS坐标，提供一个下拉菜单来选择当前的坐标系。还可以添加一个转换按钮，点击按钮后触发转换操作，并在界面上显示转换后的经纬度结果。 总之，一个GIS坐标转换为经纬度的在线小程序可以帮助用户方便地进行坐标转换操作，无需手动计算和查表。用户只需输入GIS坐标，并选择当前坐标系，点击转换即可得到对应的经纬度值。这个小程序的实现需要用到一些基本的算法和数学知识，但用户可以通过简单的操作就能完成坐标转换。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值