GMT画图中图片的网格分散问题

问题描述和解决方案

在GMT画关于地学方面的图片时，出现了下图中的网格分散问题：

GMT画图的网格分散问题.png

查找代码，以及与人交流，最后定位到此行代码问题：

gmt xyz2grd xyz.xyz -G$GRD2 -I0.0833333333333/0.0833333333333 -R$R

通过查看gmt的xyz2grd函数，发现是网格分辨率的问题：-I，以及-R。

对于xyz.xyz数据文件，使用gmtinfo命令查看其属性：

➜  gmtinfo xyz.xyz
  xyz.xyz: N = 9331200  <-179.958328/179.958328>    <-89.958336/89.958336>  <-10513/6573>

• 对应的x值为：-179.958328/179.958328，即为经度范围，间隔为：0.083328；
• 对应的y值为：-89.958336/89.958336，即为纬度范围，间隔为：0.083328；
• 对应的z值为：-10513/6573，即为数据范围。

故此，问题代码行改为：

gmt xyz2grd xyz_BED2014.xyz -G$GRD2 -I0.083328/0.083328 -R$R -V

正确的图片就出来了：

GMT画图的网格分散问题解决后的图片.png

至此，此问题就解决了。

总结一般规律

• 使用gmt xyz2grd 的时候一定要加上-V 选项，输出的信息可以帮助判断自己的转换是否正确；
• 不仅要考虑到数据xy的间隔，即为调整-I的分辨率；
• 而且xy各自的最大最小值范围，即为-R，也需要正确；

Reference

• xyz2grd — GMT模块手册
• XYZ2GRD — 官网