](a7a40474820eae152c96bb1d1e16703b.png)
参考URL
https://qiita.com/HidKamiya/items/31de6e04bcd19ed29919
うまくいけば添付した図の、全国ではなく神奈川県バージョンが出力されるはずなのですが…
すみませんが初心者なので、細かい部分はわかりません。。。
読み込むcsvは以下のような感じです。
import numpy as np import matplotlib.pyplot as plt from mpl_toolkits.basemap import Basemap def loaddata(): data = open('Kanagawa_precip.csv','r').readlines() data = [l.split(',') for l in data ] lat = [np.float(l[2]) for l in data[1:] ] lon = [np.float(l[3]) for l in data[1:] ] Jan = [np.float(l[4]) for l in data[1:] ] Feb = [np.float(l[5]) for l in data[1:] ] Mar = [np.float(l[6]) for l in data[1:] ] Apl = [np.float(l[7]) for l in data[1:] ] May = [np.float(l[8])for l in data[1:] ] Jun = [np.float(l[9]) for l in data[1:] ] Jul = [np.float(l[10]) for l in data[1:] ] Aug = [np.float(l[11]) for l in data[1:] ] Sep = [np.float(l[12]) for l in data[1:] ] Oct = [np.float(l[13]) for l in data[1:] ] Nov = [np.float(l[14]) for l in data[1:] ] Dec = [np.float(l[15]) for l in data[1:] ] Ann = [np.float(l[16]) for l in data[1:] ] #return{'lats':lats,'lons':lons,'Jan':Jan,'Feb':Feb,'Mar':Mar,'Apl':Apl,'May':May,'Jun':Jun,'Jul':Jul,'Aug':Aug, #'Sep':Sep,'Oct':Oct,'Nov':Nov,'Dec':Dec,'Ann':Ann} lonmin = 138 lonmax = 140 latmin = 35 latmax = 35.5 gspan = 5 lw = 0.5 glw = 0.3 ms = 3 ticsize = 16 labsize = 18 titsize = 20 def station_map(): data = loaddata() xlab = r'longitude($^\circ$)' # '$^\circ$':〇を上付き文字にして角度の単位 ylab = r'latitude ($^\circ$)' tit = 'Meteorological stations' # 描画画面の設定 fig = plt.figure() ax = fig.add_subplot(111) # Basemapの設定(領域指定、海岸線描画、緯度経度グリッド線描画) m = Basemap(resolution='l',llcrnrlat=latmin,llcrnrlon=lonmin, urcrnrlat=latmax,urcrnrlon=lonmax,ax=ax) m.drawcoastlines(linewidth=lw,color='k') m.drawparallels(np.arange(latmin,latmax+1,gspan),linewidth=glw) m.drawmeridians(np.arange(lonmin,lonmax+1,gspan),linewidth=glw) # プロットの描画(ks:黒い(blac'k')四角('s'quare)) m.plot(data['lons'],data['lats'],'ks',ms=ms) # 軸ラベル・表題の設定 ax.set_xticks(np.arange(lonmin,lonmax+1,gspan)) ax.set_xticklabels(map(str,np.arange(lonmin,lonmax+1,gspan)),fontsize=ticsize) ax.set_yticks(np.arange(latmin,latmax+1,gspan)) ax.set_yticklabels(map(str,np.arange(latmin,latmax+1,gspan)),fontsize=ticsize) ax.set_xlabel(xlab,fontsize=labsize) ax.set_ylabel(ylab,fontsize=labsize) ax.set_title(tit,fontsize=titsize) plt.show()
エラー
KeyError: 'PROJ_LIB'
回答1件
あなたの回答
tips
プレビュー