3

我有以下代码：

from mpl_toolkits.basemap import Basemap
map = Basemap(projection='merc', lat_0=50, lon_0=4,
    resolution = 'l', area_thresh = 0.1,
    llcrnrlon=4, llcrnrlat=50,
    urcrnrlon=40, urcrnrlat=60)

map.drawcoastlines(linewidth=0.15)
map.drawcountries(linewidth=0.15)
map.fillcontinents(color='brown',lake_color='white')
map.drawmapboundary(fill_color='white')

在此处输入图像描述

在这张地图的顶部，我想显示一个仅包含一个多边形的 shapefile。多边形定义了一个封闭区域。我找到了有关如何手动添加多边形或从 shapefile 绘制多个多边形的不同教程，但我无法针对我的情况执行此操作。shapefile 属性表仅由两个字段组成：“ c ”和“ area ”。

现在我已经到了这个

import shapefile

s = shapefile.Reader(filepath,'c',drawbounds=False)
shapes = s.shapes()
records = s.records()
for record, shape in zip(records,shapes):
    lons,lats = zip(*shape.points)
    data = np.array(map(lons, lats)).T
x, y =map(lons,lats)
4

1 回答 1

1

有同样的问题，但它是如此简单，以至于你从来没有想过用网上的许多教程和模块以及类似的问题来做这件事：

map.readshapefile('luthuania', 'any_name_you_like', drawbounds=True)

所以对于你的例子：

from mpl_toolkits.basemap import Basemap
map = Basemap(projection='merc', lat_0=50, lon_0=4,
    resolution = 'l', area_thresh = 0.1,
    llcrnrlon=4, llcrnrlat=50,
    urcrnrlon=40, urcrnrlat=60)

map.readshapefile('luthuania', 'any_name_you_like', drawbounds=True, linewidth=2, color='b')

map.drawcoastlines(linewidth=0.15)
map.drawcountries(linewidth=0.15)
map.fillcontinents(color='brown',lake_color='white')
map.drawmapboundary(fill_color='white')

这使

立陶宛

顺便说一句，Basemap 在后台使用模块 shapefile：请参阅 C:\Python33\Lib\site-packages\mpl_toolkits\basemap\shapefile.py

于 2013-12-11T07:03:02.070 回答