python - 用cartopy画一张特定国家的地图？

Question

我在这里试过这个例子，效果很好。但是我如何专注于另一个国家，即只显示德国？

这个例子的来源

http://scitools.org.uk/cartopy/docs/latest/examples/hurricane_katrina.html

我在 extend() 方法上尝试了一些坐标，但我没有设法让它看起来像美国地图。还是我必须修改形状文件？

Answer 1

使用http://www.gadm.org/country上的全球行政区数据集，只需下载德国数据集并使用 cartopy 的 shapereader（与链接示例中的操作方式相同）。

一个简短的自包含示例：

import cartopy.crs as ccrs
import cartopy.io.shapereader as shpreader
import matplotlib.pyplot as plt

# Downloaded from http://biogeo.ucdavis.edu/data/gadm2/shp/DEU_adm.zip
fname = '/downloads/DEU/DEU_adm1.shp'

adm1_shapes = list(shpreader.Reader(fname).geometries())

ax = plt.axes(projection=ccrs.PlateCarree())

plt.title('Deutschland')
ax.coastlines(resolution='10m')

ax.add_geometries(adm1_shapes, ccrs.PlateCarree(),
                  edgecolor='black', facecolor='gray', alpha=0.5)

ax.set_extent([4, 16, 47, 56], ccrs.PlateCarree())

plt.show()

高温高压

Answer 2

让我举个例子，使用来自naturalearthdata的数据。因此，可以将其扩展到任何国家。

from cartopy.io import shapereader
import numpy as np
import geopandas
import matplotlib.pyplot as plt

import cartopy.crs as ccrs

# get natural earth data (http://www.naturalearthdata.com/)

# get country borders
resolution = '10m'
category = 'cultural'
name = 'admin_0_countries'

shpfilename = shapereader.natural_earth(resolution, category, name)

# read the shapefile using geopandas
df = geopandas.read_file(shpfilename)

# read the german borders
poly = df.loc[df['ADMIN'] == 'Germany']['geometry'].values[0]

ax = plt.axes(projection=ccrs.PlateCarree())

ax.add_geometries(poly, crs=ccrs.PlateCarree(), facecolor='none', 
                  edgecolor='0.5')

ax.set_extent([5, 16, 46.5, 56], crs=ccrs.PlateCarree())

这会产生下图：