Java坐标转换的多元实现路径：在线调用、百度与高德地图API集成及纯Java代码实现——纯Java代码实现与数学模型深度剖析

最新推荐文章于 2025-10-04 14:23:41 发布

原创最新推荐文章于 2025-10-04 14:23:41 发布 · 422 阅读

5 ·

CC 4.0 BY-SA版权

文章标签：

#java #开发语言

关键词：Java坐标转换、纯Java实现、WGS84、GCJ-02、BD-09、数学模型、工程级代码

一、为何需要纯Java实现？

离线环境：政务内网、军工系统、航空器
高并发：API限流、QPS瓶颈
数据保密：坐标不能出境、不能外传
成本控制：海量坐标转换，API费用高昂

二、核心数学模型：WGS84 → GCJ-02 → BD-09

1. WGS84 → GCJ-02（国测局加密）

加密逻辑：

private static final double PI = 3.14159265358979324;
private static final double A = 6378245.0;     // 长半轴
private static final double EE = 0.00669342162296594323; // 扁率

public static double[] wgs84ToGcj02(double lng, double lat) {
    if (outOfChina(lng, lat)) return new double[]{lng, lat};

    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) / ((A * (1 - EE)) / (magic * sqrtMagic) * PI);
    dLng = (dLng * 180.0) / (A / sqrtMagic * Math.cos(radLat) * PI);

    double mgLat = lat + dLat;
    double mgLng = lng + dLng;

    return new double[]{mgLng, mgLat};
}

辅助函数：

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;
}

private static boolean outOfChina(double lng, double lat) {
    return lng < 72.004 || lng > 137.8347 || lat < 0.8293 || lat > 55.8271;
}

该算法为国测局非公开加密算法，通过多项式+三角函数拟合偏移量，精度达1米级。

2. GCJ-02 → BD-09（百度二次加密）

public static double[] gcj02ToBd09(double lng, double lat) {
    double z = Math.sqrt(lng * lng + lat * lat) + 0.00002 * Math.sin(lat * PI * 3000.0 / 180.0);
    double theta = Math.atan2(lat, lng) + 0.000003 * Math.cos(lng * PI * 3000.0 / 180.0);
    double bdLng = z * Math.cos(theta) + 0.0065;
    double bdLat = z * Math.sin(theta) + 0.006;
    return new double[]{bdLng, bdLat};
}

3. BD-09 → GCJ-02（逆向解密）

public static double[] bd09ToGcj02(double bdLng, double bdLat) {
    double x = bdLng - 0.0065;
    double y = bdLat - 0.006;
    double z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * PI * 3000.0 / 180.0);
    double theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * PI * 3000.0 / 180.0);
    double gcjLng = z * Math.cos(theta);
    double gcjLat = z * Math.sin(theta);
    return new double[]{gcjLng, gcjLat};
}

三、工程级封装：CoordinateTransformUtil

public final class CoordinateTransformUtil {

    public enum CoordType {
        WGS84, GCJ02, BD09
    }

    public static double[] transform(double lng, double lat, CoordType from, CoordType to) {
        if (from == to) return new double[]{lng, lat};

        // WGS84 → GCJ02
        if (from == CoordType.WGS84 && to == CoordType.GCJ02) {
            return wgs84ToGcj02(lng, lat);
        }
        // GCJ02 → BD09
        if (from == CoordType.GCJ02 && to == CoordType.BD09) {
            return gcj02ToBd09(lng, lat);
        }
        // BD09 → GCJ02
        if (from == CoordType.BD09 && to == CoordType.GCJ02) {
            return bd09ToGcj02(lng, lat);
        }
        // WGS84 → BD09
        if (from == CoordType.WGS84 && to == CoordType.BD09) {
            double[] gcj = wgs84ToGcj02(lng, lat);
            return gcj02ToBd09(gcj[0], gcj[1]);
        }
        // BD09 → WGS84
        if (from == CoordType.BD09 && to == CoordType.WGS84) {
            double[] gcj = bd09ToGcj02(lng, lat);
            return gcj02ToWgs84(gcj[0], gcj[1]);
        }
        throw new IllegalArgumentException("Unsupported transform");
    }
}