摘要: 1.数据准备 (1)到http://www.gadm.org/download 下载中国shapefile格式的资料,下载后的文件名为CHN_adm_shp.zip ,解压后如图: 主要用到的文件是CHN_adm1.shp 2.绘图 import time...
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: