gis 从WGS84转百度

var x_pi = 3.14159265358979324 * 3000.0 / 180.0;
var pi = 3.1415926535897932384626;
var a = 6378245.0;
var ee = 0.00669342162296594323;

//判断是否在国内，不在国内不做偏移
function out_of_china(lng, lat) {
return(lng > 73.66 && lng < 135.05 && lat > 3.86 && lat < 53.55) ? true : false;
}

function _transformlng(lng, lat) {
ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng +
0.1 * lng * lat + 0.1 * 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(lng * pi) + 40.0 *
Math.sin(lng / 3.0 * pi)) * 2.0 / 3.0
ret += (150.0 * Math.sin(lng / 12.0 * pi) + 300.0 *
Math.sin(lng / 30.0 * pi)) * 2.0 / 3.0
return ret;
}

function wgs84_to_gcj02(lng, lat) {
if(!out_of_china(lng, lat)) {
return [lng, lat];
} else {
var dlat = _transformlat(lng - 105.0, lat - 35.0);
var dlng = _transformlng(lng - 105.0, lat - 35.0);
var radlat = lat / 180.0 * pi;
var magic = Math.sin(radlat);
magic = 1 - ee * magic * magic;
var sqrtmagic = Math.sqrt(magic);
dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi);
dlng = (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * pi);
var mglat =accAdd(lat,dlat) ;
var mglng = accAdd(lng,dlng);
console.log(mglng + "------" + mglat)
return [mglng, mglat];
}
}

//ji算俩值相加高精确

function accAdd(data1, data2) {
var r1, r2, m, c;
try {
r1 = data1.toString().split(".")[1].length;
} catch(e) {
r1 = 0;
}
try {
r2 = data2.toString().split(".")[1].length;
} catch(e) {
r2 = 0;
}
c = Math.abs(r1 - r2);
m = Math.pow(10, Math.max(r1, r2));
if(c > 0) {
var cm = Math.pow(10, c);
if(r1 > r2) {
data1 = Number(data1.toString().replace(".", ""));
data2 = Number(data2.toString().replace(".", "")) * cm;
} else {
data1 = Number(data1.toString().replace(".", "")) * cm;
data2 = Number(data2.toString().replace(".", ""));
}
} else {
data1 = Number(data1.toString().replace(".", ""));
data2 = Number(data2.toString().replace(".", ""));
}
return(data1 + data2) / m;
}

//火星坐标转换百度坐标
function gcj02tobd(lng, lat) {
console.log(lng + '+++++' + lat);
var pts = wgs84_to_gcj02(lng, lat);
lng = pts[0];
lat = pts[1];
console.log(lng + '*****' + lat);
// var x_PI = 3.14159265358979324 * 3000.0 / 180.0;
// var PI = 3.1415926535897932384626;
// var a = 6378245.0;
// var ee = 0.00669342162296594323;
var z = Math.sqrt(lng * lng + lat * lat) + 0.00002 * Math.sin(lat * x_pi);
var theta = Math.atan2(lat, lng) + 0.000003 * Math.cos(lng * x_pi);
var bd_lng = z * Math.cos(theta) + 0.0065;
var bd_lat = z * Math.sin(theta) + 0.006;
return [bd_lng, bd_lat]
}

找了一堆资料终于对上了。。记录一下