Shapely: ops.split लाइनस्ट्रिंग्स और मल्टीलाइनस्ट्रिंग्स के साथ असंगत व्यवहार

ठीक। मैं यहां समस्या देखता हूं।

नेत्रहीन, वे लिनेस्ट्रिंग इस तरह दिखते हैं:

जब आप क्षैतिज ( ए ) को शीर्ष पर बड़े लंबवत ( बी ) से विभाजित करते हैं जो केवल ए को छूता है , कुछ भी नहीं होता है। लेकिन A को C से विभाजित करने से काम चल C, A को पार कर जाता है ।

अब, जब आप A को MultiLineString([B, C]) विभाजित करते हैं, तो यह तीन लिनेस्ट्रिंग में विभाजित हो जाता है, दो नहीं जैसा कि आप उम्मीद करेंगे। मैं मानता हूं कि यह असंगत लगता है।

इसका कारण बनने वाला कोड वास्तव में बहुत सुलभ है ( Permalink )।

    <strong i="26">@staticmethod</strong>
    def _split_line_with_line(line, splitter):
        """Split a LineString with another (Multi)LineString or (Multi)Polygon"""

        # if splitter is a polygon, pick it's boundary
        if splitter.type in ('Polygon', 'MultiPolygon'):
            splitter = splitter.boundary

        assert(isinstance(line, LineString))
        assert(isinstance(splitter, LineString) or isinstance(splitter, MultiLineString))

        if splitter.crosses(line):
            # The lines cross --> return multilinestring from the split
            return line.difference(splitter)
        elif splitter.relate_pattern(line, '1********'):
            # The lines overlap at some segment (linear intersection of interiors)
            raise ValueError('Input geometry segment overlaps with the splitter.')
        else:
            # The lines do not cross --> return collection with identity line
            return [line]

आप देख सकते हैं कि किसी भी विभाजन को शुरू करने के लिए splitter को line को पार करना आवश्यक है। लेकिन एक बार जब वह चेकपॉइंट पास हो जाता है, तो विभाजित ज्यामिति उत्पन्न करने के लिए difference एक साधारण कॉल की जाती है। लेकिन difference पार करने और केवल छूने के बीच भेदभाव नहीं करता है! इसलिए यदि किसी भी बिंदु पर splitter line पार कर जाता है, तो सभी केवल स्पर्श करने वाले बिंदुओं को विभाजित कर दिया जाएगा।

@ यिर्मयाह-इंग्लैंड इसमें खुदाई करने के लिए धन्यवाद! ऐसा लगता है, कम से कम लाइनों के लिए, हम क्रॉस टेस्ट को टच टेस्ट में बदल सकते हैं और कम आश्चर्यजनक परिणाम प्राप्त कर सकते हैं? तुम उसके बारे में क्या सोचते हो?

मुझे लगता है कि यह ठीक काम करेगा, हाँ। और व्यवहार मुझे समझ में आता है।

थोड़ा के लिए, मैंने सोचा था कि हम एक MultiLineString में goems फ़िल्टर कर सकते हैं splitter पार करने वालों के लिए line और केवल उन्हीं के अंतर से लेते हैं। लेकिन यह उस मामले को हल नहीं करेगा जहां एक लाइनस्ट्रिंग splitter line को छूने के लिए वापस हुक करता है जो इसे दूसरे बिंदु पर पार करता है। मुझे लगता है कि भारी भारोत्तोलन के लिए difference का उपयोग करते समय लगातार व्यवहार प्राप्त करने का एकमात्र तरीका स्पर्श रैखिक ज्यामिति को विभाजित करना है।

मुझे लगता है कि हो सकता है crosses इतना समारोह जल्दी लौट आए जब केवल की सीमा इस्तेमाल किया गया था line छू लेती splitter । दरअसल इस मामले की जांच की जा रही है । लेकिन हमारे वर्तमान मामले में, जहां की सीमा splitter छूता line , परीक्षण में नहीं है, और एक ब्लाइंड स्पॉट हो सकता है। हालांकि, मैं अन्य साधारण फीचर पुस्तकालयों के साथ अधिक अनुभव के बिना एक सुडौल मूल निवासी हूं। और मुझे नहीं पता कि केवल क्रॉसिंग करते समय विभाजन रैखिक ज्यामिति के लिए मानक व्यवहार है या नहीं।

यह 5 साल पहले की बात है, लेकिन अगर @georgeouzou को याद है तो शायद वह पुष्टि कर सकता है या समझा सकता है कि crosses का उपयोग क्यों किया गया था?

हम अपना केक ले सकते हैं और इसे खा भी सकते हैं यदि हम परीक्षण किए गए मामले में जल्दी लौटना चाहते हैं, लेकिन इस मुद्दे के मामले में बंटवारे के साथ जारी रखें। इसके लिए सीधे relate का उपयोग करना होगा। लेकिन मुझे संदेह है कि किसी भी तरह से एक बड़ा प्रदर्शन अंतर होगा।

@ यिर्मयाह-इंग्लैंड मेरी याददाश्त इस पर थोड़ी धुंधली है, लेकिन कोड को फिर से देखकर मैं पुष्टि कर सकता हूं कि यह एक अंधा स्थान था। वर्तमान मामले को पूरा करने के लिए कोड को भी परीक्षण करना चाहिए। क्या आप परिवर्तनों के साथ जारी रख सकते हैं?

अरे @georgeouzou , इसे देखने के लिए समय निकालने के लिए धन्यवाद!

और, हाँ, अगर यह सभी के साथ ठीक है, तो मैं इस मामले को सुलझाना चाहता हूँ और किसी बिंदु पर एक पीआर जमा करना चाहता हूँ।