python - Cartopy 轮廓缺乏投影

Question

我正在尝试通过 cartopy 将一些数据放到等高线图上。但是，在绘制数据之后，投影似乎仍然关闭。surface_temp.X 和 surface_temp.Y 是纬度/经度，而 masked_fill 是实际数据值。这似乎在底图中有效，但我不确定为什么它不在 cartopy 中。

卡托普：

fig = plt.figure(figsize=(12,4.76), dpi=100)
fig.clf()
ax = plt.axes(projection=ccrs.Mercator())
ax.coastlines()
ax.contourf(surface_temp.X, surface_temp.Y, surface_temp.masked_fill[:], latlon = 'true', transform = ccrs.Mercator())
plt.show()

底图：

fig = plt.figure(figsize=(15,4.76), dpi=100)
        fig.clf()
        plt.axes([0,0,1,1], frameon=False)
        plt.title(title)
        m = Basemap(projection='merc',llcrnrlat=-80,urcrnrlat=80, llcrnrlon=0,urcrnrlon=360,lat_ts=20,resolution='c')
m.contourf(surface_temp.X, surface_temp.Y, surface_temp.masked_fill[:], latlon = 'true')

底图结果：

Cartopy 结果（轮廓注释掉）：

Cartopoy 结果（轮廓）

标签： pythonmatplotlibmatplotlib-basemapcartopy

Answer 1

cartopy 的范例似乎总是在纬度/经度坐标上工作。这意味着，您不应根据投影转换数据，而应保持纬度/经度。

因此，而不是

ax.contourf(..., transform = ccrs.Mercator())

你需要

ax.contourf(..., transform = ccrs.PlateCarree())

一个完整的例子：

import matplotlib.pyplot as plt
import cartopy.crs as ccrs
from cartopy.examples.waves import sample_data

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

lons, lats, data = sample_data(shape=(20, 40))

ax.contourf(lons, lats, data, transform=ccrs.PlateCarree())

ax.coastlines()
ax.gridlines()

plt.show()