wgstobd09java_GPS经纬度转百度经纬度JAVA

第一种&＃xff1a;JAVA写法(考虑double精度问题)package com.thinkgem.jeesite.common.utils;

import com.alibaba.dubbo.common.json.JSONObject;

import java.math.BigDecimal;

/**

* 百度坐标(BD09)、国测局坐标(火星坐标&＃xff0c;GCJ02)、和WGS84坐标系之间的转换的工具

*/

public class CoordinateTransformUtil {

static String BAIDU_GEOCONV_KEY &＃61; "XXXXXXXXXXXXXXXXXXXXXXXX";

static String baiduUrl &＃61; "http://api.map.baidu.com/geocoder/v2";

static int B_DIV_SCALE &＃61; 15;

static BigDecimal x_pi &＃61; new BigDecimal("3.14159265358979324").multiply(new BigDecimal("3000")).divide(new BigDecimal("180"), B_DIV_SCALE, BigDecimal.ROUND_HALF_UP);

// π

static BigDecimal pi &＃61; new BigDecimal("3.1415926535897932384626");

// 长半轴

static BigDecimal a &＃61; new BigDecimal("6378245");

// 扁率

static BigDecimal ee &＃61; new BigDecimal("0.00669342162296594323");

static BigDecimal b_1 &＃61; new BigDecimal("-1");

static BigDecimal b_2 &＃61; new BigDecimal("10");

static BigDecimal b_3 &＃61; new BigDecimal("2");

static BigDecimal b_4 &＃61; new BigDecimal("105");

static BigDecimal b_5 &＃61; new BigDecimal("3");

static BigDecimal b_6 &＃61; new BigDecimal("35");

static BigDecimal b_7 &＃61; new BigDecimal("160");

static BigDecimal b_8 &＃61; new BigDecimal("1");

static BigDecimal b_9 &＃61; new BigDecimal("150");

static BigDecimal b_10 &＃61; new BigDecimal("180");

static BigDecimal b_11 &＃61; new BigDecimal("0.00002");

static BigDecimal b_12 &＃61; new BigDecimal("0.000003");

static BigDecimal b_13 &＃61; new BigDecimal("0.0065");

static BigDecimal b_14 &＃61; new BigDecimal("0.006");

/**

* WGS坐标转百度坐标系(BD-09)

*

* &＃64;param lng WGS84坐标系的经度

* &＃64;param lat WGS84坐标系的纬度

* &＃64;return 百度坐标数组

*/

public static double[] wgs84tobd09(double lng, double lat) {

double[] gcj &＃61; wgs84togcj02(lng, lat);

double[] bd09 &＃61; gcj02tobd09(gcj[0], gcj[1]);

return bd09;

}

/**

* 火星坐标系(GCJ-02)转百度坐标系(BD-09)

*

* &＃64;param lng 火星坐标经度

* &＃64;param lat 火星坐标纬度

* &＃64;return 百度坐标数组

* &＃64;see 谷歌、高德——>百度

*/

public static double[] gcj02tobd09(double lng, double lat) {

BigDecimal b_lng &＃61; new BigDecimal(lng);

BigDecimal b_lat &＃61; new BigDecimal(lat);

BigDecimal z &＃61; new BigDecimal(Math.sqrt(b_lng.multiply(b_lng).add(b_lat.multiply(b_lat)).doubleValue())).add(b_11.multiply(new BigDecimal(Math.sin(b_lat.multiply(x_pi).doubleValue()))));

BigDecimal theta &＃61; new BigDecimal(Math.atan2(b_lat.doubleValue(), b_lng.doubleValue())).add(b_12.multiply(new BigDecimal(Math.cos(b_lng.multiply(x_pi).doubleValue()))));

BigDecimal bd_lng &＃61; z.multiply(new BigDecimal(Math.cos(theta.doubleValue()))).add(b_13);

// z * Math.cos(theta) &＃43; 0.0065;

BigDecimal bd_lat &＃61; z.multiply(new BigDecimal(Math.sin(theta.doubleValue()))).add(b_14);

return new double[]{bd_lng.doubleValue(), bd_lat.doubleValue()};

}

/**

* WGS84转GCJ02(火星坐标系)

*

* &＃64;param lng WGS84坐标系的经度

* &＃64;param lat WGS84坐标系的纬度

* &＃64;return 火星坐标数组

*/

public static double[] wgs84togcj02(double lng, double lat) {

BigDecimal b_lng &＃61; new BigDecimal(lng &＃43; "");

BigDecimal b_lat &＃61; new BigDecimal(lat &＃43; "");

if (out_of_china(b_lng, b_lat)) {

return new double[]{b_lng.doubleValue(), b_lat.doubleValue()};

}

BigDecimal dlat &＃61; transformlat(b_lng.subtract(b_4), b_lat.subtract(b_6));

BigDecimal dlng &＃61; transformlng(b_lng.subtract(b_4), b_lat.subtract(b_6));

BigDecimal radlat &＃61; b_lat.divide(b_10, B_DIV_SCALE, BigDecimal.ROUND_HALF_UP).multiply(pi);

BigDecimal magic &＃61; new BigDecimal(Math.sin(radlat.doubleValue()));

magic &＃61; b_8.subtract(ee.multiply(magic).multiply(magic));

BigDecimal sqrtmagic &＃61; new BigDecimal(Math.sqrt(magic.doubleValue()));

dlat &＃61; (dlat.multiply(b_10)).divide((a.multiply(b_8.subtract(ee))).divide((magic.multiply(sqrtmagic)), B_DIV_SCALE, BigDecimal.ROUND_HALF_UP).multiply(pi), B_DIV_SCALE, BigDecimal.ROUND_HALF_UP);

dlng &＃61; (dlng.multiply(b_10)).divide((a.divide(sqrtmagic, B_DIV_SCALE, BigDecimal.ROUND_HALF_UP).multiply(new BigDecimal(Math.cos(radlat.doubleValue())).multiply(pi))), B_DIV_SCALE, BigDecimal.ROUND_HALF_UP);

BigDecimal mglat &＃61; b_lat.add(dlat);

BigDecimal mglng &＃61; b_lng.add(dlng);

return new double[]{mglng.doubleValue(), mglat.doubleValue()};

}

/**

* 纬度转换

*

* &＃64;param lng

* &＃64;param lat

* &＃64;return

*/

public static BigDecimal transformlat(BigDecimal lng, BigDecimal lat) {

BigDecimal bet &＃61; b_1.multiply(b_2).multiply(b_2)

.add(b_3.multiply(lng))

.add(b_5.multiply(lat))

.add(b_3.divide(b_2, B_DIV_SCALE, BigDecimal.ROUND_HALF_UP).multiply(lat).multiply(lat))

.add(b_8.divide(b_2, B_DIV_SCALE, BigDecimal.ROUND_HALF_UP).multiply(lng).multiply(lat))

.add(b_3.divide(b_2, B_DIV_SCALE, BigDecimal.ROUND_HALF_UP).multiply(new BigDecimal(Math.sqrt(Math.abs(lng.doubleValue())))));

bet &＃61; bet.add((b_3.multiply(b_2).multiply(new BigDecimal(Math.sin(b_3.multiply(b_5).multiply(lng).multiply(pi).doubleValue())))

.add(b_3.multiply(b_2).multiply(new BigDecimal(Math.sin(b_3.multiply(lng).multiply(pi).doubleValue()))))).multiply(b_3).divide(b_5, B_DIV_SCALE, BigDecimal.ROUND_HALF_UP));

bet &＃61; bet.add((b_3.multiply(b_2).multiply(new BigDecimal(Math.sin(lat.multiply(pi).doubleValue())))

.add(b_3.multiply(b_3).multiply(b_2).multiply(new BigDecimal(Math.sin(lat.divide(b_5, B_DIV_SCALE, BigDecimal.ROUND_HALF_UP).multiply(pi).doubleValue()))))).multiply(b_3).divide(b_5, B_DIV_SCALE, BigDecimal.ROUND_HALF_UP));

bet &＃61; bet.add((b_7.multiply(new BigDecimal(Math.sin(lat.divide(b_5.multiply(b_3).multiply(b_3), B_DIV_SCALE, BigDecimal.ROUND_HALF_UP).multiply(pi).doubleValue())))

.add(b_7.multiply(b_3).multiply(new BigDecimal(Math.sin(lat.multiply(pi).divide(b_5.multiply(b_2), B_DIV_SCALE, BigDecimal.ROUND_HALF_UP).doubleValue()))))).multiply(b_3).divide(b_5, B_DIV_SCALE, BigDecimal.ROUND_HALF_UP));

return bet;

}

/**

* 经度转换

*

* &＃64;param lng

* &＃64;param lat

* &＃64;return

*/

public static BigDecimal transformlng(BigDecimal lng, BigDecimal lat) {

BigDecimal bet &＃61; b_2.multiply(b_2).multiply(b_5).add(lng).add(b_3.multiply(lat)).add(b_8.divide(b_2, B_DIV_SCALE, BigDecimal.ROUND_HALF_UP).multiply(lng).multiply(lng))

.add(b_8.divide(b_2, B_DIV_SCALE, BigDecimal.ROUND_HALF_UP).multiply(lng).multiply(lat)).add(b_8.divide(b_2, B_DIV_SCALE, BigDecimal.ROUND_HALF_UP).multiply(new BigDecimal(Math.sqrt(Math.abs(lng.doubleValue())))));

bet &＃61; bet.add((b_3.multiply(b_2).multiply(new BigDecimal(Math.sin(b_3.multiply(b_5).multiply(lng).multiply(pi).doubleValue())))

.add(b_3.multiply(b_2).multiply(new BigDecimal(Math.sin(b_3.multiply(lng).multiply(pi).doubleValue()))))).multiply(b_3).divide(b_5, B_DIV_SCALE, BigDecimal.ROUND_HALF_UP));

bet &＃61; bet.add((b_3.multiply(b_2).multiply(new BigDecimal(Math.sin(lng.multiply(pi).doubleValue())))

.add(b_3.multiply(b_3).multiply(b_2).multiply(new BigDecimal(Math.sin(lng.divide(b_5, B_DIV_SCALE, BigDecimal.ROUND_HALF_UP).multiply(pi).doubleValue()))))).multiply(b_3).divide(b_5, B_DIV_SCALE, BigDecimal.ROUND_HALF_UP));

bet &＃61; bet.add((b_9.multiply(new BigDecimal(Math.sin(lng.divide(b_3.multiply(b_3).multiply(b_5), B_DIV_SCALE, BigDecimal.ROUND_HALF_UP).multiply(pi).doubleValue())))

.add(b_9.multiply(b_3).multiply(new BigDecimal(Math.sin(lng.divide(b_5.multiply(b_2), B_DIV_SCALE, BigDecimal.ROUND_HALF_UP).multiply(pi).doubleValue())))).multiply(b_3).divide(b_5, B_DIV_SCALE, BigDecimal.ROUND_HALF_UP)));

return bet;

}

/**

* 判断是否在国内&＃xff0c;不在国内不做偏移

*

* &＃64;param lng

* &＃64;param lat

* &＃64;return

*/

public static boolean out_of_china(BigDecimal lng, BigDecimal lat) {

if (lng.compareTo(new BigDecimal("73.66")) <0 || lng.compareTo(new BigDecimal("135.05")) > 0) {

return true;

} else if (lat.compareTo(new BigDecimal("3.86")) <0 || lat.compareTo(new BigDecimal("53.55")) > 0) {

return true;

}

return false;

}