层行列和经纬度坐标之间的相互转化方法（谷歌地图）

太过爱你忘了你带给我的痛 2021-09-29 01:28 448阅读 0赞

## **微信搜索：“二十同学” 公众号，欢迎关注一条不一样的成长之路** ##

谷歌地图层行列是以**左上角为原点**，瓦片编号是向右和向下，没有负数。

谷歌是从左上角经纬度(-180,90)度开始计算瓦片，在第0级时谷歌将世界地图分为1块，第1级的时候分为4块。

**层zoom(z)，列col(x)，行row(y) **

![å¾1-4][1-4]

如何在已知层行列情况计算出该瓦片对应的经纬度坐标，和已知经纬度情况下计算该点在不同层级下的瓦片位置，以左上角为原点。

/**
         * 谷歌下转换经纬度对应的层行列
         *
         * @param lon  经度
         * @param lat  维度
         * @param zoom 在第zoom层进行转换
         * @return
         */
        public static int[] GoogleLonLatToXYZ(double lon, double lat, int zoom) {
    
            double n = Math.pow(2, zoom);
            double tileX = ((lon + 180) / 360) * n;
            double tileY = (1 - (Math.log(Math.tan(Math.toRadians(lat)) + (1 / Math.cos(Math.toRadians(lat)))) / Math.PI)) / 2 * n;
    
            int[] xy = new int[2];
    
            xy[0] = (int) Math.floor(tileX);
            xy[1] = (int) Math.floor(tileY);
    
            return xy;
        }
    
        /**
         * 层行列转经纬度
         *
         * @param x
         * @param y
         * @param z
         * @return
         */
        public static double[] XYZtoLonlat(int z, int x, int y) {
    
            double n = Math.pow(2, z);
            double lon = x / n * 360.0 - 180.0;
            double lat = Math.atan(Math.sinh(Math.PI * (1 - 2 * y / n)));
            lat = lat * 180.0 / Math.PI;
            double[] lonlat = new double[2];
            lonlat[0] = lon;
            lonlat[1] = lat;
            return lonlat;
        }

如果想已**左下角为原点**的xyz,只需要将y轴反转即可:

int ext = (int) Math.pow(2, z);
    int changeRow = ext - row - 1;

以北京天安门的地理坐标点(116.407945，39.91349)计算出谷歌各层瓦片编号。

[1-4]: /images/20210920/914ccee6fde247048521584e5f4fbe46.png