百度墨卡托坐标转化笔记

已于 2023-04-20 17:33:27 修改
阅读量1.1k 收藏 3
点赞数 1
分类专栏： 工作经验 文章标签： java
于 2023-04-10 10:07:18 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ywf008/article/details/130053597
版权
一、墨卡托坐标转化

调研了python和java多种实现方式的转换，发现有的不符合需求，原因还没找到。

我是用百度地图返回的poi边界（返回的是墨卡托坐标）

转换的原理没有深入研究，直接拿来用的，测试可行，非常感谢https://blog.csdn.net/qq_33734315/article/details/121071650

挺好奇MCBAND,MC2LL 这两个变量是从哪里找到的，还是自己研究出来的，太厉害了。
下面两种实现转换后的经纬度均为百度坐标系
调研的java代码实现

package utils;

import java.io.*;
import java.util.ArrayList;
import java.util.HashMap;
import java.util.List;
import java.util.Map;

/**
 * 百度geo返回的结果转换
 * 墨卡托坐标转成百度坐标
 */
public class MocatorConvert {
    private static final Double[] MCBAND = {12890594.86, 8362377.87, 5591021d, 3481989.83, 1678043.12, 0d};
    private static final Double[][] MC2LL = {{1.410526172116255e-8, 0.00000898305509648872, -1.9939833816331, 200.9824383106796, -187.2403703815547, 91.6087516669843, -23.38765649603339, 2.57121317296198, -0.03801003308653, 17337981.2}, {-7.435856389565537e-9, 0.000008983055097726239, -0.78625201886289, 96.32687599759846, -1.85204757529826, -59.36935905485877, 47.40033549296737, -16.50741931063887, 2.28786674699375, 10260144.86}, {-3.030883460898826e-8, 0.00000898305509983578, 0.30071316287616, 59.74293618442277, 7.357984074871, -25.38371002664745, 13.45380521110908, -3.29883767235584, 0.32710905363475, 6856817.37}, {-1.981981304930552e-8, 0.000008983055099779535, 0.03278182852591, 40.31678527705744, 0.65659298677277, -4.44255534477492, 0.85341911805263, 0.12923347998204, -0.04625736007561, 4482777.06}, {3.09191371068437e-9, 0.000008983055096812155, 0.00006995724062, 23.10934304144901, -0.00023663490511, -0.6321817810242, -0.00663494467273, 0.03430082397953, -0.00466043876332, 2555164.4}, {2.890871144776878e-9, 0.000008983055095805407, -3.068298e-8, 7.47137025468032, -0.00000353937994, -0.02145144861037, -0.00001234426596, 0.00010322952773, -0.00000323890364, 826088.5}};

    public static void main(String[] args) {
        String geo = "4|13139481.25,2796369.48;13139944.08,2796651.95|1-13139481.25,2796439.55,13139533.24,2796369.48,13139572.31,2796398.44,13139925.51,2796431.39,13139944.08,2796444.84,13139801.3,2796651.95,13139481.25,2796439.55;";
        String ss = mc2jw(geo);
        System.out.println(ss);
        Double x = 13139533.24;
        Double y = 2796369.48;
        Map<String, Double> rs = convertMC2LL(x, y);
        System.out.println(rs.values());

    }

    /**
     * 墨卡托坐标转成百度坐标 返回的百度坐标
     * @param geo
     * @return
     */
    public static String mc2jw(String geo) {
        //String geo = "4|13517816.2,3597097.8;13527673.6,3607970.0|1-13520242.7,3597248.8,13520242.7,3597399.8,13519636.1,3597399.8,13519636.1,3598607.8,13519484.4,3598607.8,13519484.4,3598758.8,13519029.5,3598758.8,13519029.5,3599060.8,13518574.5,3599060.8,13518574.5,3599211.8,13518271.2,3599211.8,13518271.2,3599815.9,13517816.2,3599815.9,13517816.2,3600268.9,13518119.6,3600268.9,13518119.6,3600570.9,13518271.2,3600570.9,13518271.2,3600872.9,13518574.5,3600872.9,13518574.5,3601023.9,13519181.1,3601023.9,13519181.1,3602533.9,13518877.8,3602533.9,13518877.8,3602986.9,13519332.8,3602986.9,13519332.8,3605100.9,13518574.5,3605100.9,13518574.5,3606611.0,13518422.8,3606611.0,13518422.8,3607064.0,13518574.5,3607064.0,13518574.5,3607366.0,13518877.8,3607366.0,13518877.8,3607819.0,13519332.8,3607819.0,13519332.8,3607970.0,13519787.7,3607970.0,13519787.7,3607064.0,13521152.6,3607064.0,13521152.6,3606460.0,13522062.5,3606460.0,13522062.5,3605856.0,13521910.8,3605856.0,13521910.8,3604345.9,13522669.1,3604345.9,13522669.1,3603288.9,13523730.7,3603288.9,13523730.7,3602533.9,13524185.6,3602533.9,13524185.6,3601023.9,13524337.3,3601023.9,13524337.3,3600721.9,13526157.1,3600721.9,13526157.1,3601174.9,13527218.6,3601174.9,13527218.6,3600419.9,13527370.3,3600419.9,13527370.3,3599513.9,13527673.6,3599513.9,13527673.6,3598758.8,13527218.6,3598758.8,13527218.6,3598305.8,13526612.0,3598305.8,13526612.0,3597097.8,13524943.9,3597097.8,13524943.9,3597550.8,13524792.2,3597550.8,13524792.2,3598154.8,13522517.4,3598154.8,13522517.4,3597701.8,13522214.1,3597701.8,13522214.1,3597550.8,13520697.6,3597550.8,13520697.6,3597097.8,13520242.7,3597097.8;";
        List<String> mocatorList = new ArrayList<String>();
        if (geo.length() <= 71) {
            //没有边界，解析该点坐标
            mocatorList = parseJeoPoint(geo);
        } else {
            //有边界
            mocatorList = parseJeo(geo);
        }
        // 封装板块边界
        StringBuilder sb = new StringBuilder();
        for (int i = 0; i < mocatorList.size(); i++) {
            String[] coordinate = mocatorList.get(i).split("\\#");
            if (coordinate.length<2){
                continue;
            }
            // 墨卡托坐标转换为百度经纬度坐标
            Map<String, Double> location = convertMC2LL(Double.parseDouble(coordinate[0]), Double.parseDouble(coordinate[1]));
            Double lng = location.get("lng");
            Double lat = location.get("lat");
            String coord = lng + "," + lat;
            sb.append(coord);
            if (i < mocatorList.size() - 1) {
                sb.append(";");
            }
        }
        return sb.toString();
    }

    /**
     * 解析Jeo数据
     *
     * @param mocator
     */
    public static List<String> parseJeo(String mocator) {
        List<String> mocatorList = new ArrayList<String>();

        if (null == mocator) return null;
        /* 拆分数据 */
        String[] geos = mocator.split("\\|");
        int n = Integer.parseInt(geos[0]);
        String center = geos[1];
        String polylineMoca = geos[2]; //墨卡托坐标
       String[] plm = null;
        if (n == 2) {
            polylineMoca = "1-" + polylineMoca;
            plm = polylineMoca.split("\\;");
            String plmstr = "";
            for (String str : plm) {
                plmstr += str+",";
            }
            plm = new String[1];
            plm[0] = plmstr;
        } else {
            plm = polylineMoca.split("\\;");
        }
        /* 获取墨卡托边界 */
        String geo = null;
        if (n == 4 || n == 2) {
            for (int i = 0; i < plm.length; i++) {
                String[] geoPaths = plm[i].split("\\-");
                if (geoPaths[0].equals("1")) {
                    geo = geoPaths[1];
                }
            }
        }
        // 墨卡托坐标解析
        String[] geoPolyline = new String[9999];
        try {
            geoPolyline = geo.split("\\,");
        } catch (Exception e) {
            System.out.println(e.getMessage());
        }
        for (int i = 0; i < geoPolyline.length-1; i += 2) {
            mocatorList.add(geoPolyline[i] + "#" + geoPolyline[i + 1]);
        }
        return mocatorList;
    }

    /**
     * 解析Jeo坐标点的数据
     *
     * @param mocatorpoint
     */
    public static List<String> parseJeoPoint(String mocatorpoint) {
        List<String> mocatorList = new ArrayList<String>();

        if (null == mocatorpoint) return null;
        String str = mocatorpoint.substring(0, mocatorpoint.lastIndexOf("|"));
        String str1 = mocatorpoint.substring(str.length() + 1, mocatorpoint.length() - 1);
        /* 拆分数据 */
        // 墨卡托坐标解析
        String[] geoPolyline = str1.split("\\,");
        for (int i = 0; i < geoPolyline.length; i += 2) {
            mocatorList.add(geoPolyline[i] + "#" + geoPolyline[i + 1]);
        }
        return mocatorList;
    }

    /**
     * 墨卡托坐标转经纬度坐标
     *
     * @param x
     * @param y
     * @return
     */
    public static Map<String, Double> convertMC2LL(Double x, Double y) {
        Double[] cF = null;
        x = Math.abs(x);
        y = Math.abs(y);
        for (int cE = 0; cE < MCBAND.length; cE++) {
            if (y >= MCBAND[cE]) {
                cF = MC2LL[cE];
                break;
            }
        }
        Map<String, Double> location = converter(x, y, cF);
        location.put("lng", location.get("x"));
        location.remove("x");
        location.put("lat", location.get("y"));
        location.remove("y");
        return location;
    }

    private static Map<String, Double> converter(Double x, Double y, Double[] cE) {
        Double xTemp = cE[0] + cE[1] * Math.abs(x);
        Double cC = Math.abs(y) / cE[9];
        Double yTemp = cE[2] + cE[3] * cC + cE[4] * cC * cC + cE[5] * cC * cC * cC + cE[6] * cC * cC * cC * cC + cE[7] * cC * cC * cC * cC * cC + cE[8] * cC * cC * cC * cC * cC * cC;
        xTemp *= (x < 0 ? -1 : 1);
        yTemp *= (y < 0 ? -1 : 1);
        Map<String, Double> location = new HashMap<String, Double>();
        location.put("x", xTemp);
        location.put("y", yTemp);
        return location;
    }


}

调研的python实现

MCBAND = [12890594.86, 8362377.87, 5591021, 3481989.83, 1678043.12, 0]

LLBAND = [75, 60, 45, 30, 15, 0]

MC2LL = [[1.410526172116255e-8, 0.00000898305509648872, -1.9939833816331, 200.9824383106796, -187.2403703815547,
          91.6087516669843, -23.38765649603339, 2.57121317296198, -0.03801003308653, 17337981.2],
         [-7.435856389565537e-9, 0.000008983055097726239, -0.78625201886289, 96.32687599759846, -1.85204757529826,
          -59.36935905485877, 47.40033549296737, -16.50741931063887, 2.28786674699375, 10260144.86],
         [-3.030883460898826e-8, 0.00000898305509983578, 0.30071316287616, 59.74293618442277, 7.357984074871,
          -25.38371002664745, 13.45380521110908, -3.29883767235584, 0.32710905363475, 6856817.37],
         [-1.981981304930552e-8, 0.000008983055099779535, 0.03278182852591, 40.31678527705744, 0.65659298677277,
          -4.44255534477492, 0.85341911805263, 0.12923347998204, -0.04625736007561, 4482777.06],
         [3.09191371068437e-9, 0.000008983055096812155, 0.00006995724062, 23.10934304144901, -0.00023663490511,
          -0.6321817810242, -0.00663494467273, 0.03430082397953, -0.00466043876332, 2555164.4],
         [2.890871144776878e-9, 0.000008983055095805407, -3.068298e-8, 7.47137025468032, -0.00000353937994,
          -0.02145144861037, -0.00001234426596, 0.00010322952773, -0.00000323890364, 826088.5]]

LL2MC = [
    [-0.0015702102444, 111320.7020616939, 1704480524535203, -10338987376042340, 26112667856603880, -35149669176653700,
     26595700718403920, -10725012454188240, 1800819912950474, 82.5],
    [0.0008277824516172526, 111320.7020463578, 647795574.6671607, -4082003173.641316, 10774905663.51142,
     -15171875531.51559, 12053065338.62167, -5124939663.577472, 913311935.9512032, 67.5],
    [0.00337398766765, 111320.7020202162, 4481351.045890365, -23393751.19931662, 79682215.47186455, -115964993.2797253,
     97236711.15602145, -43661946.33752821, 8477230.501135234, 52.5],
    [0.00220636496208, 111320.7020209128, 51751.86112841131, 3796837.749470245, 992013.7397791013, -1221952.21711287,
     1340652.697009075, -620943.6990984312, 144416.9293806241, 37.5],
    [-0.0003441963504368392, 111320.7020576856, 278.2353980772752, 2485758.690035394, 6070.750963243378,
     54821.18345352118, 9540.606633304236, -2710.55326746645, 1405.483844121726, 22.5],
    [-0.0003218135878613132, 111320.7020701615, 0.00369383431289, 823725.6402795718, 0.46104986909093,
     2351.343141331292, 1.58060784298199, 8.77738589078284, 0.37238884252424, 7.45]]


def location_mokatuo(lng, lat):  # 百度坐标转百度墨卡托坐标
    cE = None
    lng = get_loop(lng, -180, 180)
    lat = get_range(lat, -74, 74)
    for i in range(len(LLBAND)):
        if lat >= LLBAND[i]:
            cE = LL2MC[i]
            break
    if cE:
        for i in range(len(LLBAND) - 1, 0, -1):
            if lat <= -LLBAND[i]:
                cE = LL2MC[i]
                break
    return converter(lng, lat, cE)


def converter(lng, lat, cE):
    lngTemp = cE[0] + cE[1] * abs(lng)
    cC = abs(lat) / cE[9]
    latTemp = cE[2] + cE[3] * cC + cE[4] * cC * cC + cE[5] * cC * cC * cC + cE[6] * cC * cC * cC * cC + cE[
        7] * cC * cC * cC * cC * cC + cE[8] * cC * cC * cC * cC * cC * cC
    lngTemp *= (-1 if lng < 0 else 1)
    latTemp *= (-1 if lat < 0 else 1)
    return [lngTemp, latTemp]


def get_loop(lng, min1, max1):
    while lng > max1:
        lng -= max1 - min1
    while lng < min1:
        lng += max1 - min1
    return lng


def get_range(lat, min1, max1):
    if min1:
        lat = max(lat, min1)
    if max1:
        lat = min(lat, max1)
    return lat


def mokatuo_location(lng,lat):  # 百度墨卡托转换百度坐标
    cF = None
    lng = abs(lng)
    lat = abs(lat)
    for i in range(len(MCBAND)):
        if lat >= MCBAND[i]:
            cF = MC2LL[i]
            break
    return converter(lng, lat, cF)


if __name__ == '__main__':
    print(location_mokatuo(118.03315104440664, 24.498307986743058))  # 经纬度转墨卡托
    print(mokatuo_location(13139533.24,2796369.48))  # 墨卡托转经纬度
二、坐标测试用例

墨卡托坐标: 13139533.24,2796369.48
百度: 118.03315104440664, 24.498307986743058
GPS: 118.02167877581007, 24.495393563196547
高德: 118.02657021322973, 24.492638370577403

三、调研资料

http://www.site-digger.com/tools/mct2latlng.html

这个工具里面测试了一下返回的经纬度坐标是百度的坐标系

https://www.cnblogs.com/feffery/p/11023673.html

这个实现用上面的测试用例测试没有通过，不知道是不是墨卡托用到的标准不同导致的？

https://blog.csdn.net/qq_33734315/article/details/121071650

主要是用到了这个大佬的思路和实现，非常感谢

默默前行的蜗牛
关注
  • 1
    点赞
  • 踩
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
百度墨卡托坐标转百度经纬度离线实现（源码）
04-19
百度墨卡托坐标转百度经纬度离线实现（源码） https://blog.csdn.net/funkstill/article/details/89400809
百度坐标，百度墨卡托坐标转换.py
09-07
百度经纬度转换百度墨卡托通常需要使用百度api进行转换，比较麻烦。经过一番搜集资料查看js代码终于实现百度经纬度与百度墨卡托离线互转，使用python代码实现
墨卡托坐标以及和经纬度坐标的转换
最新发布
qq_49973861的博客
02-19 709
提示：以下是本篇文章正文内容，下面案例可供参考。
百度摩卡托坐标、墨卡托坐标 转为 百度经纬度坐标 JavaScript
tian6632的博客
11-18 620
摩卡托坐标，墨卡托坐标转百度经纬度坐标
墨卡托坐标与经纬度坐标互换
10-30
在实际使用过程中我们用的最多的还是平面坐标，并且单位最好与常用的长度单位（米）一致。所以就产生了从曲面到平面的转换，这个过程也叫做投影，转换的结果也就是投影平面坐标系。
Mercator(墨卡托)坐标和WGS84(经纬度)互转
04-01
在WebGIS的开发中经常用到的地图投影为Web墨卡托和WGS84。该方法可以实现这两种坐标的相互转换。
convertMC2LL:将百度地图平面投影坐标（墨卡托）转换为经纬度坐标
06-19
convertMC2LL 将百度地图平面投影坐标（墨卡托）转换为经纬度坐标 用在什么场景 使用PHP来调用百度地图API，某些情况下得到的是投影坐标，需要转换为经纬度坐标 为什么要写这个东西 首先，百度地图最好用 第二，百度接口很不准，搜个西城区的地址给你定位到门头沟去 第三，百度有些隐藏接口非常好用 第四，隐藏接口坐标是投影坐标，需要转换成经纬度坐标 第五，只有百度知道怎么转换 第六，只有百度地图JS文件能看到源码 所以就写了这么个程序来将投影坐标转换为经纬度坐标。当然是根据百度地图js原来来改造的。 使用方法 非常简单 include 'src/Baidumap.php'; $baidu = new Baidumap(); $point = new b4(1294830854, 484857493); $ll = $baidu->convertMC2LL($point); var_du
墨卡托投影法将经纬度转换为平面坐标
weixin_42845306的博客
11-29 9485
地球经纬度转换为平面坐标，想必是数学建模里面经常会遇到的问题。 由于地球是一个近似椭圆，因此用经纬度确定的坐标，无法直接使用平面几何的计算公式计算距离等数据。使用墨卡托投影法将经纬度坐标投影为平面坐标。 墨卡托投影 墨卡托投影，又称正轴等角圆柱投影，圆柱投影的一种，由荷兰地图学家墨卡托 (G.Mercator) 于 1569 年提出, 为地图投影方法中影响最大的投影。假设地球被围在一中空的圆柱里，其基准纬线（赤道）与圆柱相切，然后再假想地球中心有一盏灯，把球面上的图形投影到圆柱面上，再把圆柱面展开，这就是一
将墨卡托投影（Transverse Mercator Projection）转换为地理坐标（经纬度）
weixin_44866291的博客
08-21 1092
格式的地图数据【其他格式的数据也可以转换，只要存在地理参考（geoReference）信息即可】，可以通过。
地理知识：墨卡托坐标系
qq_40206371的博客
11-30 1848
假设地球被围在一个中空的圆柱里，其赤道与圆柱相接触，然后再假想地球中心有一盏灯，把球面上的图形投影到圆柱体上，再把圆柱体展开，这就是一幅标准纬线为零度（即赤道）的“墨卡托投影”绘制出的世界地图。经纬线于任何位置皆垂直相交——>
墨卡托坐标以及 墨卡托坐标转经纬度
透明
06-27 3520
墨卡托坐标以及 墨卡托坐标转经纬度
实习一 WGS84经纬度坐标与Web墨卡托投影坐标相互转换1
08-08
运行WGS84经纬度坐标转Web墨卡托投影坐标程序时，显示如下提示：按回车后，把输入的经纬度坐标转换成墨卡托投影坐标并输出，形式如下（保留小数后面2位）：运行W
openlayers学习笔记：
Handsome2013的博客
01-24 2354
使用overlay <!DOCTYPE html> <html> <head> <title>学习overlay</title> <script src="..\bin\ol.js"></script> <link href="..\bin\ol.css"> <script src="https://cdn.staticfile.org/jquery/1.10.2/jq
利用ArcGIS进行坐标转换和纠偏
weixin_30268921的博客
07-08 987
由于测绘局对互联网地图的加密规定以及各个互联网地图处于自身利益的考虑，各个互联网地图公司的坐标系统、地图瓦片组织方式、切片规则等都不各尽相同，因此导致一套坐标数据很难在另一套地图中进行展示，如目前咱们项目组正在实施的“易行杭州”项目需要将杭州公交集团的公交自行车数据在百度地图上进行叠加展示与查询就遇到此类问题。 杭州公交公司的坐标系统未知、但可以从其页面上将所有的自行车空间...
Web地图呈现原理
腾讯位置服务
12-19 825
以下内容转载自木的树的文章《Web地图呈现原理》 作者：木的树 链接：https://www.cnblogs.com/dojo-lzz/p/9250637.html 来源：博客园 著作权归作者所有。商业转载请联系作者获得授权，非商业转载请注明出处。 腾讯位置服务致力于为各行各业提供全方位的位置服务产品。与微信、手机QQ、王者荣耀、 京东、滴滴出行等多个在各自行业具有领先地位的产品开展深度合作。...
3857墨卡托坐标系转换为4326 (WGS84)经纬度坐标
WangYan的博客
07-05 5090
地图经纬度坐标转换
火星坐标、百度坐标、WGS-84坐标相互转换及墨卡托投影坐标转经纬度JavaScript版
lmrylll的博客
07-07 286
【代码】火星坐标、百度坐标、WGS-84坐标相互转换及墨卡托投影坐标转经纬度JavaScript版。
百度地图批量坐标转换（以百度墨卡托坐标转换为百度经纬度坐标为例）
yuishelly的博客
02-21 1373
百度地图批量坐标转换（以百度墨卡托坐标转换为百度经纬度坐标为例）
墨卡托投影转换坐标C++
04-05
墨卡托投影转换坐标C是指将经纬度坐标系下的坐标转换为墨卡托投影下的坐标。具体的转换公式如下： x = R * (lon - lon0) * cos(lat0) y = R * ln(tan(π/4 + lat/2) / tan(π/4 + lat0/2)) 其中，R为地球半径，lon...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值