最近看大家说地震要频繁发生了,就挖了点数据来倒腾,这是美国地震局从1973年至今大于6级地震的历史数据:
http://dl.dropbox.com/u/3824701/earthquacks.csv
[attachment=203392,161]
数据来自美国地震局(USGS),历年按月世界范围内大于6级的地震频数热度图。颜色越深代表次数越多。
PS: 去USGS都是要翻墙的了~~[s:15]
[attachment=203392,162]
图上看出1990年后全球范围内大于6级的地震明显增多,且在其后的二十年里基本是稳定的。据USGS官方给出的原因是:咱的设备牛逼了,能探测到得地震也多了
数据里面有精度和纬度,谁能按,大陆海洋以及国家分个类不胜感激[s:11]
把代码贴出来,大家交流啊
<br />
quacks<-read.table("d:/data/earthquacks.csv",sep=",",header=TRUE)<br />
#year vs month<br />
freq.q<-table(quacks$year,quacks$month)<br />
bitmap("d:/heapmap.png",width=6,height=6,pointsize=12,res=300,family="Helvetica")<br />
heatmap.q<- heatmap(freq.q, Rowv=NA, Colv=NA,<br />
col = rainbow(256,start=.5,end=.6,alpha=1), scale="column", margins=c(5,5))<br />
dev.off();<br />
#year<br />
total.q<-data.frame(table(quacks$year))<br />
count<-total.q$Freq[1:37]<br />
#bitmap("d:/line.total.png",width=6,height=6,pointsize=12,res=300,family="Helvetica")<br />
plot(1973:2009,count,type="l",col="gold",xlab="Year",ylab="Count")<br />
title("1973-2009 6级以上地震频数图")<br />
#dev.off();<br />
</p>