经纬度坐标转地图金字塔所在块

地图金字塔从0层开始

#define LON_SPAN 360.0
// 开始经度(最左端)
#define LAT_SPAN 180.0
#define GLOBAL_LEFT -180.0
// 开始纬度(最上端)
#define GLOBAL_TOP 90.0
#define GLOBAL_RIGHT 180.0
#define GLOBAL_BOTTOM -90.0
// 地球的纬度跨度(180-(-180)):360
#define LON_SPAN (GLOBAL_RIGHT - GLOBAL_LEFT)
// 地球的纬度跨度(90-(-90)):180
#define LAT_SPAN (GLOBAL_TOP - GLOBAL_BOTTOM)

	/// <summary>
	/// 根据经纬度计算金字塔对应层的所在块的行列号
	/// </summary>
	/// <param name="lng">经度</param>
	/// <param name="lat">纬度</param>
	/// <param name="storey">层数(start from 0)</param>
	/// <param name="x">块所在的列</param>
	/// <param name="y">块所在的行</param>
	static void lonlat_to_pyramid_block(float lng, float lat, int storey, int& x, int& y)
	{
			// 计算第storey的行列(2的storey次方)
	    int n = pow(2, storey);
	    // 
    	float d_x = LON_SPAN / n; // 每块经度跨度
    	float d_y = LAT_SPAN / n; // 每块纬度跨度
    	x = (lng + LAT_SPAN) / d_x;
    	y = (GLOBAL_TOP - lat) / d_y;
	}