DEM高程數據
**
DEM高程數據包括兩個部分:ASTER GDEM30米分辨率高程數據和SRTM90米分辨率高程數據。ASTER GDEM數據來源於NASA,數據覆蓋範圍爲北緯83°到南緯83°之間的所有陸地區域,時間範圍爲2000年前後;SRTM數據來源於CIAT,數據覆蓋範圍爲北緯60°至南緯60°之間的所有陸地區域,時間範圍爲2000年前後。
ASTER GDEM 30米分辨率高程數據
本數據集利用ASTER GDEM第一版本(V1)的數據進行加工得來,是全球空間分辨率爲30米的數字高程數據產品。由於雲覆蓋,邊界堆疊產生的直線,坑,隆起,大壩或其他異常等的影響,ASTER GDEM第一版本原始數據局部地區存在異常,所以由ASTER GDEMV1加工的數字高程數據產品存在個別區域的數據異常現象,可以和全球90米分辨率數字高程數據產品互相補充使用。ASTER GDEM數據採用UTM/WGS84投影,數據格式爲IMG柵格影像,數據的值域範圍爲-152-8806米之間,比例尺爲1:25萬,其垂直精度20米,水平精度30米。
數據命名規則:ASTER GDEM基本的單元按1度X1度分片。每個GDEM數據包有兩個文件,一個數據高程文件和一個質量評估(QA)文件。每個文件的命名是根據影像幾何中心左下角的經緯度產生。例如,ASTGTM_N29E091代表左下角座標是北緯29度,東經91度。ASTGTM_N29E091_dem和ASTGTM_N29E091_num對應的分別是高程數據和質量控制數據。
SRTM 90米分辨率高程數據
SRTM(ShuttleRadarTopographyMi
SRTM數據採用WGS84橢球投影,使用16位的數值表示高程數值的( -32767米),空數據用-32726表示。數據格式爲IMG柵格影像,數據的值域範圍爲-12000-9000米之間,水平精度20米,高程精度16米。
SRTM的數據組織方式爲:每5度經緯度方格劃分一個文件,共分爲24行(-60至60度)和72列(-180至180度)。文件命名規則爲srtm_XX_YY.zip,XX表示列數(01-72),YY表示行數(01-24)。示意圖如下:
高程數據處理方法
第一次使用DEM高程數據的朋友常常遇到這個問題,IMG是壓縮包麼?怎麼不能解壓呢?爲什麼我打開之後數據是灰色的呢?明明是平原地區,爲什麼顯示的高程範圍卻在-32767-32767之間呢?爲什麼展示圖裏是五顏六色的,而我打開的卻是灰色影像呢?
首先IMG不是壓縮包,“.img”作爲一種柵格影像格式,可以直接在ArcMap、ENVI、ERDAS等遙感軟件中打開使用,無需解壓。
其次,怎麼去除高程影像中的空值(如 -32767),讓它在一個正常的範圍內顯示呢?小編這裏以TIF格式的DEM高程影像爲例(IMG的處理方式同樣),一步步帶大家來操作。
1. 在ARCMAP裏打開一幅DEM高程數據(ADD DATA),可以從左邊看到其顯示的數據範圍是-32767-32726,右側爲灰色影像。
2. 在ArcMap裏打開Spatial analyze工具,選擇Raster Calculator,設置DEM高程數據值爲0並進行計算(點擊Evaluate按鈕),頁面如下:
3. 計算之後結果顯示爲Calculation,頁面如下圖:
4. 再次打開Spatial analyze工具,選擇Reclassify,選中chang missing values to Nodata打鉤, 點擊OK:
5. 然後打開Raster Calculator面板,將上一步計算出來的[Reclass of Calculation]與[原始的DEM 高程數據]相乘, 點擊Evaluate,界面如下:
6. 計算結果,如圖:
7. 鼠標右鍵點擊Calculation2,打開Properties, 選擇Symbology選項,在左側選中Stretched, 點擊確定。如下圖:
8. 在這裏,我們看到,經過處理,高程影像的值域範圍已經變爲正常的海拔範圍(105-4336)了:
9. 我們還可以給它設一個漂亮的顯示顏色,簡單的雙擊左側的灰度條即可,到這一步,DEM數據的預處理就做完了。
DEM高程數據的應用範圍
DEM數據包含了豐富的地形、地貌、水文信息,它能夠直觀展現一個地區的地形、地貌特徵,通過DEM可以提取大量的地表形態信息,如坡度、坡向、水系等等,廣泛地應用於各種地形特徵的定量分析和三維地圖的繪製中。