在cartopy地图顶部绘制shapefile城市边界
在cartopy地图顶部绘制shapefile城市边界
提问于 2019-07-01 03:08:32
我正在尝试使用here获取的shapefile并遵循this example，在一张地形图上绘制湾区城市/城镇边界的轮廓。由于某些原因，边框没有显示，即使我通过zorder指定边框在顶部。我是不是遗漏了什么？
# import functions
import matplotlib.pyplot as plt
import cartopy.io.img_tiles as cimgt
import cartopy.crs as ccrs
from cartopy.io.shapereader import Reader
from cartopy.feature import ShapelyFeature
# Create a Stamen terrain background instance
stamen_terrain = cimgt.Stamen('terrain-background')
fig = plt.figure(figsize = (10, 10))
ax = fig.add_subplot(1, 1, 1, projection=stamen_terrain.crs)
# Set range of map, stipulate zoom level
ax.set_extent([-122.7, -121.5, 37.15, 38.15], crs=ccrs.Geodetic())
ax.add_image(stamen_terrain, 12, zorder = 0)
# Add city borders - not working
filename = r'./shapefile/ba_cities.shp' # from https://earthworks.stanford.edu/catalog/stanford-vj593xs7263
shape_feature = ShapelyFeature(Reader(filename).geometries(), ccrs.PlateCarree(), edgecolor='black')
ax.add_feature(shape_feature, zorder = 1)
plt.show()
回答 1
Stack Overflow用户
回答已采纳
发布于 2019-07-03 06:14:34
正如@ImportanceOfBeingErnest和@swatchai建议的那样，ShapelyFeature cartopy.feature.ShapelyFeature()中的CRS (坐标参考系)参数不正确。
适当的EPSG (欧洲石油勘测组？)代码可以在shapefile附带的.xml文件中找到：
<gco:CharacterString>26910</gco:CharacterString>
</code>
<codeSpace>
<gco:CharacterString>EPSG</gco:CharacterString>将此参数作为第二个参数传递到ShapelyFeature()中，就可以让shapefile正确绘制城市边界：
# Add city borders
filename = r'./shapefile/ba_cities.shp'
shape_feature = ShapelyFeature(Reader(filename).geometries(), ccrs.epsg(26910),
linewidth = 1, facecolor = (1, 1, 1, 0),
edgecolor = (0.5, 0.5, 0.5, 1))
ax.add_feature(shape_feature)
plt.show()
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/56827373
复制相关文章
相似问题
腾讯云开发者
