worldwind技術(shù)手冊(cè)

1、寫在前面雖說(shuō)是技術(shù)手冊(cè)，但是沒(méi)有源碼的研究就算是一個(gè)教程吧。由于本人能力有限，文章中有很多不當(dāng)之處，請(qǐng)各位原諒。 Bigheader 聯(lián)系：gdsqz2worldwind技術(shù)手冊(cè)目錄：1. WorldWind軟件架構(gòu) WW功能介紹WW與XMLClient與ServerWW使用鏈接本地緩存路徑本地功能配置http請(qǐng)求與aspWMS、WFS2. WorldWind金字塔體系瓦片金字塔詳解及其對(duì)應(yīng)公式NASA World Wind Tile Structure NASA World Wind Map Tile System 瓦片請(qǐng)求調(diào)度策略 2.2.1 當(dāng)前顯示請(qǐng)求及其顯示（客戶端）2.2.2目標(biāo)瓦

2、片快速搜索算法2.2.3瓦片數(shù)據(jù)的請(qǐng)求預(yù)測(cè)（服務(wù)器）2.3 可視化的地球空間數(shù)學(xué)模型2.3.1 地球三維LOD模型3. Geospatial image processing 3.1步驟：（準(zhǔn)備知識(shí)：LOTD：0層瓦片大??；TileSize： ） 3.2 NLT Landsat處理技術(shù)范例 3.3 dstile howto4. Creat a new world 4.1World Wind本地文件說(shuō)明4.1.1. Worlds4.1.2. 在安裝目錄下4.2 創(chuàng)建新的世界模型5. Making layer 5.1 添加新圖層(更詳細(xì)可參考)5.2 Point5.3 Lines5.4 Polyg

3、on5.5 Model Feature5.6 Effects and Multi-Texturing5.7 Shapefile5.8 KML/KMZ6. ddons與pluginsaddons與plugins區(qū)別插件總匯及其功能7. script腳本8. 編譯調(diào)試9. 讀書籍和論文10. 1.4.1版本發(fā)展一、 WW軟件結(jié)構(gòu)1.1 功能介紹World Wind(簡(jiǎn)稱WW，中文有人直翻譯為世界風(fēng))，是NASA發(fā)布的一個(gè)開(kāi)放源代碼(Open Source)的地理科普軟件(由NASA Research開(kāi)發(fā)，由NASA Learning Technologies來(lái)發(fā)展)，它是一個(gè)可視化地球儀，將NAS

4、A、USGS以及其它WMS服務(wù)商提供的圖像通過(guò)一個(gè)三維的地球模型展現(xiàn)，近期還包含了月球、金星、火星、天文星系等的展現(xiàn)。用戶可在所觀察的行星上隨意地旋轉(zhuǎn)、放大、縮小, 同時(shí)可以看到地名和行政區(qū)劃.軟件還包含了一個(gè)軟件包，能夠?yàn)g覽地圖及其它由因特網(wǎng)上的OpenGIS Web Mapping Service提供的圖像。目前軟件所使用的數(shù)據(jù)：低分辯率的Blue marble數(shù)據(jù)現(xiàn)在包含的初始安裝內(nèi),當(dāng)用戶放大到特定區(qū)域時(shí)，附加的高分辯率數(shù)據(jù)將會(huì)自動(dòng)從NASA服務(wù)器上被下載。 A 靜態(tài)數(shù)據(jù)層:en:Blue marble 圖像 en:Landsat 圖像 en:USGS 圖像 en:SRTM 地形數(shù)據(jù)

5、B 動(dòng)畫數(shù)據(jù)層:en:Animated Earth en:MODIS en:GLOBE1.2 WW與XMLWorldWind是一個(gè)優(yōu)秀的客戶端框架引擎，XML實(shí)現(xiàn)數(shù)據(jù)描述和軟件設(shè)置，通過(guò)WW的WorldModel、圖層、插件、http和WMS請(qǐng)求、三維渲染等實(shí)現(xiàn)交互式瀏覽。也就是說(shuō)XML是數(shù)據(jù)。WorldWind則是數(shù)據(jù)執(zhí)行的引擎，WW可擴(kuò)展性和開(kāi)發(fā)性非常強(qiáng)大。1.3 Clinet與Server1.3.1 WW使用鏈接（附表1）1.3.2本地緩存路徑（附表1）1.3.3本地功能配置（附表1）1.3.4 http請(qǐng)求與aspBlue marble、:Landsat、SRTM使用直接的http請(qǐng)求

6、。Why not use WMS? NASA Learning Technologies and the NWW user community has tried in the pastto use WMS to serve World Wind tiles. The immense demand of the NWW user base has crippled any WMS server weve tried to date. Instead of using an overly complex format(for serving predefined chunks that is)

7、we opted to go with pre-rendered, pre-defined tiles that are just stored in the file system, and can be served to the client with no further processing server side. This reduces the strain on the server incredibly. If you believe you can setup a WMS server that would be able to withstand the NWW use

8、r base, feel free,wed be happy to see that. In the mean time we are working on creating a WMS server that can read the World Wind tile structure (for compatibility with other programs).1.3. WMS、WFSWMSweb Map server 網(wǎng)絡(luò)地圖服務(wù) WFSweb Fundation server 網(wǎng)絡(luò)要素服務(wù)USGS和所有動(dòng)畫數(shù)據(jù)層使用符合OGC協(xié)議的和實(shí)現(xiàn)數(shù)據(jù)網(wǎng)絡(luò)共享。二、 WorldWind金子塔體

9、系2.1 瓦片金字塔詳解及其對(duì)應(yīng)公式 2.1.1 NASA World Wind Tile Structure采取笛卡爾坐標(biāo)，原點(diǎn)（X = 0, Y = 0.）在投影坐標(biāo)左下方，即南極點(diǎn)（-90 -180 in latitude and longitude），如圖： 圖1 圖2如圖2，World Wind使用“Level Zero Tile Size”來(lái)決定每一個(gè)瓦片寬和高的大?。ㄋ械耐咂际钦叫?。標(biāo)準(zhǔn)的level zero tile size 還沒(méi)有制定，但是它必須滿足能被180整除。The level zero tile size (以后稱作 lzts)是層與層間轉(zhuǎn)換的最簡(jiǎn)單的距離。在

10、 NLTLandsat 7 中，lzts被默認(rèn)設(shè)置為2.25度?？捎靡韵鹿接?jì)算第N層的tile size：size=lzts2N次方。也就是說(shuō)下一層將上一層一分為四。：圖3 顯示坐標(biāo)軸在X、Y方向值的增量： 2.1.2 NASA World Wind Map Tile System坐標(biāo)與圖片編號(hào)對(duì)應(yīng)公式：（1）利用下面公式,求某點(diǎn)坐標(biāo)X,Y(緯,經(jīng)),在某層N的文件號(hào),圖片號(hào)(X+90)/值*2的N次方=文件夾號(hào)(取整數(shù)如256.65即為256)(Y+180)/值*2的N次方=圖片后綴(取整數(shù))（2）利用求得的圖片編號(hào)A_B反求這張圖片的左下角坐標(biāo)A*值/2的N次方-90=緯度 B*值/2的

11、N次方-180=經(jīng)度 這樣求得圖片起始分割點(diǎn)坐標(biāo),然后以值/2的N次的間隔分割圖片,把整幅圖分成256*256 or 512*512BMNG（500m）LandSatGeocover(30m/15m)OnEarth(15m)USGS(1m)SRTM(90m/3m0(US)36 36。2.252.25。2.02.0。0.80.8。2020。2.2 瓦片請(qǐng)求調(diào)度策略2.2.1 當(dāng)前顯示請(qǐng)求及其顯示（客戶端）(1)瓦片金字塔模型構(gòu)建金字塔是一種多分辨率層次模型。在地形場(chǎng)景繪制時(shí)，在保證顯示精度的前提下為提高顯示速度，不同區(qū)域通常需要不同分辨率的數(shù)字高程模型數(shù)據(jù)和紋理影像數(shù)據(jù)。數(shù)字高程模型金字塔和影像

12、金字塔則可以直接提供這些數(shù)據(jù)而無(wú)需進(jìn)行實(shí)時(shí)重采樣。盡管金字塔模型增加了數(shù)據(jù)的存儲(chǔ)空間，但能夠減少完成地形繪制所需的總機(jī)時(shí)。分塊的瓦片金塔模型還能夠進(jìn)一步減少數(shù)據(jù)訪問(wèn)量，提高系統(tǒng)的輸入輸出執(zhí)行效率，從而提升系統(tǒng)的整體性能。當(dāng)?shù)匦物@示窗口大小固定時(shí)，采用瓦片金字塔模型可以使數(shù)據(jù)訪問(wèn)量基本保持不變。瓦片金字塔模型的這一特性對(duì)海量地形實(shí)時(shí)可視化是非常重要的。在構(gòu)建地形金字塔時(shí)，首先把原始地形數(shù)據(jù)作為金字塔的底層，即第0 層，并對(duì)其進(jìn)行分塊，形成第0 層瓦片矩陣。在第0 層的基礎(chǔ)上，按每2 2個(gè)像素合成為一個(gè)像素的方法生成第1 層，并對(duì)其進(jìn)行分塊，形成第1 層瓦片矩陣。如此下去，構(gòu)成整個(gè)瓦片金字塔。以影

13、像為例，設(shè)第l 層的像素矩陣大小為irlicl，分辨率為resl，瓦片大小為isis，則瓦片矩陣的大小trl tcl 為：trl = irl/istcl = icl/is其中“ ”為向下取整符，下同。按每22 個(gè)像素合成為1 個(gè)像素后生成的第l+1 層的像素矩陣大小irl+1icl+1 為：irl+1 = irl/2icl+1 = icl/2其分辨率resl+1 為：resl+1 = resl 2不失一般性，我們規(guī)定像素合成從像素矩陣的左下角開(kāi)始，從左至右從下到上依次進(jìn)行。同時(shí)規(guī)定瓦片分塊也從左下角開(kāi)始，從左至右從下到上依次進(jìn)行。在上述規(guī)定的約束下，影像與其瓦片金字塔模型是互逆的。同時(shí)，影像的

14、瓦片金字塔模型也便于轉(zhuǎn)換成具有更明確拓?fù)潢P(guān)系的四叉樹(shù)結(jié)構(gòu)。(2)線性四叉樹(shù)瓦片索引（WW中應(yīng)用QST，即quadrat size tree）四叉樹(shù)是一種每個(gè)非葉子節(jié)點(diǎn)最多只有四個(gè)分支的樹(shù)型結(jié)構(gòu)，也是一種層次數(shù)據(jù)結(jié)構(gòu)，其特性是能夠?qū)崿F(xiàn)空間遞歸分解。圖2 是瓦片金字塔模型的四叉樹(shù)結(jié)構(gòu)示意圖，其中矩形符號(hào)代表葉子節(jié)點(diǎn)，圓形符號(hào)代表非葉子節(jié)點(diǎn)。本文采用四叉樹(shù)來(lái)構(gòu)建瓦片索引和管理瓦片數(shù)據(jù)。在瓦片金字塔基礎(chǔ)上構(gòu)建線性四叉樹(shù)瓦片索引分三步：即邏輯分塊、節(jié)點(diǎn)編碼和物理分塊。邏輯分塊與構(gòu)建瓦片金字塔對(duì)應(yīng)，規(guī)定塊劃分從地形數(shù)據(jù)左下角開(kāi)始，從左至右，從下到上依次進(jìn)行。同時(shí)規(guī)定四叉樹(shù)的層編碼與金字塔的層編碼保持一致，

15、即四叉樹(shù)的底層對(duì)應(yīng)金字塔的底層。設(shè)(ix,iy)為像素坐標(biāo)，is 為瓦片大小，io 為相鄰?fù)咂丿B度，以像素為單位；(tx,ty)為瓦片坐標(biāo)，以塊為單位；l為層號(hào)。若瓦片坐標(biāo)(tx,ty)已知，則瓦片左下角的像素坐標(biāo)(ixlb,iylb)為：ixlb = tx isiylb = ty is瓦片右上角的像素坐標(biāo)(ixrt,iyrt)為：ixrt = (tx + 1) is + io 1iyrt = (ty + 1) is + io 1如果像素坐標(biāo)(ix,iy)已知，則像素所屬瓦片的坐標(biāo)為：tx = ix/isty = iy/is由像素矩陣行數(shù)和列數(shù)以及瓦片大小，可以計(jì)算出瓦片矩陣的行數(shù)和列數(shù)，然

16、后按從左至右，從下到上的順序依次生成邏輯瓦片，邏輯瓦片由(ixlb,iylb), (ixrt,iyrt), (tx,ty), l)唯一標(biāo)識(shí)。節(jié)點(diǎn)編碼假定用一維數(shù)組來(lái)存儲(chǔ)瓦片索引，瓦片排序從底層開(kāi)始，按從左至右，從下到上的順序依次進(jìn)行，瓦片在數(shù)組中的偏移量即為節(jié)點(diǎn)編碼。為了提取瓦片(tx,ty,l)，必須計(jì)算出其偏移量。我們采用一個(gè)一維數(shù)組來(lái)存儲(chǔ)每層瓦片的起始偏移量，設(shè)為osl。若第l 層瓦片矩陣的列數(shù)為tcl，則瓦片(tx,ty,l)的偏移量offset 為：offset = ty tcl + tx + osl 物理分塊在邏輯分塊的基礎(chǔ)上對(duì)地形數(shù)據(jù)進(jìn)行物理分塊，生成地形數(shù)據(jù)子塊。對(duì)上邊界和右邊

17、界瓦片中的多余部分用無(wú)效像素值填充。物理分塊完畢，按瓦片編號(hào)順序存儲(chǔ)。(3) 瓦片拓?fù)潢P(guān)系瓦片拓?fù)潢P(guān)系包括同一層內(nèi)鄰接關(guān)系和上下層之間的雙親與孩子關(guān)系兩個(gè)方面2。鄰接關(guān)系分別為東(E)、西(W)、南(S)、北(N)四個(gè)鄰接瓦片，如圖3(a)所示；與下層四個(gè)孩子的關(guān)系分別為西南(SW)、東南(SE)、西北(NW)、東北(NE)四個(gè)孩子瓦片，如圖3(b)所示；與上層雙親的關(guān)系是一個(gè)雙親瓦片，如圖3(c)所示。若已知瓦片坐標(biāo)為(tx,ty,l)，則該瓦片相關(guān)的拓?fù)潢P(guān)系可表示為：1) 東、西、南、北四個(gè)鄰接瓦片的坐標(biāo)分別為：(tx+1,ty,l)、(tx-1,ty,l)、(tx,ty-1,l)、(tx

18、,ty+1,l)；2) 西南、東南、西北、東北四個(gè)孩子瓦片的坐標(biāo)分別為(2tx,2ty,l-1) 、(2tx+1,2ty,l-1) 、(2tx,2ty+1,l-1) 、(2tx+1,2ty+1,l-1)；3) 雙親瓦片的坐標(biāo)為(tx/2,ty/2,l+1)。 (a) 鄰接瓦片 (b) 孩子瓦片 (c) 雙親瓦片(4)地形場(chǎng)景數(shù)據(jù)的動(dòng)態(tài)管理可見(jiàn)區(qū)域及最佳地形分辨率計(jì)算根據(jù)視點(diǎn)的坐標(biāo)和視線的方向，就可以計(jì)算出視景體與地形平均水平面相交的平面區(qū)域范圍，即地形可見(jiàn)區(qū)域范圍。圖4 為地形可見(jiàn)區(qū)域示意圖，圖中XOY 為地形平均水平面，E 為視點(diǎn)，視線EM 與地形平均水平面的交點(diǎn)為M，視點(diǎn)在XOY 上的投影

19、為M0。視景體E_ABCD 與平面XOY的四個(gè)交點(diǎn)分別為A、B、C和D，則地形可見(jiàn)區(qū)域范圍即為四邊形ABCD。有關(guān)點(diǎn)A、B、C和D的坐標(biāo)的詳細(xì)計(jì)算方法，可參考文獻(xiàn)3。地形可見(jiàn)區(qū)域的表示若用于地形顯示的屏幕窗口寬和高分別為Xw 和Yw(以像素為單位)，視景體的水平視場(chǎng)角和垂直視場(chǎng)角分別為FovX 和FovY。設(shè)當(dāng)可見(jiàn)區(qū)域ABCD恰好充滿顯示窗口時(shí)，窗口水平方向和垂直方向平均每個(gè)像素所對(duì)應(yīng)的地面距離分別為Dx 和Dy，則有：Dx = EM tan(FovX / 2) 2.0 / XwDy = EM tan(FovY / 2) 2.0 / Yw即當(dāng)?shù)匦畏直媛蔢 方向低于Dx 或Y 方向低于Dy時(shí)，地

20、形顯示的精度將會(huì)降低；反之，將會(huì)產(chǎn)生不必要的數(shù)據(jù)冗余，影響繪制效率。也就是說(shuō)，此時(shí)的Dx和Dy即為理論上地形繪制所需的最佳分辨率。由于通常地形X 方向和Y 方向的采樣間距相同，因此我們?cè)趯?shí)際應(yīng)用中取Dx 和Dy 中的小值作為最佳地形分辨率，用作后續(xù)瓦片搜索的重要依據(jù)。需要指出的是，通常該最佳分辨率被用作窗口中心瓦片所對(duì)應(yīng)的分辨率，而窗口其它位置所對(duì)應(yīng)的瓦片分辨率則根據(jù)瓦片中心到視點(diǎn)的距離作適當(dāng)?shù)慕档驼{(diào)整，因?yàn)檫@不僅符合人眼的視覺(jué)規(guī)律，而且還可以減少用于地形繪制的三角形數(shù)量。 World Wind當(dāng)前視角高度值與圖層調(diào)度間關(guān)系探討How does altitude relate to a lev

21、el in a layer being displayed?在WW中，從N-1st 轉(zhuǎn)換到 Nth level的高度值是： ( Level0TileSize * 390km / 2N ) X（因子）例如對(duì)于NLT LS7圖層，它的L0TS=2.25deg。它的0 level從390*2.25 = 877km開(kāi)始可見(jiàn)，level 1在小于439km可見(jiàn)等；（但這只是從WW的行為分析得到而不是從源代碼）對(duì)于較小的高程（小于2000km），X因子近似為1。然而，高于2000km因子就減少(可能是arctan(x)/x，或者其他) 但是，為什么剛好390km？為什么非線性非常量的因子？肯能只有WW開(kāi)發(fā)

22、者才知！A Small Analyse:(假定屏幕默認(rèn)尺寸為wWidth=992px 和 wHeight=653px, 南北方向的視角范圍為FOV=45deg) 打開(kāi)OnEarth layer（L0TS=2.0deg）并飛行到高度Alt=195km, 這時(shí)level2的瓦片開(kāi)始顯示。那么圖片區(qū)域的南北范圍尺寸（(N-S)dimension- nsDim）是多少？ nsDim = Alt * 2*tan(FOV/2)，所以在這里nsDim = 195km * 2*tan(22.5deg) = 161.5km. 將它轉(zhuǎn)換為度數(shù)(1deg=111km)，得到nsDim = 1.454deg。每個(gè)le

23、vel2的瓦片覆蓋0.5 0.5 degree的區(qū)域，因此“照相機(jī)”（可視屏幕范圍）相應(yīng)有1.454/0.5 = 2.91個(gè)瓦片(南北方向上)，或者2.91*512 = 1489 個(gè)像素。數(shù)值1489/653=2.28可以稱之為線性重采樣因子（LOF- Linear Oversampling factor）；平方后，得到區(qū)域重采樣因子（AOF-real oversampling factor）。在我們這例子為5.20。這意味著什么？5.2個(gè)瓦片圖片的像素在屏幕上顯示為1個(gè)像素。在沒(méi)有進(jìn)入下一圖層時(shí)，這些參數(shù)能減少到什么程度？在例子中，下一層的開(kāi)端為 Alt=97.5km，即195km的一半。當(dāng)

24、我們達(dá)這里時(shí)，the level2 tiles減少到只有L.O.F.=1.14或者A.O.F.=1.30。Now lets generalize it:When are these LOF and AOF becoming big (giving evidence about vasting)? When the user has - reduced the WW windows height, - increased the FOV (aka zoomed out), - tilted significantly, - gone near the Pole, -.what else have

25、 I forgot? 在這些情況下，AOF可能變得非常巨大，例如值為幾十，是否意味著100-or-so megabyte of tiles are needed for rendering of about-1 Mpix picture. 還有很多問(wèn)題例如 (e.g. the When do you know when to stop? (ie. how do you know you are at 1:1) question) 2.2.2目標(biāo)瓦片快速搜索算法目標(biāo)瓦片搜索實(shí)際上伴隨著分辨率搜索。如果采用一種分辨率，則只涉及層搜索和區(qū)域搜索。對(duì)于采用視點(diǎn)相關(guān)LOD 技術(shù)的地形可視化系統(tǒng)，要求不同

26、的區(qū)域?qū)?yīng)不同分辨率的地形數(shù)據(jù)，因此在瓦片搜索過(guò)程中需要進(jìn)行分辨率測(cè)試，圖5 為地形瓦片分辨率需求示意圖，圖中E 為視點(diǎn)，白色矩形表示與視景體相交且需要參與地形繪制的瓦片數(shù)據(jù)。矩形越小表示瓦片分辨率越高。本文目標(biāo)瓦片搜索采用由粗到細(xì)、層層推進(jìn)的策略。首先判斷可見(jiàn)區(qū)域是否在地形數(shù)據(jù)所表示的范圍內(nèi)，測(cè)試與可見(jiàn)區(qū)域交疊的頂層瓦片是否滿足分辨率要求，如果不滿足，則繼續(xù)測(cè)試瓦片的分支，直到搜索出所有滿足分辨率要求的瓦片。完整的算法描述如下：第1 步: 判斷可見(jiàn)區(qū)域是否與地形數(shù)據(jù)交疊。如果交疊，轉(zhuǎn)第2 步；否則，算法結(jié)束。第2 步: 逐個(gè)判斷頂層瓦片是否滿足分辨率要求。如果滿足，將其插入到繪制鏈表；否則，

27、將其4 個(gè)分支瓦片插入到測(cè)試鏈表1。第3 步: 判斷測(cè)試鏈表1 是否為空。如果為空，算法結(jié)束；否則，轉(zhuǎn)第4 步。第4 步: 逐個(gè)判斷測(cè)試鏈表1 中的瓦片是否與可見(jiàn)區(qū)域交疊。如果交疊，將其插入到測(cè)試鏈表2；否則，將其舍棄。第5 步: 清空測(cè)試鏈表1，轉(zhuǎn)第6 步。第6 步: 逐個(gè)判斷測(cè)試鏈表2 中的瓦片是否滿足分辨率要求。如果滿足，將其插入到繪制鏈表；否則，將其4 個(gè)分支瓦片插入到測(cè)試鏈表1。第7 步: 清空測(cè)試鏈表2，轉(zhuǎn)第3 步。交疊測(cè)試時(shí)，只要瓦片四個(gè)角點(diǎn)的任何一點(diǎn)在可見(jiàn)區(qū)域內(nèi)，即表示該瓦片與可見(jiàn)區(qū)域交疊。對(duì)于分辨率測(cè)試，首先計(jì)算瓦片四個(gè)角點(diǎn)的目標(biāo)分辨率，如果四個(gè)角點(diǎn)的目標(biāo)分辨率數(shù)值都大于或等

28、于瓦片所處金字塔層的分辨率數(shù)值，則表示該瓦片滿足分辨率要求，否則需要測(cè)試其四個(gè)分支。如圖6 所示，假設(shè)被測(cè)試瓦片的分辨率為2，由于瓦片左上角點(diǎn)目標(biāo)分辨率數(shù)值小于2，則必須對(duì)瓦片進(jìn)行細(xì)分，然后再進(jìn)行交疊測(cè)試和分辨率測(cè)試。上述目標(biāo)瓦片搜索算法是一個(gè)遞歸過(guò)程。在進(jìn)行分辨率判斷時(shí)只需對(duì)瓦片的四個(gè)角點(diǎn)進(jìn)行判斷，因而效率很高。另外，從瓦片索引數(shù)組中提取瓦片的分支無(wú)需進(jìn)行空間位置判斷，只需計(jì)算出瓦片偏移量即可，這同樣能夠顯著地提高算法的性能。2.2.3瓦片數(shù)據(jù)的請(qǐng)求預(yù)測(cè)（服務(wù)器）(在高分辨率遙感衛(wèi)星應(yīng)用中，本文作者有遙感數(shù)據(jù)庫(kù)的詳細(xì)理論)(1) 預(yù)讀?。汉线m的數(shù)據(jù)緩沖機(jī)制可以提高系統(tǒng)的性能，海量地形實(shí)時(shí)可視

29、化也不例外。當(dāng)?shù)匦螆?chǎng)景連續(xù)運(yùn)動(dòng)時(shí)，如果應(yīng)用程序能夠提前預(yù)測(cè)系統(tǒng)即將需要的數(shù)據(jù)，并能提前從磁盤把這些數(shù)據(jù)讀到緩沖區(qū)，那么當(dāng)應(yīng)用程序真正需要這些數(shù)據(jù)的時(shí)候，便可以直接從緩沖區(qū)讀取，從而省去了從磁盤讀取這些數(shù)據(jù)的時(shí)間。由于在海量地形漫游時(shí)，數(shù)據(jù)預(yù)測(cè)調(diào)度和場(chǎng)景繪制需要同時(shí)進(jìn)行，因此，我們利用操作系統(tǒng)的多線程能力，將數(shù)據(jù)預(yù)測(cè)調(diào)度和場(chǎng)景繪制分別采用兩個(gè)獨(dú)立的線程來(lái)完成。數(shù)據(jù)緩沖通過(guò)預(yù)測(cè)瓦片請(qǐng)求來(lái)實(shí)現(xiàn)。我們將瓦片請(qǐng)求劃分為高、中、低三種優(yōu)先級(jí)。高優(yōu)先級(jí)瓦片請(qǐng)求被立即處理并參與地形繪制；中優(yōu)先級(jí)瓦片請(qǐng)求在沒(méi)有高優(yōu)先級(jí)瓦片請(qǐng)求的情況下被讀取到緩沖區(qū)，但不參與地形繪制，即所謂數(shù)據(jù)緩沖；低優(yōu)先級(jí)瓦片請(qǐng)求只是在系統(tǒng)“

30、空閑”時(shí)才被讀取到緩沖區(qū)，并且也不參與地形繪制，同樣是數(shù)據(jù)緩沖。瓦片請(qǐng)求優(yōu)先級(jí)的確定方法如下：為表述方便，假設(shè)連續(xù)的三個(gè)場(chǎng)景窗口分別為窗口1、窗口2 和窗口3。為了實(shí)現(xiàn)瓦片數(shù)據(jù)緩沖，系統(tǒng)在發(fā)送窗口1 的瓦片請(qǐng)求時(shí)應(yīng)附帶發(fā)送窗口2 和窗口3 的瓦片請(qǐng)求，只是瓦片請(qǐng)求的優(yōu)先級(jí)有所不同，這便是瓦片請(qǐng)求預(yù)測(cè)。顯然，窗口1 中的瓦片請(qǐng)求應(yīng)該設(shè)為高優(yōu)先級(jí)，而窗口2 中的瓦片請(qǐng)求除窗口1 的公共部分外被設(shè)為中優(yōu)先級(jí)，窗口3 中的瓦片請(qǐng)求除與窗口1 及窗口2 的公共部分外被設(shè)為低優(yōu)先級(jí)，如圖7 所示。這樣，當(dāng)場(chǎng)景運(yùn)動(dòng)到窗口2 時(shí)，系統(tǒng)在發(fā)送窗口2 的瓦片請(qǐng)求后，數(shù)據(jù)預(yù)測(cè)調(diào)度線程就可以直接從緩沖區(qū)提取數(shù)據(jù)參與地

31、形繪制，而無(wú)需即時(shí)從磁盤讀取數(shù)據(jù)。設(shè)場(chǎng)景運(yùn)動(dòng)速度為v 瓦片/秒，運(yùn)動(dòng)方向與瓦片x 方向的夾角為，刷新頻率為f 幀/秒，則中、低優(yōu)先級(jí)瓦片請(qǐng)求的擴(kuò)展幅度tx、ty 為：tx = v cos / fty = v sin / f例如，當(dāng)v = 33.94 瓦片/秒， = 45，f = 24 幀/秒時(shí)，中、低優(yōu)先級(jí)瓦片請(qǐng)求的擴(kuò)展幅度tx、ty 都為1 個(gè)瓦片，即在當(dāng)前場(chǎng)景窗口的基礎(chǔ)上沿運(yùn)動(dòng)方向的x、y 分解方向各擴(kuò)展1 個(gè)瓦片為中優(yōu)先級(jí)瓦片請(qǐng)求，再擴(kuò)展1 個(gè)瓦片為低優(yōu)先級(jí)瓦片請(qǐng)求。（2）異構(gòu)磁盤系統(tǒng)的瓦片分布方式（高分辨率遙感衛(wèi)星應(yīng)用P205）（3）應(yīng)用實(shí)例本文的研究成果已成功應(yīng)用于多個(gè)地區(qū)的海量地形

32、數(shù)據(jù)實(shí)時(shí)可視化。以某地區(qū)地形三維可視化為例，該地區(qū)所采用的規(guī)則格網(wǎng)數(shù)字高程模型共有2376923536 個(gè)高程點(diǎn)，分辨率為60 米，彩色紋理影像大小為4753647070，分辨率為30 米。數(shù)字高程模型數(shù)據(jù)為2.13GB，紋理影像數(shù)據(jù)為6.40GB。分別對(duì)數(shù)字高程模型數(shù)據(jù)和紋理影像數(shù)據(jù)構(gòu)建9層和10 層瓦片金字塔和相關(guān)索引文件，瓦片大小均為128128，相鄰數(shù)字高程模型瓦片間重疊為1，相鄰影像瓦片間沒(méi)有重疊，生成的數(shù)字高程模型瓦片金字塔數(shù)據(jù)為2.94GB，紋理影像瓦片金字塔數(shù)據(jù)為8.77GB。實(shí)驗(yàn)計(jì)算機(jī)配置：P4 CPU 1.5GHz，512MB 內(nèi)存，GeForce4 64MB 顯卡，操作系

33、統(tǒng)為Windows2000 SP4，開(kāi)發(fā)語(yǔ)言為Visual C+ 6.0，三維圖形標(biāo)準(zhǔn)為OpenGL。利用本文研究的算法實(shí)現(xiàn)了實(shí)驗(yàn)區(qū)真實(shí)紋理映射地形的三維實(shí)時(shí)可視化，并可實(shí)現(xiàn)地形的實(shí)時(shí)漫游及縮放旋轉(zhuǎn)等交互操作。當(dāng)顯示窗口為800600 像素，屏幕顯示誤差為1 個(gè)像素時(shí)，不采用瓦片請(qǐng)求預(yù)測(cè)時(shí)平均幀頻為18.3 幀/秒，采用瓦片請(qǐng)求預(yù)測(cè)后的平均幀頻為25.6 幀/秒，緩沖區(qū)中的瓦片總數(shù)維持在70 到120 之間；屏幕顯示誤差為2 個(gè)像素時(shí)，不采用瓦片請(qǐng)求預(yù)測(cè)時(shí)平均幀頻為23.2 幀/秒，采用瓦片請(qǐng)求預(yù)測(cè)后的平均幀頻為30.1 幀/秒，緩沖區(qū)中的瓦片總數(shù)維持在50 到90 之間。圖8 為部分實(shí)驗(yàn)結(jié)果

34、。 （4）結(jié)論實(shí)驗(yàn)表明，本文提出的瓦片金字塔模型和線性四叉樹(shù)索引相結(jié)合的數(shù)據(jù)管理模式，能夠滿足海量地形數(shù)據(jù)實(shí)時(shí)可視化的需要，并且在實(shí)現(xiàn)海量地形幾何數(shù)據(jù)實(shí)時(shí)繪制的同時(shí)完成了海量紋理數(shù)據(jù)的實(shí)時(shí)映射；通過(guò)對(duì)視景體可見(jiàn)區(qū)域外地形數(shù)據(jù)的裁剪和基于分辨率測(cè)試的目標(biāo)瓦片快速搜索算法，大大減少了地形繪制的數(shù)據(jù)量，提高了系統(tǒng)的執(zhí)行效率；采用基于高、中、低優(yōu)先級(jí)的地形瓦片請(qǐng)求預(yù)測(cè)方法，進(jìn)一步提高了三維地形交互漫游的速度。本文的研究成果已擴(kuò)展應(yīng)用于基于網(wǎng)絡(luò)的分布式海量地形數(shù)據(jù)三維可視化。本文研究所采用的地形數(shù)據(jù)為規(guī)則格網(wǎng)數(shù)據(jù)，今后將對(duì)用不規(guī)則三角網(wǎng)表示的海量地形數(shù)據(jù)可視化進(jìn)行進(jìn)一步研究。此外，在海量目標(biāo)數(shù)據(jù)實(shí)時(shí)三維

35、可視化、基于全球框架的空間數(shù)據(jù)可視化等方面，作者還將作進(jìn)一步的研究。1 David Luebke, Martin Reddy, Jonathan D. Cohen, et al. Level of Detail for 3D Graphics M. San Francisco: Morgan Kaufmann Publishers, 2002.2 鄧雪清. 柵格型空間數(shù)據(jù)服務(wù)體系結(jié)構(gòu)與算法研究D. 鄭州: 解放軍信息工程大學(xué), 2003.3 王源, 劉建永, 江南, 等. 視點(diǎn)相關(guān)實(shí)時(shí)LoD 地形模型動(dòng)態(tài)構(gòu)網(wǎng)算法J. 測(cè)繪學(xué)報(bào), 2003, 32(1): 47-52.4 王永君, 龔健雅. 用

36、于大范圍虛擬地形環(huán)境的動(dòng)態(tài)場(chǎng)景管理方法 J. 測(cè)繪學(xué)院學(xué)報(bào), 2003, 20(3): 231-234.2.3 可視化的地球空間數(shù)學(xué)模型一般來(lái)說(shuō)，遙感數(shù)據(jù)圖片都是使用在WGS84坐標(biāo)下的UTM投影。對(duì)于WorldWind來(lái)說(shuō)，其地球模型為圓形，半徑與WGS84坐標(biāo)的長(zhǎng)半徑相等（也可以另設(shè)），在實(shí)時(shí)顯示時(shí)，必須將平面的高程數(shù)據(jù)、紋理數(shù)據(jù)分塊并映射到給定的球面上。2.3.1 地球三維LOD模型 (1)地球模型的經(jīng)緯網(wǎng)分割由于地球橢球模型可以在球的基礎(chǔ)上沿自轉(zhuǎn)軸方向壓縮得到(利用三維圖形庫(kù)中的縮放函數(shù)很容易做到)，所以我們先將橢球模型簡(jiǎn)化為球面模型。球面可以看作由分布均勻的經(jīng)緯網(wǎng)組成，經(jīng)緯網(wǎng)越密，地

37、球模型的精度就越高。要實(shí)現(xiàn)視點(diǎn)相關(guān)球面LOD最簡(jiǎn)單有效的方法就是根據(jù)視點(diǎn)到球面的距離來(lái)動(dòng)態(tài)調(diào)整經(jīng)緯網(wǎng)的密度，這樣就在保證視覺(jué)效果的前提下實(shí)現(xiàn)球面模型簡(jiǎn)化，圖2-5-3顯示了在不同視點(diǎn)距離下，球面模型的簡(jiǎn)化程度。 (2)球面經(jīng)緯網(wǎng)格的可見(jiàn)性判斷在近地空間環(huán)境中漫游時(shí)，由于一個(gè)時(shí)刻我們所能看見(jiàn)的都是地球表面的一部分，因此，我們沒(méi)有必要繪制看不見(jiàn)的部分。本文利用視景體裁剪的方法來(lái)判斷經(jīng)緯網(wǎng)格是否可見(jiàn)。 (3)利用凹凸貼圖(Bump Map)表現(xiàn)地球表面起伏情況基于球面的三維地球繪制方法不能用幾何形狀來(lái)表達(dá)地面起伏，但我們可以采用凹凸貼圖的方法，在一定程度上彌補(bǔ)這個(gè)缺陷。在凹凸貼圖實(shí)現(xiàn)過(guò)程中，為每一個(gè)

38、頂點(diǎn)計(jì)算光線向量是一個(gè)非常費(fèi)時(shí)的工作，為了加快凹凸貼圖的速度，可以采用基于象素著色的硬件加速方法實(shí)現(xiàn)。（2）海量數(shù)據(jù)的全球地形LOD模型(SROAM)我們將ROAMZ6 ( Real-time Optimal Adaptive Meshes，實(shí)時(shí)優(yōu)化適應(yīng)性網(wǎng)格)算法引入到具有海量地形與影像數(shù)據(jù)的全球地形三維建模中。ROAM算法是目前應(yīng)用最廣的一種視點(diǎn)相關(guān)地形實(shí)時(shí)LOD算法。其基本思想是:在對(duì)地形進(jìn)行三維顯示時(shí)，依據(jù)視點(diǎn)的位置和視線的方向等多種因素，對(duì)于表示地形表面的三角形片元進(jìn)行一系列的基于二叉樹(shù)的分裂與合并，最終形成和原始表面近似且無(wú)縫無(wú)疊的簡(jiǎn)化連續(xù)網(wǎng)格。本文在ROAM算法的基礎(chǔ)上，提出連續(xù)

39、球面地形LOD算法SROAM ( Spherical ROAM簡(jiǎn)稱SROAM )，突破了經(jīng)典的ROAM算法只適合于平面地形數(shù)據(jù)的限制。SROAM算法基本思想 由于經(jīng)典的ROAM是從分裂一個(gè)正方形網(wǎng)格開(kāi)始不斷細(xì)分地形的，在SROAM中，我們采用正方體作為球面地形的頂層，如圖2-5-4(a)。正方體的6個(gè)面分別作為ROAM算法分裂的最頂層正方形。這些正方形的公共邊設(shè)置好拓?fù)潢P(guān)系，以保證不同ROAM面在分裂、合并操作時(shí)邊界上網(wǎng)格的連續(xù)性。 當(dāng)SROAM的頂層正方體分裂成12個(gè)三角形，并設(shè)置好三角形之間正確的拓?fù)潢P(guān)系后，ROAM的原理和方法就可以完全適用了。事實(shí)上，SROAM算法比經(jīng)典的ROAD更容易

40、實(shí)現(xiàn)，由于球面沒(méi)有邊界，因此不需要判斷地形網(wǎng)格的邊界條件，可以獲得更修的繪制速度。唯一不同的是，我們需要將平面的高程數(shù)據(jù)、紋理數(shù)據(jù)分塊并映射到給定的球面上。數(shù)據(jù)模型的簡(jiǎn)化272829 要建立全球的地LOD模型，其數(shù)據(jù)從覆蓋全球的分辨率為幾分或幾秒的DEM數(shù)拒和分辨率為幾百米的衛(wèi)星影像，到局部區(qū)域高分辨率的DEM數(shù)據(jù)和正射航空影像，數(shù)據(jù)量是巨大的。對(duì)于如此大的數(shù)據(jù)，要想在中低檔微機(jī)上達(dá)到實(shí)時(shí)交互的目的，必須對(duì)數(shù)據(jù)分級(jí)、分塊管理z7，如圖2-5-5。對(duì)于地形的分級(jí)，本文采用小波變換對(duì)地形進(jìn)行簡(jiǎn)化，并分塊構(gòu)建金字塔。 小波理論是本世紀(jì)八十年代后期發(fā)展起來(lái)的一種新的信號(hào)處理工具，在時(shí)空域和頻率域都具有

41、良好的局部化性質(zhì)，可以根據(jù)不同的尺度要求，由粗到細(xì)地觀察信號(hào)，被譽(yù)為信號(hào)分析的數(shù)學(xué)顯微鏡，在圖形圖像領(lǐng)域得到了廣泛的應(yīng)用，如:基于小波變換的圖象壓縮、圖象邊緣檢測(cè)、圖象拼接、立體視覺(jué)與匹配、曲面插值等，而基于小波變換的圖象壓縮是其目前最成功的應(yīng)用之一。由于規(guī)則格網(wǎng)實(shí)際上為一二維數(shù)組表示的高度值，完全可以將高度域看成圖像的灰度，將圖像處理的算法引入格網(wǎng)簡(jiǎn)化中。本文采用多進(jìn)制小波簡(jiǎn)化DEM數(shù)據(jù)。M進(jìn)制小波變換公式為: 其中，aj,k,l為小波的低頻成分，j為任意整數(shù)。k,1 = 0,1,2,.n0、n1、n2為正整數(shù)，C為常量。根據(jù)四進(jìn)制小波將全球數(shù)據(jù)簡(jiǎn)化為金字塔數(shù)據(jù)結(jié)構(gòu)29;然后，根據(jù)計(jì)算機(jī)硬件

42、的性能把每一層分為若千塊;最后，將組織好的數(shù)據(jù)寫入文件或數(shù)據(jù)庫(kù)。 近年來(lái)對(duì)于大數(shù)據(jù)量紋理數(shù)據(jù)的處理方法研究大紋理數(shù)據(jù)調(diào)度、映射的文獻(xiàn)并不少見(jiàn)。如基于原始紋理圖像分割的傳統(tǒng)方法30和MP-Grid方法32，都是采用預(yù)先把大紋理分割成小的紋理片，使每個(gè)小紋理片生成的Mipmap圖適應(yīng)硬件存儲(chǔ)，然后根據(jù)視野的覆蓋范圍調(diào)度相關(guān)紋理片的Mipmap圖組合生成要顯示的數(shù)據(jù)。韓慧健等人提出了使用Clipmap(裁剪Mipmap)實(shí)現(xiàn)大紋理映射31 ，他的思路是先為原始大紋理生成Mipmap圖存儲(chǔ)在硬盤，再分割Mipmap各層，調(diào)入各層與視野有關(guān)的分割片數(shù)據(jù)到有限內(nèi)存，相當(dāng)于把視域覆蓋的紋理部分裁剪下來(lái)，這樣

43、可以把整個(gè)紋理指定在單個(gè)坐標(biāo)系中?？傊?，解決大紋理貼圖的問(wèn)題目前都采用分級(jí)、分塊的方法，該方法可以較好的解決硬件的存儲(chǔ)和帶寬瓶頸問(wèn)題。本文也采用分級(jí)建立影象金字塔、在每一級(jí)上分塊隨視點(diǎn)動(dòng)態(tài)調(diào)度的方法實(shí)現(xiàn)。 高程數(shù)據(jù)球面映射 對(duì)于某一塊DEM數(shù)據(jù)來(lái)說(shuō)，通常表示為一個(gè)二維數(shù)組fDEMnRownCol,nRow和nCol分別表示地形在緯度、經(jīng)度方向的采樣點(diǎn)數(shù)。由于ROAM算法的三角形網(wǎng)格的分裂、合并操作可以看作是在二維平面上進(jìn)行，三角形頂點(diǎn)的高程用頂點(diǎn)的平面位置(x,y)可以很方便的在fDEM數(shù)組中查找得到。而球面LOD算法則不然，因?yàn)槿切畏至押喜⒉僮魇窃谌S球面上進(jìn)行的。假設(shè)要繪制的地球的中心在

44、原點(diǎn)，在某個(gè)三角形分裂過(guò)程中，網(wǎng)格中新生成的某一個(gè)頂點(diǎn)計(jì)算坐標(biāo)(即被分裂的三角形斜邊中點(diǎn)坐標(biāo))為P (Px，Py，Pz )，根據(jù)球面經(jīng)緯度投影的原理，向量P的方向決定了該點(diǎn)在球面上的經(jīng)緯度坐標(biāo)。因此，我們首先建立P (Px，Py，PZ)與經(jīng)緯度平面坐標(biāo)(，)之間的映射關(guān)系(如圖2-5-6 ) 根據(jù)球面幾何關(guān)系，經(jīng)緯度(，)與P的關(guān)系為： 其中，坐標(biāo)，弧度表示。再分別將,，規(guī)劃到0-1區(qū)間內(nèi)后，我們?nèi)菀椎玫浇?jīng)緯度與fDEM坐標(biāo)之間的關(guān)系。 綜合(2-5-2)(2-5-3 )式，得到位置P (x, y,z)與DEM高程數(shù)據(jù)fDEM的坐標(biāo)映射關(guān)系：最終我們得到P點(diǎn)對(duì)應(yīng)的的實(shí)際坐標(biāo)為:其中R為地球半徑

45、。|P|為向量P的長(zhǎng)度。地圖投影模型 當(dāng)前大地坐標(biāo)系是統(tǒng)一的, 能準(zhǔn)確、唯一地描述地球上任一點(diǎn)的位置, 并能動(dòng)態(tài)適應(yīng)歷史的、現(xiàn)代的越來(lái)越精密的各種參考橢球體, 它可作為地球數(shù)字產(chǎn)品和大型G IS 最為適當(dāng)?shù)淖鴺?biāo)系, 即采用(B, L ,H) 描述地理實(shí)體幾何特征點(diǎn)。這時(shí)二維場(chǎng)(B,L ) 是三維地理信息空間沿法線的一個(gè)二維投影, 可以用地圖投影把它展平到XO Y 平面上, 也可不展, 它本身為規(guī)則的三維橢球面。由于其數(shù)學(xué)規(guī)則, 采用二維的曲線族(B, L ) 來(lái)表示其上坐標(biāo)顯得簡(jiǎn)便得多, 人們也更易理解, 它是地理信息的水平坐標(biāo)系, 而H可視為垂直坐標(biāo)。大型G IS 作用空間由小范圍局部擴(kuò)展到

46、地學(xué)大范圍, 甚至全球范圍, 地球數(shù)字產(chǎn)品的覆蓋區(qū)域也遍布全球, 在這樣大的區(qū)域內(nèi)實(shí)現(xiàn)信息共享和互操作, 除統(tǒng)一采用上述橢球坐標(biāo)系外,“地圖投影”的標(biāo)準(zhǔn)化也是十分重要的。(詳細(xì)描述見(jiàn)文章大型GIS 與數(shù)字地球的空間數(shù)學(xué)基礎(chǔ)研究)誤差分析（后加）視見(jiàn)體的計(jì)算 在視點(diǎn)漫游時(shí)，ROAM等平面LOD算法通常要對(duì)視景體范圍內(nèi)所有網(wǎng)格點(diǎn)進(jìn)行動(dòng)態(tài)誤差計(jì)算。而球面LOD算法有其特殊性，因?yàn)?，?dāng)視點(diǎn)在球面上漫游時(shí)，地平線以下的網(wǎng)格是看不見(jiàn)的。如圖2-5-80 視點(diǎn)E距離地球高度為H(先不考慮地表起伏)，地球平均半徑為R, EC為視線方向，那么地平線到視點(diǎn)的距離為:對(duì)于在球面上三角形網(wǎng)格誤差計(jì)算，我們采用的是屏幕

47、誤差P，綜合考慮了網(wǎng)格靜態(tài)誤差、三角形到視點(diǎn)距離、視線方向等多個(gè)方面的因素。計(jì)算公式為: 其中，K為比例系數(shù)，s(c;)為要判斷的三角形的靜態(tài)誤差，D為三角形到視點(diǎn)距離。 采用上述動(dòng)態(tài)誤差計(jì)算公式可以代替視景體裁減，因?yàn)?，在地平面以下的三角形滿足D幾條件，動(dòng)態(tài)誤差p =0，繪制時(shí)，將被忽略。 關(guān)于在每一塊地形上SROAM的三角形分裂、合并以及繪制優(yōu)化方法，與ROAM26算法完全相同，本文不再重復(fù)。實(shí)驗(yàn)結(jié)果與分析 本節(jié)介紹了基于橢球的地球LOD建模，提出了考慮地形的地球LOD模型(SROAM )。為了驗(yàn)證這些方法的可行性與實(shí)際效果，我們?cè)贒ELL 4550(配置同前)的微機(jī)上實(shí)現(xiàn)了這些算法，軟件

48、平臺(tái)為WindowsNT/2000, Visual C十+6. 0和三維圖形標(biāo)準(zhǔn)OnenGL，結(jié)果如下:基于橢球的地球LOD建模結(jié)果 圖2-5-9為視點(diǎn)到地球不同距離時(shí)的地球網(wǎng)格，從圖中可以看出，當(dāng)?shù)厍蜻h(yuǎn)離視點(diǎn)的時(shí)候，網(wǎng)格較稀疏，表示精度變低:當(dāng)視點(diǎn)靠近地球的時(shí)候，網(wǎng)格變密，表示精度提局。為了研究使用凹凸貼圖表現(xiàn)地面起伏實(shí)際效果，我們分別做了不進(jìn)行凹凸貼圖的地球繪制和進(jìn)行凹凸貼圖的地球繪制兩個(gè)試驗(yàn)，結(jié)果如圖2-5-10所示。很明顯，使用凹凸紋理后，地面質(zhì)感明顯增強(qiáng)(如圖2-5-10b )。地球SROAM建模結(jié)果 圖2-5-11為SROAM算法繪制帶有地形的地球表面的結(jié)果，左圖為網(wǎng)格顯示，右圖為

49、填充效果。在該圖中，視點(diǎn)在P位置，視線方向如圖箭頭方向，扇形區(qū)域?yàn)榭梢?jiàn)范圍。從圖中可以看出，在球面上視點(diǎn)區(qū)域模型繪制時(shí)具有最高的細(xì)節(jié)層次，而離視點(diǎn)較遠(yuǎn)時(shí)細(xì)節(jié)層次較低。SROAM算法平均每幀繪制20, 000個(gè)三角形，速度平均30fps，基本可以滿足實(shí)時(shí)系統(tǒng)的需要。由于受數(shù)據(jù)來(lái)源的限制，該算法對(duì)大數(shù)據(jù)量的DEM的實(shí)際效果還未進(jìn)行驗(yàn)證。三、 Geospatial image processing3.1步驟：（準(zhǔn)備知識(shí)：LOTD：0層瓦片大??；TileSize： ）1. 獲得帶有地理參考的geotiff 并將他重新投影到WGS84坐標(biāo)2. 找出geotiff的邊界3. 從geotiff左下角坐標(biāo)開(kāi)始

50、，把整個(gè)的經(jīng)緯度范圍切割成最接近LOTD的倍數(shù)。4. 在#3處理的的geotiff在左下角切下一塊L0TD大小的圖塊5. 將#4的圖塊重采樣為瓦片（512512大小或其他，），轉(zhuǎn)換圖像格式6. 將瓦片放到相應(yīng)的文件夾下，并以WW的命名規(guī)則命名瓦片的名字7. 重復(fù)切下L0TD大小的圖塊直到將TIFF右上角邊界8. 將原圖片以N=（L0TD2）大小繼續(xù)切割，從#3開(kāi)始重復(fù)9. 直到變成1：1比例，即切割大小N=(小于或等于)；3.2 NLT Landsat處理技術(shù)范例(這是NASA生產(chǎn)第一個(gè)版本Landsat 7 Mosaic（可視波段）的方法。- By Chris Maxwell (Lead D

51、eveloper, NASA World Wind) 首先從from the Global Land Coverage Facility下載Landsat7的全球衛(wèi)片Remarks: The scenes come as a set of grayscale GeoTiff files for each row/path/date scene 然后使用波段1 (Blue),2 (Green),3 (Red)合成RGB的GeoTiff image。 Remarks: 由于我不不想編程實(shí)現(xiàn)，我使用了PixelSense，但是PixelSense 提供了一個(gè) linear clip contrast

52、 stretch，因此必須 replace the image data after PixelSense outputted the RGB GeoTiff。再對(duì)RGB Landsat GeoTiff使用顏色增強(qiáng)算法Remarks: The basic idea for the algorithm I used is to divide the image into a grid with each grid square being n by n pixels, where n is about 25 for a 30meter/pixel Landsat scene. For each

53、square, find the least-squares for both the high and low, and then save this information in order to later apply a linear stretch of the pixels in that square, weighting the R, G, B bands individually to the desired tastes. I found that generally, the source Landsat scenes have too much blue and red, and I weight the stretch to compensate for this for a more pleasing image. Next, use a cubic convolution to smooth out the stretch information in the grid squares. This is only app