WGS84坐标互转百度坐标优化

于 2025-03-17 10:01:12 发布
阅读量508 收藏
点赞数 8
文章标签： java
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_56144522/article/details/146307990
版权

  1. 坐标系的加密特性
    WGS84到BD-09需经过两次加密（WGS84 → GCJ-02 → BD-09），而逆向转换需解密两次（BD-09 → GCJ-02 → WGS84）。GCJ-02的加密算法包含非线性偏移，直接逆向公式可能导致误差。使用迭代法逆向计算WGS84坐标。

  2. 高精度计算
    使用BigDecimal类型存储坐标值，避免中间过程的精度截断。确保所有数学运算（如三角函数）使用高精度库（采用maven中的Apfloat库）。

  3. 可靠的转换算法
    使用经过广泛验证的开源实现，而非自行推导公式，以降低算法误差。

  4. 测试与验证
    通过大规模坐标回环测试（WGS84 → BD-09 → WGS84）验证误差范围，确保精度可接受。

package org.example;

import java.math.BigDecimal;
import java.math.MathContext;
import java.math.RoundingMode;
import org.apfloat.Apfloat;
import org.apfloat.ApfloatMath;
import org.apfloat.ApfloatRuntimeException;

public class HighPrecisionCoordinateConverter {
    // 精度配置
    private static final MathContext MC = MathContext.DECIMAL128;
    private static final  int APFLOAT_PRECISION = 50;

    // 常量定义（使用字符串初始化）
    private static final BigDecimal PI = new BigDecimal("3.14159265358979323846264338327950288419716939937510");
    private static final BigDecimal X_PI = PI.multiply(new BigDecimal("3000")).divide(new BigDecimal("180"), MC);
    private static final BigDecimal A = new BigDecimal("6378245.0");
    private static final BigDecimal EE = new BigDecimal("0.00669342162296594323");
    private static final BigDecimal MAGIC_NUMBER = new BigDecimal("0.00002");
    private static final BigDecimal BD_OFFSET = new BigDecimal("0.0065");

    // 国内坐标范围
    private static final BigDecimal CHINA_LNG_MIN = new BigDecimal("72.004");
    private static final BigDecimal CHINA_LNG_MAX = new BigDecimal("137.8347");
    private static final BigDecimal CHINA_LAT_MIN = new BigDecimal("0.8293");
    private static final BigDecimal CHINA_LAT_MAX = new BigDecimal("55.8271");

    //======================== 公有接口 ========================
    public static double[] bd09ToWgs84(double lng, double lat) {
        BigDecimal[] gcj = bd09ToGcj02(new BigDecimal(lng), new BigDecimal(lat));
        BigDecimal[] wgs = gcj02ToWgs84(gcj[0], gcj[1]);
        return convertToDoubleArray(wgs);
    }

    public static double[] wgs84ToBd09(double lng, double lat) {
        BigDecimal[] gcj = wgs84ToGcj02(new BigDecimal(lng), new BigDecimal(lat));
        BigDecimal[] bd = gcj02ToBd09(gcj[0], gcj[1]);
        return convertToDoubleArray(bd);
    }

    //======================== 核心算法 ========================
    private static BigDecimal[] bd09ToGcj02(BigDecimal bdLng, BigDecimal bdLat) {
        BigDecimal x = bdLng.subtract(BD_OFFSET, MC);
        BigDecimal y = bdLat.subtract(new BigDecimal("0.006"), MC);

        BigDecimal z = sqrt(x.pow(2).add(y.pow(2)))
                .subtract(sin(y.multiply(X_PI)).multiply(MAGIC_NUMBER), MC);

        BigDecimal theta = atan2(y, x).subtract(
                cos(x.multiply(X_PI)).multiply(new BigDecimal("0.000003")), MC);

        return new BigDecimal[]{
                z.multiply(cos(theta), MC),
                z.multiply(sin(theta), MC)
        };
    }

    private static BigDecimal[] gcj02ToBd09(BigDecimal gcjLng, BigDecimal gcjLat) {
        BigDecimal z = sqrt(gcjLng.pow(2).add(gcjLat.pow(2)))
                .add(sin(gcjLat.multiply(X_PI)).multiply(MAGIC_NUMBER), MC);

        BigDecimal theta = atan2(gcjLat, gcjLng)
                .add(cos(gcjLng.multiply(X_PI)).multiply(new BigDecimal("0.000003")), MC);

        return new BigDecimal[]{
                z.multiply(cos(theta), MC).add(BD_OFFSET, MC),
                z.multiply(sin(theta), MC).add(new BigDecimal("0.006"), MC)
        };
    }

    private static BigDecimal[] wgs84ToGcj02(BigDecimal lng, BigDecimal lat) {
        if (isOutsideChina(lng, lat)) {
            return new BigDecimal[]{lng, lat};
        }
        BigDecimal[] delta = calculateDelta(lng, lat);
        return new BigDecimal[]{
                lng.add(delta[0], MC),
                lat.add(delta[1], MC)
        };
    }

    private static BigDecimal[] gcj02ToWgs84(BigDecimal lng, BigDecimal lat) {
        if (isOutsideChina(lng, lat)) {
            return new BigDecimal[]{lng, lat};
        }
        BigDecimal[] delta = calculateDelta(lng, lat);
        return new BigDecimal[]{
                lng.subtract(delta[0], MC),
                lat.subtract(delta[1], MC)
        };
    }

    //======================== 数学工具方法 ========================
    private static BigDecimal[] calculateDelta(BigDecimal lng, BigDecimal lat) {
        BigDecimal baseLng = lng.subtract(new BigDecimal("105"), MC);
        BigDecimal baseLat = lat.subtract(new BigDecimal("35"), MC);

        BigDecimal radLat = radians(lat);
        BigDecimal magic = sin(radLat).pow(2).multiply(EE, MC);
        magic = BigDecimal.ONE.subtract(magic, MC);
        BigDecimal sqrtMagic = sqrt(magic);

        BigDecimal dLat = transformLat(baseLng, baseLat)
                .multiply(new BigDecimal("180"), MC)
                .divide(A.multiply(BigDecimal.ONE.subtract(EE))
                        .divide(magic.multiply(sqrtMagic), MC)
                        .divide(PI, MC),MC);

        BigDecimal dLng = transformLng(baseLng, baseLat)
                .multiply(new BigDecimal("180"), MC)
                .divide(A.multiply(cos(radLat)), MC)
                .divide(sqrtMagic.multiply(PI), MC);

        return new BigDecimal[]{dLng, dLat};
    }

    private static BigDecimal transformLat(BigDecimal x, BigDecimal y) {
        MathContext mc = MathContext.DECIMAL128;

        // 第一项：-100 + 2x + 3y
        BigDecimal term1 = new BigDecimal("-100")
                .add(x.multiply(TWO, mc))
                .add(y.multiply(THREE, mc), mc);

        // 第二项：0.2y²
        BigDecimal term2 = y.pow(2, mc)
                .multiply(new BigDecimal("0.2"), mc);

        // 第三项：0.1xy
        BigDecimal term3 = x.multiply(y, mc)
                .multiply(ONE_TENTH, mc);

        // 第四项：0.2√|x|
        BigDecimal term4 = sqrt(x.abs(), APFLOAT_PRECISION)
                .multiply(new BigDecimal("0.2"), mc);

        // 合并基础项
        BigDecimal result = term1.add(term2, mc)
                .add(term3, mc)
                .add(term4, mc);

        // 第五项：正弦波动项1 (20sin6πx + 20sin2πx) * 2/3
        BigDecimal sin6x = sin(x.multiply(SIX).multiply(PI), APFLOAT_PRECISION);
        BigDecimal sin2x = sin(x.multiply(TWO).multiply(PI), APFLOAT_PRECISION);
        BigDecimal wave1 = TWENTY.multiply(sin6x.add(sin2x), mc)
                .multiply(TWO_THIRDS, mc);

        // 第六项：正弦波动项2 (20sinπy + 40sin(πy/3)) * 2/3
        BigDecimal sinY = sin(y.multiply(PI), APFLOAT_PRECISION);
        BigDecimal sinY3 = sin(y.divide(THREE, mc).multiply(PI), APFLOAT_PRECISION);
        BigDecimal wave2 = TWENTY.multiply(sinY, mc)
                .add(FORTY.multiply(sinY3, mc), mc)
                .multiply(TWO_THIRDS, mc);

        // 第七项：长周期波动项 (160sin(πy/12) + 320sin(πy/30)) * 2/3
        BigDecimal sinY12 = sin(y.divide(new BigDecimal("12"), mc).multiply(PI), APFLOAT_PRECISION);
        BigDecimal sinY30 = sin(y.divide(new BigDecimal("30"), mc).multiply(PI), APFLOAT_PRECISION);
        BigDecimal wave3 = new BigDecimal("160").multiply(sinY12, mc)
                .add(new BigDecimal("320").multiply(sinY30, mc), mc)
                .multiply(TWO_THIRDS, mc);

        // 累计所有分量
        return result.add(wave1, mc)
                .add(wave2, mc)
                .add(wave3, mc);
    }

    /**
     * 高精度经度转换
     */
    private static BigDecimal transformLng(BigDecimal x, BigDecimal y) {
        MathContext mc = MathContext.DECIMAL128;

        // 基础项：300 + x + 2y
        BigDecimal base = new BigDecimal("300")
                .add(x, mc)
                .add(y.multiply(TWO, mc), mc);

        // 二次项：0.1x² + 0.1xy
        BigDecimal quadTerm = x.pow(2, mc).multiply(ONE_TENTH, mc)
                .add(x.multiply(y, mc).multiply(ONE_TENTH, mc), mc);

        // 根号项：0.1√|x|
        BigDecimal sqrtTerm = sqrt(x.abs(), APFLOAT_PRECISION)
                .multiply(ONE_TENTH, mc);

        BigDecimal result = base.add(quadTerm, mc)
                .add(sqrtTerm, mc);

        // 高频波动项1 (20sin6πx + 20sin2πx) * 2/3
        BigDecimal sin6x = sin(x.multiply(SIX).multiply(PI), APFLOAT_PRECISION);
        BigDecimal sin2x = sin(x.multiply(TWO).multiply(PI), APFLOAT_PRECISION);
        BigDecimal wave1 = TWENTY.multiply(sin6x.add(sin2x), mc)
                .multiply(TWO_THIRDS, mc);

        // 中频波动项2 (20sinπx + 40sin(πx/3)) * 2/3
        BigDecimal sinX = sin(x.multiply(PI), APFLOAT_PRECISION);
        BigDecimal sinX3 = sin(x.divide(THREE, mc).multiply(PI), APFLOAT_PRECISION);
        BigDecimal wave2 = TWENTY.multiply(sinX, mc)
                .add(FORTY.multiply(sinX3, mc), mc)
                .multiply(TWO_THIRDS, mc);

        // 低频波动项3 (150sin(πx/12) + 300sin(πx/30)) * 2/3
        BigDecimal sinX12 = sin(x.divide(new BigDecimal("12"), mc).multiply(PI), APFLOAT_PRECISION);
        BigDecimal sinX30 = sin(x.divide(new BigDecimal("30"), mc).multiply(PI), APFLOAT_PRECISION);
        BigDecimal wave3 = new BigDecimal("150").multiply(sinX12, mc)
                .add(new BigDecimal("300").multiply(sinX30, mc), mc)
                .multiply(TWO_THIRDS, mc);

        return result.add(wave1, mc)
                .add(wave2, mc)
                .add(wave3, mc);
    }

    //======================== 高精度数学函数 ========================
    private static BigDecimal sin(BigDecimal x) {
        try {
            return new BigDecimal(String.valueOf(ApfloatMath.sin(new Apfloat(x))));
        } catch (ApfloatRuntimeException e) {
            return new BigDecimal(Math.sin(x.doubleValue()));
        }
    }

    private static BigDecimal cos(BigDecimal x) {
        try {
            return new BigDecimal(String.valueOf(ApfloatMath.cos(new Apfloat(x))));
        } catch (ApfloatRuntimeException e) {
            return new BigDecimal(Math.cos(x.doubleValue()));
        }
    }
//
    private static BigDecimal atan2(BigDecimal y, BigDecimal x) {
        try {
            Apfloat apfY = new Apfloat(y.toString(), APFLOAT_PRECISION);
            Apfloat apfX = new Apfloat(x.toString(), APFLOAT_PRECISION);
            return new BigDecimal(ApfloatMath.atan2(apfY, apfX).toString());
        } catch (ApfloatRuntimeException e) {
            return new BigDecimal(Math.atan2(y.doubleValue(), x.doubleValue()));
        }
    }

    private static BigDecimal sqrt(BigDecimal x) {
        try {
            return new BigDecimal(
                    ApfloatMath.sqrt(new Apfloat(x.toString(), APFLOAT_PRECISION))
                            .toString());
        } catch (ApfloatRuntimeException e) {
            return BigDecimal.valueOf(Math.sqrt(x.doubleValue()));
        }
    }

    private static BigDecimal sin(BigDecimal x, int precision) {
        try {
            Apfloat radians = new Apfloat(x.toString(), precision);
            return new BigDecimal(ApfloatMath.sin(radians).toString());
        } catch (ApfloatRuntimeException e) {
            // 异常时降级到JDK的双精度计算
            return BigDecimal.valueOf(Math.sin(x.doubleValue()));
        }
    }

    /**
     * 高精度平方根 (使用Apfloat实现)
     */
    private static BigDecimal sqrt(BigDecimal x, int precision) {
        try {
            Apfloat num = new Apfloat(x.toString(), precision);
            return new BigDecimal(ApfloatMath.sqrt(num).toString());
        } catch (ApfloatRuntimeException e) {
            return BigDecimal.valueOf(Math.sqrt(x.doubleValue()));
        }
    }

    //======================== 工具方法 ========================
    private static boolean isOutsideChina(BigDecimal lng, BigDecimal lat) {
        return lng.compareTo(CHINA_LNG_MIN) < 0 || lng.compareTo(CHINA_LNG_MAX) > 0 ||
                lat.compareTo(CHINA_LAT_MIN) < 0 || lat.compareTo(CHINA_LAT_MAX) > 0;
    }

    private static BigDecimal radians(BigDecimal degrees) {
        return degrees.multiply(PI).divide(new BigDecimal("180"), MC);
    }



    private static double[] convertToDoubleArray(BigDecimal[] arr) {
        return new double[]{
                arr[0].setScale(12, RoundingMode.HALF_UP).doubleValue(),
                arr[1].setScale(12, RoundingMode.HALF_UP).doubleValue()
        };
    }

    //======================== 常量缓存 ========================
    private static final BigDecimal TWO = new BigDecimal("2");
    private static final BigDecimal THREE = new BigDecimal("3");
    private static final BigDecimal SIX = new BigDecimal("6");
    private static final BigDecimal TWENTY = new BigDecimal("20");
    private static final BigDecimal FORTY = new BigDecimal("40");
    private static final BigDecimal ONE_TENTH = new BigDecimal("0.1");
    private static final BigDecimal TWO_THIRDS = new BigDecimal("2").divide(new BigDecimal("3"), 34, RoundingMode.HALF_UP);

}

wgs84坐标转百度坐标
04-04
wgs84坐标转百度坐标
火星百度坐标转WGS84坐标小工具
10-10
"火星百度坐标转WGS84坐标小工具"是一个专为解决中国地区坐标转换问题而设计的实用程序。在本篇文章中，我们将深入探讨这个工具的工作原理、涉及的技术以及其在实际应用中的重要性。 首先，我们来了解一下标题中的...
84坐标系、02坐标系、百度坐标之间相互转换算法
热门推荐
日拱一卒无有尽，功不唐捐终入海
09-11 1万+
提供了百度坐标（BD09）、国测局坐标（火星坐标，GCJ02）、和WGS84坐标系之间的转换。包含84坐标系、02坐标系、百度地图、高德地图、腾讯地图坐标之间相互转换的算法
WGS84坐标转换百度坐标
dashen110的博客
10-11 1298
目前常用的地图坐标系有多种标准，几家主流地图(如百度地图，高德地图和QQ地图等)使用的坐标系标准也各不相同。假如您有GPS坐标，想在百度地图上显示；或者有百度地图的坐标，但想在微信地图上显示，这时就需要使用地图坐标系工具来转换坐标，在目标地图上正确显示。本站提供各主流地图之间的坐标系转换工具，供您使用。
如何实现wgs84(ep4326)坐标系与百度坐标系，火星坐标系的互转，java与JavaScript都有
weixin_42759398的博客
04-03 1332
国家规定：国内出版的各种地图系统(包括电子形式)，必须至少采用GCJ-02对地理位置进行首次加密(高德地图、腾讯地图等)；工作中用到百度地图和高德地图服务，必须实现转换就成了难题了，现在它来了，java后台和纯js方法都有，直接上代码。众所周知国内的地址服务不让直接用wgs84，百度使用的是百度坐标系，高德地图使用的是火星坐标系。国际地图提供商使用的坐标系(谷歌地图国外、osm、mapbox)。百度标准，百度 SDK，百度地图，百度GeoCoding 使用；2.国测局坐标系(GCJ-02、火星坐标系)
WGS84坐标转BD-09坐标（wgs84转百度）
maohangfei.cn
05-29 5406
WGS84坐标转BD-09坐标（wgs84转百度）
地图--WGS84 高德坐标 百度坐标的相互转换（好用推荐）
VIP_CR的博客
09-21 1万+
地图 WGS84 高德坐标 百度坐标的相互转换
CGCS2000空间直角坐标与WGS84坐标互转源码
最新发布
01-13
博主在一次项目开发过程中，将CGCS2000坐标的dwg设计图转换为WGS84的kmz数据，并加载至cesium球上，当时业主提出，是否可以在cesium球上结合加载的设计图数据，对设计图数据进行指导。  因业主非专业的测绘和GIS...
（下载链接在博客）ArcGIS工具箱实现百度、火星与WGS84坐标互转（Excel版）
07-25
百度坐标、火星坐标、WGS84坐标转换工具全部升级为脱离ArcGIS环境，可在本人博客文末找到链接试用，链接呈上https://blog.csdn.net/yhiki/article/details/86476276
js中火星坐标、百度坐标、WGS84坐标转换实现方法示例
10-15
主要介绍了js中火星坐标、百度坐标、WGS84坐标转换实现方法,涉及JavaScript数值计算相关操作技巧,需要的朋友可以参考下
wgs84转百度坐标
01-13
项目中使用Echarts 地图，该地图百度已经停止维护了。由于项目需求，新增了许多行政区划区， 比如：高新区、经济开发区等 都需要显示在Echarts地图上，但Echarts地图使用的是百度的坐标系， 而我们只有wgs84的投影坐标文件。如果不加以转换则会造成严重的失准，故需要进行转化。
坐标转换（支持百度坐标系，火星坐标系，WGS84之间的互转）
02-27
坐标转换小工具，支持百度坐标系，火星坐标系，WGS84之间的互转，批量转换，非常实用的小工具，对于GIS互联网行业地图很有帮助
百度坐标转WGS84坐标程序
11-29
提供百度坐标与WGS84，高德、谷歌坐标与WGS84坐标的转换，提供百度、高德、谷歌与WGS84、火星之间的转换
百度、搞得、WGS84坐标转换 js java
09-17
百度、高得、0WGS84坐标转换 js java, 使用网上找的js的代码，自己添加了一个java的代码，方便以后使用
java 地心坐标系（ECEF）和WGS-84坐标系（WGS84）互转的实现
08-25
本文主要探讨了如何在Java中实现地心坐标系（ECEF）与WGS-84坐标系之间的转换。这两种坐标系是地球坐标系统中广泛使用的标准。 地心坐标系（ECEF，Earth-Centered Earth-Fixed）是一种笛卡尔坐标系统，其中原点位于...
wgs84转百度坐标系
weixin_49409651的博客
02-02 1744
【代码】wgs84转百度坐标系。
wgs84转百度经纬度
zmyalh的博客
07-24 652
wgs84转百度经纬度
wgs84 转百度地图坐标
lizhihua0625的博客
01-16 1249
<?php /* $lng = "112°47′23.16″"; $lat = "36°49′47.31″"; // 度分秒转十进制度 $lng = CoordinateTransform::degreeMinuteSecondToDecimal($lng); $lat = CoordinateTransform::degreeMinuteSecondToDecimal($lat); // var_dump(array($lng,$lat)); $lng='111.655390'; $lat=
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值