/ **
*百度坐标(BD09)，国测局坐标(火星坐标，GCJ02)和WGS84坐标系之间的转换的工具
*
*参考/wandergis/coordtransform实现的Java版本
* @作者feihong
* /
公共类CoordinateTransformUtil {
static double x_pi = 3.14159265358979324 * 3000.0 / 180.0 ;
// π
static double pi = 3.1415926535897932384626 ;
//长半轴
静态double a = 6378245.0 ;
//扁率
静态双 ee = 0.00669342162296594323 ;
/ **
*百度坐标系(BD-09)转WGS坐标
*
* @param lng百度坐标纬度
* @param lat百度坐标经度
* @返回 WGS84坐标数组
* /
public static double [] bd09towgs84(double lng，double lat){
double [] gcj = bd09togcj02(lng，lat);
double [] wgs84 = gcj02towgs84(gcj [ 0 ]，gcj [ 1 ]);
返回 wgs84;
}
/ **
* WGS坐标转百度坐标系(BD-09)
*
* @param lng WGS84坐标系的经度
* @param lat WGS84坐标系的纬度
* @return百度坐标数组
* /
public static double [] wgs84tobd09(double lng，double lat){
double [] gcj = wgs84togcj02(lng，lat);
double [] bd09 = gcj02tobd09(gcj [ 0 ]，gcj [ 1 ]);
返回 bd09;
}
/ **
*火星坐标系(GCJ-02)转百度坐标系(BD-09)
*
*谷歌，高德 - >百度
* @param lng火星坐标经度
* @param lat火星坐标纬度
* @return百度坐标数组
* /
public static double [] gcj02tobd09(double lng，double lat){
双 z = 数学。sqrt(lng * lng + lat * lat)+ 0.00002 * Math 。sin(lat * x_pi);
double theta = 数学。atan2(lat，lng)+ 0.000003 * Math 。cos(lng * x_pi);
double bd_lng = z * Math 。cos(θ)+ 0.0065 ;
double bd_lat = z * Math 。罪(THETA)+ 0.006 ;
返回新的double [] {bd_lng，bd_lat};
}
/ **
*百度坐标系(BD-09)转火星坐标系(GCJ-02)
*
*百度 - >谷歌，高德
* @param bd_lon百度坐标纬度
* @param bd_lat百度坐标经度
* @return火星坐标数组
* /
public static double [] bd09togcj02(double bd_lon，double bd_lat){
双 X = bd_lon - 0.0065 ;
double y = bd_lat - 0.006 ;
双 z = 数学。sqrt(x * x + y * y)- 0.00002 * Math 。sin(y * x_pi);
double theta = 数学。atan2(y，x)- 0.000003 * Math 。cos(x * x_pi);
double gg_lng = z * Math 。COS(THETA);
double gg_lat = z * Math 。罪(THETA);
返回新的double [] {gg_lng，gg_lat};
}
/ **
* WGS84转GCJ02(火星坐标系)
*
* @param lng WGS84坐标系的经度
* @param lat WGS84坐标系的纬度
* @return火星坐标数组
* /
public static double [] wgs84togcj02(double lng，double lat){
if(out_of_china(lng，lat)){
返回新的double [] {lng，lat};
}
double dlat = transformlat(lng - 105.0，lat - 35.0);
double dlng = transformlng(lng - 105.0，lat - 35.0);
double radlat = lat / 180.0 * pi;
双重魔法= 数学。罪(radlat);
魔法= 1 - EE *魔*魔;
double sqrtmagic = Math 。SQRT(魔术);
dlat =(dlat * 180.0)/((a *(1 - ee))/(magic * sqrtmagic)* pi);
DLNG =(DLNG * 180.0)/(A / sqrtmagic * 数学。 COS(radlat)* PI);
double mglat = lat + dlat;
双 mglng = lng + dlng;
返回新的double [] {mglng，mglat};
}
/ **
* GCJ02(火星坐标系)转GPS84
*
* @param lng火星坐标系的经度
* @param lat火星坐标系纬度
* @返回 WGS84坐标数组
* /
public static double [] gcj02towgs84(double lng，double lat){
if(out_of_china(lng，lat)){
返回新的double [] {lng，lat};
}
double dlat = transformlat(lng - 105.0，lat - 35.0);
double dlng = transformlng(lng - 105.0，lat - 35.0);
double radlat = lat / 180.0 * pi;
双重魔法= 数学。罪(radlat);
魔法= 1 - EE *魔*魔;
double sqrtmagic = Math 。SQRT(魔术);
dlat =(dlat * 180.0)/((a *(1 - ee))/(magic * sqrtmagic)* pi);
DLNG =(DLNG * 180.0)/(A / sqrtmagic * 数学。 COS(radlat)* PI);
double mglat = lat + dlat;
双 mglng = lng + dlng;
返回新的double [] {lng * 2 - mglng，lat * 2 - mglat};
}
/ **
*纬度转换
*
* @param lng
* @param lat
* @返回
* /
public static double transformlat(double lng，double lat){
double ret = - 100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + 0.1 * lng * lat + 0.2 * Math 。SQRT(数学。 ABS(LNG));
RET + =(20.0 * 数学。 SIN(6.0 * LNG * PI)+ 20.0 * 数学。 SIN(2.0 * LNG * PI))* 2.0 / 3.0 ;
RET + =(20.0 * 数学。罪(LAT * PI)+ 40.0 * 数学。罪(LAT / 3.0 * PI))* 2.0 / 3.0 ;
RET + =(160.0 * 数学。罪(LAT / 12.0 * PI)+ 320 * 数学。罪(LAT * PI / 30.0))* 2.0 / 3.0 ;
返回 ret;
}
/ **
*经度转换
*
* @param lng
* @param lat
* @返回
* /
public static double transformlng(double lng，double lat){
double ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + 0.1 * lng * lat + 0.1 * Math 。SQRT(数学。 ABS(LNG));
RET + =(20.0 * 数学。 SIN(6.0 * LNG * PI)+ 20.0 * 数学。 SIN(2.0 * LNG * PI))* 2.0 / 3.0 ;
RET + =(20.0 * 数学。罪(LNG * PI)+ 40.0 * 数学。罪(LNG / 3.0 * PI))* 2.0 / 3.0 ;
RET + =(150.0 * 数学。罪(LNG / 12.0 * PI)+ 300.0 * 数学。罪(LNG / 30.0 * PI))* 2.0 / 3.0 ;
返回 ret;
}
/ **
*判断是否在国内，不在国内不做偏移
*
* @param lng
* @param lat
* @返回
* /
public static boolean out_of_china(double lng，double lat){
如果(lng < 72.004 || lng > 137.8347){
返回true ;
} else if(lat < 0.8293 || lat > 55.8271){
返回true ;
}
返回false ;
}
}
