最近做项目需要百度地图坐标转换到高德地图坐标,高德官方也给出了转换接口(百度地图也给出了转换接口)
http://lbs.amap.com/api/javascript-api/reference/lnglat-to-address#t1
但是还想想直接通过后台一次性处理好坐标,因此需要找到坐标转换算法,封装起来调用。
1 百度官方对百度坐标为何有偏移的解释
国际经纬度坐标标准为WGS-84,国内必须至少使用国测局制定的GCJ-02,对地理位置进行首次加密。百度坐标在此基础上,进行了BD-09二次加密措施,更加保护了个人隐私。百度对外接口的坐标系并不是GPS采集的真实经纬度,需要通过坐标转换接口进行转换。
2 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换算法
GCJ-02(火星坐标) 和 BD-09 (百度坐标)
算法代码如下,其中
封装经纬度:
/** 经纬度点封装
* Created by 明明如月 on 2017-03-22.
*/
public class LngLat {
private double longitude;//经度
private double lantitude;//维度
public LngLat() {
}
public LngLat(double longitude, double lantitude) {
this.longitude = longitude;
this.lantitude = lantitude;
}
public double getLongitude() {
return longitude;
}
public void setLongitude(double longitude) {
this.longitude = longitude;
}
public double getLantitude() {
return lantitude;
}
public void setLantitude(double lantitude) {
this.lantitude = lantitude;
}
@Override
public String toString() {
return "LngLat{" +
"longitude=" + longitude +
", lantitude=" + lantitude +
'}';
}
}
/** 百度地图坐标和火星坐标转换
* Created by 明明如月 on 2017-03-22.
*/
public class CoodinateCovertor {
private static double x_pi = 3.14159265358979324 * 3000.0 / 180.0;
/**
* 对double类型数据保留小数点后多少位
* 高德地图转码返回的就是 小数点后6位,为了统一封装一下
* @param digit 位数
* @param in 输入
* @return 保留小数位后的数
*/
static double dataDigit(int digit,double in){
return new BigDecimal(in).setScale(6, BigDecimal.ROUND_HALF_UP).doubleValue();
}
/**
* 将火星坐标转变成百度坐标
* @param lngLat_gd 火星坐标(高德、腾讯地图坐标等)
* @return 百度坐标
*/
public static LngLat bd_encrypt(LngLat lngLat_gd)
{
double x = lngLat_gd.getLongitude(), y = lngLat_gd.getLantitude();
double z = sqrt(x * x + y * y) + 0.00002 * sin(y * x_pi);
double theta = atan2(y, x) + 0.000003 * cos(x * x_pi);
return new LngLat(dataDigit(6,z * cos(theta) + 0.0065),dataDigit(6,z * sin(theta) + 0.006));
}
/**
* 将百度坐标转变成火星坐标
* @param lngLat_bd 百度坐标(百度地图坐标)
* @return 火星坐标(高德、腾讯地图等)
*/
static LngLat bd_decrypt(LngLat lngLat_bd)
{
double x = lngLat_bd.getLongitude() - 0.0065, y = lngLat_bd.getLantitude() - 0.006;
double z = sqrt(x * x + y * y) - 0.00002 * sin(y * x_pi);
double theta = atan2(y, x) - 0.000003 * cos(x * x_pi);
return new LngLat( dataDigit(6,z * cos(theta)),dataDigit(6,z * sin(theta)));
}
//测试代码
public static void main(String[] args) {
LngLat lngLat_bd = new LngLat(120.153192,30.25897);
System.out.println(bd_decrypt(lngLat_bd));
}
}
经验证和高德地图自带接口转换出效果一致。