用GMT画M-t图
1 时间类型的数据格式必须为 yyyy-mm-ddT::,中间的那个T表示为时间格式的数据,必需的;2 在利用-R参数指定绘图区间时,同样也要用到这样的时间格式;
3 在利用-B参数指定边框刻度时,注意需要根据所绘的时间长度来指定边框刻度的标识,常用的有Y(年)、O(月)、d(日)、H(时)、M(分)、C(秒)
4 一般利用JX投影方式
#!/usr/bin/perl -w
open(SRC,"wenchuan.eqt");
@aa=<SRC>;
chomp @aa;
close SRC;
$start_t="2008-05-01T"; #开始时间
$end_t="2011-10-30T"; #结束时间
$start_m=0; #最小震级
$end_m=8; #最大震级
$region="$start_t/$end_t/$start_m/$end_m";
`psbasemap -JX7i/4i -R$region -Ba1Yf4O:"time":/1:"magnitude":WSne -K -P > mt.ps`; #每年标注一次,每4个月画一个刻度
open(REC,">tmp");
foreach $a(@aa){
$yr=substr($a,1,4);
$mm=substr($a,5,2);
$dd=substr($a,7,2);
$HH=substr($a,9,2);
$MM=substr($a,11,2);
$SS=substr($a,13,2);
$mag=substr($a,28,3);
$line1=$yr."-".$mm."-".$dd."T".$HH.":".$MM.":".$SS." $mag";
$line2=$yr."-".$mm."-".$dd."T".$HH.":".$MM.":".$SS." $start_m";
print REC "$line1\n$line2\n>\n";
}
close REC;
`psxy tmp -R -JX -m -K -P -O >> mt.ps`;
页:
[1]