GMT修改nc网格文件的Z值并制图

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

---

前言

本节主要激励利用GMT修改nc网格文件的Z值并制图，如有错误，请多多指教。

一、先上结果

二、代码

set GMT_SESSION_NAME=12648
gmt begin figure2 png
    REM 读取网格文件的相关信息
    gmt grdinfo SRTM15_V2.4.nc 
    REM 对网格文件中的Z值进行截断
    gmt grdclip SRTM15_V2.4.nc -Gnew_SRTM15.nc -R110/120/10/20 -Sa0/NaN
    REM 对新网格文件进行制图  
    gmt grdimage new_SRTM15.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 

三、代码注释

grdclip SRTM15_V2.4.nc -Gnew_SRTM15.nc -R110/120/10/20 -Sa0/NaN

grdclip:对网格文件中的Z值进行截断;
SRTM15_V2.4.nc:原来的网格文件的名称
-G：输出新的网格文件名；
new_SRTM15.nc：新的网格文件名
-R：截取范围，110/120/10/20 经纬度；
-Sa0/NaN：将所有大于 0的值设置为 NaN;