百度坐标转化成墨卡托坐标
根据网上其他博客的转化规则测试了一下,还是有问题,我根据自己项目所用的地方加以修改和测试了,最终在3d地图上展示了具体坐标位置。
下面具体介绍一下
火星坐标:火星坐标的真实名称应该是GCJ-02坐标。
百度坐标:火星坐标是在国际标准坐标WGS-84上进行的一次加密,由于国内的电子地图都要至少使用火星坐标进行一次加密,百度直接就任性一些,直接自己又研究了一套加密算法,来了个 二次加密,这就是我们所熟知的百度坐标(BD-09)。
WGS-84坐标:WGS-84坐标是一个国际的标准,一般卫星导航,原始的GPS设备中的数据都是采用这一坐标系。国外的Google地图、OSM等采用的都是这一坐标。
所以第一步:百度坐标转火星坐标(GCJ-02坐标)
第二步:火星坐标(GCJ-02坐标) 转 WGS-84坐标
中间做了一个判断是否为国内坐标,不是不做纠偏
第三步:WGS-84坐标转墨卡托坐标
代码如下:
import java.util.HashMap;
import java.util.Map;
public class BaiduToMkt {
/**
* 百度坐标转火星坐标(GCJ-02坐标)
*/
public static Map<String,Double> baiduTomars(Double lon, Double lat) {
double x_pi=3.14159265358979324 * 3000.0 / 180.0;
Map<String,Double> map = new HashMap<>();
double x = lon-0.0065;
double y = lat-0.006;
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 mars_lon=z * Math.cos(theta);
double mars_lat=z * Math.sin(theta);
map = GcjToWgs(mars_lon,mars_lat);
return map;
}
/**
* 火星坐标(GCJ-02坐标) 转 WGS-84坐标
*/
public static Map<String,Double> GcjToWgs(Double lng,Double lat){
double PI = 3.1415926535897932384626;
double a = 6378245.0;
double ee = 0.00669342162296594323;
Map<String,Double> map = new HashMap<>();
if (outOfChina(lng, lat)) {
map.put("lng",lng);
map.put("lat",lat);
return map;
} else {
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;
map.put("lng",lng * 2 - mglng);
map.put("lat",lat * 2 - mglat);
return map;
}
}
//是否在国内 不在国内不做偏移
public static Boolean outOfChina(Double lon, Double lat){
if ((lon < 72.004 || lon > 137.8347)&&(lat < 0.8293 || lat > 55.8271)){
return true;
}else {
return false;
}
}
// 纠偏Lat
public static Double transformlat(Double lng,Double lat){
double pi = 3.1415926535897932384626;
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;
}
// 纠偏Lon
public static Double transformlng(Double lng,Double lat){
double pi = 3.1415926535897932384626;
double 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;
}
// 转化墨卡托
public static Map<String,Double> LonlattoWebMercator(Double lon,Double lat) {
Map<String,Double> map = new HashMap<>();
map = baiduTomars(lon,lat);
double x = map.get("lng") / 180.0 * 20037508.34;
double y = Math.log(Math.tan((90 + map.get("lat")) * Math.PI / 360)) / (Math.PI / 180);
y = y * 20037508.34 / 180;
Map<String,Double> map1 = new HashMap<>();
map1.put("lon",x);
map1.put("lat",y);
return map1;
}
public static void main(String[] args){
double lon = 116.180928000000000;
double lat = 36.357202000000000;
System.out.println(LonlattoWebMercator(lon,lat));
}
}