GMT画矢量和椭圆笔记

GMT画矢量和椭圆笔记

plot是GMT最常用的画图模块之一，输入数据的格式是

x坐标 y坐标 方位角 长度 #画矢量时-SV选项对应的输入数据
x坐标 y坐标 长轴方位角 长轴长度 短轴长度 #画椭圆时-SE选项对应的输入数据

今天在画矢量和椭圆的过程中发现他们的长度并不匹配。于是想到他们的长度单位是什么？

画矢量时单位是Jx范围/Rx范围

例如画横线

echo 0 0 90 0.5 | gmt plot -JX2c/2c -R-1/1/-10/10 -SV1c+jc -BWSne -Baf -png tmp

图片长宽各为2cm，x方向上数据范围是2（-1到1）y方向上的数据范围是20（-10到10），矢量方位角是90度，矢量长度0.5的单位是2cm/2=1cm

画椭圆时长轴长度单位是Jx范围/Rx范围，短轴长度单位是Jy范围/Ry范围

例如画椭圆（实际效果不是圆而是一条横线）

echo 0 0 90 1 1 | gmt plot -JX2c/4c -R-1/1/-20/20 -SE -Baf -BWSne -png tmp 

图片x方向长为2cm，x方向上数据范围是2 （-1 到1 ），椭圆长轴长度1的单位是2cm/2=1cm
图片y方向长为4cm，y方向上数据范围是40（-20到20），椭圆短轴长度1单位是是4cm/40=0.1cm
所以画出来是一条线。

一个例子

#!/bin/bash
gmt begin tmp png 
gmt basemap -JX4c/4c -R-1/1/-20/20 -Baf -BWSne 
echo 0 0 45 1 10 | gmt plot -SE      -W0.1p,blue 
echo 0 0 45 2    | gmt plot -SV1c+jc -W0.1p,red
gmt end

生成图片效果如下

更多资料请访问GMT中文社区：https://gmt-china.org/