杜靈瑀，馬秋禾，賁 進(jìn)，王 蕊

信息工程大學(xué)地理空間信息學(xué)院，河南 鄭州 450001

格網(wǎng)系統(tǒng)是一種多分辨率柵格層次數(shù)據(jù)結(jié)構(gòu)，本質(zhì)是通過(guò)規(guī)則格網(wǎng)單元的分解與聚合將連續(xù)平面映射到多層次格網(wǎng)中，且每個(gè)單元具有包含層次關(guān)系的唯一編碼[1]。格網(wǎng)系統(tǒng)的離散性使其更適合計(jì)算機(jī)處理，有助于多源、異構(gòu)、多尺度空間數(shù)據(jù)的融合及多要素的空間疊加和復(fù)合分析[2]。三角形格網(wǎng)系統(tǒng)廣泛應(yīng)用于地形建模[3-5]、空間數(shù)據(jù)索引與可視化[6]、地圖綜合模型[7]；矩形(四邊形)格網(wǎng)系統(tǒng)廣泛應(yīng)用于資源環(huán)境綜合科學(xué)調(diào)查[8]、遙感影像數(shù)據(jù)存儲(chǔ)與管理[9]及地圖多級(jí)顯示[10]。

相比于上述格網(wǎng)，六邊形格網(wǎng)具有形態(tài)緊湊、平面覆蓋率和角度分辨率高及一致相鄰等特性[11]，更有利于空間數(shù)據(jù)的組織、處理和分析。然而，六邊形不具備自相似性，格網(wǎng)層次關(guān)系描述及計(jì)算是研究難點(diǎn)之一[12]。文獻(xiàn)[13—14]提出“廣義平衡三進(jìn)制”(generalized balanced ternary，GBT)，將高分辨率六邊形單元聚合為低分辨率單元，但其相鄰層次單元的方向連續(xù)旋轉(zhuǎn)，導(dǎo)致相關(guān)應(yīng)用較復(fù)雜[15]。文獻(xiàn)[1，16—17]分別設(shè)計(jì)了三孔六邊形格網(wǎng)系統(tǒng)編碼方案，并借助改進(jìn)的GBT實(shí)現(xiàn)不同層次單元快速索引，其編碼運(yùn)算無(wú)法完全用二進(jìn)制優(yōu)化，效率不甚理想[18]。文獻(xiàn)[12]以“格元—格點(diǎn)—格邊—格心”概念模型[19]為基礎(chǔ)，建立了三孔六邊形格網(wǎng)系統(tǒng)編碼方案，并揭示其數(shù)學(xué)本質(zhì)。文獻(xiàn)[18]針對(duì)四孔六邊形格網(wǎng)系統(tǒng)設(shè)計(jì)了“六邊形四元平衡結(jié)構(gòu)”(hexagonal quaternary balanced structure，HQBS)，但在編碼運(yùn)算過(guò)程中頻繁出現(xiàn)“正則化”失敗回退重算的情況[12]。針對(duì)這一缺陷，文獻(xiàn)[20]設(shè)計(jì)了“六邊形格點(diǎn)四叉樹(shù)”(hexagon lattice quad tree,HLQT)，其編碼運(yùn)算效率優(yōu)于HQBS及PYXIS，但其碼元分布不對(duì)稱，不利于鄰近單元查詢。文獻(xiàn)[21]采用菱形塊四叉樹(shù)組織球面六邊形格網(wǎng)層次結(jié)構(gòu)，并提出格網(wǎng)編碼方案。在實(shí)踐應(yīng)用方面，ESRI公司在ArcGIS GeoAnalytics Server中采用六邊形格網(wǎng)系統(tǒng)分析、展示數(shù)據(jù)(http:∥www.esrichina.com.cn/openclass-2017.html)。UBER公司在其流處理系統(tǒng)中采用六邊形格網(wǎng)系統(tǒng)實(shí)時(shí)采集、分析車輛信息(http:∥www.infoq.com/cn/presentations/uber-stream-processing-system-and-practice)。歐空局的SMOS衛(wèi)星也采用六邊形格網(wǎng)存儲(chǔ)、處理影像數(shù)據(jù)[22]。

綜上所述，當(dāng)前六邊形格網(wǎng)系統(tǒng)的研究已取得可喜進(jìn)展，但現(xiàn)有成果仍存在編碼方案不完善、運(yùn)算效率待提升的問(wèn)題。本文引入復(fù)進(jìn)制數(shù)理論,依據(jù)間隔層次單元隸屬關(guān)系，建立平面四孔六邊形格網(wǎng)系統(tǒng)數(shù)學(xué)模型，在此基礎(chǔ)上提出具有結(jié)構(gòu)對(duì)稱性的等效單元編碼方案、定義編碼運(yùn)算并歸納運(yùn)算規(guī)則、設(shè)計(jì)編碼索引及編碼與笛卡兒坐標(biāo)互換算法，嘗試提高編碼操作效率。

1 復(fù)進(jìn)制數(shù)建模

復(fù)進(jìn)制定位計(jì)數(shù)系統(tǒng)(complex radix positional system)由計(jì)算機(jī)科學(xué)泰斗Donald E.Knuth于1960年提出[23]，它的基數(shù)(進(jìn)制)、數(shù)字位集合元素可以是復(fù)數(shù)，本質(zhì)上是二進(jìn)制、十進(jìn)制等常見(jiàn)定位計(jì)數(shù)系統(tǒng)在復(fù)數(shù)域上的推廣。復(fù)數(shù)在復(fù)進(jìn)制定位計(jì)數(shù)系統(tǒng)中的表示形式稱為復(fù)進(jìn)制數(shù)(complex radix number)。在全球離散格網(wǎng)的相關(guān)研究中，通常采用笛卡兒坐標(biāo)描述單元位置[24]，該方式雖然直觀，但與編碼關(guān)聯(lián)性不強(qiáng)。研究表明，以復(fù)進(jìn)制數(shù)形式表示單元位置，更有助于編碼設(shè)計(jì)及其運(yùn)算規(guī)律推導(dǎo)[25-26]。

作為位置、對(duì)象及數(shù)據(jù)與格網(wǎng)單元的關(guān)聯(lián)點(diǎn)，通常利用格心代替格網(wǎng)單元[19]。四孔六邊形格網(wǎng)的剖分規(guī)律如圖1所示，第n層格心由第n-1層格心及單元各邊中心組成，相鄰層次單元邊長(zhǎng)為2倍關(guān)系，且單元方向不隨層次改變。

圖1 四孔六邊形格網(wǎng)系統(tǒng)示意圖Fig.1 Diagram of aperture 4 hexagon grid system

式中，“∪”、“+”是集合間的并運(yùn)算、加法運(yùn)算符。

證明：設(shè)第n-1層格網(wǎng)單元外接圓半徑為rn-1，第n層格心集合為L(zhǎng)n，根據(jù)格網(wǎng)系統(tǒng)生成規(guī)律，Ln可表示為

所以

此時(shí)有

(1)

不同于HLQT、HQBS等四孔六邊形格網(wǎng)系統(tǒng)逐層編碼方案以及PYXIS三孔六邊形格網(wǎng)系統(tǒng)奇、偶層分別編碼方案[1]，這里的定位計(jì)數(shù)系統(tǒng)建立在層次遞推關(guān)系的基礎(chǔ)上。隨層次增加，格心呈現(xiàn)向內(nèi)凹陷加密的趨勢(shì)。如圖4所示，用不同大小和顏色的點(diǎn)表示L1—L5中的格點(diǎn)，圖中空缺部分可由首層相鄰單元生成的子格點(diǎn)填補(bǔ)。

2 等效編碼方案

由上述規(guī)則直接給出L1與L2的編碼，如圖5所示。以L1、L2和集合D的編碼為基礎(chǔ)，依據(jù)層次遞推關(guān)系建立Ln(n≥3)編碼，由Ln-1得到Ln中繼承子單元的過(guò)程編碼層次增加但格心位置不變，因此規(guī)定Ln(n≥3)中繼承子單元編碼由Ln-1編碼末尾添加“00”得到，如圖5L3編碼中黑色單元所示。剖分子單元的格心由第n-2層格網(wǎng)單元依據(jù)集合D剖分得到，但第n層與第n-2層并非相鄰層次，因此在Ln-2編碼末尾添加“00”表示跨層后再添加集合D對(duì)應(yīng)的等效碼元，如圖5L3編碼中紅色單元所示。

根據(jù)上述編碼規(guī)則，每層碼元由2位字符構(gòu)成，由于單元?dú)w屬明確，單元編碼具有唯一性；復(fù)進(jìn)制數(shù)集合D中的元素關(guān)于原點(diǎn)對(duì)稱，編碼具有對(duì)稱性，有利于編碼索引操作。

圖2 基于層次遞推的格網(wǎng)系統(tǒng)層次關(guān)系Fig.2 Hierarchical relation of grid system based on hierarchical recursion

圖4 L1—L5在復(fù)平面上的分布Fig.4 Distribution of L1—L5 in the complex plane

圖3 集合D在復(fù)平面上的表示Fig.3 Representation of set D in the complex plane

3 編碼運(yùn)算規(guī)則和索引算法

編碼記錄了格網(wǎng)單元相對(duì)原點(diǎn)的距離和方位，可通過(guò)一維編碼運(yùn)算實(shí)現(xiàn)二維復(fù)平面內(nèi)的幾何操作，達(dá)到降維目的。相較于浮點(diǎn)數(shù)，編碼更適合計(jì)算機(jī)存儲(chǔ)與處理，有利于提高運(yùn)算效率。

3.1 編碼加法運(yùn)算規(guī)則

與矢量加法類似，可利用平行四邊形法則定義編碼加法運(yùn)算[27]，以符號(hào)⊕表示。按照層次關(guān)系，Ln(n≥3)與Ln-1和Ln-2均相關(guān)，在編碼加法中，以4位(即每?jī)蓪?編碼為基礎(chǔ)碼元

圖5 編碼示意圖Fig.5 Diagram of code

{0000,0001,0002,0003,0004,0005,0006}

{0010,0020,0030,0040,0050,0060}

{0100,0200,0300,0400,0500,0600}

{1000,2000,3000,4000,5000,6000}

依照平行四邊形法則形成25×25加法查找表，見(jiàn)表1。將加法表結(jié)果統(tǒng)一記為8位，則加法運(yùn)算的進(jìn)位操作僅存在于相鄰基礎(chǔ)碼元。

表1 編碼加法查找略表

圖6 錯(cuò)位相加法示意圖Fig.6 Diagram of staggered addition

對(duì)于第5層編碼加法0000010003⊕0000000003，錯(cuò)位相加法得到的結(jié)果為0000010300，而實(shí)際結(jié)果為0000001000。分析可知，錯(cuò)位相加法以單元0000010003為中心加上0000000003得到最終結(jié)果，而實(shí)際編碼0000001000由第4層繼承得到，其生成中心為原點(diǎn)單元，中心單元的不一致導(dǎo)致計(jì)算結(jié)果的不一致。出現(xiàn)該問(wèn)題的編碼存在連續(xù)非“00”編碼的明顯特征，以編碼α=α1α2…α2n-1α2n為例，αk-3αk-2αk-1αk是出現(xiàn)連續(xù)非“00”編碼的位置，將α1α2…αk轉(zhuǎn)換為α1α2…αk-3αk-200⊕00…00αk-1αk的結(jié)果，并在末尾加上αk后的編碼αk+1…α2n可將編碼α轉(zhuǎn)化為正確編碼，避免“回退”操作。

3.2 編碼乘法運(yùn)算規(guī)則

編碼乘法對(duì)應(yīng)旋轉(zhuǎn)和縮放操作[26]，以符號(hào)?表示。集合{01,02,03,04,05,06}中的碼元分別對(duì)應(yīng)旋轉(zhuǎn)0°、逆時(shí)針旋轉(zhuǎn)60°、逆時(shí)針旋轉(zhuǎn)120°、旋轉(zhuǎn)180°、順時(shí)針旋轉(zhuǎn)120°、順時(shí)針旋轉(zhuǎn)60°。據(jù)此可得編碼乘法運(yùn)算查找表，見(jiàn)表2。

表2 編碼乘法運(yùn)算查找表

乘法運(yùn)算的碼元及結(jié)果均為2位編碼，將編碼按2位字符分組，逐組運(yùn)算并記錄結(jié)果即可。以0000060002?02為例，依乘法表，00?02=00，02?02=03，02?06=01，則結(jié)果編碼為0000010003，與逆時(shí)針旋轉(zhuǎn)60°結(jié)果一致。

3.3 編碼索引算法

層次單元查找分為子單元和父單元查找，按照第3部分繼承子單元和剖分子單元的編碼規(guī)律即可實(shí)現(xiàn)子單元的查找。父單元查找也分為鄰層和隔層兩種情況，對(duì)于第n(n≥3)層編碼α=α1α2…α2n-1α2n，若α2n-1α2n=00，則該單元繼承于第n-1層，將末尾編碼“00”去掉可得到鄰層父單元。若α2n-1α2n≠00，則該單元由第n-2層的單元剖分得到，將末尾4位編碼去掉可得到隔層父單元。

4 平面笛卡兒坐標(biāo)與編碼轉(zhuǎn)換

4.1 格網(wǎng)編碼到平面坐標(biāo)的轉(zhuǎn)換

4.2 平面坐標(biāo)到格網(wǎng)編碼的轉(zhuǎn)換

5 試驗(yàn)與分析

本文編碼方案可借助十六進(jìn)制數(shù)實(shí)現(xiàn)。碼元{00,01,02,03,04,05,06,10,20,30,40,50,60}與十六進(jìn)制數(shù){0,1,2,3,4,5,6,A,B,C,D,E,F}對(duì)應(yīng),可將字符編碼轉(zhuǎn)化為更適合計(jì)算機(jī)存儲(chǔ)的十六進(jìn)制整數(shù)，并借助二進(jìn)制位運(yùn)算提高計(jì)算效率。相比于其他現(xiàn)有六邊形格網(wǎng)編碼方案，HQBS和HLQT在編碼結(jié)構(gòu)和編碼運(yùn)算效率方面具有較大優(yōu)勢(shì)[18,20]。因此，選擇上述兩編碼方案進(jìn)行對(duì)比試驗(yàn)。本文采用Visual Studio 2012開(kāi)發(fā)工具實(shí)現(xiàn)本文編碼方案對(duì)應(yīng)的編碼加法、鄰近單元查詢、坐標(biāo)與編碼互換算法，并與HQBS及HLQT相應(yīng)算法比較。全部代碼均編譯為Release版本，并在相同臺(tái)式機(jī)(硬件配置：Intel Core i5-4590M CPU @ 3.30 GHz，8 GB RAM；操作系統(tǒng)：Windows 7 x64旗艦版)上測(cè)試。4—11層中逐層隨機(jī)生成2500個(gè)單元，不重復(fù)相加，得到3 123 750組加法運(yùn)算實(shí)例；隨機(jī)生成1 000 000組坐標(biāo)作為笛卡兒坐標(biāo)到編碼轉(zhuǎn)換試驗(yàn)數(shù)據(jù)；第4層隨機(jī)選取256個(gè)單元，5—11層以4倍遞增確定隨機(jī)選取單元個(gè)數(shù)作為編碼到坐標(biāo)轉(zhuǎn)換及鄰近單元查詢?cè)囼?yàn)數(shù)據(jù)。以單位時(shí)間內(nèi)完成操作的次數(shù)為效率指標(biāo)，統(tǒng)計(jì)10次測(cè)試的平均值，部分試驗(yàn)結(jié)果如圖9—圖12所示。

圖7 第5層編碼加法示意圖Fig.7 Diagram of code addition in the 5th level

圖8 笛卡兒坐標(biāo)到編碼轉(zhuǎn)換示意圖Fig.8 Diagram of Cartesian coordinates to code

圖9 編碼加法效率對(duì)比Fig.9 Efficiency comparison of code addition

分析試驗(yàn)結(jié)果，得出以下結(jié)論：

(1) 本文方案編碼加法的平均效率約為HQBS的10.11倍，HLQT的2.00倍。HQBS格心與頂點(diǎn)混合編碼，運(yùn)算規(guī)律復(fù)雜且易出現(xiàn)編碼“正則化”失敗回退重算的情況。HLQT逐層進(jìn)行碼元相加，計(jì)算非首位碼元相加時(shí)，除了得到當(dāng)前位計(jì)算結(jié)果碼元，還會(huì)在左側(cè)所有位產(chǎn)生進(jìn)位碼元。雖然相比于HQBS和HLQT，本文方案的加法表規(guī)模較大。但任何編碼方案中的加法表均以二維數(shù)組形式存儲(chǔ)，通過(guò)碼元在二維數(shù)組中的位置，可直接定位相應(yīng)結(jié)果，因此加法表的規(guī)模幾乎不會(huì)影響其運(yùn)算效率。同時(shí)，本文方案以每?jī)蓪訛榛居?jì)算單元，且進(jìn)位只發(fā)生在相鄰基本單元之間，相比之下計(jì)算量更小，因而效率更高。

圖10 笛卡兒坐標(biāo)到編碼轉(zhuǎn)換效率對(duì)比Fig.10 Efficiency comparison of Cartesian coordinates to code

圖11 編碼到笛卡兒坐標(biāo)轉(zhuǎn)換效率對(duì)比Fig.11 Efficiency comparison of code to Cartesian coordinates

圖12 鄰近單元查詢效率對(duì)比Fig.12 Efficiency comparison of search for neighboring cells

(2) 本文方案笛卡兒坐標(biāo)轉(zhuǎn)換為編碼的平均效率約為HQBS的4.33倍、HLQT的0.80倍。本文方案與HLQT在i、j軸上的編碼分布均具有二進(jìn)制數(shù)特征，可采用位運(yùn)算加速，因此兩者效率均高于HQBS。相較于HLQT，由于本文方案的編碼具有對(duì)稱性，在i、j軸上的部分編碼與二進(jìn)制數(shù)的轉(zhuǎn)換略復(fù)雜，導(dǎo)致轉(zhuǎn)換效率略低。筆者認(rèn)為，與對(duì)稱編碼在測(cè)試(1)、(3)、(4)中帶來(lái)的諸多優(yōu)勢(shì)相比，該轉(zhuǎn)換算法約20%的效率損失可以接受。

(3) 本文方案編碼轉(zhuǎn)換為笛卡兒坐標(biāo)的平均效率約為HQBS的6.29倍、HLQT的2.36倍。HQBS算法涉及較多的浮點(diǎn)數(shù)運(yùn)算，本文方案與HLQT均采用整數(shù)代替浮點(diǎn)數(shù)運(yùn)算，因此兩者效率均高于HQBS。相較于HLQT，由于本文方案的編碼以“00”和非“00”碼元交替形式存在，計(jì)算中只需考慮非“00”碼元，因而效率更高。

6 總結(jié)與展望

本文結(jié)合平面四孔六邊形格網(wǎng)結(jié)構(gòu)特點(diǎn)，引入復(fù)進(jìn)制數(shù)理論，通過(guò)間隔層次格網(wǎng)單元隸屬關(guān)系建立數(shù)學(xué)模型，以此為基礎(chǔ)提出等效編碼方案、定義編碼運(yùn)算并歸納運(yùn)算規(guī)則、設(shè)計(jì)編碼索引及編碼與坐標(biāo)相互轉(zhuǎn)換算法。相較于同類成果，本文編碼方案具有結(jié)構(gòu)對(duì)稱性，在編碼加法、鄰近單元查詢及本方案編碼轉(zhuǎn)換為笛卡兒坐標(biāo)等方面效率優(yōu)勢(shì)明顯，有助于各類格網(wǎng)處理、分析算法的實(shí)現(xiàn)，具有較大的應(yīng)用潛力。

本文僅在平面內(nèi)建立了編碼方案并實(shí)現(xiàn)了相關(guān)操作，借助多面體及合適的映射關(guān)系可將本文編碼方案擴(kuò)展到球(參考橢球)面，實(shí)現(xiàn)全球四孔六邊形離散格網(wǎng)系統(tǒng)的構(gòu)建。在此過(guò)程中涉及的相關(guān)問(wèn)題，將是筆者下一步研究的重點(diǎn)內(nèi)容。