41I 2 3 6 7 8 本文件按照GB/T1.1—2020《標(biāo)準(zhǔn)化工作導(dǎo)則第1部分：標(biāo)準(zhǔn)化文件的結(jié)構(gòu)和起草規(guī)則》的規(guī)定請(qǐng)注意本文件的某些內(nèi)容可能涉及專利。本文件的發(fā)布機(jī)構(gòu)不承擔(dān)識(shí)別專利的責(zé)任。本文件由河南省自然資源標(biāo)準(zhǔn)化技術(shù)委員會(huì)（HN/TC1DB41/T2917—網(wǎng)格由兩組或多種曲線(面)集所包絡(luò)的空間區(qū)域。地質(zhì)空間剖分將地質(zhì)空間劃分為形狀近似、尺度連續(xù)、無縫測(cè)地線三棱臺(tái)基于測(cè)地線中分和地球半徑二分規(guī)則，將(橢)球體遞歸剖分為具有多分辨率、三棱臺(tái)特征的離散以符合時(shí)空大數(shù)據(jù)和云計(jì)算技術(shù)應(yīng)用和管理為目標(biāo)，按照三棱臺(tái)八叉樹網(wǎng)格體系的特征進(jìn)行剖分2DB41/T2917—2025在同一剖分系統(tǒng)內(nèi)，各網(wǎng)格位置編碼是唯一的，并采用統(tǒng)一的編碼模型。4.3穩(wěn)定性編碼模型是對(duì)地質(zhì)離散空間位置特征的邏輯抽象，采用一致的遞歸編碼機(jī)理，具有穩(wěn)定性。4.4擴(kuò)展性將八分體球面三角形每一條邊的中點(diǎn)用測(cè)地線連接，形成4個(gè)球面三角形，依此進(jìn)行遞歸剖分。剖分層次n=1,n=2時(shí)的情5.4.1徑向等長剖分。從球面網(wǎng)格的每個(gè)節(jié)點(diǎn)向球心引直線(球的半徑R),以R/2n將半徑等分，在各分層處以球面進(jìn)行剖分(以球心為原點(diǎn)，以K×R/2n為半徑畫球，(K為1,2,...,2"-1的整數(shù)序?yàn)槔钠史诌^程，如圖2所示。3a)n=1圖2三棱臺(tái)八叉樹網(wǎng)格徑向剖分效果(徑向等長)5.4.2徑向變長剖分。從球面網(wǎng)格的每個(gè)節(jié)點(diǎn)向球心引直線(球的半徑R)。以一定比例對(duì)半徑分割，RRRRRAR0圖3三棱臺(tái)八叉樹網(wǎng)格剖分效果(徑向變長)5.5三棱臺(tái)八叉樹網(wǎng)格特征三棱臺(tái)八叉樹網(wǎng)格(徑向等長剖分方式)特征及數(shù)量見附錄A。6三棱臺(tái)八叉樹網(wǎng)格編碼6.1編碼結(jié)構(gòu)與原則6.1.1編碼由卦限標(biāo)識(shí)碼、球面位置碼、徑向類型碼、徑向深度碼和拓展編碼構(gòu)成，拓展編碼可根據(jù)46.1.2用四個(gè)分隔符“_”將球體三棱臺(tái)八叉樹網(wǎng)格碼分為五個(gè)碼段，如圖4所示，其中卦限標(biāo)識(shí)碼1位，球面位置碼24位，徑向類型碼1位，徑向深度碼24位，拓展編碼根據(jù)應(yīng)用需求自行編碼不固定長度?！卣咕幋a——分隔符(英文下劃線)——徑向深度碼(二進(jìn)制)——分隔符(英文下劃線)——徑向類型碼(二進(jìn)制)——分隔符(英文下劃線)——球面位置碼(四進(jìn)制)——分隔符(英文下劃線)——卦限標(biāo)識(shí)碼(八進(jìn)制)圖4三棱臺(tái)八叉樹網(wǎng)格編碼規(guī)則51467“3”。編碼取值范圍為：0-333333333333333333333333,333333333333333333333333為球面剖分第2411023a)上三角形編碼b)下三角形編碼圖6固定方向編碼示意圖56.4徑向類型碼球體沿半徑方向剖分，采用徑向等長剖分方式時(shí)為0,采用徑向變長剖分方式時(shí)為1。6.5徑向深度碼6.5.1球體沿半徑方向剖分時(shí)，等長剖分方式和變長剖分方式，分別將徑向二分并遞歸，按照從地心到地表的方向，賦予順序編號(hào)。以等長為例，剖分3次后，P點(diǎn)的徑向深度編碼為110,如圖7所示。01PP十剖分2次剖分3次圖7徑向剖分示意圖(以等長剖分為例)6.5.2徑向深度碼的取值范圍為：0-111111111111111111111111。111111111111111111111111為地球半徑剖分第24層次時(shí)二進(jìn)制碼的最大值。6.6編碼示例與方法三棱臺(tái)八叉樹網(wǎng)格編碼示例見附錄B,編解碼方法見附錄C。6素的平均幾何特征與數(shù)量如表A.1所示。這里采用徑向等長剖分方式，示例的地球半徑取CGCS2000橢球長半軸、短半軸的平均值6367.444657km。1234567897三棱臺(tái)八叉樹網(wǎng)格編碼示例見表B.1，這里的編碼100球面與徑向等長方式同260面與徑向等長方式同時(shí)剖310面剖分13次，徑向等長方441剖分14次得到的網(wǎng)格（球511位于鄭州“鄭州二七罷在球面剖分24次、徑向變長方式剖分23次，開孔所610位于鄭州“鄭州二七罷在球面剖分24次、徑向等長方式剖分23次，終孔1000米深度（假設(shè)鉆孔完全按照地球半徑方向鉆8網(wǎng)格的編碼是指已知點(diǎn)的大地坐標(biāo)求其網(wǎng)格編碼，解碼是指已知網(wǎng)格編碼求其網(wǎng)格中心的大地坐根據(jù)八分體跨經(jīng)緯90°的剖分機(jī)理，卦限碼的計(jì)算公式見C.1：λ——為經(jīng)度（0°~360°,注意與東西經(jīng)的轉(zhuǎn)換）；φ——為緯度，北緯為正，南緯為負(fù)。如：點(diǎn)（E113.6775°,N34.7500°)的卦限碼為1。三棱臺(tái)八叉樹剖分的球面三角形的位置碼表示為：=123??,1到是個(gè)四分碼，根里以“行列逼近法”（趙學(xué)勝等，2003，QTM地址碼與經(jīng)緯度坐標(biāo)的快速轉(zhuǎn)換算法）為例說C.3.2.1根據(jù)三角形的層次求出其最大的行列數(shù)(,)，見公式C.3。C.3.2.2確定經(jīng)緯度坐標(biāo)(,)在網(wǎng)格地址碼的行列數(shù)(,)，見公式C.4。C.3.2.3根據(jù)所求位置在網(wǎng)格的行列數(shù)(,)和網(wǎng)格最大列數(shù)9b)如果小于或等于/2，分兩種情況：c)重復(fù)a）、b）步驟，求下一個(gè)地址示例：設(shè)剖分層次為24層，點(diǎn)（E113.6775°,N34.7500°)的球面位置碼為：020230212003022223131131。C.3.3.1根據(jù)三角形的層次求出其最大的行列數(shù)(,)C.3.3.2根據(jù)網(wǎng)格的最大行列數(shù)(,),把地址碼根據(jù)固定的方向按次序進(jìn)行逐層次地在格網(wǎng)的行數(shù)（注意：由于地址碼為“1”時(shí)，三角形可為頂或底，這對(duì)計(jì)算行數(shù)有很大影響。所以，d)重復(fù)步驟a）、b）、c直到最后求出三角形C.3.3.3根據(jù)行數(shù)和格網(wǎng)最大行列數(shù)(,)，把地址碼按次序進(jìn)行逐層次地解算，求出位置在格網(wǎng)的e)重復(fù)步驟（1234直到最后求出三角形所在的列數(shù)。C.3.3.4根據(jù)網(wǎng)格地址碼的行列數(shù)(,)確定三角形中點(diǎn)的經(jīng)緯度坐標(biāo)(,)，見公式C.6。分?jǐn)?shù)，取其分子的整數(shù)部分加上1（若為整數(shù)則不加1），然后再將其轉(zhuǎn)化為位的二進(jìn)制深度為00101000001101000101101。變長剖分的原理與之相同，但剖分的比例