本文糅合
https://www.cnblogs.com/xiaozhi_5638/p/4748186.html - [只提供了C语言版本]
https://www.jianshu.com/p/063bee79507a - [无法处理负值经纬度]
贴出改版墨卡托转换代码
1 public class BDMercatorToLonLatAll { 2 //以下是根据百度地图JavaScript API破解得到 百度坐标<->墨卡托坐标 转换算法 3 private static double[] array1 = {75, 60, 45, 30, 15, 0}; 4 private static double[] array3 = {12890594.86, 8362377.87, 5591021, 3481989.83, 1678043.12, 0}; 5 private static double[][] array2 = {new double[]{-0.0015702102444, 111320.7020616939, 1704480524535203D, -10338987376042340D, 26112667856603880D, -35149669176653700D, 26595700718403920D, -10725012454188240D, 1800819912950474D, 82.5} 6 , new double[]{0.0008277824516172526, 111320.7020463578, 647795574.6671607, -4082003173.641316, 10774905663.51142, -15171875531.51559, 12053065338.62167, -5124939663.577472, 913311935.9512032, 67.5} 7 , new double[]{0.00337398766765, 111320.7020202162, 4481351.045890365, -23393751.19931662, 79682215.47186455, -115964993.2797253, 97236711.15602145, -43661946.33752821, 8477230.501135234, 52.5} 8 , new double[]{0.00220636496208, 111320.7020209128, 51751.86112841131, 3796837.749470245, 992013.7397791013, -1221952.21711287, 1340652.697009075, -620943.6990984312, 144416.9293806241, 37.5} 9 , new double[]{-0.0003441963504368392, 111320.7020576856, 278.2353980772752, 2485758.690035394, 6070.750963243378, 54821.18345352118, 9540.606633304236, -2710.55326746645, 1405.483844121726, 22.5} 10 , new double[]{-0.0003218135878613132, 111320.7020701615, 0.00369383431289, 823725.6402795718, 0.46104986909093, 2351.343141331292, 1.58060784298199, 8.77738589078284, 0.37238884252424, 7.45}}; 11 private static double[][] array4 = {new double[]{1.410526172116255e-8, 0.00000898305509648872, -1.9939833816331, 200.9824383106796, -187.2403703815547, 91.6087516669843, -23.38765649603339, 2.57121317296198, -0.03801003308653, 17337981.2} 12 , new double[]{-7.435856389565537e-9, 0.000008983055097726239, -0.78625201886289, 96.32687599759846, -1.85204757529826, -59.36935905485877, 47.40033549296737, -16.50741931063887, 2.28786674699375, 10260144.86} 13 , new double[]{-3.030883460898826e-8, 0.00000898305509983578, 0.30071316287616, 59.74293618442277, 7.357984074871, -25.38371002664745, 13.45380521110908, -3.29883767235584, 0.32710905363475, 6856817.37} 14 , new double[]{-1.981981304930552e-8, 0.000008983055099779535, 0.03278182852591, 40.31678527705744, 0.65659298677277, -4.44255534477492, 0.85341911805263, 0.12923347998204, -0.04625736007561, 4482777.06} 15 , new double[]{3.09191371068437e-9, 0.000008983055096812155, 0.00006995724062, 23.10934304144901, -0.00023663490511, -0.6321817810242, -0.00663494467273, 0.03430082397953, -0.00466043876332, 2555164.4} 16 , new double[]{2.890871144776878e-9, 0.000008983055095805407, -3.068298e-8, 7.47137025468032, -0.00000353937994, -0.02145144861037, -0.00001234426596, 0.00010322952773, -0.00000323890364, 826088.5}}; 17 18 public static Map<String, Double> convertLonLatToMC (Double lon, Double lat) {//百度坐标转墨卡托 19 double[] arr = null; 20 double n_lat = lat > 74 ? 74 : lat; 21 n_lat = n_lat < -74 ? -74 : n_lat; 22 for (int i = 0; i < array1.length; i++) { 23 if (lat >= array1[i]) { 24 arr = array2[i]; 25 break; 26 } 27 } 28 if (arr == null) { 29 for (int i = array1.length - 1; i >= 0; i--) { 30 if (lat <= -array1[i]) { 31 arr = array2[i]; 32 break; 33 } 34 } 35 } 36 return converter(lon, lat, arr); 37 } 38 39 public static Map<String, Double> convertMCToLonLat (Double x, Double y) {//墨卡托坐标转百度 - Mercator 40 double[] arr = null; 41 for (int i = 0; i < array3.length; i++) { 42 if (Math.abs(y) >= array3[i]) { 43 arr = array4[i]; 44 break; 45 } 46 } 47 Map<String,Double> location = converter(Math.abs(x), Math.abs(y), arr); 48 location.put("lng",location.get("x")); 49 location.remove("x"); 50 location.put("lat",location.get("y")); 51 location.remove("y"); 52 return location; 53 } 54 55 private static Map<String, Double> converter(double x, double y, double[] param) { 56 Double T = param[0] + param[1] * Math.abs(x); 57 Double cC = Math.abs(y) / param[9]; 58 Double cF = param[2] + param[3] * cC + param[4] * cC * cC + param[5] * cC * cC * cC + param[6] * cC * cC * cC * cC + param[7] * cC * cC * cC * cC * cC + param[8] * cC * cC * cC * cC * cC * cC; 59 T *= (x < 0 ? -1 : 1); 60 cF *= (y < 0 ? -1 : 1); 61 Map<String, Double> location = new HashMap<String, Double>(); 62 location.put("x", T); 63 location.put("y", cF); 64 return location; 65 } 66 }
后期传参过程中，将-值改为M即可
至于说 https://www.cnblogs.com/xiaozhi_5638/p/4748186.html 第三段话
测试代码如下
1 double lat = -34.7715d; 2 double lon = -113.7276d; 3 4 System.out.println("经纬度坐标x--->"+lon); 5 System.out.println("经纬度坐标y--->"+lat); 6 7 Map<String, Double> stringDoubleMap = BDMercatorToLonLatAll.convertLonLatToMC(lon, lat); 8 //System.out.println("墨卡托坐标--->"+stringDoubleMap); 9 10 Double xMC = stringDoubleMap.get("x"); 11 Double yMC = stringDoubleMap.get("y"); 12 System.out.println("墨卡托坐标x--->"+xMC); 13 System.out.println("墨卡托坐标y--->"+yMC); 14 15 //根据百度经纬度计算出来墨卡托坐标后，将结果除以地图分辨率Math.Pow(2,18-zoom)即可得到平面像素坐标，然后将像素坐标除以256分别得到瓦片的行列号。 16 17 int zoom = 13; 18 System.out.println("地图层级--->"+zoom); 19 20 double dpi = Math.pow(2, 18 - zoom); 21 System.out.println("地图分辨率--->"+dpi); 22 23 double xPlanePixelCoordinates = xMC / dpi; 24 double yPlanePixelCoordinates = yMC / dpi; 25 System.out.println("平面像素坐标x--->"+xPlanePixelCoordinates); 26 System.out.println("平面像素坐标y--->"+yPlanePixelCoordinates); 27 28 int tileImgDpi = 256; 29 System.out.println("瓦片大小--->"+tileImgDpi); 30 31 double xTile = xPlanePixelCoordinates / tileImgDpi; 32 double yTile = yPlanePixelCoordinates / tileImgDpi; 33 //System.out.println("瓦片行列号x--->"+xTile); 34 //System.out.println("瓦片行列号y--->"+yTile); 35 36 int xTileI = (int) xTile; 37 int yTileI = (int) yTile; 38 System.out.println("瓦片行列号x--->"+xTileI); 39 System.out.println("瓦片行列号y--->"+yTileI);
测试通过。