ArcMap में लाइन शेपफाइल की "कनेक्टिविटी" की जाँच की जा रही है?

Question

मुझे ArcMap में एक लाइन शेपफाइल विरासत में मिली है जिसका उपयोग सड़क नेटवर्क का प्रतिनिधित्व करने के लिए किया जाता है। समस्या यह है कि सड़क नेटवर्क बहुत बड़ा है और यह मेरे लिए असंभव है:

1. यह बताएं कि क्या प्रत्येक रेखा का अंत (सड़क खंड) एक किनारे, शीर्ष या दूसरी पंक्ति (सड़क खंड) के अंत में "जुड़ा हुआ" है; तथा
2. यह बताएं कि क्या पूरा रोड नेटवर्क "कनेक्टेड" है - यानी, क्या मैं रोड नेटवर्क पर किसी भी बिंदु से सड़क नेटवर्क पर किसी अन्य बिंदु पर एक मार्ग का पता लगा सकता हूं?

वहाँ ArcMap में यह पूरा करने के लिए एक रास्ता है? यही है, क्या ऐसा कोई फ़ंक्शन है जो मुझे लाइन सेगमेंट को उजागर करने की अनुमति देगा जो अन्य लाइन सेगमेंट से कनेक्ट नहीं होता है, या कुछ इसी तरह से जो मुझे ऊपर दिए गए दो प्रश्नों का उत्तर देने की अनुमति देगा?

Answer 1

हाँ, लेकिन की तरह। आर्कगिस में अब लाइन-नोड टोपोलॉजी नहीं है जो उपयोगकर्ता को यह बताने में सक्षम करता है कि उनके छोर (नोड्स) पर कितने आर्क्स (लाइनें) जुड़े हुए हैं।

जांच करना एक बात है, लेकिन इसके बजाय कैसे ठीक करना है? यदि आप आर्कपेज़ में फीचर क्लास खोलते हैं और फिर प्लेनरीज़ लाइन्स (एक सहिष्णुता दें) का उपयोग करते हैं और लाइनों को काट दिया जाएगा और चौराहे पर विभाजित किया जाएगा - बहुत काम बचाता है। यदि आप लाइनों को विभाजित नहीं करना चाहते हैं, तो उपकरण को एकीकृत करें पर विचार करें, लेकिन बहुत छोटी सहिष्णुता का उपयोग करने के लिए सावधान रहें, यह एक साथ छोरों को स्नैप करेगा लेकिन लाइनों को एक साथ स्नैप भी कर सकता है। इंटीग्रेट का उपयोग करने से पहले एक बैकअप रखें क्योंकि यह आपके डेटा को नष्ट कर सकता है!

अब, डिस्कनेक्ट किए गए सिरों को खोजने के लिए पॉइंट पॉइंट्स को पाने के लिए फ़ीचर वर्टिक्स का उपयोग करें और फिर ईवेंट्स लीजिए, जो आपको एंड पॉइंट पॉइंट्स की संख्या के साथ एक फीचर क्लास देगा, इस स्तर पर 1 की कोई भी घटना संदिग्ध है इसलिए आपको अलग करने की आवश्यकता होगी ये बाहर।

यदि यह कनेक्ट होना चाहिए, तो यह काम करना अगला काम है, जेनरेट नियर टेबल (फिर से उपयुक्त सहिष्णुता के साथ) और निकटतम के विकल्प का उपयोग करें = सभी मूल लाइनों के खिलाफ 1 की गिनती के साथ घटनाओं का उपयोग करके, फिर सारांश सांख्यिकी का उपयोग करके आप पा सकते हैं प्रत्येक बिंदु के लिए एक केस फील्ड के रूप में IN_FID और NEAR_FID के रूप में "COUNT" के सांख्यिकीय प्रकार के साथ डेटा फ़ील्ड के रूप में रिकॉर्ड की गणना।

पास की मेज से आसान अर्क की तरह बनाने के लिए तालिका चयन का उपयोग करके 0 से अधिक दूरी के साथ रिकॉर्ड । प्रत्येक ईवेंट को वह लाइन मिल जाएगी जिसने उसे उत्पन्न किया था, लेकिन दूरी 0 होगी, अगर इसे दूसरी लाइन से ठीक से जोड़ा जाए (एक शीर्ष पर) तो दूरी भी 0 होगी, इसलिए अब कोई भी घटना जिसके पास तालिका में रिकॉर्ड शेष है संभवतः असंतुष्ट हैं, लेकिन इन्हें मैन्युअल रूप से देखने की आवश्यकता होगी।

Comments

धन्यवाद। जब मैं कल काम करने जाऊंगा और वापस रिपोर्ट करूंगा, तो मैं इन फंक्शन की जांच करूंगा।

— derNincompoop

"फीचर वर्टिस टू पॉइंट" में दंगल पॉइंट को आउटपुट करने का विकल्प है। मुझे लगता है कि यह आपके वर्कफ़्लो को सरल करेगा।

— क्लेविस

वाह @klewis, वह कब दिखाई दिया? यह निश्चित रूप से 9.3 में नहीं था। अच्छी तरह से देखा गया!

— माइकल स्टिम्सन

Answer 2

एक और दृष्टिकोण एमएपी टोपोलॉजी का उपयोग करना है। मैंने VBA कोड के इस बिट को खटखटाया जो कि झूलने वाले किनारों की पहचान करेगा। यदि आप नेटवर्क के अपेक्षित छोरों के बजाय नेटवर्क के भीतर झूलते हुए किनारों को देखते हैं तो डिस्कनेक्ट होना चाहिए।

कोड आपके द्वारा VBA स्थापित होने, संपादन मोड में होने पर निर्भर करता है और पॉलीलाइन लेयर को मैप टोपोलॉजी में जोड़ा जाता है।

Public Sub SelectDanglingPolylines()
 ' Description: Takes a polyline dataset and select all dangling polylines.
 '
 ' Requirements: You need to be in edit mode and have added the layer to a MAP TOPOLOGY,
 ' also polyline layer must be first in TOC.
 '
 ' Limitations: Large datasets take a long time to build the cache and may even fail.
 '
 ' Author: Duncan Hornby
 ' Created: 11/12/2011
 '

 ' Get map and then first layer, must be of polyline type
 Dim pMXDocument As IMxDocument
 Set pMXDocument = ThisDocument
 Dim pMap As IMap
 Set pMap = pMXDocument.FocusMap
 Dim pLayer As ILayer
 Set pLayer = pMap.Layer(0)
 Dim pFeatureLayer As IFeatureLayer
 Set pFeatureLayer = pLayer
 Dim pFeatureClass As IFeatureClass
 Set pFeatureClass = pFeatureLayer.FeatureClass
 If pFeatureClass.ShapeType <> esriGeometryPolyline Then
     MsgBox "This code works only with polylines!", vbExclamation, "Wrong data type at layer 0"
     Exit Sub
 End If

 ' Get editor and topology extension
 Dim pEditor As IEditor
 Dim pID As New UID
 Dim pTopologyExtension As ITopologyExtension
 Dim pTEID As New UID
 pID = "esriEditor.editor"
 Set pEditor = Application.FindExtensionByCLSID(pID)
 pTEID = "esriEditorExt.TopologyExtension"
 Set pTopologyExtension = Application.FindExtensionByCLSID(pTEID)
 If pTopologyExtension.CurrentTopology Is Nothing Then Exit Sub

 ' Get a MAP topology not a geodatabase topology
 Dim pMapTopology As IMapTopology
 If TypeOf pTopologyExtension.CurrentTopology Is IMapTopology Then
     Set pMapTopology = pTopologyExtension.MapTopology
 Else
     ' Not a Map Topology
     Exit Sub
 End If

 ' This is the colection that FID are added to
 Dim aColl As New Collection

 ' Build cache
 Application.StatusBar.Message(0) = "Building MAP TOPOLOGY cache, this can take a long time on large datasets..."
 DoEvents
 Dim pGeoDataset As IGeoDataset
 Set pGeoDataset = pFeatureClass
 Dim pEnvelope As IEnvelope
 Set pEnvelope = pGeoDataset.Extent
 pMapTopology.Cache.Build pEnvelope, False

 ' Identify dangling nodes and add polyline FID to collection
 Application.StatusBar.Message(0) = "Identifying dangling nodes..."
 DoEvents
 Dim pEnumTopologyParent As IEnumTopologyParent
 Dim pTopologyNode As ITopologyNode
 Dim pEnumTopologyNode As IEnumTopologyNode
 Set pEnumTopologyNode = pMapTopology.Cache.Nodes
 pEnumTopologyNode.Reset
 Set pTopologyNode = pEnumTopologyNode.Next
 While Not pTopologyNode Is Nothing
 If pTopologyNode.Degree = 1 Then
     ' As this has 1 degree it has only 1 parent polyline
     Set pEnumTopologyParent = pTopologyNode.Parents
     pEnumTopologyParent.Reset
     aColl.Add (pEnumTopologyParent.Next.m_FID) 'Adds polyline FID to collection
 End If
 Set pTopologyNode = pEnumTopologyNode.Next
 Wend

 ' Use collection to select polylines
 Application.StatusBar.Message(0) = "Selecting polylines..."
 DoEvents
 Dim pFeatureSelection As IFeatureSelection
 Set pFeatureSelection = pFeatureLayer
 Dim X As Variant
 For Each X In aColl
     pFeatureSelection.SelectionSet.Add CLng(X)
 Next
 pMXDocument.ActiveView.PartialRefresh esriViewGeoSelection, Nothing, pEnvelope
 Application.StatusBar.Message(0) = ""
 DoEvents
End Sub

Comments

मुझे वह कोड पसंद है, अगर आपको कोई आपत्ति नहीं है तो मैं उसके कुछ अंशों की नकल करूँगा; इसे इस तरह से करने का लाभ यह है कि यह लाइव है और जैसा कि त्रुटियों को ठीक किया जाता है वे गायब हो जाएंगे, खराब हिस्सा यह है कि यदि कोई रेखा लटक रही है तो यह दिखाया गया है - क्या एक छोर जुड़ा हुआ है या लाइन पूरी तरह से डिस्कनेक्ट हो गई है। मेरा सुझाव है कि आप कोड को बदल कर नोड "वेलेंस" (कितनी कनेक्टेड लाइन्स) या सिर्फ दंगल नोड्स प्रदर्शित करते हैं (यह मेरा वर्कस्टेशन फिर से आ रहा है!) ... वह भयानक होगा ... लेकिन केवल अगर ओपी! VBA / VB.net क्षमता है।

— माइकल स्टिम्सन

कोई भी कोड को मोड़ सकता है और किनारों को 1 की वैधता प्राप्त कर सकता है, जो पूरी तरह से डिस्कनेक्ट की गई रेखा को खोजने के लिए है, या 1 के अंत में 1 की तरह विविधताएं और दूसरे पर 2, मेरे द्वारा अनुमान किए गए नेटवर्क पर निर्भर करता है?

— हॉर्नबीड अग्ली

Answer 3

यह एक पुरानी पोस्ट है, लेकिन मुझे लगता है कि सबसे सरल उपाय यह है:

1. अपनी पॉलीलाइन सुविधा को भंग करें
2. दंगल विकल्प के साथ पॉइंट्स के लिए फ़ीचर वर्टिक्स का उपयोग करें
3. स्थानिक द्वारा जुड़ें मूल पॉलीलाइन सुविधा जिसके परिणामस्वरूप बिंदु परत है। "इसके द्वारा संचालित" विकल्प का उपयोग करें।

परिणाम में आपकी परत में प्रत्येक पंक्ति के लिए "काउंट" फ़ील्ड होगा। यदि गणना 1 से अधिक है, तो रेखा बाकी लाइनों से "कनेक्ट" नहीं है।

वैचारिक रूप से: चरण 2 यहां एकल कनेक्टेड एज (एक लाइन "में जा रहा है", शून्य "बाहर जा रहा है") के साथ कोने पर बिंदु बनाता है। चूँकि "कनेक्टेड" नेटवर्क के भीतर प्रत्येक लाइन में अधिकतम 1 ऐसा वर्टेक्स होगा, जिसमें 1 से अधिक वाली कोई भी लाइन नेटवर्क का हिस्सा नहीं है और इसलिए वह "कनेक्टेड" नहीं है।

Comments

यह विधि कई झूठी सकारात्मकता पैदा करती है। कहते हैं एक सड़क एक टी जंक्शन की ओर जाती है। T के शीर्ष भाग में दो लटकते हुए कोने होंगे। यदि आप दो की गिनती के साथ सुविधाओं को लेते हैं और फिर उन्हें सड़क नेटवर्क के साथ फिर से जोड़ते हैं, तो स्पर्श नहीं करने वाले द्वीप हैं।

— woot

दरअसल, मैंने जिस दूसरे इंटरसेक्ट के बारे में लिखा था, वह केवल अलग-थलग सड़कों को ही ढूंढता है। यह पता नहीं मिल रहा होता अगर नेटवर्क कई भागों में है

— woot

Answer 4

यहाँ एक तरीका है, जो मैंने कुछ दोस्तों की मदद से मॉडल बिल्डर और गेफी के साथ किया था। चरण 1 आर्कमॉडल लिंक / किनारों की तालिका बना रहा है (यदि वांछित है तो प्रत्येक पंक्ति पर नोड्स जोड़ना) चरण 2 लिंक / किनारों के Gephi आयात फिर घटक आईडी जोड़ना चरण 3 आर्कमॉडल घटक आईडी को मूल पंक्ति में वापस जोड़ता है

चरण 1 सभी इनपुट सुविधाओं को चौराहों पर विभाजित करता है ताकि यह सुनिश्चित किया जा सके कि नेटवर्क नोड मौजूद है और Gephi में आयात होने के लिए एक तालिका बनाता है। ये चरण हैं: प्रक्रिया: फ़ीचर टू लाइन (कई इनपुट ले सकते हैं) प्रक्रिया: मल्टीपार्ट टू सिंगलपार्ट प्रक्रिया: रिपेयर ज्योमेट्री प्रोसेस: एड स्टार्ट एंड कोऑर्डर्स (ज्योमेट्री विशेषताएँ जोड़ें) प्रक्रिया: फ़ील्ड जोड़ें "सोर्स" प्रोसेस: फील्ड "टारगेट" जोड़ें प्रक्रिया: "सोर्स" (स्टार्ट एक्सएंडवाईवाई के रूप में) की गणना करें प्रक्रिया: "टारगेट" (एक्स एंड एंडवाई के रूप में) की गणना करें प्रक्रिया: फ़ील्ड्स (भ्रम से बचने के लिए अतिरिक्त फ़ील्ड्स को हटाएं) प्रक्रिया: जीडीबी टेबल से सीएसवी

चरण 2 Gephi प्रक्रिया (मुफ्त डाउनलोड) -Sport और CSV आउटपुट स्रोत और लक्ष्य के साथ नोड फ़ील्ड्स लिंक के रूप में -अनुकूलित घटक के रूप में अप्रत्यक्ष (सांख्यिकीय उपकरण के तहत) के रूप में नोड्स के लिए डेटा प्रयोगशाला से निर्यात CSV (नोड आईडी और घटक आईडी शामिल)

चरण 3: Gephi आउटपुट लेता है और मूल लाइनों के लिए घटक / नेटवर्क विशेषता जोड़ता है प्रक्रिया: तालिका से तालिका (सबसे अनोखी ऑब्जेक्ट को संरचना करने के लिए Geodatabase में Gephi आउटपुट) प्रक्रिया: फ़ील्ड से जुड़ें (स्रोत नोड मान लाइनों को घटक संख्या को पॉप करने के लिए Gephi आउटपुट के साथ शामिल हुए) घटक द्वारा प्रतीक

इस बिंदु के बाद सफाई संभवतः जांच की एक मैन्युअल प्रक्रिया होगी जहां डिस्कनेक्ट लाइनों के साथ हो रहे हैं और यदि डिस्कनेक्ट एक वैध वास्तविक-दुनिया अलगाव या सिर्फ एक डेटा दोष है।