GMT 添加断层名

在 GMT 中绘制断层很简单,只要有断层的坐标数据,直接用 psxy 即可绘制。

有时不仅要绘制断层,还需要在断层线的旁边加上断层名,最直接的想法是用 pstext 将断层名写上去。使用 pstext 添加断层名时会遇到两个问题:

  1. text 的位置;
  2. text 的角度(一般情况下,断层名与断层线平行);

如果只是几条断层的话,可以手动指定并微调 text 的位置和角度,以达到想要的效果。如果要标记的断层名有点多,或者想要批量绘制不同区域的多张断层图,手动微调的方法就有些麻烦了。当然,如果真有需求的话,可以手动微调,建立一个断层名的数据库,以后用起来就方便多了。

GMT 绘制图中图

英文一般称为 inset map 或 insert map,中文还不确定有没有专门的名词,一般常称为图中图、小图、嵌入图。

所谓 insert map,即在绘制小区域地图时,在图上的某个角落里绘制一张更大范围的地图,并在大范围地图上绘制一个方框,以更清楚地表明小区域的地理位置。

若是了解 GMT 中图层的概念以及图层叠加的原理,就可以很容易地把 insert map 画出来。由于 GMT4 和 GMT5 语法上的差异,这里同时给出 GMT4 和 GMT5 两个版本的绘图脚本。二者原理是差不多的,只是语法上有些差异,因而下面只解释 GMT5 脚本中的细节,对 GMT4 脚本不做过多介绍。

GMT 区域填色

提出问题

GMT 中可以用 psxy 绘制多段数据,画出多个闭合的多边形。一个示例命令如下:

gmt psxy input.dat -R0/360/-90/90 -JX15c/10c -Bx60 -By30 -W1p -L > test.ps