北斗网格码作为前沿主流技术,在低空空域的数字化管理中发挥着重要作用。 北斗网格位置码相关描述见国标: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)根据经纬度点、高度及级别绘制点所在网格 外链:https://zhuanlan.zhihu.com/p/8428345687
标签:const,level,--,32,网格,空域,toString,binaryStr,Cesium From: https://www.cnblogs.com/bimgis/p/18589121