平常使用GMT画图,如果需要用到栅格数据的话,则需要用“.grd”的数据,早前都是使用matlab将二进制的数据转换成xyz的三列的文本,再用gmt 内置的 xyz2grd进行转换,这样的转换不仅效率低,而且极其占用硬盘空间。
而GDAL工具中的gdal_translate工具则很好地解决了这个问题,只需要将数据转换成tif格式,再利用gdal_translate就可以把栅格数据转换成GMT所需的grd数据了。
为了得到tif数据,可以使用GAMMA软件中自带的"data2geotiff"工具,最后用gdal_translate就可以了,在gdal_translate中,需要将输出的格式设置为GSBG,例如:
gdal_translate -of GSBG SouthCA.tif SouthCA.grd
对于其他支持格式,可以参考:
gdal_translate --formats