使用Python实现绘制地图的示例详解

import cartopy.crs as ccrs
import matplotlib.pyplot as plt
ax = plt.axes(projection=ccrs.PlateCarree())
ax.coastlines()
plt.show()

type(ax)
# <class 'cartopy.mpl.geoaxes.GeoAxes'>

ax = plt.axes(projection=ccrs.PlateCarree(180))
ax.stock_img()
plt.show()

import numpy as np
import cartopy.crs as ccrs
def sample_data(shape=(20, 30)):
    crs = ccrs.RotatedPole(pole_longitude=177.5, pole_latitude=37.5)
    x = np.linspace(311.9, 391.1, shape[1])
    y = np.linspace(-23.6, 24.8, shape[0])
    X, Y = np.deg2rad(np.meshgrid(x, y))
    u = 10 * (2 * np.cos(2 * X + 3 * (Y + 30)) ** 2)
    v = 20 * np.cos(6 * np.deg2rad(X))
    return x, y, u, v, crs

import matplotlib.pyplot as plt

fig = plt.figure(figsize=(10, 5))
ax = fig.add_subplot(1, 1, 1, projection=ccrs.PlateCarree())
ax.set_extent([-90, 80, 10, 85], crs=ccrs.PlateCarree())
ax.stock_img()
ax.coastlines()

x, y, u, v, vector_crs = sample_data(shape=(10, 14))
ax.barbs(x, y, u, v, length=5,
    sizes=dict(emptybarb=0.25, spacing=0.2, height=0.5),
    lw=0.95, transform=vector_crs)
plt.show()