ｖ
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
from mpl_toolkits.basemap import Basemap
from matplotlib.collections import PolyCollection
import numpy as np
import matplotlib as mpl
fig=plt.figure(figsize=(5,4),dpi=200)
mpl.rcParams[‘font.family’] = ‘IPAGothic’
map =Basemap(llcrnrlon=125,llcrnrlat=25,urcrnrlon=146,urcrnrlat=48,resolution=’l’)
ax = Axes3D(fig)
ax.set_axis_off() #座標軸、面を非
ax.azim = -50 #方位角
ax.dist = 7
ax.add_collection3d(map.drawcoastlines(linewidth=0.2))
lons = np.array([141.347899,140.740593,141.152667,140.872103,140.102334,
140.363634,140.467521,140.446793,139.883565,139.060156,
139.648933,140.123308,139.691704,139.642514,139.023221,
137.211338,136.625573,136.221642,138.568449,138.181224,
136.722291,138.383054,136.906565,136.508591,135.86859,
135.755608,135.519711,135.183025,135.832744,135.167506,
134.237672,133.050499,133.934675,132.459622,131.4705,
134.559303,134.043444,132.765362,133.53108,130.418314,
130.298822,129.873756,130.741667,131.612591,131.423855,130.557981,127.680932])
#緯度y
lats = np.array([43.063968,40.824623,39.703531,38.268839,39.7186,
38.240437,37.750299,36.341813,36.565725,36.391208,
35.857428,35.605058,35.689521,35.447753,37.902418,
36.69529,36.594682,36.065219,35.664158,36.651289,
35.391227,34.976978,35.180188,34.730283,35.004531,
35.021004,34.686316,34.691279,34.685333,34.226034,
35.503869,35.472297,34.661772,34.39656,34.186121,
34.06577,34.340149,33.84166,33.559705,33.606785,
33.249367,32.744839,32.789828,33.238194,31.91109,31.560148,26.212401])
Prefectures_name = np.array([“北海道”, 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,
“京都府”,26,”大阪府”,27,”兵庫県”,28,”奈良県”,29,”和歌山県”,30,
“鳥取県”,31,”島根県”,32,”岡山県”,33,”広島県”,34,”山口県”,35,
“徳島県”,36,”香川県”,37,”愛媛県”,38,”高知県”,39,”福岡県”,40,
“佐賀県”,41,”長崎県”,42,”熊本県”,43,”大分県”,44,”宮崎県”,45,”鹿児島県”,46,”沖縄県”,47])
dd=[
[“Tokyo”,13-1,51],
[“Hokkaido”, 1-1,70],
[“Aichi”,23-1,28],
[“Kanagawa”,14-1,21],
[“Wakayama”,30-1,13],
[“Chiba”,12-1,14],
[“Kumamoto”,43-1,5],
[“Ishikawa”,17-1,7],
[“Okinawa”,47-1,3],
[“Kyouto”,26-1,2],
[“Tpchigi”, 9-1,1],
[“Saitama”,11-1,2],
[“Fukuoka”,40-1,2],
[“Nagano”,20-1,2],
[“Mie”,24-1,1],
[“Nara”,29-1,1],
[“Osaka”,27-1,5],
[“Tokushima”,36-1,1],
[“Gifu”,21-1,2],
[“Miyagi”,4-1,1]
]
l=len(dd)
population=np.zeros(len(lons))
for i in range(l):
population[dd[i][1]]=dd[i][2]/60
x, y = map(lons, lats)
ax.bar3d(x, y, np.zeros(len(x)), .4, .4, population, color= ‘r’, alpha=0.9) #
for i in range(l):
ax.text(x[dd[i][1]]-0.9,y[dd[i][1]]-0.9,0.015*dd[i][2]+0.1,dd[i][0]+” “+str(dd[i][2]))
plt.show()
Yukki7's Blog 