地图坐标系转换

已于 2024-06-20 16:17:50 修改
阅读量963 收藏 8
点赞数 27
文章标签： java javascript
于 2024-06-19 18:30:59 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_33538365/article/details/139810301
版权

目录

一、坐标系定义

地理坐标系（Geographic Coordinate System）是用于在地球表面定位的系统。它基于一个三维的球体（地球）来定义位置，通常使用经度（longitude）和纬度（latitude）这两个主要参数来描述一个点的位置。

经度（Longitude）

定义：从格林威治子午线（0°经线）起，向东或向西量度的地球表面上的弧长。
范围：从-180°（西经）到+180°（东经）。
表示方法：E（东经）或W（西经）后跟度数。

纬度（Latitude）

定义：从赤道（0°纬线）起，向北或向南量度的地球表面上的弧长。
范围：从-90°（南纬）到+90°（北纬）。
表示方法：N（北纬）或S（南纬）后跟度数。

除了经度和纬度，地理坐标系还可以包括海拔（Elevation）或高程（Altitude），这通常用于描述地面上的点相对于平均海平面的高度。

二、常见坐标系

①WGS-84坐标系：

定义：地心坐标系，GPS原始坐标体系。
特性：在全球定位系统（GPS）中广泛使用，是全球通用的坐标系统。

②GCJ-02坐标系（火星坐标系）：

定义：国测局坐标，一种对经纬度数据的加密算法，即加入随机的偏差。
特性：在中国，互联网地图必须至少使用GCJ-02进行首次加密，不允许直接使用WGS-84坐标下的地理数据。
使用情况：是国内最广泛使用的坐标体系，高德、腾讯、Google中国地图都使用它。

③CGCS2000坐标系（2000国家大地坐标系）：

定义：通过中国GPS连续运行基准站、空间大地控制网以及天文大地网与空间地网联合平差建立的地心大地坐标系统。
特性：属于地心坐标系，是我国当前使用的主要大地坐标系之一。

④BD-09坐标系：

定义：百度中国地图所采用的坐标系，由GCJ-02进行进一步的偏移算法得到。
特性：百度地图特有的坐标系统。

⑤搜狗坐标系：

定义：搜狗地图所采用的坐标系，同样基于GCJ-02进行偏移。
特性：搜狗地图专用的坐标系统。

⑥图吧坐标系：

定义：图吧地图所采用的坐标系，也是基于GCJ-02的偏移算法。
特性：图吧地图特有的坐标系统。

⑦1954北京坐标系：

定义：基于克拉索夫斯基椭球的一种参心大地坐标系。
特性：是我国早期使用的大地坐标系之一。

⑧1980西安坐标系：

定义：基于1975年国际大地测量与地球物理联合会（IUGG）第十六届大会推荐的地球椭球参数的参心坐标系。
特性：在我国地质矿产勘查领域广泛使用。

⑨高斯-克吕格平面直角坐标系：

定义：一种将椭球面上的点投影到平面上，并保持投影后角度不变的投影坐标系。
特性：常用于地形图制作和测量工程。

⑩其他特定行业坐标系：

在某些特定行业或应用中，可能还存在其他特定的地理坐标系，如UTM坐标系（通用横轴墨卡托投影坐标系）等，这些坐标系通常根据特定需求和应用场景进行定义和使用。

三、坐标系转换代码

①java

/**
 * @author LiChenHao
 * @date 2023/8/25 13:20
 * 坐标转换工具类
 */
public class GPSUtil {

    public static double pi = 3.1415926535897932384626;
    public static double x_pi = 3.14159265358979324 * 3000.0 / 180.0;

    public static double a = 6378245.0;
    public static double ee = 0.00669342162296594323;

    public static double transformLat(double x, double y) {
        double ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y
                + 0.2 * Math.sqrt(Math.abs(x));
        ret += (20.0 * Math.sin(6.0 * x * pi) + 20.0 * Math.sin(2.0 * x * pi)) * 2.0 / 3.0;
        ret += (20.0 * Math.sin(y * pi) + 40.0 * Math.sin(y / 3.0 * pi)) * 2.0 / 3.0;
        ret += (160.0 * Math.sin(y / 12.0 * pi) + 320 * Math.sin(y * pi / 30.0)) * 2.0 / 3.0;
        return ret;
    }

    public static double transformLon(double x, double y) {
        double ret = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1
                * Math.sqrt(Math.abs(x));
        ret += (20.0 * Math.sin(6.0 * x * pi) + 20.0 * Math.sin(2.0 * x * pi)) * 2.0 / 3.0;
        ret += (20.0 * Math.sin(x * pi) + 40.0 * Math.sin(x / 3.0 * pi)) * 2.0 / 3.0;
        ret += (150.0 * Math.sin(x / 12.0 * pi) + 300.0 * Math.sin(x / 30.0
                * pi)) * 2.0 / 3.0;
        return ret;
    }

    public static double[] transform(double lat, double lon) {
        if (outOfChina(lat, lon)) {
            return new double[]{lat, lon};
        }
        double dLat = transformLat(lon - 105.0, lat - 35.0);
        double dLon = transformLon(lon - 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) / ((a * (1 - ee)) / (magic * sqrtMagic) * pi);
        dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * pi);
        double mgLat = lat + dLat;
        double mgLon = lon + dLon;
        return new double[]{mgLat, mgLon};
    }

    public static boolean outOfChina(double lat, double lon) {
        if (lon < 72.004 || lon > 137.8347)
            return true;
        if (lat < 0.8293 || lat > 55.8271)
            return true;
        return false;
    }

    /**
     * 84 to 火星坐标系 (GCJ-02) World Geodetic System ==> Mars Geodetic System
     *
     * @param lat
     * @param lon
     * @return
     */
    public static double[] gps84_To_Gcj02(double lat, double lon) {
        if (outOfChina(lat, lon)) {
            return new double[]{lat, lon};
        }
        double dLat = transformLat(lon - 105.0, lat - 35.0);
        double dLon = transformLon(lon - 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) / ((a * (1 - ee)) / (magic * sqrtMagic) * pi);
        dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * pi);
        double mgLat = lat + dLat;
        double mgLon = lon + dLon;
        return new double[]{mgLat, mgLon};
    }

    /**
     * * 火星坐标系 (GCJ-02) to 84 * * @param lon * @param lat * @return
     */
    public static double[] gcj02_To_Wgs84(double lat, double lon) {
        double[] gps = transform(lat, lon);
        double lontitude = lon * 2 - gps[1];
        double latitude = lat * 2 - gps[0];
        return new double[]{latitude, lontitude};
    }

    /**
     * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换算法 将 GCJ-02 坐标转换成 BD-09 坐标
     *
     * @param lat
     * @param lon
     */
    public static double[] gcj02_To_Bd09(double lat, double lon) {
        double x = lon, y = lat;
        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 tempLon = z * Math.cos(theta) + 0.0065;
        double tempLat = z * Math.sin(theta) + 0.006;
        double[] gps = {tempLat, tempLon};
        return gps;
    }

    /**
     * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换算法 * * 将 BD-09 坐标转换成GCJ-02 坐标
     *  @param lat
     *  @param lon
     *  @return
     */
    public static double[] bd09_To_Gcj02(double lat, double lon) {
        double x = lon - 0.0065, 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 tempLon = z * Math.cos(theta);
        double tempLat = z * Math.sin(theta);
        double[] gps = {tempLat, tempLon};
        return gps;
    }

    /**
     * 将WGS84(GPS坐标系)转为bd09
     *
     * @param lat
     * @param lon
     * @return
     */
    public static double[] wgs84_To_Bd09(double lat, double lon) {
        double[] gcj02 = gps84_To_Gcj02(lat, lon);
        double[] bd09 = gcj02_To_Bd09(gcj02[0], gcj02[1]);
        return bd09;
    }

    public static double[] bd09_To_Wgs84(double lat, double lon) {
        double[] gcj02 = bd09_To_Gcj02(lat, lon);
        double[] gps84 = gcj02_To_Wgs84(gcj02[0], gcj02[1]);
        //保留小数点后六位
        gps84[0] = retain7(gps84[0]);
        gps84[1] = retain7(gps84[1]);
        return gps84;
    }

    /**
     *  两点间距离  1° = 111 km (or 60 nautical miles) 0.1° = 11.1 km.
     *  单位：km
     */
    public static double getDistance(double lat1, double lon1, double lat2, double lon2){
        double latitudeDifference = lat1 - lat2;
        double longitudeDifference = lon1 - lon2;
        double sqrt = Math.sqrt((latitudeDifference * latitudeDifference) + (longitudeDifference * longitudeDifference));
        return Math.ceil(sqrt * 111);
    }

    /**
     * 保留小数点后六位
     *
     * @param num
     * @return
     */
    private static double retain7(double num) {
        String result = String.format("%.7f", num);
        return Double.valueOf(result);
    }

	//test
    public static void main(String[] args) {
        double lat1 = 31.026679;
        double lon1 = 121.611927;

        double lat2 = 31.026679;
        double lon2 = 121.911927;
        System.out.println(getDistance(lat1, lon1, lat2, lon2));
    }
}

②js

//定义一些常量
var x_PI = 3.14159265358979324 * 3000.0 / 180.0;
var PI = 3.1415926535897932384626;
var a = 6378245.0;
var ee = 0.00669342162296594323;


/**
 * WGS84 与百度坐标系 (BD-09) 的转换
 * @param lng
 * @param lat
 * @returns {*[]}
 */
export function wgs84tobd09(lng, lat) {
  let pointArr = wgs84togcj02(lng,lat);
  return gcj02tobd09(pointArr[0],pointArr[1]);
}

/**
 * 百度坐标系 (BD-09) 与 WGS84的转换
 * @param lng
 * @param lat
 * @returns {*[]}
 */
export function bd09towgs84(lng, lat) {
  let pointArr = bd09togcj02(lng,lat);
  return gcj02towgs84(pointArr[0],pointArr[1]);
}

/**
 * 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换
 * 百度 转 谷歌、高德
 * @param lng
 * @param lat
 * @returns {*[]}
 */
export function bd09togcj02(bd_lon, bd_lat) {
  var x_pi = 3.14159265358979324 * 3000.0 / 180.0;
  var x = bd_lon - 0.0065;
  var y = bd_lat - 0.006;
  var z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * x_pi);
  var theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * x_pi);
  var gg_lng = z * Math.cos(theta);
  var gg_lat = z * Math.sin(theta);
  return [gg_lng, gg_lat]
}

/**
 * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换
 * 即谷歌、高德 转 百度
 * @param lng
 * @param lat
 * @returns {*[]}
 */
export function gcj02tobd09(lng, lat) {
  var z = Math.sqrt(lng * lng + lat * lat) + 0.00002 * Math.sin(lat * x_PI);
  var theta = Math.atan2(lat, lng) + 0.000003 * Math.cos(lng * x_PI);
  var bd_lng = z * Math.cos(theta) + 0.0065;
  var bd_lat = z * Math.sin(theta) + 0.006;
  return [bd_lng, bd_lat]
}

/**
 * WGS84转GCj02
 * @param lng
 * @param lat
 * @returns {*[]}
 */
export function wgs84togcj02(lng, lat) {
  if (out_of_china(lng, lat)) {
    return [lng, lat]
  }
  else {
    var dlat = transformlat(lng - 105.0, lat - 35.0);
    var dlng = transformlng(lng - 105.0, lat - 35.0);
    var radlat = lat / 180.0 * PI;
    var magic = Math.sin(radlat);
    magic = 1 - ee * magic * magic;
    var sqrtmagic = Math.sqrt(magic);
    dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * PI);
    dlng = (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * PI);
    var mglat = lat + dlat;
    var mglng = lng + dlng;
    return [mglng, mglat]
  }
}

/**
 * GCJ02 转换为 WGS84
 * @param lng
 * @param lat
 * @returns {*[]}
 */
export function gcj02towgs84(lng, lat) {
  if (out_of_china(lng, lat)) {
    return [lng, lat]
  }
  else {
    var dlat = transformlat(lng - 105.0, lat - 35.0);
    var dlng = transformlng(lng - 105.0, lat - 35.0);
    var radlat = lat / 180.0 * PI;
    var magic = Math.sin(radlat);
    magic = 1 - ee * magic * magic;
    var sqrtmagic = Math.sqrt(magic);
    dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * PI);
    dlng = (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * PI);
    var mglat = lat + dlat;
    var mglng = lng + dlng;
    return [lng * 2 - mglng, lat * 2 - mglat]
  }
}

export function transformlat(lng, lat) {
  var 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
}

export function transformlng(lng, lat){
  var 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
}
export function out_of_china(lng, lat) {
  return (lng < 72.004 || lng > 137.8347) || ((lat < 0.8293 || lat > 55.8271) || false)
}

/**
 *
 * @param area 后台查询到的区域树dtArea
 * @param names 需要匹配的名称集合【上海市，上海市，浦东新区】
 * @returns codes
 */
export function provinceCityCountyNameToCode(area, names) {
  if(!area || !names){
    return null;
  }
  var codes=[];
  area.forEach(province=>{
    if(province.areaName ===names[0] || province.areaName+'市' === names[0]){
      codes.push(province.id);
      province.children.forEach(city=>{
        if(city.areaName === names[1]){
          codes.push(city.id);
          city.children.forEach(county=>{
            if(county.areaName === names[2]){
              codes.push(county.id);
            }
          })
        }
      })
    }
  })
  return codes;
}
李晨豪 Lch
关注
  • 27
    点赞
  • 踩
  • 8
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
如何批量轻松实现地图坐标系转换
tool_lu_的博客
01-04 1488
实现WGS84（国际通用）、GCJ02（高德\QQ）、BD09（百度地图）互相转换，支持Excel批量导入、批量转换、Excel批量导出。
不同地图坐标系转换工具.zip
07-17
本压缩包“不同地图坐标系转换工具.zip”可能包含一系列软件或脚本，用于在不同的坐标系统间进行转换。下面将详细阐述地图坐标系的基本知识以及转换工具的重要性和使用方法。 地图坐标系是为了精确地定位地球上任何...
地图坐标系转换简介
leokingszx的博客
08-28 1493
地图坐标系转换简介及问题分享 背景介绍 从android系统定位获取到的经纬度，到谷歌地图查询，发现定位偏移很大，以为是定位误差或者定位策略问题导致，实际上并非如此，详情以下分解。 各坐标系简介 WGS84坐标系 即地球坐标系，国际上通用的坐标系。 设备一般包含GPS芯片或者北斗芯片获取的经纬度为WGS84地理坐标系。谷歌地图采用的是WGS84地理坐标系（中国范围内谷歌中国地图采用的是GCJ02地理坐标系)，android手机定位使用 GCJ02坐标系 即火星坐标系，WGS84坐标系经加
百度地图坐标系转换
fantudan00的博客
02-02 2532
做一个包含地图的项目，地图切片用的是百度的。实际使用时发现，地图上数据的点位有误差。这是因为坐标系不一致导致。正常手机定位采用的是GPS，上报的点位都是WGS84坐标系，而百度地图采用的是BD09坐标系。
常见地图坐标系以及转换方法、转换工具
泛、不系之舟
06-08 4923
因为国内对地理坐标有特殊的政策，GPS获取的原始坐标不能在互联网上对外提供，所有电子地图必须对位置做偏移处理，这就导致从底层接口得到的坐标展示在地图上会有偏移。不同的坐标系之间可能有几十到几百米的偏移，所以在地图上绘制矢量图形、做地理数据可视化时，需要修正不同坐标系之间的偏差。 包含常见地图坐标系以及转换方法、转换工具......
易语言-地图坐标系转换模块
06-25
网上翻了好久，没找到易语言坐标系转换模块，正好最近需要，就参考网上java源码写了个易语言的模块，分享给大家，用户Google国内地图、高德、QQ地图转地球坐标系，国际通用坐标系，。有不对的地方也欢迎指正。 注：...
地图坐标系转换模块新增转换规则-易语言
06-11
在IT行业中，地图坐标系转换是一项重要的技术，尤其在地理信息系统（GIS）和导航软件开发中扮演着核心角色。易语言是一种中国本土的编程语言，它以简洁的语句和直观的语法为特点，便于初学者上手。在这个场景中，...
《JavaSE》---6.＜基础语法（Java三大程序控制结构）＞
m0_73456341的博客
07-14 823
本篇博客主要讲解Java基础语法中的三大结构，一种顺序结构、两大分支结构i（if-else、swich-case）、四大循环结构（while、do while、fot、foreach）
Spring Security 授权
persistence_PSH的博客
07-14 895
在 loadUserByUsername 方法，在查询数据库用户信息的时候，同时查询出用户的权限，这里以角色名代指权限。在 loadUserByUsername 方法，在查询数据库用户信息的时候，同时查询出用户的权限，这里以角色名代指权限。在 loadUserByUsername 方法，在查询数据库用户信息的时候，同时查询出用户的权限，这里以角色名代指权限。通过 RBAC 获取到用户的具体权限后，再通过 Security 的 用户-权限-资源 来进行权限控制。HttpSecurity 权限配置。
【Go系列】 Sync并发控制
u013379032的博客
07-14 1021
在上一篇文章中，我们介绍了goroutine和channel，理论上，通过channel可以实现并发控制，但是其他语言的开发者可能更习惯一些原子操作的库。当然Go语言也会提供这样的库，所以我们今天了解一下sync库。
day04：Redis和店铺营业状态设置
m0_69266818的博客
07-15 825
介绍了Redis命令和用java操作redis还有营业额状态接口的书写
【JAVA】相关程序组件
qq_44077764的博客
07-15 156
【JAVA】相关程序组件
【Java算法】前缀和 下
最新发布
2302_79806056的博客
07-18 686
这段代码实现了一个寻找数组中具有相等数量的0和1的最长子数组长度的算法。具体来说，它使用了前缀和（prefix sum）的概念以及哈希映射（HashMap）来优化查找过程。初始化哈希映射：计算前缀和：检查相等数量的0和1：更新哈希映射：返回结果：这种方法的时间复杂度为O(n)，因为我们只遍历数组一次，并且对每个元素执行常数时间的操作。空间复杂度也是O(n)，最坏情况下需要存储n个前缀和值。 初始化阶段:开始遍历数组:继续
Java 中有哪几种基本数据类型？请分别列出它们并简述每种数据类型的特点及其在内存中的占用空间？
liangzai215的专栏
07-15 497
在Java的世界里，数据是构建应用程序的基石。为了高效地处理这些数据，Java设计了一系列基础数据类型，它们直接映射到计算机硬件上，因此在性能和内存使用上更为高效。我们常说的Java八大基本数据类型，涵盖了整数、浮点数、字符和布尔值，下面我将一一介绍它们的特点以及在内存中的占用空间，并通过简单的代码实例来加深理解。
maven项目容器化运行之1-基于1Panel软件将docker镜像构建能力分享给局域网
qq_37372909的博客
07-15 512
docker核心功能包括镜像的构建和容器的运行，它可以开放端口将镜像构建的能力暴露。docker daeson就是docker的守护进程，开放能力是通过指定docker守护进程监听端口来实现的。后面开发人员就可以通过maven中docker插件来远程调用docker镜像构建能力了。
深入剖析 Java 中 Spring Bean 的生命周期
weixin_42545951的博客
07-12 704
在 Java 开发中，Spring 框架是广泛使用的企业级应用开发框架。理解 Spring Bean 的生命周期对于开发者来说至关重要，这也是面试中经常会被问到的重要知识点。
从零手写实现 nginx-27-return 指令
07-14 668
大家好，我是老马。很高兴遇到你。我们为 java 开发者实现了 java 版本的 nginx如果你想知道 servlet 如何处理的，可以参考我的另一个项目：手写从零实现简易版 tomcatminicat。
地图投影与地理坐标系
地图投影是将三维的地球表面转换到二维平面上的过程，而坐标系则是定位地球上任何位置的基础。 1、地理空间 地理空间是指地球表层，包括陆地和海洋，是人类活动的主要区域。在GIS中，地理空间分为绝对空间和相对...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

打赏作者

李晨豪 Lch

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

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

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

打赏作者

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

抵扣说明：

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

余额充值