gmt绘制江苏省高程异常图

最新推荐文章于 2025-03-12 19:18:53 发布

snowflakeeeee

最新推荐文章于 2025-03-12 19:18:53 发布

阅读量2.9k

点赞数 5

本文链接：https://blog.csdn.net/qq_44749244/article/details/105446875

版权

本文详细介绍如何使用GMT工具进行地图绘制，包括设置边框格式、经纬度标注格式、定义图幅范围和投影格式。通过具体步骤展示如何将txt数据转换为grd格式，利用gmtsurface生成网格数据，并使用gmtpsbasemap、gmtmakecpt、gmtgrdimage等命令完成地图底层绘制、颜色表建立、高程异常绘制及色条显示。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

gmt set MAP_FRAME_TYPE fancy+ //设置边框格式为黑白框，+号使得边框四角为圆角
gmt set FORMAT_GEO_MAP ddd:mmF //设置经纬度标注格式
set R=116.0/122.0/30.0/36.0 //定义图幅范围
set j=M112/23/6i //定义投影格式

绘图过程中需要使用江苏省高程异常网格文件。gmt支持将txt转换成grd文件。

gmt surface %data% -G%grdfile% -R116.0/122.0/30.0/36.0 -I1m
//%data%处是txt文件，%grdfile%处为希望生成的grd文件名，文件会在同级目录下生成

gmt psbasemap -R%R% -J%J% -B1.5NEWS -K >> JS.ps 
//绘制底层
gmt makecpt -T22/42/2 -Z -Crainbow >mycpt
//建立自己的颜色表，从22到42每隔2用一种颜色表示，-Z表示颜色连续
gmt psclip jiangsu.txt -R%R% -J%J% -K -O >>JS.ps
//剪切，只在jiangsu.txt范围内画图
gmt grdimage heigh.grd -R%R% -J%J% -Cmycpt -BWSEN -K -O>>JS.ps
//绘制高程异常
gmt psclip -C -O -K>>JS.ps
//剪切停止，之后的图层不再只在jiangsu.txt范围内绘制
gmt psscale -D18/9/16/0.7 -Cmycpt -By+l"m" -S -O -K>>JS.ps
//显示色条