पोस्टजीआईएस में अन्य लाइनों के लिए नोड्स की शुरुआत और अंत में तड़क-भड़क


9

ऐसे बहुत से उदाहरण हैं जो बताते हैं कि लाइनों को बिंदुओं पर कैसे स्नैप किया जाता है, लेकिन मैं अन्य लाइनों के नोड्स के लिए लाइन स्ट्रिंग्स के प्रारंभ और अंत नोड्स को स्नैप करने का कोई भी (तेज़!) तरीका नहीं खोज पाया।

अनिवार्य रूप से मैं पोस्टगीस (2.0) में अपनी परत को "साफ" करना चाहता हूं, लगभग समान बिंदुओं को एक साथ ले जाना और लाइन स्ट्रिंग्स के बीच छोटे उद्घाटन को सिलाई करना।

इससे कोई फर्क नहीं पड़ता कि मैं एक और नोड जोड़ता हूं, या तो लाइन के पहले / अंतिम नोड को स्थानांतरित करता हूं, या दोनों बिंदुओं को केंद्र में ले जाता हूं।

मुझे दो विकल्प मिले हैं, लेकिन मुझे यकीन नहीं है कि दोनों में से कैसे शुरू करें:

दूसरा विकल्प संभव लगता है, लेकिन इस विधि का पालन करने के बारे में कोई मदद बहुत सराहना की जाएगी।

जवाबों:


6

मैंने इसका उल्लेख करने के लिए प्रबंधित किया है, बिना वर्णित GRASS टूल्स या टोपोलॉजिकल फ़ंक्शन का उपयोग किए।

मूल रूप से मैं सभी स्टार्ट-एंड एंडोड्स लेती हूं, उन्हें एक नई, अस्थायी तालिका में डालती हूं, उनके चारों ओर एक बफर डालती हूं, बफर ऑब्जेक्ट्स को यूनियन करती हूं, और प्रत्येक बफर में सभी पाए गए नोड्स को बफर के सेंट्रोइड पर ले जाती हूं।

जब मैंने ऐसा किया है तो मैं नए स्थान पर मूल आरंभ और अंत बिंदुओं को स्थानांतरित करता हूं।

उम्मीद की तुलना में आसान है, और अभी भी तेज है, लेकिन मुझे उम्मीद थी कि PostGIS के पास इसके लिए कुछ अंतर्निहित फ़ंक्शन होंगे - जो कि और भी तेज़ होगा।

संपादित करें: समुदाय को वापस देने के हित में, यह अब के लिए मेरा (काफी भद्दा) कोड है।

drop table if exists nodes;
drop table if exists nodes2;
drop table if exists buffers;

-- Get Start and End nodes
select ST_StartPoint(wkb_geometry) startnode,  ST_EndPoint(wkb_geometry) endnode,    ogc_fid into nodes  from sourceTable;
-- Combine all nodes into one table for easier queries
select startnode node, ogc_fid into nodes2 from nodes;
insert into nodes2 select endnode node, ogc_fid from nodes;

-- Some indexes to speed everything up
CREATE INDEX nodesstart_idx ON nodes USING gist  (startnode);
CREATE INDEX nodesend_idx ON nodes USING gist  (endnode);
CREATE INDEX nodes2_idx ON nodes2 USING gist  (node);
CREATE INDEX nodes_ogcfid_idx ON nodes USING btree (ogc_fid ASC NULLS LAST);

-- Create buffers, combine them, split combined objects again
select (ST_Dump(ST_Union(ST_Buffer(node, 1)))).geom geom into buffers from nodes2;
CREATE INDEX buffers_idx ON buffers USING gist  (geom);

-- Update start/end nodes table
UPDATE nodes SET startnode = ST_Centroid((select geom from buffers WHERE geom && startnode));
UPDATE nodes SET endnode = ST_Centroid((select geom from buffers WHERE geom && endnode));
-- Update original points
update sourceTable set wkb_geometry = ST_SetPoint(
ST_SetPoint(wkb_geometry, 0, (select startnode from nodes where ogc_fid=sourceTable.ogc_fid)), 
ST_NumPoints(wkb_geometry) - 1, (select endnode from nodes where ogc_fid=sourceTable.ogc_fid));

DROP TABLE nodes;
DROP TABLE nodes2;
DROP TABLE buffers;

यह जवाब मेरे जवाब से "गैर-टोपोलॉजिकल" सुझाव जैसा लगता है। यह अच्छा होगा यदि आपने उत्तर दिया या उत्तर का चयन किया। वे हैं जो यहाँ समुदाय को खिलाते हैं :)
katahdin

आप सही हे। मैंने आपके उत्तर को बदल दिया है, और अपनी प्रतिक्रिया को अपने कोड को शामिल करने के लिए संपादित करूंगा।
जेलर बैस

4

यहाँ तीन विकल्प हैं। उम्मीद है कि कोई मदद करेगा।

v.clean

QGIS में GRASS टूल्स का उपयोग करके आप एक स्थानिक वस्तु की टोपोलॉजी को साफ कर सकते हैं। उपयोगकर्ता @RK एक अलग प्रश्न के उत्तर में यह करने के लिए निर्देशों का एक अच्छा सेट देता है । GRASS को जो फायदा होता है, वह है शेपफाइल के टोपोलॉजी को। आपकी स्थिति के लिए नुकसान यह है कि आपका डेटा एक आकार में नहीं है। बेशक आप Postgres से डेटा को "AddGGIS लेयर" टूल का उपयोग करके एक आकृति में निर्यात कर सकते हैं, लेकिन यह एक अतिरिक्त कदम है।

गैर-सामयिक पोस्टगिस फ़ंक्शन

PostGIS में आप एक linestring के लिए अंत और प्रारंभ बिंदु प्राप्त करने के लिए ST_EndPoint और ST_StartPoint फ़ंक्शन का उपयोग कर सकते हैं । फिर, ST_DWithi n और ST_Distance के संयोजन का उपयोग करके , आप नजदीकी ज्यामिति पर निकटतम प्रारंभ या समाप्ति बिंदु पा सकते हैं। यदि आपके पास बहुत सारे बिंदु हैं तो ST_DWithin क्वेरी को बहुत तेज़ी से बढ़ाएगा - यह मानते हुए कि आपके पास एक सूचकांक है। वहां से आपको एक नियम स्थापित करना होगा जो परिभाषित करता है कि किन बिंदुओं को संशोधित किया जाना है, और जो तय हो गए हैं।

यहां लाभ यह है कि आपको अपना डेटा सफाई के लिए GRASS को नहीं भेजना है, लेकिन बाहर देखने के लिए कुछ नुकसान हैं।

सामयिक पोस्टगिस कार्य

प्रश्न PostGIS सामयिक कार्यों को संदर्भित करता है। ये महान काम करते हैं, लेकिन, जैसा कि विकी का वर्णन है , आपको किनारों, नोड्स और चेहरों को स्पष्ट रूप से परिभाषित करना होगा। स्पष्ट रूप से यह आपके डेटा सेट के लिए एक समस्या होगी क्योंकि आप टोपोलॉजी के साथ समस्याओं को जानते हैं।


1

PostGIS में तड़क-भड़क वाले कार्य होते हैं .. शायद वे मदद करेंगे?

ST_Snap: स्नैप सेगमेंट और इनपुट ज्यामिति के कोने एक संदर्भ ज्यामिति के कोने तक।

ST_SnapToGrid: इनपुट ज्यामिति के सभी बिंदुओं को एक नियमित ग्रिड में स्नैप करें।


1
धन्यवाद, मैं इन कार्यों से अवगत हूँ। ST_Snap सभी नोड्स को स्नैप करता है , मैं केवल प्रारंभ और समाप्ति नोड चाहता हूं। ST_SnapToGrid वास्तव में अनुकूल नहीं है क्योंकि यह सभी मौजूदा ज्यामिति को संशोधित करता है, और उन नोड्स को स्थानांतरित करने का एक मौका होता है जो आगे दूर हैं, क्योंकि वे अभी मुश्किल से दूसरे खंड में आते हैं।
जैलमेर बास
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.