[转]GMT 以固定直径画圆方法

小二

GMT 以固定直径画圆方法

雷达高度计的足迹直径大概是5Km,那么为了在图中表示点的覆盖范围，以及同一轨道不同周期点的足迹偏差，那么有必要以实际半径在图中绘制足迹，而不是以相对于纸张大小的半径绘制足迹位置。

关键就在使用psxy的-SE密令，这是个绘制椭圆的命令，只要设置长半轴和短轴相等即可画圆，且单位是Km。

下面绘制TP雷达GDR和SGR 1HZ(10个原始采样数据平均)不同周期的数据足迹范围，可见TP的足迹覆盖是相当准确的。

代码如下：

gawk "{print $2,$1,$3}" ExportAsciitpgdr.txt > ExportAsciitpgdr.d

gawk "{print $3,$2,$1}" tpSdr.txt> tpSdr.d

psbasemap -P -K -R87/87.5/25/26 -JM2i  -BNEWSf0.1a0.2g0.1/f0.1a0.2g0.1:."Ground points of T/P GDR(LB) and SDR(LY)":  --PLOT_DEGREE_FORMAT=ddd:mmF > Gtrack.ps

pscoast -R -J -W1/0.8p  -N1/0.5p -A100l  -Df -V -O -K -P>>Gtrack.ps

gawk "{print $1, $2, 0, 5, 5}" tpSdr.d | psxy -R -Glightyellow -J -O -K -SE -Wthinnest >> Gtrack.ps