全球数字高程模型GTOPO30 下载

GTOPO30是USGS公布的全球数字高程模型，其采样间隔为30弧秒。


下载


使用wget批量下载所有数据文件:
wget -r -c -np -nd ftp://edcftp.cr.usgs.gov/data/gtopo30/global/




解压


for i in *.tar.gz;
do
    tar -zxvf $i;
done




说明


在README里有详细说明，重点归纳如下

• 
  将全球分为33个区块，如下图所示，北纬90度到南纬60度分为3份，每份纬向50度，经向40度，共计3*9=27个区块；南极附近为南纬60度到90度，经向60度，计6个区块，共计33个区块。
  
• 
  采样间隔30弧秒，即0.008333333333333度，全球分为21600行、43200列。
• 
  坐标参考为WGS84；
• 
  其值为垂直方向高出海平面的距离，即海拔，单位为m；
• 
  海拔范围-407到8752，这里不包含海洋深度信息，负值应该是大陆架的海拔；
• 
  海洋处海拔标记为-9999，海岸线以为大陆海拔至少为1；
• 
  小于1平方千米的岛不考虑；
• 
  文件以该区块的左上角坐标为文件名；
• 
  文件后缀：DEM为数据高程模型；其他的文件都是一些辅助信息；
• 
  纬度50度*经度40度文件大小：57600000 bytes
• 
  纬度30度*经度60度文件大小：28800000 bytes

复制


将所有DEM文件拷贝至$GMTHOME/share/dbase


修改grdraster.info


# GTOPO30 N90-N40
11 "GTOPO30 W180N90"            "m"     -R-180/-140/40/90       -I0.5m          PG i 1           0    -9999      W180N90.DEM B
12 "GTOPO30 W140N90"            "m"     -R-140/-100/40/90       -I0.5m          PG i 1           0    -9999      W140N90.DEM B
13 "GTOPO30 W100N90"            "m"     -R-100/-60/40/90        -I0.5m          PG i 1           0    -9999      W100N90.DEM B
14 "GTOPO30 W060N90"            "m"     -R-60/-20/40/90         -I0.5m          PG i 1           0    -9999      W060N90.DEM B
15 "GTOPO30 W020N90"            "m"     -R-20/20/40/90          -I0.5m          PG i 1           0    -9999      W020N90.DEM B
16 "GTOPO30 E020N90"            "m"     -R20/60/40/90           -I0.5m          PG i 1           0    -9999      E020N90.DEM B
17 "GTOPO30 E060N90"            "m"     -R60/100/40/90          -I0.5m          PG i 1           0    -9999      E060N90.DEM B
18 "GTOPO30 E100N90"            "m"     -R100/140/40/90         -I0.5m          PG i 1           0    -9999      E100N90.DEM B
19 "GTOPO30 E140N90"            "m"     -R140/180/40/90         -I0.5m          PG i 1           0    -9999      E140N90.DEM B

# GTOPO30 N40-S10
20 "GTOPO30 W180N40"            "m"     -R-180/-140/-10/40      -I0.5m          PG i 1           0    -9999      W180N40.DEM B
21 "GTOPO30 W140N40"            "m"     -R-140/-100/-10/40      -I0.5m          PG i 1           0    -9999      W140N40.DEM B
22 "GTOPO30 W100N40"            "m"     -R-100/-60/-10/40       -I0.5m          PG i 1           0    -9999      W100N40.DEM B
23 "GTOPO30 W060N40"            "m"     -R-60/-20/-10/40        -I0.5m          PG i 1           0    -9999      W060N40.DEM B
24 "GTOPO30 W020N40"            "m"     -R-20/20/-10/40         -I0.5m          PG i 1           0    -9999      W020N40.DEM B
25 "GTOPO30 E020N40"            "m"     -R20/60/-10/40          -I0.5m          PG i 1           0    -9999      E020N40.DEM B
26 "GTOPO30 E060N40"            "m"     -R60/100/-10/40         -I0.5m          PG i 1           0    -9999      E060N40.DEM B
27 "GTOPO30 E100N40"            "m"     -R100/140/-10/40        -I0.5m          PG i 1           0    -9999      E100N40.DEM B
28 "GTOPO30 E140N40"            "m"     -R140/180/-10/40        -I0.5m          PG i 1           0    -9999      E140N40.DEM B

# GTOPO30 S10-S60
29 "GTOPO30 W180S10"            "m"     -R-180/-140/-60/-10     -I0.5m          PG i 1           0    -9999      W180S10.DEM B
30 "GTOPO30 W140S10"            "m"     -R-140/-100/-60/-10     -I0.5m          PG i 1           0    -9999      W140S10.DEM B
31 "GTOPO30 W100S10"            "m"     -R-100/-60/-60/-10      -I0.5m          PG i 1           0    -9999      W100S10.DEM B
32 "GTOPO30 W060S10"            "m"     -R-60/-20/-60/-10       -I0.5m          PG i 1           0    -9999      W060S10.DEM B
33 "GTOPO30 W020S10"            "m"     -R-20/20/-60/-10        -I0.5m          PG i 1           0    -9999      W020S10.DEM B
34 "GTOPO30 E020S10"            "m"     -R20/60/-60/-10         -I0.5m          PG i 1           0    -9999      E020S10.DEM B
35 "GTOPO30 E060S10"            "m"     -R60/100/-60/-10        -I0.5m          PG i 1           0    -9999      E060S10.DEM B
36 "GTOPO30 E100S10"            "m"     -R100/140/-60/-10       -I0.5m          PG i 1           0    -9999      E100S10.DEM B
37 "GTOPO30 E140S10"            "m"     -R140/180/-60/-10       -I0.5m          PG i 1           0    -9999      E140S10.DEM B

# GTOPO30 S60-S90
38 "GTOPO30 W180S60"            "m"     -R-180/-120/-90/-60     -I0.5m          PG i 1           0    -9999      W180S60.DEM B
39 "GTOPO30 W120S60"            "m"     -R-120/-60/-90/-60      -I0.5m          PG i 1           0    -9999      W120S60.DEM B
40 "GTOPO30 W060S60"            "m"     -R-60/0/-90/-60         -I0.5m          PG i 1           0    -9999      W060S60.DEM B
41 "GTOPO30 W000S60"            "m"     -R0/60/-90/-60          -I0.5m          PG i 1           0    -9999      W000S60.DEM B
42 "GTOPO30 E060S60"            "m"     -R60/120/-90/-60        -I0.5m          PG i 1           0    -9999      E060S60.DEM B
43 "GTOPO30 E120S60"            "m"     -R120/180/-90/-60       -I0.5m          PG i 1           0    -9999      E120S60.DEM B




其他说明

• 这些DEM文件为Pixel Registration格式
• 数据值存储为有符号16位整型
• 这些数据都是Big-endian的

绘图测试

1  2  3  4  5  6  7  8  9 10 11 12 13

#!/bin/bash ID=11 Range=-180/-140/40/90 B=a10g10 Inc=0.5m PS=${ID}.ps grdraster $ID -R$Range -I$Inc -Gout.grd makecpt -Ctopo -T1/8000/1000 -Z > colors.cpt grdimage out.grd -B$B -R$Range -JQ15c -Ccolors.cpt > $PS rm out.grd colors.cpt

这里只绘制一个第一个区块的地形，若需要同时绘制多个区块，需要多次调用grdraster从不同的网格文件中提取网格，然后用grdpaste粘贴成一个网格再绘图。


南极区域


在南极区域还有一个特殊的高程数据ANTARCPS.DEM，其包含了S60-S90区域的六个区块的数据，但是是经过极坐标方位立体投影的。因而横坐标范围为（-2700，2700），纵坐标范围为（-2700，2700）。这个处理起来要麻烦一点。