原始GPS座標（WGS-84）轉成百度的經緯度，零誤差

前言：在博客上找了很多，也看了很多瀏覽度高的，但是一測試，還是有幾百米或者幾十米的誤差

於是我就借鑑幾個博客做了調整，經度度與百度手機定位幾乎一致。

借鑑博客地址：https://blog.csdn.net/coolypf/article/details/8569813

                         還有一個博客：找不到了，看到很多博客，都是轉載或者抄襲借鑑，暫時沒找到原文，找到了再放上來

 

一、84 to 火星座標系 (GCJ-02) to 百度地圖座標

/**
 * @Author zr
 * @Date 2019-08-29 17:03
 * @Version 2.0
 */
public class BdMapUtil {
    /** 圓周率 */
    public static double pi = 3.1415926535897932384626;
    public static double a = 6378245.0;
    public static double ee = 0.00669342162296594323;
    private static double x_pi = 3.14159265358979324 * 3000.0 / 180.0;

    /**
     * 得到百度地圖經緯度
     *
     * @author zr
     * @date 2019-09-09
     * @param lat WGS-84 gps的緯度
     * @param lon WGS-84 gps的經度
     * @return com.construn.vehicle.message.util.Gps
     */
    public static Gps getBd09(double lat, double lon){
        //得到火星座標經緯度
        Gps gps = gps84ToGcj02(lat, lon);
        if(gps ==null){
            return null;
        }
        //得到百度座標經緯度
        return gcj02ToBd09(gps.getWgLat(), gps.getWgLon());
    }
    /**
     * 火星座標系 (GCJ-02)轉百度地圖座標系
     *
     * @author zr
     * @date 2019-09-09
     * @param  * @param gg_lat
     * @param gg_lon
     * @return void
     */
    private static Gps gcj02ToBd09(double gg_lat, double gg_lon) {
        double x = gg_lon, y = gg_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 bd_lon = z * Math.cos(theta) + 0.0065;
        double bd_lat = z * Math.sin(theta) + 0.006;
        return new Gps(bd_lat,bd_lon);
    }
    /**
      * 84 to 火星座標系 (GCJ-02)
      *
      * @author zr
      * @date 2019-09-09
      * @param  * @param lat
      * @param lon
      * @return void
      */
    private static Gps gps84ToGcj02(double lat, double lon) {
        if (outOfChina(lat, lon)) {
            return null;
        }
        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 Gps(mgLat,mgLon);
    }

    /**
     * 判斷有沒有超過中國區
     *
     * @author zr
     * @date 2019-09-09
     * @param  * @param lat
     * @param lon
     * @return boolean
     */
    private 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;
    }

    /**
     * 處理緯度
     *
     * @author zr
     * @date 2019-09-09
     * @param  * @param x
     * @param y
     * @return double
     */
    private 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;
    }

    /**
     * 處理經度
     *
     * @author zr
     * @date 2019-09-09
     * @param  * @param x
     * @param y
     * @return double
     */
    private 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;
    }
}

Gps類

import lombok.Getter;
import lombok.Setter;

/**
 * gps 經緯度座標類
 *
 * @Author zr
 * @Date 2019-09-09 19:12
 * @Version 2.0
 */
@Getter
@Setter
public class Gps {
    private double wgLat;
    private double wgLon;

    public Gps() {
    }
    public Gps(double wgLat, double wgLon) {
        setWgLat(wgLat);
        setWgLon(wgLon);
    }
    @Override
    public String toString() {
        return "Gps{" +
                "wgLat=" + wgLat +
                ", wgLon=" + wgLon +
                '}';
    }
}

二、說明下、百度的地圖加密原始gps地理座標是這樣的：百度地圖加密的是 火星座標系 (GCJ-02)，而GCJ-02加密的是原始gps座標系。所以從原始gps座標，到百度地圖座標，順序是WGS-84 --> GCJ-02 --> BD-09