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/TC1地質(zhì)空間三棱臺(tái)八叉樹(shù)網(wǎng)格剖分編碼規(guī)范3.13.23.33.43.53.6以符合時(shí)空大數(shù)據(jù)和云計(jì)算技術(shù)應(yīng)用和管理為目標(biāo)，按照三棱臺(tái)八叉樹(shù)網(wǎng)格體系的特征進(jìn)行剖分和編碼，實(shí)現(xiàn)本編碼到大地坐標(biāo)系相互轉(zhuǎn)換，面向2在同一剖分系統(tǒng)內(nèi)，各網(wǎng)格位置編碼是唯一的，并采用統(tǒng)一的編碼模4.3穩(wěn)定性編碼模型是對(duì)地質(zhì)離散空間位置特征的邏輯抽象4,4擴(kuò)展性5三棱臺(tái)八叉樹(shù)網(wǎng)格剖分5.2八分體剖分將八分體球面三角形每一條邊的中點(diǎn)用測(cè)地線連接，形成4個(gè)球面三角形，依此進(jìn)行遞歸剖分。剖5.4.1徑向等長(zhǎng)剖分。從球面網(wǎng)格的每個(gè)節(jié)點(diǎn)向球心引直線（球的半徑R），以R/2n將半徑等分，在3三棱臺(tái)八叉樹(shù)網(wǎng)格（徑向等長(zhǎng)剖分方式）特征及數(shù)量見(jiàn)附6.1.1編碼由卦限標(biāo)識(shí)碼、球面位置碼、徑向類型碼、徑向深度碼和拓展編碼構(gòu)成，拓展編碼可根據(jù)46.2卦限標(biāo)識(shí)碼根據(jù)八分體所在卦限進(jìn)行編號(hào)，對(duì)第0卦限八分體進(jìn)行球面位置碼設(shè)計(jì)和徑向碼設(shè)計(jì)；依據(jù)卦限對(duì)稱性質(zhì)，實(shí)現(xiàn)與其它7個(gè)八分體編碼的相互轉(zhuǎn)換。卦限標(biāo)識(shí)碼編碼如圖6.3球面位置碼三角形編碼為“1”，按照上三角形逆時(shí)針的順序、下三角形順時(shí)針“3”。編碼取值范圍為：0-333333333333333333333333，333333333333333333333333為球面剖分第24層次時(shí)四進(jìn)制碼的最大值。56.4徑向類型碼球體沿半徑方向剖分，采用徑向等長(zhǎng)剖分方式時(shí)為0，采用徑向變長(zhǎng)剖分方式6.5徑向深度碼6.5.1球體沿半徑方向剖分時(shí)，等長(zhǎng)剖分方式和變長(zhǎng)剖分方式，分別將徑向二分并遞歸，按照從地心6.6編碼示例與方法三棱臺(tái)八叉樹(shù)網(wǎng)格編碼示例見(jiàn)附錄B，編解碼方法6素的平均幾何特征與數(shù)量如表A.1所示。這里采用徑向等長(zhǎng)剖分方式，示例的地球半徑取CGCS2000橢球長(zhǎng)半軸、短半軸的平均值6367.444657km。1234567897三棱臺(tái)八叉樹(shù)網(wǎng)格編碼示例見(jiàn)表B.1，這里的編碼100球面與徑向等長(zhǎng)方式同260面與徑向等長(zhǎng)方式同時(shí)剖310面剖分13次，徑向等長(zhǎng)方441剖分14次得到的網(wǎng)格（球511位于鄭州“鄭州二七罷在球面剖分24次、徑向變長(zhǎng)方式剖分23次，開(kāi)孔所610位于鄭州“鄭州二七罷在球面剖分24次、徑向等長(zhǎng)方式剖分23次，終孔1000米深度（假設(shè)鉆孔完全按照地球半徑方向鉆8網(wǎng)格的編碼是指已知點(diǎn)的大地坐標(biāo)求其網(wǎng)格編碼，解碼是指已知網(wǎng)格編碼求其網(wǎng)格中心的大地坐根據(jù)八分體跨經(jīng)緯90°的剖分機(jī)理，卦限碼Q的計(jì)算公式見(jiàn)C.1：Q=λ\90°-4×((φ-90°)\90°)………………（C.1）λ——為經(jīng)度（0°~360°,注意與東西經(jīng)的轉(zhuǎn)換）；φ——為緯度，北緯為正，南緯為負(fù)。如：點(diǎn)（E113.6775°,N34.7500°)的卦限碼為1。 (C.2）三棱臺(tái)八叉樹(shù)剖分的球面三角形A的位置碼表示為：A=a1a2a3??ak,a1到ak是k個(gè)四分碼，根據(jù)本省地質(zhì)應(yīng)用需求，本文件規(guī)定該碼最多為24位，k為剖分層次。固定方向編碼的解碼方法很多，這里以“行列逼近法”（趙學(xué)勝等，2003，QTM地址碼與經(jīng)緯度坐標(biāo)的快速轉(zhuǎn)換算法）為例說(shuō)C.3.2.1根據(jù)三角形的層次K求出其最大的行列數(shù)(I,J)，見(jiàn)公式C.3。I=2k,J=2×I-1………（C.3）C.3.2.2確定經(jīng)緯度坐標(biāo)(λ,φ)在網(wǎng)格地址碼的行列數(shù)(i,j)，見(jiàn)公式C.4。C.3.2.3根據(jù)所求位置在網(wǎng)格的行列數(shù)(i,j)和網(wǎng)格最大列數(shù)(I,J)，按一定的方向逐層次地進(jìn)行遞歸9a)如果i大于n/2，地址碼ak為“1”，同時(shí)n=n/2,I=(i-n)；b)如果i小于或等于n/2，分兩種情況：1)如果j小于地址碼ak為“2”，同時(shí)n=n/2；2)否則，若j大于n，則地址碼ak為“3”，同時(shí)n=n/2,j=(j-2n)。若j小于或等于n，則地址碼ak為“0”，同時(shí)j=(j-n+2i-1)，I=(-i+1),n=n/2。c)重復(fù)a）、b）步驟，求下一個(gè)地址編碼ak+1。示例：設(shè)剖分層次為24層，點(diǎn)（E113.6775°,N34.7500°)的球面位置碼為：020230212003022223131131。C.3.3.1根據(jù)三角形的層次K求出其最大的行列數(shù)(I,J)I=2k,J=2×I-1……………（C.5）C.3.3.2根據(jù)網(wǎng)格的最大行列數(shù)(I,J),把地址碼根據(jù)固定的方向按次序進(jìn)行逐層次地解算，求出位置在格網(wǎng)的行數(shù)i（注意：由于地址碼為“1”時(shí)，三角形可為頂或底，這對(duì)計(jì)算行數(shù)有很大影響。所以，a)如果地址碼ak=“1”：若p=1，則i=i+；b)如果地址碼ak=“2”或“3”：若p=-1，則；若p=1,則n=c)如果地址碼ak=“0”：若p=-1，則i；若p=1,則n=n/2，且p=-p；d)重復(fù)步驟a）、b）、c直到最后求出三角形所在的行數(shù)i。C.3.3.3根據(jù)行數(shù)i和格網(wǎng)最大行列數(shù)(I,J)，把地址碼按次序進(jìn)行逐層次地解算，求出位置在格網(wǎng)的a)如果地址碼ak=“1”：則n=n/2，i=i-n；d)如果地址碼ak=“0”：則n=n/2，i=n-i+1,j=j+(2×(n-i)+1)；e)重復(fù)步驟（1234直到最后求出三角形所在的列數(shù)j。C.3.3.4根據(jù)網(wǎng)格地址碼的行列數(shù)(i,j)確定三角形中點(diǎn)的經(jīng)緯度坐標(biāo)(λ,φ)，見(jiàn)公式C.6。分?jǐn)?shù)，取其分子的整數(shù)部分加上1（若為整數(shù)則不加1），然后再將其轉(zhuǎn)化為k位的二進(jìn)深度為00101000001101000101101。變長(zhǎng)剖分的原理與之相同，但剖分的比例不同。上述例子的變