2-8-grd文件转换-grdconvert

17 Jul 2018

命令:grdconvert

解决的问题:
.grd后缀的文件实际上有很多种格式,标准的格式是netcdf,但是有时候输入文件并不是标准的netcdf,而是其他写入格式的,比如surfer 6 的TXT格式。这种格式在GMT绘图时会出现种种问题。

问题

当R选项的范围比实际的范围小的时候,绘制出来的图会与边框有很多缝隙,或者直接会发生偏移,如图:

nanhai_surfer.jpg

而正确的结果应当是:

nanhai_netcdf.jpg

而且,就算用grdcut把文件区域裁剪下来还是可能会出现这些问题。


解决办法

最好的解决办法就是把其他格式的grd文件转化成netcdf标准格式
命令:grdconvert
用法:

gmt grdconvert inputgrd=id -Goutgrd=id

有时候加入-N,不输出头段数据,在之后使用grdraster时有用 id表示grd文件格式的代码
那么如何知道文件格式代码呢

gmt grdinfo nanhai_surfer.grd

nanhai_surfer.grd: Title: Grid imported via GDAL nanhai_surfer.grd: Command: nanhai_surfer.grd: Remark: nanhai_surfer.grd: Gridline node registration used [Cartesian grid] nanhai_surfer.grd: Grid file format: gd = Import/export through GDAL nanhai_surfer.grd: x_min: 105.0083 x_max: 120.0083 x_inc: 0.202702702703 name: x n_columns: 75 nanhai_surfer.grd: y_min: 5.002 y_max: 24.9926 y_inc: 0.201925252525 name: y n_rows: 100 nanhai_surfer.grd: z_min: -128.49 z_max: 328.43 name: z nanhai_surfer.grd: scale_factor: 1 add_offset: 0

第五行:grid file format: gd= import/export through GDAL
表示这个文件的格式代号是gd
所以利用grdconvert把grd文件转换格式后,再用GMT绘图就没有问题了。

通常grdconvert可以自动识别grd的格式,完成转换,若不能识别,再使用a.grd=id手动加入id网格文件格式。 GMT netCDF format (32-bit float, COARDS, CF-1.5)