分享一个java处理坐标转换问题的工具类【天地图、 高德地图、百度地图】

/**
 * @FileName: CoordinatesTransformControllerUtil
 * @Description: 坐标转换工具类 天地图使用的是地球坐标系(WGS84); 高德地图使用的是火星坐标系(GCJ02);
 *               百度地图使用的是百度坐标(bd09II)
 */
public class CoordinateTransformUtil {

	private static final double X_PI = 3.14159265358979324 * 3000.0 / 180.0;
	private static final double PI = 3.14159265358979324;
	// 卫星椭球坐标投影到平面地图坐标系的投影因子。 地球长半径
	private static final double EARTH_MAJOR_RADIUS = 6378245.0;
	// 椭球的偏心率。
	private static final double ECCENTRICITY_RATIO = 0.00669342162296594323;

	/**
	 * @Discription: 百度坐标(bd09II)转火星坐标(GCJ02)
	 * @Param lat 百度坐标纬度
	 * @Param lng 百度坐标经度
	 * @Created on: 2020/12/23 15:40
	 * @author muyuanpei
	 */
	public static Map<String, Double> baiduTomars(double lat, double lng) {
		Map<String, Double> mars_point = new HashMap<>();
		mars_point.put("lng", 0D);
		mars_point.put("lat", 0D);

		double x = lng - 0.0065;
		double y = lat - 0.006;
		double z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * X_PI);
		double theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * X_PI);
		mars_point.put("lng", z * Math.cos(theta));
		mars_point.put("lat", z * Math.sin(theta));
		return mars_point;
	}

	/**
	 * @Discription: 火星坐标系(GCJ02)转百度坐标系(bd09II)
	 * @Param gcjLat 火星坐标纬度
	 * @Param gcjlng 火星坐标经度
	 * @Created on: 2020/12/23 15:42
	 * @author muyuanpei
	 */
	public static Map<String, Double> marsTobaidu(double gcjLat, double gcjlng) {
		HashMap<String, Double> baidu_point = new HashMap<>();
		baidu_point.put("lng", 0D);
		baidu_point.put("lat", 0D);

		double x = gcjlng;
		double y = gcjLat;
		double z = Math.sqrt(x * x + y * y) + 0.00002 * Math.sin(y * X_PI);
		double theta = Math.atan2(y, x) + 0.000003 * Math.cos(x * X_PI);
		baidu_point.put("lng", z * Math.cos(theta) + 0.0065);
		baidu_point.put("lat", z * Math.sin(theta) + 0.006);
		return baidu_point;
	}

	/**
	 * @Discription: 火星坐标系(GCJ02)转地球坐标系(WGS84)
	 * @Param lat 火星坐标纬度
	 * @Param lng 火星坐标经度
	 * @Created on: 2020/12/23 15:42
	 * @author muyuanpei
	 */
	public static Map<String, Double> transformGCJ2WGS(double gcjLat, double gcjlng) {
		Map<String, Double> map = delta(gcjLat, gcjlng);
		double lat = map.get("lat");
		double lng = map.get("lng");
		map.put("lat", gcjLat - lat);
		map.put("lng", gcjlng - lng);
		return map;
	}

	private static Map delta(double lat, double lng) {
		double dLat = transformLat(lng - 105.0, lat - 35.0, PI);
		double dlng = transformlng(lng - 105.0, lat - 35.0, PI);
		double radLat = lat / 180.0 * PI;
		double magic = Math.sin(radLat);
		magic = 1 - ECCENTRICITY_RATIO * magic * magic;
		double sqrtMagic = Math.sqrt(magic);
		dLat = (dLat * 180.0) / ((EARTH_MAJOR_RADIUS * (1 - ECCENTRICITY_RATIO)) / (magic * sqrtMagic) * PI);
		dlng = (dlng * 180.0) / (EARTH_MAJOR_RADIUS / sqrtMagic * Math.cos(radLat) * PI);

		Map<String, Double> map = new HashMap<>();
		map.put("lat", dLat);
		map.put("lng", dlng);
		return map;
	}

	/**
	 * @Discription: 地球坐标系(WGS - 84)转火星坐标系(GCJ)
	 * @Param wgLat 地球坐标纬度
	 * @Param wglng 地球坐标经度
	 * @Created on: 2020/12/23 15:42
	 * @author muyuanpei
	 */
	public static Map<String, Double> transform(double wgLat, double wglng) {

		HashMap<String, Double> mars_point = new HashMap<>();
		mars_point.put("lng", 0D);
		mars_point.put("lat", 0D);
		// 如果是国外坐标点，则直接返回传进来的坐标点
		if (outOfChina(wglng, wgLat)) {
			mars_point.put("lng", wglng);
			mars_point.put("lat", wgLat);
			return mars_point;
		}
		// 如果为国内坐标点，则计算偏移量
		double dLat = transformLat(wglng - 105.0, wgLat - 35.0, PI);
		double dlng = transformlng(wglng - 105.0, wgLat - 35.0, PI);
		double radLat = wgLat / 180.0 * PI;
		double magic = Math.sin(radLat);
		magic = 1 - ECCENTRICITY_RATIO * magic * magic;
		double sqrtMagic = Math.sqrt(magic);
		dLat = (dLat * 180.0) / ((EARTH_MAJOR_RADIUS * (1 - ECCENTRICITY_RATIO)) / (magic * sqrtMagic) * PI);
		dlng = (dlng * 180.0) / (EARTH_MAJOR_RADIUS / sqrtMagic * Math.cos(radLat) * PI);
		mars_point.put("lng", wglng + dlng);
		mars_point.put("lat", wgLat + dLat);

		return mars_point;
	}

	// 公共方法
	/* 判断是否在国内，不在国内则不做偏移 */
	private static boolean outOfChina(double lng, double lat) {
		if ((lng < 72.004 || lng > 137.8347) && (lat < 0.8293 || lat > 55.8271)) {
			return true;
		} else {
			return false;
		}
	}

	private static double transformLat(double x, double y, double pi) {
		double ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * Math.sqrt(Math.abs(x));
		ret += (20.0 * Math.sin(6.0 * x * pi) + 20.0 * Math.sin(2.0 * x * pi)) * 2.0 / 3.0;
		ret += (20.0 * Math.sin(y * pi) + 40.0 * Math.sin(y / 3.0 * pi)) * 2.0 / 3.0;
		ret += (160.0 * Math.sin(y / 12.0 * pi) + 320 * Math.sin(y * pi / 30.0)) * 2.0 / 3.0;
		return ret;
	}

	private static double transformlng(double x, double y, double pi) {
		double ret = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1 * Math.sqrt(Math.abs(x));
		ret += (20.0 * Math.sin(6.0 * x * pi) + 20.0 * Math.sin(2.0 * x * pi)) * 2.0 / 3.0;
		ret += (20.0 * Math.sin(x * pi) + 40.0 * Math.sin(x / 3.0 * pi)) * 2.0 / 3.0;
		ret += (150.0 * Math.sin(x / 12.0 * pi) + 300.0 * Math.sin(x / 30.0 * pi)) * 2.0 / 3.0;
		return ret;
	}

	public static void main(String[] args) {

		// 测试坐标使用的都是"北京市西城区南礼士路66号建威大厦的坐标"

		Map<String, Double> map = new HashMap<>();

		double bdLat = 39.914366; // bd09II 纬度lat
		double bdlng = 116.359607; // bd09II 经度lng
		map = CoordinateTransformUtil.baiduTomars(bdLat, bdlng);

		System.out.println("火星坐标---->     " + map.get("lng") + "," + map.get("lat"));

		double gcjLat = 34.604665; // gcj02 纬度lat
		double gcjlng = 119.231148; // gcj02 经度lng
		map = CoordinateTransformUtil.marsTobaidu(gcjLat, gcjlng);

		System.out.println("百度坐标---->    " + map.get("lng") + "," + map.get("lat"));

		map = CoordinateTransformUtil.transformGCJ2WGS(gcjLat, gcjlng);

		System.out.println("地球坐标---->     " + map.get("lng") + "," + map.get("lat"));

		double wgLat = 39.9073; // wgs84 纬度lat
		double wglng = 116.3470; // wgs84 经度lng
		map = CoordinateTransformUtil.transform(wgLat, wglng);

		System.out.println("火星坐标---->     " + map.get("lng") + "," + map.get("lat"));
	}

}