北斗网格位置码在低空空域管理的Cesium实现
北斗网格码作为前沿主流技术，在低空空域的数字化管理中发挥着重要作用。
北斗网格位置码相关描述见国标：GB/T 39409-2020
（1）特定区域、特定级别的网格显示
//特定区域如果跨东西半球、南北半球，分部分处理，遍历循环时比较方便
//特定区域：中国区域版图
const lon_west=72.0 //矩形左下角经度
const lon_east=138.0 //矩形左上角经度
const lat_south=16.0 //矩形左下角纬度
const lat_north=56.0 //矩形左上角纬度
分别按剖分规则遍历经纬度、高度，绘制立方体
//标准划分
const level_1 = 4; // 一级 4° 纬度
const level_1_0 = 6; // 一级 6° 经度
const level_2 = 30 /60; // 二级 30′
const level_3 = 15 /60; // 三级 15′ //经度上15/60，
const level_3_0 = 10 /60; //纬度上应该为10/60
const level_4 = 1 /60; // 四级 1′
const level_5 = 4 /3600; // 五级 4″
const level_6 = 2 /3600; // 六级 2″
const level_7 = 1 /(3600*4); // 七级 1/4″
const level_8 = 1 /(3600*32); // 八级 1/32″
const level_9 = 1 /(3600*256); // 九级 1/256″
const level_10 = 1 /(3600*2048); // 十级 1/2048″ 1.5cm
一级剖分：
五级剖分：
（2）北斗二维编码
按照国标要求编码，注意编码是不等长的，级别不同，长度不同。将编码绘制到网格的中心。
二级的二维编码如下：
（3）北斗三维编码
高度维编码获取方式
// 根据公式 C.14 逐个提取每个部分的编码
const a11 = parseInt(binaryStr.slice(32 - 3, 32), 2).toString(8); // --从低位（右侧）取第 1 位到第 3 位
const a10 = parseInt(binaryStr.slice(32-6, 32-3), 2).toString(8);// --从低位（右侧）取第 4 位到第 6 位
const a9 = parseInt(binaryStr.slice(32-9, 32-6), 2).toString(8);// --从低位（右侧）取第 7 位到第 9 位
const a8 = parseInt(binaryStr.slice(32-12, 32-9), 2).toString(8);// --从低位（右侧）取第 10 位到第 12 位
const a7 = parseInt(binaryStr[32 - 13], 2).toString(2);// --从低位（右侧）取第 13 位
const a6 = parseInt(binaryStr.slice(32-17, 32-13), 2).toString(16).toUpperCase();// --从低位（右侧）取第 14 位到第 17 位
const a5 = parseInt(binaryStr.slice(32-21, 32-17), 2).toString(16).toUpperCase();// --从低位（右侧）取第 18 位到第 21 位
const a4 = parseInt(binaryStr[32 - 22], 2).toString(2);// --从低位（右侧）取第 22 位
const a3 = parseInt(binaryStr.slice(32-25, 32-22), 2).toString(8); // --从低位（右侧）取第 23 位到第 25 位
const a21 = parseInt(binaryStr.slice(32-31, 32-25), 2).toString(10).padStart(2, '0');// --从低位（右侧）取第 26 位到第 31 位
const a0 = parseInt(binaryStr[32 - 32], 2).toString(2);// --从低位（右侧）取第 32 位
二级网格的三维编码如下：
（4）用户可视范围、用户地图比例尺动态确定区域及级别
视频见文末链接
（5）根据经纬度点、高度及级别绘制点所在网格
(6)动态隔离飞行器所在北斗网格空域
更多内容见外链：https://zhuanlan.zhihu.com/p/8428345687
