GMT读取GeoTiff绘图

GMT grdimage 命令

gdal_translate 将GeoTiff转为GMT识别的grd文件

可将apps目录价添加至环境变量

可查看支持转换格式


命令行调用gdal_translate

GMT脚本

batch 代码片.

rem makecpt 默认rainbow
rem @+2上标

gmt begin 
gmt set PS_CHAR_ENCODING Standard+
gmt set FONT_LABEL 18p,39,black
set projection="M10c"
set region="1/2/3/4"
gmt makecpt -Cwhite,black -T0/20000

    gmt figure a png E720
    gmt grdimage a.grd -J%projection% -R%region% -Baf -BWSen
    gmt colorbar -Bxaf+l"方差(m@+2)"

    gmt figure b png E720
    gmt grdimage b.grd -J%projection% -R%region% -Baf -BWSen
    gmt colorbar -Bxaf+l"方差(m@+2)"

    gmt figure c png E720
    gmt grdimage c.grd -J%projection% -R%region% -Baf -BWSen
    gmt colorbar -Bxaf+l"方差(m@+2)"

gmt end show

参考：https://docs.gmt-china.org/latest/