Basemap绘制四川轮廓

Python与开源GIS

Basemap绘制四川轮廓

2017-10-31 作者: xuzhiping 浏览: 4916 次

摘要: 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:

关注“开源集思”公众号
获取免费资源

随机推荐


Copyright © 2014-2019 OSGeo中国中心 吉ICP备05002032号

Powered by TorCMS

OSGeo 中国中心 邮件列表

问题讨论 : 要订阅或者退订列表,请点击 订阅

发言 : 请写信给: osgeo-china@lists.osgeo.org