层行列和经纬度坐标之间的相互转化方法

公式

Lon./lat. to tile numbers

n = 2 ^ zoom
xtile = n * ((lon_deg + 180) / 360)
ytile = n * (1 - (log(tan(lat_rad) + sec(lat_rad)) / π)) / 2

sec函数和cos函数互为倒数关系，即sec=1/cos，cos=1/sec。

Tile numbers to lon./lat.

n = 2 ^ zoom
lon_deg = xtile / n * 360.0 - 180.0
lat_rad = arctan(sinh(π * (1 - 2 * ytile / n)))
lat_deg = lat_rad * 180.0 / π

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

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

int long2tilex(double lon, int z) 
{ 
	return (int)(floor((lon + 180.0) / 360.0 * (1 << z))); 
}

int lat2tiley(double lat, int z)
{ 
    double latrad = lat * M_PI/180.0;
	return (int)(floor((1.0 - asinh(tan(latrad)) / M_PI) / 2.0 * (1 << z))); 
}

double tilex2long(int x, int z) 
{
	return x / (double)(1 << z) * 360.0 - 180;
}

double tiley2lat(int y, int z) 
{
	double n = M_PI - 2.0 * M_PI * y / (double)(1 << z);
	return 180.0 / M_PI * atan(0.5 * (exp(n) - exp(-n)));
}

 

缩放等级

zoom参数是介于0（缩小）和18（放大）之间的整数。18通常是最大值，但有些磁贴服务器可能会超过这个值。

（*）虽然以度为单位的宽度（经度）是恒定的，但给定缩放级别，对于所有平铺，高度不会发生这种情况。一般来说，属于同一行的瓦片在度数上具有相等的高度，但从赤道移动到极点的高度会减少。