GMT读取nc文件并制图

提示：文章写完后，目录可以自动生成，如何生成可参考右边的帮助文档

前言

新手学习GMT制图，如有错误，请多多指教。
本节主要通过GMT读取nc文件并制图。

1、先上结果图

2、代码

set GMT_SESSION_NAME=12648
gmt begin figure1 pnggmt grdimage SRTM15_V2.4.nc -JM10c -R110/120/10/20 -Baf -BWSen -Cmatlab/jet -I+d Gmt colorbar -C -DjMR+w8c/0.3c+o-1c/0c+m -Ba1000f500 -By+lm
gmt end showecho 

3、代码详解

• gridimage :读取一个grd网格文件，在每个网格节点上绘制小矩形并根据z值填充颜色；
• SRTM15_V2.4.nc:官网下载地形数据；
• -JM10c：J投影 M选择墨卡托 10c:宽度为10cm;
• -R110/120/10/20:范围 东经110到120 北纬10度到20度；
• -Baf :自动调节横轴与纵轴的刻度与标准；
• -BWSen：W S 表示西和南 有刻度与标注；en表示东和北有刻度无标注；
• -Cmatlab/jet：选择颜色；
• -I+d ：光照效果 有立体感

第二行

• colorbar：右侧图例；
• -C 按照指定的cpt,其实就是第一行的Cmatlab/jet；
• -DjMR+w8c/0.3c+o-1c/0c+m -D表示色标的尺寸和位置。j表示通过 2 个字母的对齐方式码指定锚点作为参考点，M不懂，R代表右侧；+w8c/0.3c表示 长8厘米宽0.3厘米； m表示图例的文字部分放在右侧，不然默认在左侧；
• -Ba1000f500：每隔1000标注一下，隔500有个刻度；
• -By+lm： 添加图例上面的m。