Shapely: 3차원 드롭 기능 추가

3D 지오메트리를 2D로 변환하는 것에 대한 GIS Stack Exchange에 대한 오래된 질문이 있습니다: Convert 3D WKT to 2D Shapely Geometry . 이 기능은 Shapely에 포함되어야 하므로 다음과 같이 사용할 수 있습니다.

>>> from shapely.geometry import Polygon
>>> p = Polygon([(0, 0, 0), (1, 0, 0), (1, 1, 0)])
>>> p.wkt
'POLYGON Z ((0 0 0, 1 0 0, 1 1 0, 0 0 0))'
>>> p2 = p.drop_z
>>> p2.wkt
'POLYGON ((0 0, 1 0, 1 1, 0 0))'

나는 답변 중 하나와 문서에서 3차원 이 _기하학적 분석에 영향을 미치지 않기 때문에 _이 작업이 필요하지 않음 _을 보았습니다. 그러나 분할 지오메트리의 왼쪽을 결정하기 위해 자체 기능을 구현할 때(아직 구현되지 않음: https://github.com/Toblerity/Shapely/issues/589) 이 기능이 유용할 수 있음을 확인했습니다.

from shapely.geometry import (LinearRing, 
                              LineString, 
                              Polygon)


def is_left(polygon: Polygon,
            line: LineString) -> bool:
    """
    Determines if the polygon is on the left side of the line
    according to:
    https://stackoverflow.com/questions/50393718/determine-the-left-and-right-side-of-a-split-shapely-geometry
    """
    ring = LinearRing([*line.coords, *polygon.centroid.coords])
    return ring.is_ccw

이 코드는 3D 지오메트리에 실패합니다.

p = Polygon([(0, 0, 0), (1, 0, 0), (1, 1, 0)])
l = LineString([(0, 0, 0), (1, 0, 0)])
is_left(p, l)

이 오류가 발생합니다.

AttributeError                            Traceback (most recent call last)
~/miniconda3/lib/python3.7/site-packages/shapely/speedups/_speedups.pyx in shapely.speedups._speedups.geos_linearring_from_py()

AttributeError: 'list' object has no attribute '__array_interface__'

During handling of the above exception, another exception occurred:

IndexError                                Traceback (most recent call last)
<ipython-input-55-555b9e2533fa> in <module>()
----> 1 is_left(p, l)

<ipython-input-52-7fad75b19ce3> in is_left(polygon, line)
      6     https://stackoverflow.com/questions/50393718/determine-the-left-and-right-side-of-a-split-shapely-geometry
      7     """
----> 8     ring = LinearRing([*line.coords, *polygon.centroid.coords])
      9     return ring.is_ccw

~/miniconda3/lib/python3.7/site-packages/shapely/geometry/polygon.py in __init__(self, coordinates)
     51         BaseGeometry.__init__(self)
     52         if coordinates is not None:
---> 53             self._set_coords(coordinates)
     54 
     55     <strong i="18">@property</strong>

~/miniconda3/lib/python3.7/site-packages/shapely/geometry/polygon.py in _set_coords(self, coordinates)
     66     def _set_coords(self, coordinates):
     67         self.empty()
---> 68         ret = geos_linearring_from_py(coordinates)
     69         if ret is not None:
     70             self._geom, self._ndim = ret

~/miniconda3/lib/python3.7/site-packages/shapely/speedups/_speedups.pyx in shapely.speedups._speedups.geos_linearring_from_py()

IndexError: tuple index out of range

이는 Polygon 의 중심이 항상 2D로 반환되고(https://github.com/Toblerity/Shapely/issues/554) LinearRing 는 반환될 수 없기 때문입니다. 차원 수가 다른 점으로 구성됩니다.

drop_z 메소드가 있다면 line = LineString([xy[:2] for xy in list(line.coords)]) 같은 코드로 코드를 어지럽히거나 이에 대한 함수를 구현하는 대신 ring = LinearRing([*line.drop_z.coords, *polygon.centroid.coords]) 만 작성합니다. 또는 더 나은 방법으로 최상위 수준의 파일에서 읽은 원래 상위 다각형의 0으로만 구성된 중복 3차원을 삭제하여 모든 하위 도형이 2차원만 갖도록 합니다.

Shapely 버전: 1.6.4.post1, conda에서 설치됨.