- 积分
- 736
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-2-7
- 最后登录
- 1970-1-1
|
发表于 2014-2-12 18:12:14
|
显示全部楼层
本帖最后由 阿阿飞飞 于 2014-2-12 18:19 编辑
可以使用basemap读取shapefile,这样最为标准,下方代码结果为
结果输出
,如果读取带省界的shape文件,自然可以显示省界。等值线直接用plt的contour命令或basemap的contour命令皆可,格式为plt.contour(x,y,value),同时可以再设置间隔,标注等。- mymap = Basemap(projection='cyl',llcrnrlon=100, urcrnrlon=150, \
- llcrnrlat=0, urcrnrlat=50, \
- lon_0=125, lat_0=25, resolution='l')
- #===================================================
- #设定海岸线边界为白色(这样与后面使用shapefile不重叠)
- mymap.drawcoastlines(linewidth=0.4,color='w')
- #导入shapefile文件,map目录下有个World.shp文件
- shpdir='D:\\ipython_data\\DataforPython\\map\\World'
- s=mymap.readshapefile(shpdir,'World',linewidth=0.6)
复制代码 |
|