Basemap绘制四川轮廓

1.数据准备

（1）到http://www.gadm.org/download 下载中国shapefile格式的资料，下载后的文件名为CHN_adm_shp.zip ，解压后如图：

主要用到的文件是CHN_adm1.shp

2.绘图

import time

start = time.clock()

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
from matplotlib.patches import Polygon
fig = plt.figure()
ax1 = fig.add_axes([0.1,0.1,0.8,0.8])

map = Basemap(llcrnrlon=80.33,
              llcrnrlat=3.01,
              urcrnrlon=138.16,
              urcrnrlat=56.123,
             resolution='h', projection='cass', lat_0 = 42.5,lon_0=120,ax=ax1)

shp_info = map.readshapefile("F:\\dataMining\\GpsDataMining\\data\\CHN_adm1",'states',drawbounds=True) # CHN_adm1的数据是中国各省区域

for info, shp in zip(map.states_info, map.states):
    proid = info['NAME_1']  # 可以用notepad打开CHN_adm1.csv文件，可以知道'NAME_1'代表各省的名称
    if proid == 'Sichuan':
        poly = Polygon(shp,facecolor='g',edgecolor='c', lw=3) # 绘制广东省区域
        ax1.add_patch(poly)

#map.shadedrelief() # 绘制阴暗的浮雕图

map.drawcoastlines()
map.drawcountries()
map.drawcounties()
end=time.clock()
print(end-start)
plt.show()

a:浮雕图

b: