java中WGS84坐标（ios）转换BD09坐标（百度坐标）

public class CoordinateConversion {static double x_PI &＃61; 3.14159265358979324 * 3000.0 / 180.0;static double PI &＃61; 3.1415926535897932384626;static double a &＃61; 6378245.0;static double ee &＃61; 0.00669342162296594323;/*** 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换* 即 百度 转 谷歌、高德* &＃64;param bd_lon* &＃64;param bd_lat* &＃64;returns {*[]}*/public String bd09togcj02(double bd_lon, double bd_lat){double x &＃61; bd_lon - 0.0065;double y &＃61; bd_lat - 0.006;double z &＃61; Math.sqrt(x * x &＃43; y * y) - 0.00002 * Math.sin(y * x_PI);double theta &＃61; Math.atan2(y, x) - 0.000003 * Math.cos(x * x_PI);double gg_lng &＃61; z * Math.cos(theta);double gg_lat &＃61; z * Math.sin(theta);// Point point&＃61;new Point(gg_lng, gg_lat);// return point;return gg_lng&＃43;","&＃43;gg_lat;}/*** 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换* 即谷歌、高德 转 百度* &＃64;param lng* &＃64;param lat* &＃64;returns {*[]}*/public static String gcj02tobd09(double lng, double lat){double z &＃61; Math.sqrt(lng * lng &＃43; lat * lat) &＃43; 0.00002 * Math.sin(lat * x_PI);double theta &＃61; Math.atan2(lat, lng) &＃43; 0.000003 * Math.cos(lng * x_PI);double bd_lng &＃61; z * Math.cos(theta) &＃43; 0.0065;double bd_lat &＃61; z * Math.sin(theta) &＃43; 0.006;//Point point&＃61;new Point(bd_lng, bd_lat);// return point;return bd_lng&＃43;","&＃43;bd_lat;};/*** WGS84转GCj02* &＃64;param lng* &＃64;param lat* &＃64;returns {*[]}*/public static String wgs84togcj02(double lng, double lat){double dlat &＃61; transformlat(lng - 105.0, lat - 35.0);double dlng &＃61; transformlng(lng - 105.0, lat - 35.0);double radlat &＃61; lat / 180.0 * PI;double magic &＃61; Math.sin(radlat);magic &＃61; 1 - ee * magic * magic;double sqrtmagic &＃61; Math.sqrt(magic);dlat &＃61; (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * PI);dlng &＃61; (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * PI);double mglat &＃61; lat &＃43; dlat;double mglng &＃61; lng &＃43; dlng;//Point point&＃61;new Point(mglng, mglat);// return point;return mglng&＃43;","&＃43;mglat;};/*** GCJ02 转换为 WGS84* &＃64;param lng* &＃64;param lat* &＃64;returns {*[]}*/public String gcj02towgs84(double lng, double lat){double dlat &＃61; transformlat(lng - 105.0, lat - 35.0);double dlng &＃61; transformlng(lng - 105.0, lat - 35.0);double radlat &＃61; lat / 180.0 * PI;double magic &＃61; Math.sin(radlat);magic &＃61; 1 - ee * magic * magic;double sqrtmagic &＃61; Math.sqrt(magic);dlat &＃61; (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * PI);dlng &＃61; (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * PI);double mglat &＃61; lat &＃43; dlat;double mglng &＃61; lng &＃43; dlng;// Point point&＃61;new Point(mglng, mglat);// return point;return mglat&＃43;","&＃43;mglng;};/*** WGS84 转换为 BD-09* &＃64;param lng* &＃64;param lat* &＃64;returns {*[]}* */private static String wgs84tobd09(double lng, double lat){//第一次转换double dlat &＃61; transformlat(lng - 105.0, lat - 35.0);double dlng &＃61; transformlng(lng - 105.0, lat - 35.0);double radlat &＃61; lat / 180.0 * PI;double magic &＃61; Math.sin(radlat);magic &＃61; 1 - ee * magic * magic;double sqrtmagic &＃61; Math.sqrt(magic);dlat &＃61; (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * PI);dlng &＃61; (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * PI);double mglat &＃61; lat &＃43; dlat;double mglng &＃61; lng &＃43; dlng;//第二次转换double z &＃61; Math.sqrt(mglng * mglng &＃43; mglat * mglat) &＃43; 0.00002 * Math.sin(mglat * x_PI);double theta &＃61; Math.atan2(mglat, mglng) &＃43; 0.000003 * Math.cos(mglng * x_PI);double bd_lng &＃61; z * Math.cos(theta) &＃43; 0.0065;double bd_lat &＃61; z * Math.sin(theta) &＃43; 0.006;return bd_lng&＃43;","&＃43;bd_lat;}private static double transformlat(double lng,double lat){double ret&＃61; -100.0 &＃43; 2.0 * lng &＃43; 3.0 * lat &＃43; 0.2 * lat * lat &＃43; 0.1 * lng * lat &＃43; 0.2 * Math.sqrt(Math.abs(lng));ret &＃43;&＃61; (20.0 * Math.sin(6.0 * lng * PI) &＃43; 20.0 * Math.sin(2.0 * lng * PI)) * 2.0 / 3.0;ret &＃43;&＃61; (20.0 * Math.sin(lat * PI) &＃43; 40.0 * Math.sin(lat / 3.0 * PI)) * 2.0 / 3.0;ret &＃43;&＃61; (160.0 * Math.sin(lat / 12.0 * PI) &＃43; 320 * Math.sin(lat * PI / 30.0)) * 2.0 / 3.0;return ret;}private static double transformlng(double lng,double lat){double ret &＃61; 300.0 &＃43; lng &＃43; 2.0 * lat &＃43; 0.1 * lng * lng &＃43; 0.1 * lng * lat &＃43; 0.1 * Math.sqrt(Math.abs(lng));ret &＃43;&＃61; (20.0 * Math.sin(6.0 * lng * PI) &＃43; 20.0 * Math.sin(2.0 * lng * PI)) * 2.0 / 3.0;ret &＃43;&＃61; (20.0 * Math.sin(lng * PI) &＃43; 40.0 * Math.sin(lng / 3.0 * PI)) * 2.0 / 3.0;ret &＃43;&＃61; (150.0 * Math.sin(lng / 12.0 * PI) &＃43; 300.0 * Math.sin(lng / 30.0 * PI)) * 2.0 / 3.0;return ret;}public static void main(String[] args) {// 两次谷歌转换为百度坐标// 第一次 WGS84转GCj02String lnglat&＃61;wgs84togcj02(117.20296517261839,31.841652709281103); double lng&＃61;Double.parseDouble(lnglat.split(",")[0]);double lat&＃61;Double.parseDouble(lnglat.split(",")[1]);System.out.println("第一次转换的结果:"&＃43;lng&＃43;","&＃43;lat);// 第二次 gcj02tobd09System.out.println("第二次转换的结果:"&＃43;gcj02tobd09(lng,lat));// 谷歌转百度一次转换System.out.println("谷歌转换为百度一次转换的结果:"&＃43;wgs84tobd09(117.20296517261839,31.841652709281103));}
}