最近做项目需要百度地图坐标转换到高德地图坐标，高德官方也给出了转换接口(百度地图也给出了转换接口)
/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(百度坐标)
算法代码如下，其中bd_encrypt将GCJ-02坐标转换成BD-09坐标，bd_decrypt反之。
封装经纬度：
/** 经纬度点封装
* Created by 明明如月 on -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 -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));
}
}
经验证和高德地图自带接口转换出效果一致。
