纸箱中的溢出错误

时间:2015-09-23 17:01:08

标签: python shapely cartopy

我正在使用cartopy绘制一些地图。在某些情况下,当我在我的轴上调用.set_extent()时,我收到此错误:

Traceback (most recent call last):
  File "<pyshell#315>", line 1, in <module>
    ax.set_extent([bounds.X1.min(), bounds.X2.max(), bounds.Y1.min(), bounds.Y2.max()], cartopy.crs.AlbersEqualArea())
  File "C:\FakeProgs\Python27\lib\site-packages\cartopy\mpl\geoaxes.py", line 587, in set_extent
    projected = self.projection.project_geometry(domain_in_crs, crs)
  File "C:\FakeProgs\Python27\lib\site-packages\cartopy\crs.py", line 172, in project_geometry
    return getattr(self, method_name)(geometry, src_crs)
  File "C:\FakeProgs\Python27\lib\site-packages\cartopy\crs.py", line 178, in _project_line_string
    return cartopy.trace.project_linear(geometry, src_crs, self)
  File "lib\cartopy\trace.pyx", line 109, in cartopy.trace.project_linear (lib/cartopy\trace.cpp:1135)
  File "lib\cartopy\trace.pyx", line 71, in cartopy.trace.geos_from_shapely (lib/cartopy\trace.cpp:838)
OverflowError: Python int too large to convert to C long

事情是这种行为有点随机。不是每次拨打.set_extent()都会这样做。这是一个解释器会话的摘录(bounds是一个pandas DataFrame,它保存了我打算稍后添加到轴上的各种形状的边界框的坐标。)

>>> ax = pyplot.axes(projection=cartopy.crs.AlbersEqualArea())
... ax.set_extent([bounds.X1.min(), bounds.X2.max(), bounds.Y1.min(), bounds.Y2.max()], cartopy.crs.AlbersEqualArea())
# result is exception shown above
>>> [bounds.X1.min(), bounds.X2.max(), bounds.Y1.min(), bounds.Y2.max()]
[-2218681.0391451684,
 -2103178.2838086924,
 -195096.93292225525,
 7468.2970529943705]
>>> [int(x) for x in [bounds.X1.min(), bounds.X2.max(), bounds.Y1.min(), bounds.Y2.max()]]
[-2218681, -2103178, -195096, 7468]
>>> [long(x) for x in [bounds.X1.min(), bounds.X2.max(), bounds.Y1.min(), bounds.Y2.max()]]
[-2218681L, -2103178L, -195096L, 7468L]
>>> ax = pyplot.axes(projection=cartopy.crs.AlbersEqualArea())
... ax.set_extent([bounds.X1.min(), bounds.X2.max(), bounds.Y1.min(), bounds.Y2.max()], cartopy.crs.AlbersEqualArea())
# works without problem!

相同的代码可以工作,而不会更改其间的任何变量。

trace,pyx中的这一行似乎引发了错误:

cdef ptr geos_geom = shapely_geom._geom

我做了一些搜索,发现an old commitsome mailing list上提出的类似问题有关。

我对这个问题的理解是这些Shapely对象的_geom属性存储了某种C库中某个对象的指针。如果此指针的整数值对于C long而言太大,则会引发错误。该错误无法重现,因为每次创建新的_geom时都会创建一个新的GeoAxes,而新的_geom可能会或可能不会太大。

但是,令人费解的是,我能找到的关于此错误的大部分内容(例如,上述提交中的提交消息)表明它应该只是32位系统的问题,但我正在使用所有库的64位版本的64位Python 2.7。

所以我的问题是:我对正在发生的事情是对的吗?如果是这样,为什么我仍然在64位系统上出现这些错误?有没有办法解决它?

1 个答案:

答案 0 :(得分:1)

  我正确的发生了什么事吗?

我无法断然确认你是对的,但它看起来似乎有道理。我之前从未见过这个问题,但同样我也不会经常在Windows上使用cartopy。

  

如果是这样,为什么我仍然在64位系统上出现这些错误?

您的计算机可能是64位,但是您使用的是64位的Python吗?

  

有没有办法解决它?

鉴于它似乎是随机的,解决方法可能是:

for attempt in range(10):
    try:
        ...
    except OverflowError:
        print('Failed attempt {}, retrying upto 10 times.'.format(attempt))

它当然不是很漂亮,但可能是此时解决方法的唯一方法。

显然你发现的是一个错误,所以我认为cartopy issue tracker是找到问题的长期解决方案的正确位置。我认为提供您正在使用的软件版本是个好主意,理想情况下您找到的坐标会触发问题(即使是随机的)。

HTH