9月に入ってから、何回か地震のような微震を感じ、地震だと思って天気サイトの地震情報を見ても地震は発生していないという現象を経験している。自分が揺れているだけなのか、それとも本当に自分の家の近辺でtenki.jpに載らない地震が発生しているのかを知るために、9月に入ってからの全ての地震をプロットして検証してみることにした。
import folium
import pandas as pd
df = pd.read_csv('sep_2019_earthquakes.csv',encoding='utf-8')
ロードしたデータをチラ見する。
df.head()
df.tail()
9月1日〜9月21日の3週間で10209回も地震が発生しているようである。マグニチュードがマイナスって何なのかを調べてみたらこのサイトにて以下のように説明されていた。
(梅田):そうです。マグニチュードマイナス2ですと、1.6m=160cmくらいの断層が動いたことになります。
マイナスマグニチュードは非常に微小な地震だということらしい。
df1 = df[9609:10209]
m = folium.Map(location=[35.658099,139.741357],tiles="OpenStreetMap",zoom_start=5)
for i in range(0,len(df1)):
if df1.iloc[i]['マグニチュード']>=-1 and df1.iloc[i]['マグニチュード']<1.0:
folium.CircleMarker([df1.iloc[i]['経度'],df1.iloc[i]['緯度']],
radius=10*df1.iloc[i]['マグニチュード']+1,
popup=folium.Popup('<div style="font-size: 18pt; color : black">'
+"発生日:"+df1.iloc[i]['日']+"<br>"+"発生時刻:"+df1.iloc[i]['時']+
"<br>"+"震源地:"+df1.iloc[i]['地名']+"<br>"+"マグニチュード:"+
df1.iloc[i]['マグニチュード'].astype(str)+"<br>"+"震源の深さ:"+
df1.iloc[i]['深さ'].astype(str)+"km"+'</div>',max_width=330,min_width=330,
parse_html=False),color='purple',fill_color='purple',).add_to(m)
if df1.iloc[i]['マグニチュード']>=1.0 and df1.iloc[i]['マグニチュード']<3.0:
folium.CircleMarker([df1.iloc[i]['経度'],df1.iloc[i]['緯度']],
radius=9.5*df1.iloc[i]['マグニチュード']+1,
popup=folium.Popup('<div style="font-size: 18pt; color : black">'
+"発生日:"+df1.iloc[i]['日']+"<br>"+"発生時刻:"+df1.iloc[i]['時']+
"<br>"+"震源地:"+df1.iloc[i]['地名']+"<br>"+"マグニチュード:"+
df1.iloc[i]['マグニチュード'].astype(str)+"<br>"+"震源の深さ:"+
df1.iloc[i]['深さ'].astype(str)+"km"+'</div>',max_width=330,min_width=330,
parse_html=False),color='darkyellow',fill_color='darkyellow',
).add_to(m)
if df1.iloc[i]['マグニチュード']>=3.0 and df1.iloc[i]['マグニチュード']<4.0:
folium.CircleMarker([df1.iloc[i]['経度'],df1.iloc[i]['緯度']],
radius=9*df1.iloc[i]['マグニチュード']+1,
popup=folium.Popup('<div style="font-size: 18pt; color : black">'
+"発生日:"+df1.iloc[i]['日']+"<br>"+"発生時刻:"+df1.iloc[i]['時']+
"<br>"+"震源地:"+df1.iloc[i]['地名']+"<br>"+"マグニチュード:"+
df1.iloc[i]['マグニチュード'].astype(str)+"<br>"+"震源の深さ:"+
df1.iloc[i]['深さ'].astype(str)+"km"+'</div>',max_width=330,min_width=330,
parse_html=False),color='green',fill_color='green',).add_to(m)
if df1.iloc[i]['マグニチュード']>=4.0 and df1.iloc[i]['マグニチュード']<5.0:
folium.CircleMarker([df1.iloc[i]['経度'],df1.iloc[i]['緯度']],
radius=8.5*df1.iloc[i]['マグニチュード']+1,
popup=folium.Popup('<div style="font-size: 18pt; color : black">'
+"発生日:"+df1.iloc[i]['日']+"<br>"+"発生時刻:"+df1.iloc[i]['時']+
"<br>"+"震源地:"+df1.iloc[i]['地名']+"<br>"+"マグニチュード:"+
df1.iloc[i]['マグニチュード'].astype(str)+"<br>"+"震源の深さ:"+
df1.iloc[i]['深さ'].astype(str).astype(str)+"km"+'</div>',max_width=330,min_width=330,
parse_html=False),color='violet',fill_color='violet',).add_to(m)
if df1.iloc[i]['マグニチュード']>=5.0:
folium.CircleMarker([df1.iloc[i]['経度'],df1.iloc[i]['緯度']],
radius=8*df1.iloc[i]['マグニチュード']+1,
popup=folium.Popup('<div style="font-size: 18pt; color : black">'
+"発生日:"+df1.iloc[i]['日']+"<br>"+"発生時刻:"+df1.iloc[i]['時']+
"<br>"+"震源地:"+df1.iloc[i]['地名']+"<br>"+"マグニチュード:"+
df1.iloc[i]['マグニチュード'].astype(str)+"<br>"+"震源の深さ:"+
df1.iloc[i]['深さ'].astype(str)+"km"+'</div>',max_width=330,min_width=330,
parse_html=False),color='blue',fill_color='blue',).add_to(m)
m