Arcgis栅格数据转至GMT支持的.grd格式绘图

前言

GMT画图支持“.grd”格式数据，

一、方法1（复杂）

（1）Arcgis栅格数据-----栅格转点-----添加xy坐标----属性表导出文本格式

（2）gmt xyz2grd-----grdimage实现GMT绘图

二、方法2（推荐）

GDAL工具中的gdal_translate很好地解决了这个问题，只需要将Arcgis栅格数据导出成tif格式，再利用gdal_translate就可以把栅格数据转换成GMT所需的grd数据。

类似的方法还支持其他软件如GAMMA等软件生成的结果转换。

Ubuntu/Debian 下:

sudo apt-get install gdal-bin

CentOS/Fedora 下:

sudo yum install gdal

安装之后执行

gdal_translate -of GMT test.tif test.grd

对于其他支持格式，可通过以下语法查询

gdal_translate --formats

GMT生成的grd格式数据导入Arcgis也可用其转换

gdal_translate -of GTIFF vel_asc_clip.grd vel_asc_clip.tif