[GMT学习笔记]GMT中有用的Tips备忘

最新推荐文章于 2022-07-15 18:07:12 发布

hr11asdfg

最新推荐文章于 2022-07-15 18:07:12 发布

阅读量752

点赞数

分类专栏：其他笔记文章标签： linux

本文链接：https://blog.csdn.net/weixin_42638388/article/details/108179553

版权

其他笔记专栏收录该内容

24 篇文章 13 订阅

订阅专栏

关于gmt5中-K和-O选项的用法（GMT5手册中P136）
查看网格文件.nc格式
ncdump -c ESMGFZ_NTAL_cm_v1.3_24h_2016.nc_subset16.nc
如果没有安装，则需要先运行sudo apt install netcdf-bin
关于gmt画图坐标轴的间隔设置，一直很令人头疼。在GMT5.4手册的指导下，做了几组实验，对比看一下不同-B选项画出来图的效果。为了更加直观地发现其中的差异，我们把不同的图放到同一张图上，用子图的形式来表达：

#!/bin/bash
PS=test_B.ps

J=X10c/10c
R=0/10/0/10
gmt psxy -J$J -R$R -T -K > $PS

gmt psbasemap -JX10c/10c -R0/10/0/10 -Ba2f1g1 -K -O >> $PS
gmt psbasemap -JX10c/10c -R0/10/0/10 -Bxa2 -Bya1 -X12c -K -O >> $PS
gmt psbasemap -JX10c/10c -R0/10/0/10 -Bxafg -Byafg -X12c -K -O >> $PS
gmt psbasemap -JX10c/10c -R0/10/0/10 -Ba2+1f1g1 -K -X12c -O >> $PS

gmt psxy -J$J -R$R -T -O >> $PS
gmt psconvert -A2c/2c/2c/25c -Tg -P $PS
rm -rf $PS!

rm gmt.*

生成的图如下：
在这里插入图片描述
由上图可以发现，各项的含义分别为a（标注）、f（刻度）、g（网格线），如果想让x轴和y轴不相同，则需要向上面那样分别设置。

gmt在绘图时，内置有已经搭配好的cpt颜色表文件，位于安装目录/share/cpt文件夹下。但是在绘图中，如果想自己定义cpt颜色表，就需要自己动手制作想要的颜色搭配。这里比如想画一个地表沉降的图，用格网文件读取其中的数据，然后蓝色表示沉降，红色表示抬升，类似于这样的效果：（图来自于姚朝龙, 刘立龙*, 林旭等 .0.利用GPS与环境负荷形变数据研究台风引起的垂向地表位移. 地球物理学报,(): 1-,doi: doi:10.6038/cjg2020N0116）

想制作一个这样的颜色表，查阅GMT5.4手册，发现一个小练习，位于P110。主要介绍了如何制作cpt文件以及如何将内置的cpt文件翻转、改变前后景色等，感兴趣的可以自行查阅。其生成的图如下：

可以发现，在截取的和翻转的颜色表中，有我们想要的颜色，于是打开其中生成的cpt文件，将蓝色部分和红色部分截取出来（对照这个个图把cpt文件里面相应key值对应的颜色复制出来就行，我这里舍去了中间绿色的部分），生成一个新的cpt文件，命名为my_new.cpt，注意复制过来的时候要修改key的值，让其按照从小到大的顺序。

新的cpt文件如下：如果想要直接用，可以将其保存为name.cpt放到安装目录/share/cpt文件夹下，在绘图的时候记得引用这个cpt文件，就可以达到想要的效果。

-10	8/12/21.2	-9	8/12/21.2
-9	23.6/36/63.6	-8	23.6/36/63.6
-8	38/60/106	-7	38/60/106
-7	44.4/76/127.6	-6	44.4/76/127.6
-6	50.2/91.4/149.2	-5	50.2/91.4/149.2
-5	60.6/119.8/180.4	-4	60.6/119.8/180.4
-4	75.6/157.8/215.4	-3	75.6/157.8/215.4
-3	90/185/233	-2	90/185/233
-2	130.8/205/237.8	-1	130.8/205/237.8
-1	203.4/237/248.6	0	203.4/237/248.6
 0	255/215/215	 1	255/215/215
 1	255/196/196	 2	255/196/196
 2	255/158/158	 3	255/158/158
 3	255/90/90	 4	255/90/90
 4	238/80/78	 5	238/80/78
 5	255/160/69	 6	255/160/69
 6	247/215/104	 7	247/215/104
 7	240/236/121	 8	240/236/121
 8	205/255/162	 9	205/255/162
 9	172/245/168	 10	172/245/168
 B	black
F	white
N	gray

下面看一下我绘图的效果：
在这里插入图片描述