.
.
数语记吾学,以备不时之需,若遇同仁愿得赐教。
.
.
GMT6.0绘制青藏高原某湖泊 ,GMT中文手册104页已有详细介绍,仅做部分修改。
gmt begin 1 png
gmt coast -R85/86/30.5/31.3 -JM12c -W0.2p -A500 -BWSen+t"zhari_nam_co" -Gwhite -Sazure1 --FORMAT_GEO_MAP=dddF
gmt inset begin -DjBL+w3c/3.6c+o0.1c -F+gwhite+p0.3p
gmt coast -R73/135/3.5/53.5 -JM? -N1/0.5p -A10000 -W0.5p,black
echo 85 30.5 86 32 | gmt plot -Sa0.3c+s -Gred -W0.1p,red
gmt inset end
gmt plot 2.txt -Baf -BWSen -Sc0.1c -W0.001p,black -Gred
gmt end show
inset begin 定义了小图的位置位于大图左下角(-DjBL),小图区域的宽度为 3 厘米,高度为 3.6 厘米(+w3c/3.6c),并且相对大图左下角偏移 0.1 厘米(+o0.1c)。同时还设置了小图区域的背景色为白色(+gwhite),并绘制了小图区域的边框(+p0.3p)。
在小图区域内,我们使用 coast 模块绘制了亚洲全图,小图投影参数为 -JM?,其中 ?表示根据 inset begin 中 -D 选项指定的宽度自动决定小图宽度。(同时可以使用 plot 模块的 -Sr+s 选项在小图中绘制了一个对应于大图区域的矩形框,该命令需要输入矩形区域的两个对角顶点的经纬度。)此处,我更改了 -Sr 为 -Sa 并指定半径为0.3绘制了一个五角星,-Gred填充为红色。
使用 gmt plot 模块绘制了湖面足迹。文件名为“2.txt” 共两列 第一列经度 第二列 纬度。绘制圆形(-Sc)半径0.1 线条粗细(0.001)黑色,填充红色(-Gred)
绘制的图片如下图所示:
脚本尚有不完善的地方,时间仓促,以实现绘图功能为目标,未做详细修改。
致谢:
GMT中文社区
GMT 中文手册 v6.0
陈师兄分享GMT学习资源,提供学习指导。