小二 发表于 2013-2-23 10:05:08

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

GMT 以固定直径画圆方法
雷达高度计的足迹直径大概是5Km,那么为了在图中表示点的覆盖范围,以及同一轨道不同周期点的足迹偏差,那么有必要以实际半径在图中绘制足迹,而不是以相对于纸张大小的半径绘制足迹位置。关键就在使用psxy的-SE密令,这是个绘制椭圆的命令,只要设置长半轴和短轴相等即可画圆,且单位是Km。下面绘制TP雷达GDR和SGR 1HZ(10个原始采样数据平均)不同周期的数据足迹范围,可见TP的足迹覆盖是相当准确的。

代码如下:gawk "{print $2,$1,$3}" ExportAsciitpgdr.txt > ExportAsciitpgdr.dgawk "{print $3,$2,$1}" tpSdr.txt> tpSdr.dpsbasemap -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.pspscoast -R -J -W1/0.8p-N1/0.5p -A100l-Df -V -O -K -P>>Gtrack.psgawk "{print $1, $2, 0, 5, 5}" tpSdr.d | psxy -R -Glightyellow -J -O -K -SE -Wthinnest >> Gtrack.ps

页: [1]
查看完整版本: [转]GMT 以固定直径画圆方法