R

気象庁地震カタログ(米子市周辺）

ggmap パッケージ

(参考)

気象庁 震源リスト(最近9ヶ月、日データ)
これより過去の震源については地震月報(カタログ編)

記事「気象庁地震カタログ02」の続き

ggmapパッケージを使って「米子市周辺」の震源マップを描く。

背景となる地図を得る

1
2
3
4
5
6
7
8
library(ggmap)
library(mapproj)
#center <- geocode('tottori')
#center <- geocode('kurayoshi')
center <- geocode('yonago')
#maptype = c(“terrain”, “satellite”, “roadmap”, “hybrid”, “toner”, “watercolor”)
map <- get_map(c(center$lon, center$lat), zoom = 10, maptype = 'roadmap')
#ggmap(map)

マグニチュードで色分け

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
#load("tottori1998_2014.RData")
#
subeqdata<-subset(eqdata,longitude>=132.8 & longitude<=133.8 & latitude>=35 &  latitude<=35.8)
#M>=3、微小地震、極微小地震で色を買えるので分類する
#極微小地震
subeqdata1<-subset(subeqdata,mag<1)
#微小地震
subeqdata2<-subset(subeqdata,mag>=1 & mag<3)
#M>=3(マグニチュードの昇順に並べ替え)
subeqdata3<-subset(subeqdata,mag>=3)
sortlist <- order(subeqdata3$mag)
subeqdata3<- subeqdata3[sortlist,]
#
head(subeqdata1) ; tail(subeqdata1)
head(subeqdata2) ; tail(subeqdata2)
head(subeqdata3) ; tail(subeqdata3)
#
#
#png("eqgg01.png")
alpha=0.2
ggmap(map) + 
geom_point(data=subeqdata1, aes(x=longitude, y=latitude), color=rgb(1,0,0,alpha=alpha),size=1)+ 
geom_point(data=subeqdata2, aes(x=longitude, y=latitude), color=rgb(0,1,0,alpha=alpha),size=1)+
geom_point(data=subeqdata3, aes(x=longitude, y=latitude), color=rgb(0,0,1,alpha=alpha),size=1)
#dev.off()

1998年から2014年までの震源マップ（マグニチュードで色分け）

  • mag<1:赤
  • mag>=1 & mag<3:緑
  • mag>=3:青

震源の深さで色分け

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
#depthによって色分け
#depth<20
subeqdata4<-subset(subeqdata,depth<20)
#depth>=20
subeqdata5<-subset(subeqdata,depth>=20)
#
head(subeqdata4) ; tail(subeqdata4)
head(subeqdata5) ; tail(subeqdata5)
#
#
#png("eqgg02.png")
alpha=0.2
ggmap(map) + 
geom_point(data=subeqdata5, aes(x=longitude, y=latitude), color=rgb(0,0,1,alpha=alpha),size=1)+ 
geom_point(data=subeqdata4, aes(x=longitude, y=latitude), color=rgb(1,0,0,alpha=alpha),size=1)
#dev.off()

1998年から2014年までの震源マップ(震源の深さで色分け)

  • depth<20:赤
  • depth>=20:青

プロットする順序を変えたもの

1
2
3
4
5
6
#png("eqgg03.png")
alpha=0.2
ggmap(map) + 
geom_point(data=subeqdata4, aes(x=longitude, y=latitude), color=rgb(1,0,0,alpha=alpha),size=1)+
geom_point(data=subeqdata5, aes(x=longitude, y=latitude), color=rgb(0,0,1,alpha=alpha),size=1)
#dev.off()