DB41T2917-2025地质空间三棱台八叉树网格剖分编码规范

A DB41/T2917— 发DB41/TDB41/T2917—DB41/TDB41/T2917— 前 附录A（资料性）三棱台八叉树网格径向等长剖分几何特征与数 附录B（资料性）三棱台八叉树网格编码示 附录C（资料性）三棱台八叉树网格的编解码方 本文件按照GB/T1.1—2020《标准化工作导则第1部分：标准化文件的结构和起草规则》的规定DB41/TDB41/T2917—DB41/TDB41/T2917—[来源：GB/T39409— 分层次𝑛=1，𝑛=2时的情况如图1所示。a)𝑛= b)𝑛=图1 各分层处以球面进行剖分（以球心为原点，以𝐾𝑅⁄2𝑛为半径画球，(𝐾1，2,...，2n-12 𝑛= b)𝑛=图2三棱台八叉树网格径向剖分效果（径向等长3 3 3

3（√2𝑛𝑅√2𝑛𝑅√2𝑛𝑅，√

=𝑛=23𝑛= b)𝑛=图3三棱台八叉树网格剖分效果（径向变长 图4 图5 上三角形编 b)下三角形编图6 3，P1107图7径向剖分示意图（以等长剖分为例 附录长半轴、短半轴的平均值6367.444657km。表A.1三棱台上表面边长三棱台上表面面积三棱台体积附录B.1分5（棱台式剖分5（面三棱台剖分14次得到的网格（面三棱台（E113.6775°，（E113.6775°，长方式剖分23次，终孔全按照地球半径方向钻进）附录网格的编码是指已知点的大地坐标求其网格编码，解码是指已知网格编码求其网格中心的大地坐𝑄=λ\90°−4×((φ− \λ——为经度（0°～360°，注意与东西经的转换φ——为纬度，北纬为正，南纬为负。如：点（E113.6775°，N34.7500°）的卦限码为1𝜆=45°+𝑄×90°（𝑄≤𝜆=45°+(𝑄−4)×90°（𝑄≥ 三棱台八叉树剖分的球面三角形𝐴的位置码表示为：𝐴=𝑎1𝑎2𝑎3𝑎𝑘𝑎1到𝑎𝑘是𝑘个四分码，根根据三角形的层次𝐾求出其最大的行列数(𝐼,𝐽)C.3I=2𝑘,𝐽=2×𝐼− 确定经纬度坐标(𝜆,𝜑)在网格地址码的行列数(𝑖,𝑗)C.4𝜑×𝑖=𝑖𝑛𝑡

]+𝑗=2×𝑖𝑛𝑡[𝜆×(𝐼−𝑖+1)]+ 根据所求位置在网格的行列数(𝑖,𝑗)和网格最大列数(𝐼,𝐽)如果𝑖大于𝑛/2，地址码𝑎𝑘为“1”，同时𝑛=𝑛/2𝐼=(𝑖−如果𝑗小于2𝑛𝑖+1],地址码

为“2”，同时𝑛= 否则，若𝑗大于𝑛，则地址码𝑎𝑘为“3”，同时𝑛=𝑛/2,𝑗=(𝑗2𝑛)。若𝑗小于或等于𝑛

为“0”，同时𝑗=(𝑗𝑛+2𝑖−1)，𝐼=(𝑛𝑖+1)𝑛=𝑛/2设剖分层次为层，点根据三角形的层次𝐾求出其最大的行列数(𝐼,I=2𝑘,𝐽=2×𝐼− 根据网格的最大行列数(𝐼,𝐽),把地址码根据固定的方向按次序进行逐层次地解算，求出位置在层次解算中，首先要设置一个标志位𝑃“1”为顶三角形时，𝑃1“1”为底三角形时，𝑃

=“1”：若𝑃=1，则𝑖=𝑖+𝑛；若𝑃=−1,则𝑛==“2”或“3”：若𝑃=−1，则𝑖=𝑖+𝑛；若𝑃=1,则𝑛==“0”：若𝑃=−1，则𝑖=𝑖+𝑛；若𝑃=1,则𝑛=𝑛/2，且𝑝=−根据行数𝑖和格网最大行列数(𝐼,𝐽)，把地址码按次序进行逐层次地解算，求出位置在格网的如果地址码𝑎𝑘=“1”：则𝑛=𝑛/2，𝑖=𝑖−如果地址码𝑎𝑘=“2”，则𝑛=如果地址码𝑎𝑘=“3”：则𝑗=𝑗+𝑛，𝑛=如果地址码𝑎𝑘=“0”：则𝑛=𝑛/2，𝑖=𝑛−𝑖+1𝑗=𝑗+(2(𝑛𝑖根据网格地址码的行列数(𝑖𝑗)确定三角形中点的经纬度坐标(𝜆,𝜑)，C.6𝜆=90°× 𝜑= ×(𝑖

式即可。上述变长径