画小范围高精度的地形图（SRTM和ASTERDEM）

Posted 2020-10-04 淘淘麻麻

参考seisman http://seisman.info/convert-geotiff-to-gmt.html

目的
画**地区地形图
特点小范围高精度
试验数据两种
（一）SRTM 3弧秒 90米

1 从网页下载
http://srtm.csi.cgiar.org/SRT-ZIP/SRTM_V41/SRTM_Data_GeoTiff/srtm_61_07.zip

http://srtm.csi.cgiar.org/SRT-ZIP/SRTM_V41/SRTM_Data_GeoTiff/srtm_61_07.zip 

2 GDAL grd数据拼接
gdal_merge.py -init 255 -o out.tif srtm_60_07.tif srtm_61_07.tif 
%输出的名称必须是out.tif

3 GDAL格式转换 tif2grd
gdal_translate -of GMT out.tif ningde.grd

4 gmt画图

 （二）ASTER GDEM 1弧秒 30米

ubuntu安装sudo apt-get install gdal-bin

gdal_translate -of GMT ASTGTM2_N28E121_dem .tif ASTGTM2_N28E121_dem.grd

gmt画图代码

   

 

 

参考seisman http://seisman.info/convert-geotiff-to-gmt.html

下面的命令可以将 GeoTiff 格式的数据转换为 GMT 可识别的 netCDF 格式:

gdal_translate -of GMT srtm_56_05.tif srtm_56_05.nc

gdal_translate 的语法如下:

Usage: gdal_translate [--help-general] [--long-usage]

[-ot {Byte/Int16/UInt16/UInt32/Int32/Float32/Float64/

CInt16/CInt32/CFloat32/CFloat64}] [-strict]

[-of format] [-b band] [-mask band] [-expand {gray|rgb|rgba}]

[-outsize xsize[%] ysize[%]]

[-unscale] [-scale[_bn] [src_min src_max [dst_min dst_max]]]* [-exponent[_bn] exp_val]*

[-srcwin xoff yoff xsize ysize] [-projwin ulx uly lrx lry] [-epo] [-eco]

[-a_srs srs_def] [-a_ullr ulx uly lrx lry] [-a_nodata value]

[-gcp pixel line easting northing [elevation]]*

[-mo"META-TAG=VALUE"]* [-q] [-sds]

[-co"NAME=VALUE"]* [-stats] [-norat]

src_dataset dst_dataset

-of 后接要转换的数据格式，使用 gdal_translate --long-usage 可以看到，支持的格式如下:

VRT: Virtual Raster

GTiff: GeoTIFF

NITF: National Imagery Transmission Format

HFA: Erdas Imagine Images (.img)

ELAS: ELAS

AAIGrid: Arc/Info ASCII Grid

GSAG: Golden Software ASCII Grid (.grd)

GSBG: Golden Software Binary Grid (.grd)

JPEG: JPEG JFIF

BT: VTP .bt (Binary Terrain) 1.3 Format

IDA: Image Data and Analysis

USGSDEM: USGS Optional ASCII DEM (and CDED)

ADRG: ARC Digitized Raster Graphics

BLX: Magellan topo (.blx)

Rasterlite: Rasterlite

SAGA: SAGA GIS Binary Grid (.sdat)