全球遙感數(shù)據(jù)剖分組織的GeoSOT網(wǎng)格應用

宋樹華,程承旗,濮國梁,安豐光,羅 旭

1．北京大學遙感與地理信息系統(tǒng)研究所,北京 100871;2．北京大學航空航天信息工程研究所,北京 100871;3．中煤科技集團公司,北京 100013

全球遙感數(shù)據(jù)剖分組織的GeoSOT網(wǎng)格應用

宋樹華1,程承旗2,濮國梁2,安豐光1,羅 旭3

1．北京大學遙感與地理信息系統(tǒng)研究所,北京 100871;2．北京大學航空航天信息工程研究所,北京 100871;3．中煤科技集團公司,北京 100013

針對目前不同部門按自身行業(yè)特點采用不同數(shù)據(jù)組織網(wǎng)格的問題,為尋求更適合于現(xiàn)有測繪數(shù)據(jù)組織體系兼容的遙感數(shù)據(jù)組織網(wǎng)格,提出基于GeoSOT網(wǎng)格的遙感數(shù)據(jù)組織方案,理論證明了Geo-SOT網(wǎng)格與國家地形圖圖幅和Worldwind、Google Earth、Google Maps、Bing Maps、天地圖等網(wǎng)格具有很好的同構性,有利于對傳統(tǒng)測繪數(shù)據(jù)的繼承。同時,在不改變現(xiàn)有數(shù)據(jù)組織體系的前提下,提出基于GeoSOT全球遙感數(shù)據(jù)“虛擬一張網(wǎng)”的數(shù)據(jù)組織模型和數(shù)據(jù)整合方法。通過試驗證明,基于GeoSOT遙感影像“虛擬一張網(wǎng)”的數(shù)據(jù)組織可有效提高遙感數(shù)據(jù)整合效率。

GeoSOT網(wǎng)格;剖分面片;邏輯剖分;虛擬一張網(wǎng);同構性

在實際應用中,空間數(shù)據(jù)種類繁多、關系復雜,且不同行業(yè)有不同內容的專業(yè)數(shù)據(jù),空間數(shù)據(jù)的生產(chǎn)、維護分散在不同的單位且采用的數(shù)據(jù)網(wǎng)格標準各異,造成了部門或系統(tǒng)之間、不同歷史階段之間數(shù)據(jù)產(chǎn)品規(guī)格不統(tǒng)一[1]。例如,衛(wèi)星地面數(shù)據(jù)接收中心的原始遙感影像數(shù)據(jù)按照軌道景或軌道條帶組織,而軌道景、軌道條帶主要依據(jù)格網(wǎng)參考系GRS(Grid Reference System)和世界參考系WRS(Worldwide Reference System)所形成的固定參考網(wǎng)格并利用Path/Row對軌道景編碼[2-3];測繪部門的正射影像數(shù)據(jù)產(chǎn)品,高級數(shù)據(jù)產(chǎn)品(4級以上)按照地圖圖幅標準進行組織,低級數(shù)據(jù)產(chǎn)品按照軌道景數(shù)據(jù)組織;廣泛使用的空間數(shù)據(jù)服務,如Worldwind、Google Earth、天地圖等,均采用構建影像金字塔的方式進行數(shù)據(jù)組織。這種不同部門建立各自獨立的遙感數(shù)據(jù)組織和索引方式,可很好地適應本部門內業(yè)務需要,但跨部門間數(shù)據(jù)應用時,不同的數(shù)據(jù)組織和記錄方式給用戶帶來了很大的困難,不利于空間信息的檢索、整合、共享與分發(fā)等。

為解決上述問題,本文提出基于GeoSOT (geographical coordinate subdividing grid with one dimension integer coding on 2n-tree)網(wǎng)格[4]的全球遙感影像數(shù)據(jù)剖分組織,為遙感數(shù)據(jù)高效應用奠定基礎。

1 國內外研究現(xiàn)狀

近年來,國內外學者對全球離散網(wǎng)格[5-6]有很深入的研究,取得了很多成果,具體歸納為如下幾類。第1類是基于柏拉圖立體的正多面體(4面體、立方體、8面體、12面體和20面體)網(wǎng)格[7-8],如O-QTM(octahedral-quaternary triangular mesh)[9-10]、SQT(sphere quadtree)[11-12]和STQIE(spherical triangle quadtree based on Icosahedron and ERLRP)[13-15];第2類是基于正六邊形的全球離散網(wǎng)格[8,16-18];第3類是基于Voronoi圖的自適應全球離散網(wǎng)格[19-21];第4類是基于地理坐標劃分的全球離散網(wǎng)格,如多級網(wǎng)格[22]、Worldwind網(wǎng)格[2]、Google Earth和Google Maps網(wǎng)格[2]、Bing Maps網(wǎng)格[23]、天地圖網(wǎng)格[23]和我國測繪數(shù)據(jù)網(wǎng)格[24]等。

這些全球離散網(wǎng)格具有自身的特點。第1類全球離散網(wǎng)格可使得各網(wǎng)格單元在全球范圍內的形狀相似、各向同性,但各單元為球面三角形,從多面體到球體的映射關系比較復雜[25];第2類全球離散網(wǎng)格最緊湊、具有唯一的鄰域,可用6個離散的速度向量描述連續(xù)的各向同性的流體[26],但上下層之間的網(wǎng)格單元不像三角形網(wǎng)格那樣具有層次性和對稱性,而是存在一定角度的旋轉;第3類網(wǎng)格可針對空間數(shù)據(jù)的分布情況動態(tài)生成網(wǎng)格,適合組織空間點數(shù)據(jù),但該網(wǎng)格單元不規(guī)則、無層次性,無法與比例尺和分辨率對應,且動態(tài)生成Voronoi圖計算復雜,結構繁瑣[27];第4類全球離散網(wǎng)格多是按經(jīng)緯度網(wǎng)格進行劃分,坐標轉換方便,但高緯度地區(qū)存在一定的變形[27]。針對遙感影像數(shù)據(jù)的矩形點陣結構的特點,基于三角形、六邊形和Voronoi圖等的網(wǎng)格,由于網(wǎng)格單元與影像數(shù)據(jù)塊單元不一致,在組織管理遙感影像數(shù)據(jù),需大量復雜數(shù)據(jù)轉換工作,使得這3類網(wǎng)格僅用于理論研究和少數(shù)應用[28-29],組織管理遙感影像數(shù)據(jù)均存在一定的局限性?；诮?jīng)緯度坐標的全球離散網(wǎng)格單元與影像數(shù)據(jù)點陣相似以及坐標轉換簡單等優(yōu)點,被Google Earth、Worldwind、天地圖等廣泛地應用。因此,本文主要對經(jīng)緯度坐標全球離散網(wǎng)格討論。同時,我國遙感影像數(shù)據(jù)產(chǎn)品多按地圖圖幅組織,Worldwind、Google Earth、Bing Maps和天地圖等網(wǎng)格單元大小并不是整度、整分和整秒劃分,不能無縫地拼接為地圖圖幅(具體見后面理論證明),對跨部門之間遙感數(shù)據(jù)組織網(wǎng)格的兼容性和歷史數(shù)據(jù)繼承等問題考慮較少,在制作遙感影像數(shù)據(jù)產(chǎn)品還需數(shù)據(jù)拼接、裁剪等數(shù)據(jù)處理。

在此背景下,以北京大學程承旗教授為首的研究團隊提出了GeoSOT剖分網(wǎng)格理論,期望利用GeoSOT網(wǎng)格獨特的性質,解決跨部門之間、部門內各業(yè)務階段遙感影像數(shù)據(jù)組織基準不統(tǒng)一的問題。

2 GeoSOT網(wǎng)格及特點分析

2．1 GeoSOT網(wǎng)格與編碼

為了使得全球網(wǎng)格單元(剖分面片)大小為二進制整型、便于與地圖圖幅網(wǎng)格兼容,GeoSOT網(wǎng)格通過地球表面經(jīng)緯度范圍空間經(jīng)過3次擴展后再對其進行嚴格的遞歸四叉剖分。具體GeoSOT網(wǎng)格劃分方法如下:首先將地球通過簡單投影變換到平面,將180°×360°地球表面空間擴展為512°×512°,并將該空間作為第0級剖分面片(cell),且面片中心與赤道和本初子午線的交點重合(見圖1(a)和圖1(b))。在此基礎上,對第0級剖分面片進行遞歸四叉剖分,直到1°剖分面片。為了進行整分剖分,對1°剖分面片先進行擴展,即從60′擴展為64′(見圖1(e)),再對擴展的剖分面片遞歸四叉剖分,直到1′大小的剖分面片。對于1′網(wǎng)格面片,同樣先將其擴展到64″(見圖1 (f)),然后遞歸四叉剖分直到32級(1/2048)″。這樣,GeoSOT網(wǎng)格產(chǎn)生整度、整分、整秒和秒下級等多級剖分網(wǎng)格(見圖1(c))。這樣,經(jīng)緯度坐標空間經(jīng)過3次擴展后,整個地球表面經(jīng)緯度空間在經(jīng)線方向和緯線方向可嚴格的整型二分,由此將整個地球分割為大到全球、小到厘米級的整度、整分和整秒的層次網(wǎng)格。

圖1 GeoSOT網(wǎng)格多級剖分與3次擴展示意圖Fig．1 Multi-level subdivision and three extension of GeoSOT

針對GeoSOT網(wǎng)格剖分方案,采用64位編碼對各級剖分面片進行標識。其中,第0級剖分面片標識為0,第1級剖分面片編碼分別為00、01、02和03(見圖2(a)),并從第2級開始,下一級剖分面片的編碼在上一級剖分面片編碼后追加0、1、2和3。具體編碼規(guī)則是,距赤道和本初子午線的交點最近的剖分面片為0,最遠的為3,然后按照先沿緯線方向再沿經(jīng)線方向對其他兩個剖分面片分別為1和2(見圖2(b))。

圖2 GeoSOT網(wǎng)格單元的編碼Fig．2 GeoSOT cells’codes

2．2 GeoSOT網(wǎng)格特點分析

由于GeoSOT網(wǎng)格是基于經(jīng)緯度整度、整分、整秒的遞歸四叉剖分,與Google Earth、Worldwind、Bing Maps和天地圖等格比較,GeoSOT網(wǎng)格與測繪數(shù)據(jù)地圖圖幅具有很好的同構性,證明如下。

為了方便后面的證明,首先介紹有關網(wǎng)格同構的幾個概念,且討論的全球經(jīng)緯度網(wǎng)格是將球面網(wǎng)格投影到平面后的全球經(jīng)緯度網(wǎng)格。

2．2．1 基本概念

定義1:在全球經(jīng)緯度網(wǎng)格G中,若映射f適合

式中,t、b、l、r代表網(wǎng)格單元的四邊長度;θ代表任意兩鄰邊的夾角,則稱C為網(wǎng)格G的單元。其中,t、b、l、r、θ單位為度(°)。

第i層全球經(jīng)緯度網(wǎng)格G的單元C,記作C(G,i)或Ci,表示如下

在全球經(jīng)緯度網(wǎng)格G單元Ci＝f(ti,bi,li, ri,θ)中,參數(shù)取值如下

式中,Lt、Lb、Bl、Br為整個地球或者第1層單元格的經(jīng)度和緯度長度,單位為度(°)。當Lt、Lb、Bl、Br、k1、k2、k3、k4和θ0取不同值時,即可表示不同全球經(jīng)緯度網(wǎng)格。

根據(jù)定義1,NASA Worldwind、Google Earth、Bing Maps、天地圖、我國測繪數(shù)據(jù)規(guī)格網(wǎng)格和GeoSOT等網(wǎng)格單元Ci＝f(ti,bi,li,ri, θ),其參數(shù)取值如表1所示。

定義2:對于全球經(jīng)緯度網(wǎng)格G1和G2,C (G1,i)＝f1(ti,bi,li,ri,θ1)和C(G2,j)＝f2(tj,bj,lj,rj,θ2)(i,j∈N)分別是G1第i層和G2第j層的網(wǎng)格單元。若n(n∈N)個相鄰的C(G1,i)覆蓋的空間范圍與C(G2,j)覆蓋的空間范圍相同,則稱G1可聚合成G2。特別的,G1與G2可以是同一個全球經(jīng)緯度網(wǎng)格。

若tj＝λti,bj＝λbi,lj＝μli,rj＝μri,θ1＝θ2＝θ且n＝λμ(1≤λ≤n,1≤μ≤n,λ∈N,μ∈N, n∈N)時,稱為長度聚合(如圖3),表示為

定義3:對于全球經(jīng)緯度網(wǎng)格G1和G2, C(G1,i)(i∈N)和C(G2,j)(j∈N)分別是G1和G2的單元。若有限個C(G1,i)聚合而成的空間范圍與C(G2,j)的覆蓋空間范圍不相等,但形狀相同,則稱G1弱聚合成G2。

表1 不同經(jīng)緯度全球剖分網(wǎng)格單元f(ti,bi,li,ri,θ)參數(shù)取值Tab．1 Values of f(ti,bi,li,ri,θ)in different global subdivision grids

圖3 長度聚合示意圖Fig．3 Length of aggregation

定義4:對于全球經(jīng)緯度網(wǎng)格G1、G2和G3,若G1可聚合成G2,而G1可弱聚合成G3,則稱G2與G3弱同構。

因此,對于全球經(jīng)緯度網(wǎng)格G1和G2,若G1和G2的θ相等,G1和G2必是弱同構。其證明見附錄。

定義5:對于全球經(jīng)緯度網(wǎng)格G1、G2和G3,若G1可聚合成G2,同時G1也可聚合成G3,則稱G2與G3同構。特別的,全球經(jīng)緯度網(wǎng)格G與其本身同構。

2．2．2 GeoSOT網(wǎng)格與其他數(shù)據(jù)網(wǎng)格的同構性分析

根據(jù)定義4,若要證明GeoSOT網(wǎng)格與其他全球經(jīng)緯度網(wǎng)格同構,只需證明GeoSOT網(wǎng)格可聚合成其他全球經(jīng)緯度網(wǎng)格即可,具體證明如下。

證明:對GeoSOT網(wǎng)格G1與其他全球經(jīng)緯度網(wǎng)格G2,假設C(G1,i)(i∈N)和C(G2,j)(j∈N)分別是G1第i層和G2第j層的網(wǎng)格單元。要證明G1與G2同構,也就是要證明存在正整數(shù)m(m∈N),使得G1中m個單元C(G1,i)可聚合成G2單元C(G2,j),即

根據(jù)定義2,證明式(6)成立等價于證明存在λ∈N,μ∈N,使得下面3個等式成立

也就是,當GeoSOT網(wǎng)格G1層級i≥6時,總可使得式(16)中μ為正整數(shù)。

因此,根據(jù)式(14)和式(17),對于表1中的其他全球經(jīng)緯度網(wǎng)格G2,只要GeoSOT網(wǎng)格G1層級i≥7時,總存在正整數(shù)m＝λμ使得式(6)成立,即證GeoSOT網(wǎng)格與其他全球經(jīng)緯度網(wǎng)格同構。

而Worldwind、Google Earth、Bing Maps、天地圖等網(wǎng)格與測繪數(shù)據(jù)規(guī)格網(wǎng)格具有弱同構性,可采用反證法進行證明。

證明:假設Worldwind、Google Earth、Bing Maps、天地圖等網(wǎng)格與測繪數(shù)據(jù)規(guī)格網(wǎng)格同構性。

根據(jù)定義2,只要存在λ整數(shù),使得式(18)成立,即證明Worldwind、Google Earth、Bing Maps、天地圖等網(wǎng)格與測繪數(shù)據(jù)規(guī)格網(wǎng)格同構性

由于2j始終為偶數(shù),k1＝{60,90}為2、3、5等素數(shù)的合數(shù),因此要使得λ為整數(shù),則式(20)中左邊k必須為3的倍數(shù)。而Google Earth、Worldwind、Bing Maps、天地圖等網(wǎng)格模型中k取值為{1,2,10},不是3的倍數(shù),所以式(20)中左邊結果為非整數(shù),即λ為小數(shù),與假設λ為整數(shù)矛盾,所以假設為錯誤。而Worldwind、Google Earth等網(wǎng)格與測繪數(shù)據(jù)規(guī)格網(wǎng)格的夾角θ相等,因此Worldwind、Google Earth等網(wǎng)格與測繪數(shù)據(jù)規(guī)格網(wǎng)格弱同構性。

3 基于GeoSOT網(wǎng)格的遙感數(shù)據(jù)組織

鑒于GeoSOT網(wǎng)格的多尺度性以及與我國測繪數(shù)據(jù)規(guī)格網(wǎng)格同構的特點,在不改變現(xiàn)有數(shù)據(jù)組織的前提下,本文提出通過邏輯剖分,建立以GeoSOT網(wǎng)格剖分面片為遙感影像景劃分單元,實現(xiàn)基于空間區(qū)域的遙感數(shù)據(jù)“全球一張網(wǎng)”組織,提高多源多類型遙感數(shù)據(jù)使用效率。

3．1 遙感數(shù)據(jù)“全球一張網(wǎng)”數(shù)據(jù)組織

我國經(jīng)過幾十年的發(fā)展積累了大量的遙感數(shù)據(jù),若這些數(shù)據(jù)都按GeoSOT剖分面片進行物理分景處理,必是一個耗時耗力的巨大工程,且可能對目前遙感數(shù)據(jù)生產(chǎn)流程帶來影響。因此,本文提出的基于GeoSOT遙感數(shù)據(jù)組織思路是,在不改變目前遙感數(shù)據(jù)組織模式的前提下,首先根據(jù)遙感數(shù)據(jù)的地理坐標范圍,確定該數(shù)據(jù)的最小外包面片MBC(minimum bounding cell),并將MBC面片編碼作為整幅遙感數(shù)據(jù)的剖分標識,寫入剖分索引表中,然后根據(jù)遙感數(shù)據(jù)的空間分辨率確定所在GeoSOT層,并通過GeoSOT剖分面片單元4個角點地理坐標與遙感數(shù)據(jù)的地理坐標對應,對遙感數(shù)據(jù)按GeoSOT網(wǎng)格剖分面片空間范圍進行“虛擬劃分”,即邏輯剖分,將遙感影像數(shù)據(jù)覆蓋剖分面片的編碼與剖分面片4個角點對應遙感影像數(shù)據(jù)像素行列號等信息,寫入與遙感影像文件名相同的剖分索引文件中,由此建立基于GeoSOT剖分面片編碼的遙感數(shù)據(jù)塊剖分索引,從而實現(xiàn)遙感數(shù)據(jù)按空間區(qū)域數(shù)據(jù)組織。由于GeoSOT網(wǎng)格剖分面片覆蓋全球且大小固定,與剖分面片關聯(lián)的遙感數(shù)據(jù)也虛擬覆蓋全球且遙感數(shù)據(jù)大小規(guī)格一定,形成遙感數(shù)據(jù)虛擬“全球一張網(wǎng)”的數(shù)據(jù)組織體系,如圖4所示。

圖4 遙感數(shù)據(jù)建立邏輯剖分的流程圖Fig．4 Logical subdivision of remote sensing data

此方案中,遙感數(shù)據(jù)與GeoSOT網(wǎng)格層級l的對應關系采用的如下公式確定

式中,r為遙感數(shù)據(jù)分辨率,單位為°/像素;l為GeoSOT網(wǎng)格層級;k為常數(shù),一般取2－n,n∈N。

同時,由于各級剖分面片空間范圍固定,而不同等級的影像數(shù)據(jù)產(chǎn)品幾何精度不同,導致基于剖分面片“分割”的遙感影像數(shù)據(jù)會產(chǎn)生一些問題。例如,對地理精度誤差較大的低級別(3級以下)的遙感影像,由于其存在較大的幾何精度誤差,導致本應該與剖分面片A關聯(lián)的數(shù)據(jù)塊,分割后可能與A相鄰的剖分面片B關聯(lián)。為了避免出現(xiàn)上述現(xiàn)象,在對不同地理精度誤差的遙感影像數(shù)據(jù)邏輯剖分時,分別采用冗余邏輯剖分和精確邏輯剖分。對地理精度誤差較大的低級別(3級以下)遙感影像數(shù)據(jù)邏輯剖分時,將剖分面片空間范圍向周圍擴展遙感影像地理誤差的寬度,并按擴展后的剖分面片“分割”遙感影像,即為冗余邏輯剖分,從而使得地球表面與剖分面片地理范圍對應區(qū)域的遙感影像數(shù)據(jù)始終在關聯(lián)的剖分面片內。圖5為地理精度誤差為d的遙感數(shù)據(jù)產(chǎn)品冗余邏輯剖分示意圖,實線單元格為GeoSOT剖分面片,虛線網(wǎng)格為根據(jù)精度誤差d擴展后的冗余剖分面片,冗余邏輯剖分即按虛線對影像數(shù)據(jù)進行“分割”。而對于地理精度誤差較小的高級(4級以上)影像數(shù)據(jù)產(chǎn)品,直接根據(jù)GeoSOT剖分面片的空間范圍進行邏輯剖分,建立遙感數(shù)據(jù)邏輯剖分索引,即為精確邏輯剖分,如圖6所示。

圖5 低級影像數(shù)據(jù)產(chǎn)品的冗余邏輯剖分示意圖Fig．5 Redundant logical subdivision of lowerlevel remote sensing data products

圖6 高級影像數(shù)據(jù)產(chǎn)品的精確邏輯剖分示意圖Fig．6 Precise logical subdivision of advanced remote sensing data products

3．2 基于剖分面片聚合的遙感數(shù)據(jù)整合

基于GeoSOT遙感數(shù)據(jù)“全球一張網(wǎng)”的數(shù)據(jù)組織,使得遙感數(shù)據(jù)具有共同的數(shù)據(jù)組織基礎,可避免不同部門、不同型號衛(wèi)星遙感數(shù)據(jù)采用的數(shù)據(jù)組織網(wǎng)格不一致問題。同時使得大區(qū)域的遙感數(shù)據(jù)可直接聚合得到,即當利用遙感數(shù)據(jù)制作某規(guī)格數(shù)據(jù)產(chǎn)品時,根據(jù)規(guī)格數(shù)據(jù)產(chǎn)品的空間范圍信息,例如地圖圖幅范圍,計算出規(guī)格產(chǎn)品范圍覆蓋的剖分面片,并根據(jù)剖分面片編碼提取相關的遙感數(shù)據(jù)。最后,根據(jù)GeoSOT網(wǎng)格與其他規(guī)格數(shù)據(jù)網(wǎng)格同構性,通過剖分面片聚合的方式即可快速生成規(guī)格數(shù)據(jù)產(chǎn)品(如圖7),可有效減少參與計算的遙感數(shù)據(jù)數(shù)據(jù)量,從而提高遙感數(shù)據(jù)整合效率。

圖7 基于GeoSOT網(wǎng)格剖分面片聚合的數(shù)據(jù)產(chǎn)品生產(chǎn)過程Fig．7 Process of products aggregated by GeoSOT cells

4 試驗及分析

根據(jù)GeoSOT網(wǎng)格及編碼模型,在Worldwind開源代碼的基礎上,利用Microsoft Visual studio 2008 C＃作為開發(fā)平臺,Microsoft Office Access2007作為后端數(shù)據(jù)庫系統(tǒng),在CPU主頻為2．4 GHz、內存為2 GB、硬盤轉速為7200 rpm的PC機上進行了遙感數(shù)據(jù)GeoSOT組織與數(shù)據(jù)整合的試驗驗證。試驗數(shù)據(jù)為局部地區(qū)某衛(wèi)星1A、2B和3B級約32 GB遙感影像數(shù)據(jù),分辨率分別為2 m、5 m和10 m,其格式為TIF文件。當式(21)中k取值為1/512,則上述分辨率的遙感影像分別對應GeoSOT網(wǎng)格第15層、第14層和第13層。對地理精度相對較低的1A級影像數(shù)據(jù)采用冗余邏輯剖分,而對地理精度相對較高的2B級和3B級影像數(shù)據(jù)采用精確邏輯剖分處理。3種不同級別影像建立邏輯剖分索引的耗時如表2。

表2 試驗數(shù)據(jù)統(tǒng)計表Tab．2 Experiment data in GeoSOT system

在影像邏輯剖分組織的基礎上,利用GeoSOT網(wǎng)格與測繪地形圖圖幅之間的同構關系,通過剖分面片編碼對單個剖分面片對應的影像數(shù)據(jù)以及不同區(qū)域、比例尺地圖圖幅的影像數(shù)據(jù)進行了數(shù)據(jù)檢索和提取試驗,其試驗結果如表3所示。

表3 數(shù)據(jù)整合試驗Tab．3 Data Integration in GeoSOT system

從表3可看出,由于遙感數(shù)據(jù)建立了基于GeoSOT網(wǎng)格的邏輯剖分索引,在對剖分面片或規(guī)格產(chǎn)品范圍對應的遙感數(shù)據(jù)進行數(shù)據(jù)整合具有較高的效率。

需要說明的是,本測試系統(tǒng)主要是驗證遙感數(shù)據(jù)按GeoSOT網(wǎng)格進行邏輯剖分組織的效率和基于剖分面片聚合的數(shù)據(jù)整合的效率,以此說明基于GeoSOT的遙感數(shù)據(jù)組織的可行性和遙感數(shù)據(jù)應用的高效性。

5 結 論

本文對目前常用的全球經(jīng)緯度網(wǎng)格和GeoSOT網(wǎng)格進行了分析,理論證明GeoSOT網(wǎng)格較Google Earth、Worldwind、Bing Maps和天地圖等網(wǎng)格,與現(xiàn)有的測繪數(shù)據(jù)規(guī)格網(wǎng)格具有更好的同構性。鑒于GeoSOT網(wǎng)格的多尺度性和對我國測繪數(shù)據(jù)規(guī)格網(wǎng)格同構性特點,在不改變現(xiàn)有遙感數(shù)據(jù)組織的前提下,提出了通過基于GeoSOT網(wǎng)格的邏輯剖分方法,實現(xiàn)基于空間區(qū)域的遙感數(shù)據(jù)“全球一張網(wǎng)”組織。同時,本文還提出了基于剖分面片聚合的規(guī)格數(shù)據(jù)產(chǎn)品生成方法,可有效減少參與計算的遙感數(shù)據(jù)量,從而提高遙感數(shù)據(jù)整合效率。最后,通過試驗驗證,通過對現(xiàn)有的遙感數(shù)據(jù)建立邏輯剖分索引,實施代價不大,且基于GeoSOT的統(tǒng)一數(shù)據(jù)組織有助于規(guī)格數(shù)據(jù)產(chǎn)品的數(shù)據(jù)整合。

[1] SHAO Zhenfeng,LI Deren．Study on Spatial Information Multi-grid under Grid Computing Environment[J]．Geomatics World,2005,3(2):31-35．(邵振峰,李德仁．基于網(wǎng)格計算環(huán)境下的空間信息多級格網(wǎng)研究[J]．地理信息世界,2005,3(2):31-35．)

[2] SONG Shuhua．Research on a Universal Scene Model of Remote Sensing Data and Its Applications[D]．Beijing: Peking University,2011．(宋樹華．遙感數(shù)據(jù)標準分景框架及其應用模型研究[D]．北京:北京大學,2011．)

[3] LI Jinping,LI Yonggang,LIU Xiufang．Research on the Estimate Method of Landsat WRS Grid Coordinate Position[J]．Science of Surveying and Mapping,2006,31 (3):52-53．(李錦萍,李永剛,劉秀芳．Landsat衛(wèi)星WRS格網(wǎng)坐標位置估算方法的研究[J]．測繪科學, 2006,31(3):52-53．)

[4] CHENG Chengqi,REN Fuhu,PU Guoliang,et al．An Introduce to Spatial Informaiton Subdivision Organization [M]．Beijing:Science Press,2012．(程承旗,任伏虎,濮國梁,等．空間信息剖分組織導論[M]．北京:科學出版社,2012．)

[5] SAHR K,WHITE D,KIMERLING A J．Geodesic Discrete Global Grid Systems[J]．Cartography and Geographic Information Science,2003,30(2):121-134．

[6] GOODCHILD M．Discrete Global Grids for Digital Earth [C]∥Proceedings of International Conference on Discrete Global Grids．Santa Barbara:[s．n．],2000．

[7] WHITE D．Global Grids from Recursive Diamond Subdivisions of the Surface of an Octahedron or Icosahedrons[J]．Environmental Monitoring and Assessment,2000,64 (1):93-103．

[8] GOLAY J E．Hexagonal Parallel Pattern Transformations [J]．IEEE Transactions on Computers,1969,18(8): 733-739．

[9] DUTTON G．A Hierarchical Coordinate System for Geoprocessing and Cartography[M]．Berlin:Springer-Verlag,1999．

[10] DUTTON G．Encoding and Handling Geospatial Data with Hierarchical Triangular Meshes[C]∥Proceedings of 7th International Symposium on Spatial Data Handling．Amsterdam:Taylors and Francis,1996．

[11] FEKETE G．Rendering and Managing Spherical Data with Sphere Quadtrees[C]∥Proceedings on Visualization’90．New York:ACM,1990．

[12] FEKETE G,TREINISH L．Sphere Quadtrees:A New Data Structure to Support the Visualization of Spherically Distributed Data[C]∥Extracting Meaning from Complex Data:Processing,Display,Interaction．Santa Clara: SPIE,1990,1259:242-253．

[13] YUAN Wen,CHENG Chengqi,MA Ainai,et al．L Curve for Spherical Triangle Region Quantrees[J]．Science China:Technological Science,2004,47(3):265-270．

[14] YUAN Wen．STQIE Model for GeoGrid and the Prototype System[D]．Beijing:Peking University,2004．(袁文．地理格網(wǎng)STQIE模型及原型系統(tǒng)[D]．北京:北京大學,2004．)

[15] YUAN Wen,MA Ainai,GUAN Xiaojing．A New Projection for Spherical Triangle:Equal Angle Ratio Projection (EARP)[J]．Acta Geodaetica et Cartographica Sinica,2005, 34(1):78-84．(袁文,馬藹乃,管曉靜．一種新的球面三角投影:等角比投影(EARP)[J]．測繪學報,2005,34(1): 78-84．)

[16] TONG Xiaochong,BEN Jin,QIN Zhiyuan,et al．The Subdivision of Partial Grid Based on Discrete Global Grid Systems[J]．Acta Geodaetica et Cartographica Sinica, 2009,38(6):506-513．(童曉沖,賁進,秦志遠,等．基于全球離散網(wǎng)格框架的局部網(wǎng)格劃分[J]．測繪學報,2009, 38(6):506-513．)

[17] BEN Jin,TONG Xiaochong,ZHANG Yongsheng,et al．Research on Generating Algorithm and Software Model of Discrete Global Grid Systems[J]．Acta Geodaetica et Cartographica Sinica,2007,36(2):187-191．(賁進,童曉沖,張永生,等．球面等積網(wǎng)格系統(tǒng)生成算法與軟件模型研究[J]．測繪學報,2007,36(2):187-191．)

[18] ZHANG Yongsheng,BEN Jin,TONG Xiaochong,et al．Geospatial Information Processing Method Based on Spherical Hexagon Grid System[J]．Journal of Zhengzhou Institute of Surveying and Mapping,2006,23(2):110-114．(張永生,賁進,童曉沖,等．基于球面六邊形網(wǎng)格系統(tǒng)的空間信息處理方法[J]．測繪科學技術學報,2006,23 (2):110-114．)

[19] XU Yinfeng,YANG Boting．Popularization of Voronoi Diagram [J]．Chinese Journal of Engineering&Mathematics, 1995,12(2):93-96．(徐寅峰,楊波艇．Voronoi圖在球面上的推廣[J]．工程數(shù)學學報,1995,12(2):93-96．)

[20] ZHAO X S,CHEN J,LI Z L．A Mix Global Data Structure Based on QTM and Voronoi[C]∥Proceedings of XXth ISPRS Congress．Istanbul:ISPRS,2004:791-796．

[21] ZHAO Xuesheng,CHEN Jun,WANG Jinzhuang．QTM-based Algorithm for the Generating of Voronoi Diagram for Spherical Objects[J]．Acta Geodaetica et Cartographica Sinica,2002,31(2):157-163．(趙學勝,陳軍,王金莊．基于O-QTM的球面Voronoi圖的生成算法[J]．測繪學報,2002,31(2):157-163．)

[22] LI Deren,ZHU Xinyan,GONG Jianya．From Digital Map to Spatial Information Multi-grid:A Thought of Spatial Information Multi-grid Theory[J]．Geomatics and Information Science of Wuhan University,2003,28(6):642-650．(李德仁,朱欣焰,龔健雅．從數(shù)字地圖到空間信息網(wǎng)格:空間信息多級網(wǎng)格理論思考[J]．武漢大學學報:信息科學版,2003,28(6):642-650．)

[23] LV X F,CHENG C Q,GONG J Y,et al．Review of Data Storage and Management Technologies for Massive Remote Sensing Data[J]．Science China:Technological Science,2011,54(12):3220-3232．

[24] Institute on Standardiztion of National Adminstration of Suveying, Mapping and Geoinformation．GB/T 13989-2012 Subdivision and Numbering for the National Primary Scale Topographic Maps[S]．Beijing:China Zhijian Publishing House,2012．(GB/T13989-92,(國家測繪地理信息局測繪標準化研究所．GB/T13989-2012國家基本比例尺地形圖分幅和編號[S]．北京:中國質檢出版社,2012．)

[25] LI Deren,XIAO Zhifeng,ZHU Xinyan,et al．Research on Grid Division and Encoding of Spatial Information Multigrids[J]．Acta Geodaetica et Cartographica sinica,2006,35 (1):52-56．(李德仁,肖志峰,朱欣焰,等．空間信息多級網(wǎng)格的劃分方法及編碼研究[J]．測繪學報,2006,35(1): 52-56．)

[26] FRISCH U,HASSLACHER B,POMEAU Y．Lattice-gas Automata for the Navier-Stokes Equations[J]．Physical Review Letters,1986,56(3):1505-1508．

[27] SONG Shuhua,CHENG Chengqi,GUAN Li,et al．Analysis on Global Geodata Partitioning Models[J]．Geography and Geo-Information Science,2008,24(4):11-15．(宋樹華,程承旗,關麗,等．全球空間數(shù)據(jù)剖分模型分析[J]．地理與地理信息科學,2008,24(4):11-15．)

[28] BARTHOLDI J,GOLDSMAN P．Continuous Indexing of Hierarchical Subdivisions of the Globe[J]．International Journal of Geographical Information Science,2001,15 (6):489-522．

[29] SAHR K,WHITE D,KIMMERLING A．Geodesic Discrete Global Grid Systems[J]．Cartography and Geographic Information Science,2003,30(2):121-134．

(責任編輯:叢樹平)

附 錄:

全球經(jīng)緯度網(wǎng)格G1和G2弱同構證明。

證明:對于全球經(jīng)緯度網(wǎng)格G1和G2,C(G1, i)＝f1(ti,bi,li,ri,θ1)和C(G2,j)＝f2(tj, bj,lj,rj,θ2)(i,j∈N)分別是G1第i層和G2第j層的網(wǎng)格單元,且tj＞ti＞0、lj＞li＞0、θ1＝θ2。

由于θ1＝θ2,則C(G1,i)和C(G2,j)可以進行長度聚合。而要證明G1和G2弱同構,則只要證明存在正整數(shù)n1和n2,使得n1ti＞tj, n2li＞lj即可。

取N1＝[tj/ti]＋1,N2＝[lj/li]＋1,對于任意的正整數(shù)n1和n2,當n1＞N1,n2＞N2時, n1ti＞N1ti＞titj/ti＝tj,n2li＞N2li＞lilj/li＝lj同時成立。

即證,G1與G2是弱同構。

Global Remote Sensing Data Subdivision Organization Based on GeoSOT

SONG Shuhua1,CHENG Chengqi2,PU Guoliang2,AN Fengguang1,LUO Xu3
1．Institute of Remote Sensing and GIS,Peking University,Beijing 100871,China;2．Institute for Aeronautics and Astronautics Information Engineering,Peking University,Beijing 100871,China;3．China Coal Technologies Group Corp,Beijing 100013,China

At present,there are various data grids to organize data in different department data centers．In order to seek a remote sensing image data organization grid,which is compatible with the existing surveying and mapping data,a scheme of remote sensing data organization based on GeoSOT,geographical coordinate subdividing grid with one dimension integer coding on 2n-tree,is proposed．it theoretically proves that GeoSOT has good isomorphism with National Topographic Map and other grids,such as Worldwind,Google Earth,Google Maps,Bing Maps and Mapworld,which makes GeoSOT grid inherit easily traditional surveying and mapping data and organize global remote sensing data．Under the premise of keeping the existing data organization,a virtual one global grid for global remote sensing data organization based on GeoSOT and a method of fast generating specification data products by GeoSOT cells aggregation are introduced．The test shows that it is very significantly to prove data integration efficiency with the virtual one global grid for global remote sensing data organization based on GeoSOT．

GeoSOT;subdivision cell;logical subdivision;virtual one global grid;isomorphism

SONG Shuhua(1980—),male,PhD, majors in theory and application of GIS．

P208

A

1001-1595(2014)08-0869-08

國家973計劃(61399)

2012-10-10

宋樹華(1980—),男,博士,研究方向為GIS理論與應用。

E-mail:song-shuhua＠163．com

SONG Shuhua,CHENG Chengqi,PU Guoliang,et al．Global Remote Sensing Data Subdivision Organization Based on GeoSOT [J]．Acta Geodaetica et Cartographica Sinica,2014,43(8):869-876．(宋樹華,程承旗,濮國梁,等．全球遙感數(shù)據(jù)剖分組織的GeoSOT網(wǎng)格應用[J]．測繪學報,2014,43(8):869-876．)

10．13485/j．cnki．11-2089．2014．0103

修回日期:2013-07-05