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

地图金字塔从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)

	/// 

	/// 根据经纬度计算金字塔对应层的所在块的行列号
	/// 
	/// 经度
	/// 纬度
	/// 层数(start from 0)
	/// 块所在的列
	/// 块所在的行
	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;
	}