摘要: Basemap是一个强大的绘制地图工具包，配合matplotlib，可以绘制一些漂亮的地图，通过网上搜集的数据就可以绘制关于人口分布、天气等不同因素在不同地域的分布情况。 在win10环境下使用的，安装basemap 在http://www.lfd.uci.e...
Basemap是一个强大的绘制地图工具包，配合matplotlib，可以绘制一些漂亮的地图，通过网上搜集的数据就可以绘制关于人口分布、天气等不同因素在不同地域的分布情况。
在win10环境下使用的，安装basemap 在http://www.lfd.uci.edu/~gohlke/pythonlibs/#basemap 这个网站下载.whl 格式的文件安装，我几乎所有的python package都是在这里下载安装的，使用.whl 文件安装可以避免一些直接用pip 命令安装产生的莫名其妙的错误，建议顺便把gdal安装了（同样在刚才那个网站查找gdal），因为官网的教程里有些是需要import osgeo ，要导入这个模块就需要安装gdal。
参考教程：http://basemaptutorial.readthedocs.io/en/latest/index.html
程序使用的美国人口及地理数据（含经纬度）：https://github.com/plotly/datasets/blob/master/2014_us_cities.csv 在excel打开预览如下：
由于是csv文件，可以用pandas模块的read_csv(filename) 函数方便读取，具体程序如下：
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import pandas as pd
import numpy as np
from matplotlib import cm
# 绘制基础地图，选择绘制的区域，因为是绘制美国地图，故选取如下经纬度，lat_0和lon_0是地图中心的维度和经度
map = Basemap(projection='stere',lat_0=90,lon_0=-105,\
llcrnrlat=23.41 ,urcrnrlat=45.44,\
llcrnrlon=-118.67,urcrnrlon=-64.52,\
rsphere=6371200.,resolution='l',area_thresh=10000)
# map = Basemap(projection='stere',
# lat_0=0, lon_0=280,
# llcrnrlon=73.33,
# llcrnrlat=3.51,
# urcrnrlon=112.16,
# urcrnrlat=53.123)
map.drawmapboundary() # 绘制边界
#map.fillcontinents() # 填充大陆，发现填充之后无法显示散点图，应该是被覆盖了
map.drawstates() # 绘制州
map.drawcoastlines() # 绘制海岸线
map.drawcountries() # 绘制国家
map.drawcounties() # 绘制县
parallels = np.arange(0.,90,10.)
map.drawparallels(parallels,labels=[1,0,0,0],fontsize=10) # 绘制纬线
meridians = np.arange(-110.,-60.,10.)
map.drawmeridians(meridians,labels=[0,0,0,1],fontsize=10) # 绘制经线
posi=pd.read_csv("D:\\Files\\python\\2014_us_cities.csv") # 读取数据
## 原始数据有3228组数据，我只选择了180个城市的数据
lat = np.array(posi["lat"][0:180]) # 获取维度之维度值
lon = np.array(posi["lon"][0:180]) # 获取经度值
pop = np.array(posi["pop"][0:180],dtype=float) # 获取人口数，转化为numpy浮点型
size=(pop/np.max(pop))*1000 # 绘制散点图时图形的大小，如果之前pop不转换为浮点型会没有大小不一的效果
x,y = map(lon,lat)
# plt.text(x, y, 'Lagos',fontsize=12,fontweight='bold',
# ha='left',va='bottom',color='k')
#
#
# x, y = map(lon[0], lat[0])
#
# plt.text(x, y, 'Barcelona',fontsize=12,fontweight='bold',
# ha='left',va='center',color='k',
# bbox=dict(facecolor='b', alpha=0.2))
#
# plt.scatter(x,y,s=size,cmap=cm.hsv,edgecolors=None,facecolors='c')
# plt.scatter(x,y,s=size,cmap=cm.hsv) # 使用matplotlib的散点图绘制函数
map.scatter(x,y,s=size) # 也可以使用Basemap的methord本身的scatter
plt.title('Population distribution in America')
plt.show()
最终效果：
由图可以直观看出纽约、旧金山、芝加哥、华盛顿是美国人口数量最多的城市。
