मेरे पोस्टग्रेक्यूएल डेटाबेस में मेरे पास स्ट्रीट (रोड सेंटर लाइन) और बहुभुज टेबल हैं। नमूना परिदृश्य इस प्रकार है:
मुसीबत:
मुझे गली के चारों ओर 50 मीटर बफर के चौराहे पर सड़क के साथ समानांतर रेखाओं की गणना करने की आवश्यकता है और दोनों तरफ निकटतम इमारत बहुभुज। वांछित आउटपुट परिदृश्य है:
मैंने क्या कोशिश की है:
मेरा दृष्टिकोण था:
1) Generate 50m buffer around street layer
2) get the intersection of buffer and polygons
3) Compute the distance
4) Draw offset curves (parallel lines) at both sides of street layer
5) Merge both curves to get parallel lines at the intersection
यहाँ मेरा प्रयास है:
WITH street_buffer AS (
SELECT
street.gid street_id,
street.geom street_geom,
ST_Buffer(street.geom, 50, 'endcap=square join=round') geom1,
building.geom geom2
FROM street
LEFT JOIN building on ST_DWithin(building.geom, street.geom, 50)
ORDER BY street_id
),
selected_buildings AS (
SELECT
street_id,
street_geom,
ST_Intersection(geom1, geom2) geom
FROM street_buffer
),
distance AS (
SELECT
street_id,
street_geom,
ST_Distance(street_geom, geom) as dist
FROM selected_buildings
),
curves AS (
SELECT
street_id,
ST_OffsetCurve(ST_LineMerge(street_geom), dist) as curve1,
ST_OffsetCurve(ST_LineMerge(street_geom), -dist) as curve2
FROM distance
ORDER BY street_id
)
SELECT
street_id,
ST_Union(curve1, curve2) geom
FROM curves
ORDER BY street_id
उपरोक्त कोड के साथ समस्या यह है कि यह वांछित आउटपुट के अनुसार समानांतर लाइनों को वापस नहीं करता है, अर्थात, सभी बहुभुज चौराहों पर समानांतर रेखाएं निकटतम बहुभुज के चौराहे पर उत्पन्न हो रही हैं।
EDIT_1:
उपरोक्त कोड का वास्तविक उत्पादन है:
हालांकि, उपरोक्त आउटपुट में केवल पीली समानांतर रेखाएं (सड़क के दोनों ओर निकटतम बहुभुजों के लिए वक्र घटती हैं) की आवश्यकता होती है:
क्या कोई मुझे सुझाव दे सकता है कि वांछित आउटपुट कैसे प्राप्त करें?