gdal_translate用法

1、geotiff转换为GMT所需的grd文件
gdal_translate -of GMT -a_nodata NAN input.tif output.tif
2、有时候geotiff不止一个波段，需要转成单波段的geotiff才能进一步转化为GMT的grd文件
gdal_translate -b 1 input.tif output.tif
3、更改投影系统，比如把横轴墨卡托投影转换为带有经纬度的WGS84坐标系
gdalwrap -s_srs def -t_srs def input.tif output.tif
s_srs是源投影系统，t_srs是目标投影系统有几种定义投影系统的方式，比如常用的，直接写WGS84，或者用EPSG:number来定义，WGS84的EPSGnumber为4326，Web墨卡托投影的EPSGnumber为900913（谷歌地图使用的都是该投影），也可用gdalinfo直接查看目标投影tif文件的EPSGnumber
4、由geotiff文件提取每个像元值，即txt文件
gdal_translate input.tif -of XYZ -ot float32 -r bilinear-a_nodata NAN output.tif

参考