Блог

Главная - Вопросы по программированию - Ошибка рекурсии: максимальная глубина рекурсии превышена по сравнению с работой с пандами

Ошибка рекурсии: максимальная глубина рекурсии превышена по сравнению с работой с пандами

#python #pandas #geopandas #shapely

Вопрос:

Я работаю с некоторыми геоданными. Я хочу получить точечный объект по широте и долготе, которые у меня уже есть. Мой код выглядит следующим образом:

 for i in range(len(map_polygon)):
    x = map_polygon.at[i, 'lon'] 
    y = map_polygon.at[i, 'lat']
    map_polygon.at[i, 'point'] = Point(x, y)

Однако я получаю следующую ошибку:
Ошибка рекурсии: максимальная глубина рекурсии превышена по сравнению с

Что я делаю не так? Длина фрейма данных составляет всего около 20 тысяч строк.

Спасибо!

PS Я прикрепил все сообщение об ошибке ниже:

 map_polygon.at[i, 'point'] = Point(x, y)
  File "/Library/Frameworks/Python.framework/Versions/3.8/lib/python3.8/site-packages/pandas/core/indexing.py", line 2167, in __setitem__
    return super().__setitem__(key, value)
  File "/Library/Frameworks/Python.framework/Versions/3.8/lib/python3.8/site-packages/pandas/core/indexing.py", line 2118, in __setitem__
    self.obj._set_value(*key, value=value, takeable=self._takeable)
  File "/Library/Frameworks/Python.framework/Versions/3.8/lib/python3.8/site-packages/pandas/core/frame.py", line 3278, in _set_value
    self.loc[index, col] = value
  File "/Library/Frameworks/Python.framework/Versions/3.8/lib/python3.8/site-packages/pandas/core/indexing.py", line 692, in __setitem__
    iloc._setitem_with_indexer(indexer, value, self.name)
  File "/Library/Frameworks/Python.framework/Versions/3.8/lib/python3.8/site-packages/pandas/core/indexing.py", line 1635, in _setitem_with_indexer
    self._setitem_with_indexer_split_path(indexer, value, name)
  File "/Library/Frameworks/Python.framework/Versions/3.8/lib/python3.8/site-packages/pandas/core/indexing.py", line 1720, in _setitem_with_indexer_split_path
    self._setitem_single_column(loc, value, pi)
  File "/Library/Frameworks/Python.framework/Versions/3.8/lib/python3.8/site-packages/pandas/core/indexing.py", line 1813, in _setitem_single_column
    ser._mgr = ser._mgr.setitem(indexer=(pi,), value=value)
  File "/Library/Frameworks/Python.framework/Versions/3.8/lib/python3.8/site-packages/pandas/core/internals/managers.py", line 568, in setitem
    return self.apply("setitem", indexer=indexer, value=value)
  File "/Library/Frameworks/Python.framework/Versions/3.8/lib/python3.8/site-packages/pandas/core/internals/managers.py", line 427, in apply
    applied = getattr(b, f)(**kwargs)
  File "/Library/Frameworks/Python.framework/Versions/3.8/lib/python3.8/site-packages/pandas/core/internals/blocks.py", line 967, in setitem
    return self.astype(dtype).setitem(indexer, value)
  File "/Library/Frameworks/Python.framework/Versions/3.8/lib/python3.8/site-packages/pandas/core/internals/blocks.py", line 967, in setitem
    return self.astype(dtype).setitem(indexer, value)
  File "/Library/Frameworks/Python.framework/Versions/3.8/lib/python3.8/site-packages/pandas/core/internals/blocks.py", line 967, in setitem
    return self.astype(dtype).setitem(indexer, value)
  [Previous line repeated 977 more times]
  File "/Library/Frameworks/Python.framework/Versions/3.8/lib/python3.8/site-packages/pandas/core/internals/blocks.py", line 966, in setitem
    dtype, _ = maybe_promote(np.array(value).dtype)
  File "/Library/Frameworks/Python.framework/Versions/3.8/lib/python3.8/site-packages/shapely/geometry/point.py", line 111, in array_interface
    if self.is_empty:
  File "/Library/Frameworks/Python.framework/Versions/3.8/lib/python3.8/site-packages/shapely/geometry/base.py", line 699, in is_empty
    return (self._geom is None) or bool(self.impl['is_empty'](self))
  File "/Library/Frameworks/Python.framework/Versions/3.8/lib/python3.8/site-packages/shapely/predicates.py", line 25, in __call__
    return self.fn(this._geom)
  File "/Library/Frameworks/Python.framework/Versions/3.8/lib/python3.8/site-packages/shapely/geos.py", line 583, in errcheck_predicate
    if result == 2:
RecursionError: maximum recursion depth exceeded in comparison

Ответ №1:

При добавлении столбца геометрии зацикливание и установка каждого значения по отдельности обычно является плохой идеей (медленной), поэтому вы, вероятно, захотите сделать это следующим образом:

 map_polygon['point'] = [Point(p) for p in zip(df.lon, df.lat)]

Комментарии:

1. В примере OP Point(...) принимает 2 аргумента, здесь вы дали ему один аргумент в виде кортежа. Вероятно, должно быть [Point(*p) for p in zip(df.lon, df.lat)] или [Point(lon, lat) for lon, lat in zip(df.lon, df.lat)]

2. Согласно документации шейпли, Точка принимает либо аргументы кортежа, либо позиционные аргументы: «Конструктор точек принимает значения координат положения или параметры точечного кортежа. «источник = shapely.readthedocs.io/en/stable/manual.html#points Так что и то, и другое работает, и передача кортежа означает меньше ввода

Метки: