matplotlib / matplotlib/basemap

Basemap creation fails for map covering pole

未关闭
#161 0 条评论 0 个 reaction 已指派 0 人 在 GitHub 查看

还没有人认领这个 Issue。

主要语言
Python
星标
817
派生
395
PR 合并指标
30 天内没有已合并 PR

描述

``` python
>>> Basemap(width=1000000, height=1000000, projection='stere', lat_0=-88, lon_0=-3)

Traceback (most recent call last):
File "", line 1, in
File "..\site-packages\mpl_toolkits\basemap\__init__.py", line 1085, in __init__
self._readboundarydata('gshhs',as_polygons=True)
File "..\site-packages\mpl_toolkits\basemap\__init__.py", line 1428, in _readboundarydata
if not poly.is_valid(): poly=poly.fix()
File "_geoslib.pyx", line 234, in _geoslib.BaseGeometry.fix (src/_geoslib.c:2233)
File "_geoslib.pyx", line 334, in _geoslib.Polygon.__init__ (src/_geoslib.c:3479)
IndexError: index -1 is out of bounds for axis 0 with size 0
```

Strange thing is that e.g. using `lon_0=-2` works.

I'm using basemap 1.0.8.

贡献指南

这个仓库没有索引到贡献指南

从这里开始

  1. 先读完整个 Issue,再读项目的贡献指南。
  2. 在 Issue 下留言说明你要接手 —— 这能避免两个人做同样的事。
  3. Fork 仓库,在一个分支上完成修改。
  4. 提交 Pull Request,并在描述里引用这个 Issue 编号。

调研方向

复现使用 projection='stere'、lat_0=-88 和 lon_0=-3 的确切 Basemap 调用,然后将其与正常工作的 lon_0=-2 情况进行比较。从 mpl_toolkits/basemap/__init__.py 中的 _readboundarydata 开始,跟踪该失败直到 _geoslib.pyx 中的 BaseGeometry.fix。完成标准是覆盖极点的地图不再引发 IndexError。

由索引模型根据 Issue 内容生成。

评估

技术栈
python
领域
data-visualization
Issue 类型
缺陷
难度
3/5
预计耗时
1-2 天
活跃度
停滞
描述清晰度
基本清楚
新手友好度
42/100

把新 issue 发到你的邮箱

精选适合新手参与的 GitHub issue 摘要。