आप प्रत्येक बिंदु के निकटतम पड़ोसी का पता लगाने के लिए एक पुनरावर्ती क्वेरी का उपयोग कर सकते हैं जो आपके द्वारा बनाई गई लाइनों के प्रत्येक ज्ञात छोर से शुरू होती है।
पूर्वापेक्षाएँ : अपने बिंदुओं के साथ एक पोस्टगिस लेयर तैयार करें और दूसरा आपके सड़क पर एक एकल मल्टी-लिनस्ट्रिंग ऑब्जेक्ट के साथ। दो परतें एक ही CRS पर होनी चाहिए। यहाँ मेरे द्वारा बनाए गए परीक्षण डेटा-सेट के लिए कोड है, कृपया इसे आवश्यकतानुसार संशोधित करें। (पोस्ट 9.2 और पोस्टगिस 2.1 पर परीक्षण किया गया)
WITH RECURSIVE
points as (SELECT id, st_transform((st_dump(wkb_geometry)).geom,2154) as geom, my_comment as com FROM mypoints),
roads as (SELECT st_transform(ST_union(wkb_geometry),2154) as geom from highway),
यहाँ कदम हैं :
प्रत्येक बिंदु के लिए उत्पन्न करें प्रत्येक पड़ोसी की सूची और उनकी दूरी जो तीन मानदंडों को पूरा करती है।
- दूरी उपयोगकर्ता परिभाषित सीमा से अधिक नहीं होनी चाहिए (यह अलग-अलग बिंदुओं से जुड़ने से बच जाएगी)
graph_full as (
SELECT a.id, b.id as link_id, a.com, st_makeline(a.geom,b.geom) as geom, st_distance(a.geom,b.geom) as distance
FROM points a
LEFT JOIN points b ON a.id<>b.id
WHERE st_distance(a.geom,b.geom) <= 15
),
- सीधे रास्ते पर सड़क पार नहीं करनी चाहिए
graph as (
SELECt graph_full.*
FROM graph_full RIGHT JOIN
roads ON st_intersects(graph_full.geom,roads.geom) = false
),
दूरी निकटतम पड़ोसी से एक उपयोगकर्ता परिभाषित अनुपात से अधिक नहीं होनी चाहिए (यह तय दूरी से अनियमित डिजिटलाइजेशन के लिए बेहतर समायोजित करना चाहिए) यह हिस्सा वास्तव में लागू करने के लिए बहुत कठिन था, निश्चित खोज त्रिज्या से चिपका हुआ
चलो इस तालिका को "ग्राफ" कहते हैं
ग्राफ़ में शामिल होकर लाइन पॉइंट के अंत का चयन करें और केवल उस बिंदु को रखें, जिसमें ग्राफ़ में एक प्रविष्टि हो।
eol as (
SELECT points.* FROM
points JOIN
(SELECT id, count(*) FROM graph
GROUP BY id
HAVING count(*)= 1) sel
ON points.id = sel.id),
आइए इस तालिका को "ईओएल" (लाइन का अंत)
आसान कहते हैं? यह एक शानदार ग्राफ बनाने के लिए इनाम है, लेकिन अगले कदम पर चीजें पागल हो जाएंगी
एक पुनरावर्ती क्वेरी सेट करें जो प्रत्येक ईओल से शुरू होने वाले पड़ोसियों से पड़ोसियों तक साइकिल चलाएगी
- ईओल तालिका का उपयोग करके पुनरावर्ती क्वेरी की शुरुआत करें और गहराई के लिए एक काउंटर जोड़कर, मार्ग के लिए एक एग्रीगेटर और लाइनों के निर्माण के लिए एक ज्यामिति निर्माता
- ग्राफ़ का उपयोग करके निकटतम पड़ोसी पर स्विच करके अगले पुनरावृत्ति पर जाएं और जांचें कि आप पथ का उपयोग करके कभी भी पीछे नहीं जाते हैं
- पुनरावृत्ति समाप्त होने के बाद प्रत्येक प्रारंभिक बिंदु के लिए केवल सबसे लंबा रास्ता रखें (यदि आपके डेटासेट में संभावित लाइनों के बीच संभावित प्रतिच्छेदन शामिल है जो उस हिस्से को अधिक परिस्थितियों की आवश्यकता होगी)
recurse_eol (id, link_id, depth, path, start_id, geom) AS (--initialisation
SELECT id, link_id, depth, path, start_id, geom FROM (
SELECT eol.id, graph.link_id,1 as depth,
ARRAY[eol.id, graph.link_id] as path,
eol.id as start_id,
graph.geom as geom,
(row_number() OVER (PARTITION BY eol.id ORDER BY distance asc))=1 as test
FROM eol JOIn graph ON eol.id = graph.id
) foo
WHERE test = true
UNION ALL ---here start the recursive part
SELECT id, link_id, depth, path, start_id, geom FROM (
SELECT graph.id, graph.link_id, r.depth+1 as depth,
path || graph.link_id as path,
r.start_id,
ST_union(r.geom,graph.geom) as geom,
(row_number() OVER (PARTITION BY r.id ORDER BY distance asc))=1 as test
FROM recurse_eol r JOIN graph ON r.link_id = graph.id AND NOT graph.link_id = ANY(path)) foo
WHERE test = true AND depth < 1000), --this last line is a safe guard to stop recurring after 1000 run adapt it as needed
चलो इस तालिका को "recurse_eol" कहते हैं
प्रत्येक प्रारंभ बिंदु के लिए केवल सबसे लंबी लाइन रखें और हर सटीक डुप्लिकेट पथ को हटा दें उदाहरण: पथ 1,2,3,5 और 5,3,2,1 एक ही रेखा है जिसे दो लाइन "लाइन के अंत" द्वारा खोजा गया है।
result as (SELECT start_id, path, depth, geom FROM
(SELECT *,
row_number() OVER (PARTITION BY array(SELECT * FROM unnest(path) ORDER BY 1))=1 as test_duplicate,
(max(depth) OVER (PARTITION BY start_id))=depth as test_depth
FROM recurse_eol) foo
WHERE test_depth = true AND test_duplicate = true)
SELECT * FROM result
शेष त्रुटियों (अलग-अलग बिंदुओं, अतिव्यापी लाइनों, अजीब आकार की सड़क) को मैन्युअल रूप से जांचता है
जैसा कि वादा किया गया था, मुझे अभी भी यह पता नहीं चल सका है कि कभी-कभी एक ही पंक्ति के विपरीत ईओएल से शुरू करते समय पुनरावर्ती क्वेरी सटीक परिणाम क्यों नहीं देती है, इसलिए कुछ डुप्लिकेट अब परिणाम परत में बने रह सकते हैं।
बेझिझक पूछें मुझे पूरी तरह से लगता है कि इस कोड को अधिक टिप्पणियों की आवश्यकता है। यहाँ पूरी क्वेरी है:
WITH RECURSIVE
points as (SELECT id, st_transform((st_dump(wkb_geometry)).geom,2154) as geom, my_comment as com FROM mypoints),
roads as (SELECT st_transform(ST_union(wkb_geometry),2154) as geom from highway),
graph_full as (
SELECT a.id, b.id as link_id, a.com, st_makeline(a.geom,b.geom) as geom, st_distance(a.geom,b.geom) as distance
FROM points a
LEFT JOIN points b ON a.id<>b.id
WHERE st_distance(a.geom,b.geom) <= 15
),
graph as (
SELECt graph_full.*
FROM graph_full RIGHT JOIN
roads ON st_intersects(graph_full.geom,roads.geom) = false
),
eol as (
SELECT points.* FROM
points JOIN
(SELECT id, count(*) FROM graph
GROUP BY id
HAVING count(*)= 1) sel
ON points.id = sel.id),
recurse_eol (id, link_id, depth, path, start_id, geom) AS (
SELECT id, link_id, depth, path, start_id, geom FROM (
SELECT eol.id, graph.link_id,1 as depth,
ARRAY[eol.id, graph.link_id] as path,
eol.id as start_id,
graph.geom as geom,
(row_number() OVER (PARTITION BY eol.id ORDER BY distance asc))=1 as test
FROM eol JOIn graph ON eol.id = graph.id
) foo
WHERE test = true
UNION ALL
SELECT id, link_id, depth, path, start_id, geom FROM (
SELECT graph.id, graph.link_id, r.depth+1 as depth,
path || graph.link_id as path,
r.start_id,
ST_union(r.geom,graph.geom) as geom,
(row_number() OVER (PARTITION BY r.id ORDER BY distance asc))=1 as test
FROM recurse_eol r JOIN graph ON r.link_id = graph.id AND NOT graph.link_id = ANY(path)) foo
WHERE test = true AND depth < 1000),
result as (SELECT start_id, path, depth, geom FROM
(SELECT *,
row_number() OVER (PARTITION BY array(SELECT * FROM unnest(path) ORDER BY 1))=1 as test_duplicate,
(max(depth) OVER (PARTITION BY start_id))=depth as test_depth
FROM recurse_eol) foo
WHERE test_depth = true AND test_duplicate = true)
SELECT * FROM result
