Java【代码分享 04】坐标系说明+WGS84 GCJ02 BD09坐标系转换工具+Java坐标系转换及验证源代码分享（粘贴可用

BD09 是在 GCJ02 坐标系的基础上再次加密偏移后形成的坐标系，官方解释为进一步保护用户隐私，当前只适用于百度地图。

GPS经纬度是39°54'26.2"N 116°23'28.4"E，转化为度的单位就是39.907270 116.391213（小数部分 = 分 / 60 + 秒 / 3600）

2. 坐标系转换

• 国测局规定：互联网地图在国内必须至少使用 GCJ02 进行首次加密，不允许直接使用 WGS84 坐标下的地理数据，同时任何坐标系均不可转换为 WGS84 坐标。因此不存在将 GCJ02 坐标转换为 WGS84 坐标的官方转换方法。
• 目前百度 API 提供了从其它坐标系转换为 BD09 坐标系的 API，但却没有从 BD09 坐标系转为其他坐标系的API。

2.1 WGS84转换

WGS84转换测试工具 可将 WGS84 转换成 GCJ02、BD09ll 坐标系：

2.2 高德坐标系转换

高德坐标系转换工具 能够将用户输入的非高德坐标（GPS坐标、mapbar坐标、baidu坐标）转换成高德坐标：

以上的两个工具可以用来验证Java坐标系转换结果是否正确。

3. Java坐标系转换工具

package yz.com.javautil.utilpackage;


import yz.com.javautil.vo.TransformPoint;

import java.math.BigDecimal;
import java.math.MathContext;
import java.math.RoundingMode;
import java.util.HashMap;
import java.util.Map;
import java.util.Optional;
import java.util.function.Consumer;
import java.util.function.Supplier;


/\*\*
 \* 坐标系转换工具类
 \*
 \* @author Administrator
 \*/
public class CoordinateTransferUtils {

    /\*\*
 \* 国内坐标边界
 \*/
    private static final double MIN_LON = 72.004D;
    private static final double MAX_LON = 137.8347D;
    private static final double MIN_LAT = 0.8293D;
    private static final double MAX_LAT = 55.8271D;

    /\*\*
 \* PI 圆周率
 \*/
    private static final double PI = 3.14159265358979324D;

    /\*\*
 \* A WGS 长轴半径
 \*/
    private static final double A = 6378245.0D;

    /\*\*
 \* EE WGS 偏心率的平方
 \*/
    private static final double EE = 0.00669342162296594323D;


    /\*\*
 \* WGS84转换GCJ02适用于继承的情况
 \*
 \* @param t 继承TransformPoint类的子类对象
 \*/
    public static <T extends TransformPoint> void wgs84ToGcj02(T t) {
        wgs84ToGcj02(t::getLng, t::getLat, t::setTransformLng, t::setTransformLat);
    }

    /\*\*
 \* WGS84转换GCJ02核心方法
 \*
 \* @param fromLon 转换前的经度
 \* @param fromLat 转换前的纬度
 \* @return 转换后的经纬度map对象
 \*/
    public static Map<String, Double> wgs84ToGcj02(double fromLon, double fromLat) {
        HashMap<String, Double> transformRes = new HashMap<>(2);
        // 国外坐标不用进行加密
        if (outOfChina(fromLon, fromLat)) {
            transformRes.put("lon", fromLon);
            transformRes.put("lat", fromLat);
            return transformRes;
        }
        // 计算转换后的经纬度坐标
        double dLat = transformLat(fromLon - 105.0, fromLat - 35.0);
        double dLon = transformLon(fromLon - 105.0, fromLat - 35.0);
        double radLat = fromLat / 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 = fromLat + dLat;
        double mgLon = fromLon + dLon;
        transformRes.put("lon", new BigDecimal(mgLon + "", new MathContext(9, RoundingMode.HALF_UP)).doubleValue());
        transformRes.put("lat", new BigDecimal(mgLat + "", new MathContext(9, RoundingMode.HALF_UP)).doubleValue());

        return transformRes;
    }

    /\*\*
 \* GCJ02转换WGS84
 \*
 \* @param lon 转换前的经度
 \* @param lat 转换后的纬度
 \* @return 转换后的经纬度map对象
 \*/
    public static Map<String, Double> gcj02ToWgs84(double lon, double lat) {
        Map<String, Double> transformRes = new HashMap<>(2);
        double longitude = lon \* 2 - wgs84ToGcj02(lon, lat).get("lon");
        double latitude = lat \* 2 - wgs84ToGcj02(lon, lat).get("lat");
        transformRes.put("lon", longitude);
        transformRes.put("lat", latitude);
        return transformRes;
    }

    /\*\*
 \* GCJ02转换BD09
 \*
 \* @param gcjLat GCL纬度坐标
 \* @param gcjLng GCL经度坐标
 \*/
    public static Map<String, Double> Gcj02ToBd09(double gcjLat, double gcjLng) {
        Map<String, Double> transformRes = new HashMap<>(2);
        double z = Math.sqrt(gcjLng \* gcjLng + gcjLat \* gcjLat) + 0.00002 \* Math.sin(gcjLat \* PI);
        double theta = Math.atan2(gcjLat, gcjLng) + 0.000003 \* Math.cos(gcjLng \* PI);
        transformRes.put("lon", z \* Math.cos(theta) + 0.0065);
        transformRes.put("lat", z \* Math.sin(theta) + 0.006);
        return transformRes;
    }

    /\*\*
 \* BD09转换GCJ02
 \*
 \* @param bdLat 百度纬度坐标
 \* @param bdLng 百度经度坐标
 \*/
    public static Map<String, Double> bd09ToGcj02(double bdLat, double bdLng) {
        Map<String, Double> transformRes = new HashMap<>(2);
        double x = bdLng - 0.0065, y = bdLat - 0.006;
        double z = Math.sqrt(x \* x + y \* y) - 0.00002 \* Math.sin(y \* PI);
        double theta = Math.atan2(y, x) - 0.000003 \* Math.cos(x \* PI);

        transformRes.put("lon", z \* Math.cos(theta));
        transformRes.put("lat", z \* Math.sin(theta));
        return transformRes;
    }

    /\*\*
 \* WGS84转换GCJ02优化（处理转换前经纬度坐标为null问题）
 \*
 \* @param fromLon 转换前的经度
 \* @param fromLat 转换前的纬度
 \* @param toLon 转换后的经度
 \* @param toLat 转换后的纬度
 \*/
    private static void wgs84ToGcj02(Supplier<BigDecimal> fromLon, Supplier<BigDecimal> fromLat, Consumer<BigDecimal> toLon, Consumer<BigDecimal> toLat) {
        wgs84ToGcj02(Optional.ofNullable(fromLon.get()).orElse(BigDecimal.ZERO), Optional.ofNullable(fromLat.get()).orElse(BigDecimal.ZERO), toLon, toLat);
    }

    /\*\*
 \* WGS84转换GCJ02优化（处理转换后经纬度坐标）
 \*
 \* @param fromLon 转换前的经度
 \* @param fromLat 转换前的纬度
 \* @param toLon 转换后的经度
 \* @param toLat 转换后的纬度
 \*/
    private static void wgs84ToGcj02(BigDecimal fromLon, BigDecimal fromLat, Consumer<BigDecimal> toLon, Consumer<BigDecimal> toLat) {
        final Map<String, Double> transformRes = wgs84ToGcj02(fromLon.doubleValue(), fromLat.doubleValue());
        toLon.accept(new BigDecimal(transformRes.get("lon") + ""));
        toLat.accept(new BigDecimal(transformRes.get("lat") + ""));
    }

    /\*\*
 \* 坐标是否在国外
 \*
 \* @param lon 经度
 \* @param lat 纬度
 \* @return true 国外坐标 false 国内坐标
 \*/
    private static boolean outOfChina(double lon, double lat) {
        if (lon < MIN_LON || lon > MAX_LON) {
            return true;
        }
        return lat < MIN_LAT || lat > MAX_LAT;
    }

    /\*\*
 \* 转换经度坐标
 \*
 \* @param x 偏移后的经度
 \* @param y 偏移后的纬度
 \* @return double 转换经度坐标
 \*/
    private static double transformLat(double x, double y) {