问题标签 [cartopy]

我想在用 cartopy 制作的地图上绘制等高线。如果我在以下位置运行示例：

http://scitools.org.uk/cartopy/docs/latest/matplotlib/introductory_examples/03.contours.html

使用来自 github 的最新代码结帐，我得到：

应该如何做到这一点？

（作为对原始公告的回应http://matplotlib.1069221.n5.nabble.com/ANN-Cartopy-A-new-mapping-library-using-matplotlib-td39886.html）

是否可以使用最新版本的 cartopy 在极地立体地图上绘制等高线？我想看一个例子来说明这是如何完成的，因为我自己正在努力解决这个问题！

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

谢谢

我收到了以下电子邮件，并希望确保每个人都能获得此问题的答案：

你好，

我想使用 cartopy 设置一个简单的纬度经度图，它穿过日期变更线，左侧显示东亚，右侧显示北美西部。以下谷歌地图大致是我所追求的：

https://maps.google.co.uk/?ll=56.559482,-175.253906&spn=47.333523,133.066406&t=m&z=4

这可以用 Cartopy 完成吗？

我想绘制一个北半球的极地立体图，在图的底部有 180，这样我就可以强调太平洋地区。我正在使用来自 git 的最新 cartopy，可以制作极地立体图，但我不知道如何更改图底部的经度。我尝试将经度范围设置为 [-180, 180] 但这无济于事，并且 NorthPolarStereo() 不接受任何关键字参数，如 central_longitude。目前这可能吗？

我正在使用 cartopy 生成带有立体投影的北极地图，然后在顶部绘制一条线（以显示横截面的位置）。如果我使用下面的代码，那么这条线不会沿着一条直线穿过极点，而是沿着一条纬线。

为了解决这个问题，我必须将 x 和 y 更改为：

所以在北极有两个数据点。有没有更好的解决方案？

编辑：附上图片

谢谢，

蒂姆

如果我定义一组具有给定高度和宽度的（地理）轴，我如何确保绘图将填充这些轴？

这会为地图生成一个具有合理纵横比的图，但我希望它改为填充轴，从而导致图的高度大于宽度。这只是一个示例，但在现实世界中这样做很重要。

编辑

只是为了澄清，我实际上希望结果被扭曲，所以在我的例子中，我真的想要一张全局范围大于宽度的地图。使用ax.set_aspect('auto')似乎适用于PlateCarree和NorthPolarStereographic预测，但可能不适用于所有人（OSGB例如）。

我一直在按照用户 pelson 的指示创建一个填充了国家形状的地图。（在python底图中填写国家）

现在我很想把这一步更进一步并创建一个像这样的 html 网站： http ://www.ethnologue.com/region/NEU 我不需要那些花哨的弹出窗口，而是那些链接（http://www.w3schools .com/tags/att_area_href.asp）对于每个国家来说都是非常好的。是否可以使用 cartopy 创建这些坐标列表？我正在寻找一个生成静态 html 文件的全自动脚本。

我有一些卫星图像数据想用 Cartopy 显示。我已成功遵循此处详述的图像示例。导致此代码：

此代码生成以下图像

我的问题是卫星图像数据不在 PlateCarree 投影中，而是在墨卡托投影中。

但是当我得到轴对象时

我失去了海岸线。

我看到了这里报告的问题。但

结果如下图：

数据不存在，圣地亚哥位于错误的位置。纬度/经度范围也发生了变化。我究竟做错了什么？

讨论后更新

主要问题是我没有使用该transform_points方法正确指定目标投影中的图像范围。imshow正如菲尔建议的那样，我还必须具体说明该方法中的坐标参考系。这是正确的代码：

产生这个正确的地理投影卫星图像：

我正在尝试绘制一些海冰数据图。数据在 EASE-North 网格中提供，示例文件 (HDF4) 可在以下位置下载：

我为 EASE-Grid 创建了一个自定义投影类，它似乎正在工作（海岸线与数据很好地对齐）。

当我尝试添加自然地球功能时，它返回一个空的 Matplotlib 图。

上面的脚本工作正常并产生以下结果：

但是当我取消注释时ax.add_feature(land)它会失败而没有任何错误，只返回空图。我错过了一些明显的东西吗？

这是 IPython 笔记本： http: //nbviewer.ipython.org/6779935

我的 Cartopy 版本是 Christoph Gohlke 网站上的 0.9 版（谢谢！）。

编辑：

尝试保存该图确实会引发异常：

检查“土地”cartopy.feature没有发现任何问题，所有多边形都通过了.isvalid()，所有环（ext en int）都是 4 个或更多元组。所以输入形状似乎不是问题（并且在 PlateCaree() 中工作正常）。

也许一些环（如南半球）在转换为 EASE_North 后会“损坏”？

编辑2：

当我删除内置的 NE 功能并加载相同的 shapefile（但剪裁低于 40N 的任何内容）时，它可以工作。所以这似乎是某种重投影问题。