通过js实现地理信息坐标WGS84、GCJ-02、BD-09坐标互转

时间：2024-11-16 08:36:49

WGS84坐标系

WGS84坐标系是国际上定义的统一的经纬度坐标规范，目前常见的地图应用有谷歌地图、天地图，这些地图都是用的WGS84坐标系。

GCJ-02坐标系

GCJ-02坐标系是中国国家测绘局定制的一套地理信息坐标系，它是在WGS84基础上进行了加密处理，GCJ-02坐标系也称之为火星坐标系；目前常见的地图应用有高德地图、腾讯地图，这些地图都以GCJ-02为坐标参考系。

BD-09坐标系

BD-09坐标系是百度定制的一套加密规则，BD-09在GCJ-02的基础上进行加密处理，进一步提高了地图数据的安全性，应用于百度地图。

了解了这三种坐标系，自然就能够理解它们之间的坐标区别，所以如果用一个经纬度要在各自的地图上标点就会出现坐标精度有偏差，那么就需要将该坐标转换成对应坐标系的经纬度坐标，这样绘制在地图上才不会有偏差。

WGS84转GCJ-02

function transformLatWithXY(x, y) {
    let pi = 3.14159265358979324;
    let lat = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * Math.sqrt(Math.abs(x));
    lat += (20.0 * Math.sin(6.0 * x * pi) + 20.0 * Math.sin(2.0 * x * pi)) * 2.0 / 3.0;
    lat += (20.0 * Math.sin(y * pi) + 40.0 * Math.sin(y / 3.0 * pi)) * 2.0 / 3.0;
    lat += (160.0 * Math.sin(y / 12.0 * pi) + 320 * Math.sin(y * pi / 30.0)) * 2.0 / 3.0;
    return lat;
}

function transformLonWithXY(x, y) {
    let pi = 3.14159265358979324;
    let lon = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1 * Math.sqrt(Math.abs(x));
    lon += (20.0 * Math.sin(6.0 * x * pi) + 20.0 * Math.sin(2.0 * x * pi)) * 2.0 / 3.0;
    lon += (20.0 * Math.sin(x * pi) + 40.0 * Math.sin(x / 3.0 * pi)) * 2.0 / 3.0;
    lon += (150.0 * Math.sin(x / 12.0 * pi) + 300.0 * Math.sin(x / 30.0 * pi)) * 2.0 / 3.0;
    return lon;
}


/**
 *  将WGS-84(国际标准)转为GCJ-02(火星坐标):
 */
function transformFromWGSToGCJ(longitude,latitude) {
    let lat = "";
    let lon = "";
    let ee = 0.00669342162296594323;
    let a = 6378245.0;
    let pi = 3.14159265358979324;
    let adjustLat = transformLatWithXY(longitude - 105.0, latitude - 35.0);
    let adjustLon = transformLonWithXY(longitude - 105.0, latitude - 35.0);
    let radLat = latitude / 180.0 * pi;
    let magic = Math.sin(radLat);
    magic = 1 - ee * magic * magic;
    let sqrtMagic = Math.sqrt(magic);
    adjustLat = (adjustLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * pi);
    adjustLon = (adjustLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * pi);
    lon = longitude + adjustLon;
    lat = latitude + adjustLat;
    
    return {  longitude: lon,latitude: lat };

}

GCJ-02转BD-09

/**
 *  将GCJ-02(火星坐标)转为百度坐标:
 */
function transformFromGCJToBaidu(longitude,latitude) {
    let pi = 3.14159265358979324 * 3000.0 / 180.0;

    let z = Math.sqrt(longitude * longitude + latitude * latitude) + 0.00002 * Math.sin(latitude * pi);
    let theta = Math.atan2(latitude, longitude) + 0.000003 * Math.cos(longitude * pi);
    let a_longitude = (z * Math.cos(theta) + 0.0065);
    let a_latitude = (z * Math.sin(theta) + 0.006);

    return { longitude: a_longitude, latitude: a_latitude };
}

WGS84转BD-09

/**
 *  将WGS-84转为百度:
 */

function transformWGS84toBaidu(longitude,latitude){
    if (!isNaN(latitude)&&!isNaN(longitude)){
        let lat = parseFloat(latitude);
        let lon = parseFloat(longitude);
        let res= transformFromWGSToGCJ(lon,lat);
        let resBaidu = transformFromGCJToBaidu(res.longitude,res.latitude);
        return {longitude: resBaidu.longitude,latitude: resBaidu.latitude};
    } else{
        return { longitude: longitude, latitude: latitude };
    }

}

相关文章