CesiumLab 地形數據處理

最近接連有用戶反應地形數據處理的各種問題,我也是各種測試,想想還是整理一個文檔徹底說明一下。

地形柵格數據格式,一般是tif ,也有dem或者img,但是我個人強烈建議使用tif格式,因爲cesiumlab都是用這個格式測試的。
我們先說一下地形柵格數據的幾個關鍵參數:

強烈建議處理之前,先用arcmap打開,確保以下參數正常
1,空間參考(spatial reference)

spatial reference 定義了 柵格數據的座標空間,下圖爲正常參數,如果arcmap 這裏顯示undefined,那表示空間參考未知,必然無法正常處理。

CesiumLab 地形數據處理
柵格的空間參考信息
如果tif中不含空間參考,但是空間參考數據處理人員又明確知道,那麼可以在cesiumlab中指定。反過來說如果不確定,那就不要亂嘗試,只會浪費處理時間。
CesiumLab 地形數據處理
cesiumlab中指定空間索引
2,無數據值(nodata)
nodata定義了刪格中的無效數據數據值,也就是說如果柵格里某個像素的值 = nodata,意味着該像素沒有數據。

對於nodata的處理,cesiumlab把該位置使用高程0填充。
但是原始數據中這個值很亂,完全是用戶設置的,柵格里可能沒有這個nodata信息,也可能爲任意值。一般用戶會設置一個超大的負數值,比如 -10000  ,-99999等。這裏強烈建議,不要設置數據極值,例如-1.797693e+308 ,這種極值會導致中間計算出現非數字,導致地形處理異常。

CesiumLab 地形數據處理
無數據值
如果原始柵格中沒有nodata值,只有一種情況處理地形高程會正常(數據的空間參考是wgs84,地形數據邊界爲矩形),除此之外都必須設置此值。在cesiumlab選中對應tif,如下設置。
CesiumLab 地形數據處理
這個-9999.005859是怎麼來的? 我們還是必須通過arcmap打開tif,使用indentify工具,點擊我們明確知道是無數據區域,那麼獲得他的值。
CesiumLab 地形數據處理
arcgis里人工判定nodata
3,地理範圍(extent)

柵格數據實際上是一張圖片,必須要一個轉換,把像素座標轉到地理座標,這就是extent,在gdal裏是geotransform。在arcmap裏如下圖所示
CesiumLab 地形數據處理
柵格的地理範圍
正常這裏應該顯示的地理座標,如果你發現這個extent的四個值裏有0,那麼基本意味着,你的tif的地理範圍信息缺失。

tif 一般把這個信息存儲在文件內,還有另一種方式存儲爲 同名的tfw 文件。
所以我們處理或者拷貝數據的時候,如果有tfw,請一併拷貝。

CesiumLab 地形數據處理
tfw文件
確保上面三個關鍵參數正確,我們再分析cesiumlab處理地形的異常情況原因
1, 瞬間結束      空間參考缺失  或者  地理範圍缺失
2,處理完成之後半邊黑球      nodata 過大
3,地形高程異常      nodata 缺失或者不正確
4,位置不正確或伴隨高程不正常      一般是空間參考不正確根據異常原因,請在arcmap使用相關工具處理你的柵格數據。

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章