【Python】　Basemapでプロットいっぱい　－深発地震－

目　次

データファイルから複数地点プロット　－Python－
深発地震

データファイルから複数地点プロット　－Python－

　データファイル（eq20190728.csv）の中身は、

PREF_NO,PREF,MUNI_NO,MUNI,LATITUDE,LONGITUDE,LEVEL
4,宮城県,4341,伊具郡 丸森町,37.9114293,140.7655281,4
4,宮城県,4361,亘理郡 亘理町,38.0374595,140.8525812,3
7,福島県,7204,いわき市,37.0504195,140.8876817,3
・
・
・

　データは防災科学技術研究所　J-RISQ　地震速報（www.j-risq.bosai.go.jp/report/）などから。
　　※　『気象庁発表の情報と一致しない場合がある。』

　以下、.pyファイル。
　.csvファイルの読み込みなどに
　pandasライブラリ（pandas.pydata.org/）
　を利用。

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import numpy as np
import pandas as pd

fig=plt.figure(figsize=(10,10))
map = Basemap(projection='merc',llcrnrlat = 28.0,urcrnrlat = 42.0,llcrnrlon = 126.0,urcrnrlon = 150.0,resolution = "i")

map.drawcoastlines(linewidth=0.5,color='#000000')
map.fillcontinents(color = '#9acd32',lake_color='aqua')
map.drawmapboundary(fill_color='#7fffd4')
map.drawmeridians(np.arange(125,150,5),labels=[0,0,0,1],fontsize=12)
map.drawparallels(np.arange(25,50,5),labels=[1,0,0,0],fontsize=12)

point = pd.read_csv('eq20190728.csv',index_col=0)

point4 = point[point['LEVEL']==4]
point3 = point[point['LEVEL']==3]

lon4 = point4['LONGITUDE'].tolist()
lat4 = point4['LATITUDE'].tolist()
xpt4,ypt4 = map(lon4,lat4)

map.plot(xpt4,ypt4,"ro",markersize=7)

lon3 = point3['LONGITUDE'].tolist()
lat3 = point3['LATITUDE'].tolist()
xpt3,ypt3 = map(lon3,lat3)

map.plot(xpt3,ypt3,"mo",markersize=5)

xptx,yptx = map([137.40],[33.16])

map.plot(xptx,yptx,"rx",markersize=7)

base = "plotippai.png"
plt.savefig(base)
plt.show()