ETOPO

download ETOPO files from http://www.ngdc.noaa.gov/mgg/global/

根据给出的两个版本，下载其中的grid-registered版本的binary文件

etopo1_ice_g_i2.zip

解压后得到bin文件和hdr文件，hdr决定grdraster.info如何写。

将.bin文件copy到GMT/share/dbase下

修改grdraster.info文件加入(/GMT/share/dbase)

8 "ETOPO1 Ice Surface"  "m"  -R-180/180/-90/90  -I1m  GG i 1  0 -32768 etopo1_ice_g_i2.bin L

如果下载的是netCDF格式的网格文件，需要利用如下命令将数据转换为binary格式：

grdreformat ETOPO1_Ice_g_gmt4.grd etopo1_ice_g_i2.bin=bs -N -V

后面可以依次加入更多的其他数据，如DEM数据。

10 "GTOPO30 W180N90" "m" -R-180/-140/40/90 -I0.5m PG i 1 0 -9999 W180N90.DEM B

直接用

grdraster 9 -Rg -I1m -Goutput.grd -V
makecpt -Cglobe -T-5000/8000/1000 -Z -V > color.cpt
grdimage output.grd -Ba60g30 -Rg -Yc -Xc -JN0/25c -Ccolor.cpt -K -V > etopo1.ps
pssacle -Ba2500f500::/:"m": -Ccolor.cpt -D12.5c/-2c/15c/0.35ch -O -V >>etopo1.ps
rm output.grd color.cpt

GMT地形数据有三种etopo、GTOPO30、SRTM，精度由低到高依次是

etopo5: 5 arc minutes

etopo2:2 arc minutes

etopo1:1 arc minutes

GTOPO30:30 arc second

SRTM30_PLUS: 30arc second

SRTM3:Shuttle Radar Topography Mission进一步处理为DEM数据,精度3arc second

ASTER GREM:

etopo数据包含了全球地形和海水深度数据，GTOPO30包含了全球地形数据，SRTM30_PLUS，网格精度1km，包括全球地形和海洋深度数据，大部分区域使用SRTM30中的全球地形数据，高纬度区域使用GTOPO30数据（无SRTM数据），南纬81度到北纬81度的海底地形数据来自于Smith和Sandwell的全球1分数据等;SRTM3数据包含了全球大部分区域的地形数据。ASTER GDEM数据覆盖范围为北纬83°到南纬83°之间的所有陆地区域，达到了地球陆地表面99%。

对于全球尺度或者大尺度地形图来说，可以使用etopo1数据；对于3度*3度-10度*10度的区域，可以使用GTOPO30数据；对于3度*3度以下的区域，可以使用SRTM3数据。

GTOPO30

解压复制：将所有DEM文件拷贝至/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

将全球分为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文件为Pixel Registration格式

数据值存储为有符号16位整型

这些数据都是Big-endian的

/转：http://blog.sciencenet.cn/blog-2637373-974284.html

STRM3:

SRTM，全称为 Shuttle Radar Topography Mission，该项目获取了北纬 60 度至南纬 60 度之间的雷达影像数据，进一步处理制成了数字地形高程模型（DEM），即现在的 SRTM 地形产品数据。
SRTM3 精度为 3 弧秒，即 90m 一个点，包括非洲、北美、南美、欧亚、澳大利亚以及部分岛屿；
SRTM1 精度为 1 弧秒，即 30m 一个点，仅限美国地区；

申请数据：https://harvestchoice.wufoo.com/forms/download-cgiarcsi-srtm/

wget -c -r -np --user=xxx --password=xxxxx http://data.cgiar-csi.org/srtm/tiles/ＧeoTIFF 下载

Geotiff 格式的数据需要使用 GDAL 提供的命令 gdal_translate 转换为 netCDF 格式
GeoTiff 格式转 netCDF 格式

sudo apt-get install gdal-bin
gdal_translate -of GMT srtm_56_05.tif srtm_56_05.grd

数据合并：有时需要将多个数据合并成更大区域的数据。常见的做法是先将每个 GeoTiff 文件转换成 netCDF 格式，然后利用 grdcut 和 grdpaste 进行裁剪与合并。实际操作时，合并总会遇到一些问题，可能是由于数据的转换过程中存在问题导致的。推荐的做法是先将多个 GeoTiff 文件合并起来，然将合并后的 GeoTiff 文件转换成 netCDF 格式。合并 GeoTiff 文件可以gdal_merge.py，CentOS 下要安装 gdal-python ；Ubuntu 下要安装 python-gdal 或 python3-gdal。转换命令如下:

gdal_merge.py -init 255 -o out.tif srtm_56_07.tif srtm_57_07.tif

高精度SRTM

部分网站将数据分割为 5 度 *5度单元，其相对1度*1 度单元使用起来免去了 cut 和 paste 的麻烦。在经度方向上，从 - 180 度到 180 度，划分为 72 个区块；在纬度方向上，从 60 度到 - 60 度，划分为 24 个区块，可以参考这里的图示，其文件名格式为 srtm_经度区块数_纬度区块数. zip。

或者也可以这样计算:

x = int((lon+180.0)/5.0)+1;
y = int((60-lat)/5.0)+1;

ETOPO

ETOPO

results matching ""

No results matching ""