Shapely: तीसरा आयाम छोड़ने के लिए कार्यक्षमता जोड़ें

3D ज्यामिति को 2D में परिवर्तित करने के बारे में GIS स्टैक एक्सचेंज पर एक पुराना प्रश्न है: 3D WKT को 2D शेपली ज्योमेट्री में कनवर्ट करें । मुझे लगता है कि इस कार्यक्षमता को शेपली में शामिल किया जाना चाहिए, ताकि हम इसका उपयोग कर सकें, उदाहरण के लिए, जैसे:

>>> 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))'

मैंने एक उत्तर में और डॉक्स में देखा है कि _यह ऑपरेशन आवश्यक नहीं है_ क्योंकि तीसरा आयाम _ज्यामितीय विश्लेषण पर कोई प्रभाव नहीं डालता है । लेकिन एक विभाजित ज्यामिति के बाईं ओर निर्धारित करने के लिए अपने स्वयं के कार्य को लागू करते समय (यह अभी तक लागू नहीं हुआ है: 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]) लिखता। या इससे भी बेहतर, मैं मूल माता-पिता बहुभुज से केवल शून्य से मिलकर तीसरे आयाम को छोड़ दूंगा, जिसे मैंने शीर्ष स्तर पर एक फ़ाइल से पढ़ा है, इसलिए सभी बाल ज्यामिति में केवल 2 आयाम होंगे।

सुडौल संस्करण: 1.6.4.post1, कोंडा से स्थापित।