2

我真的很喜欢 cartopy 可以在不同的地图投影中自动绘制的想法。但是，我不知道如何处理 Iris 立方体。作为一个姊妹项目，我希望我能做到。有可能做这样的事情吗？

import iris as I
import cartopy.crs as ccrs
import matplotlib.pyplot as plt

someCube = I.load('someCube.pp')
ax = plt.axes(projection=ccrs.Robinson())
I.plot.contourf(someCube, transform=ccrs.Robinson())
plt.show()

谢谢

4

1 回答 1

4

我采用了您的伪代码，并使其可以使用 Iris 的示例数据运行：

import iris
import iris.plot as iplt
import cartopy.crs as ccrs
import matplotlib.pyplot as plt


fname = iris.sample_data_path('air_temp.pp')
air_temp = iris.load_cube(fname)

ax = plt.axes(projection=ccrs.Robinson())
iplt.contourf(air_temp, transform=ccrs.Robinson(central_longitude=180))
ax.coastlines()

plt.show()

如果您运行此代码，您将收到以下异常：

Traceback (most recent call last):
  File "using_custom_projections.py", line 11, in <module>
    iris.plot.contourf(air_temp, transform=ccrs.Robinson())
  File "lib/iris/plot.py", line 452, in contourf
    result = _draw_2d_from_points('contourf', None, cube, *args, **kwargs)
  File "lib/iris/plot.py", line 263, in _draw_2d_from_points
    result = _map_common(draw_method_name, arg_func, iris.coords.POINT_MODE, cube, data, *args, **kwargs)
  File "lib/iris/plot.py", line 406, in _map_common
    assert 'transform' not in kwargs, 'Transform keyword is not allowed.'
AssertionError: Transform keyword is not allowed.

哪个试图告诉您您不需要告诉它立方体在哪个“变换”（或坐标系）中。原因是虹膜立方体应该包含有关基础数据的完整元数据：坐标系是该元数据的一部分。

因此，要使示例正常工作，您只需在contourf调用中删除 transform 关键字参数：

import iris
import iris.plot as iplt
import cartopy.crs as ccrs
import matplotlib.pyplot as plt


fname = iris.sample_data_path('air_temp.pp')
air_temp = iris.load_cube(fname)

ax = plt.axes(projection=ccrs.Robinson(central_longitude=180))
iplt.contourf(air_temp)
ax.coastlines()

plt.show()

轮廓结果

在iris 库中有一个类似的例子，特别是http://scitools.org.uk/iris/docs/latest/examples/graphics/rotated_pole_mapping.html#rotated-pole-mapping-03（例子中的最后一个情节）。

高温下，

于 2012-12-10T09:15:21.690 回答