पॉलीगॉन के भीतर आने वाले बिंदुओं के साथ डेटाबेस पंक्ति को अपडेट करें

Question

मेरे पास एक PostGIS / Postgresql DB है जिसमें दो टेबल हैं। एक बिंदु ज्यामितीय और दूसरा जो देश की सीमाओं को बहुभुज के रूप में दर्शाता है। मैं देश का नाम जोड़ना चाहूंगा कि प्रत्येक बिंदु मेरी पंक्ति में प्रत्येक पंक्ति के साथ प्रतिच्छेदन करता है। शायद एक बड़ी अपडेट क्वेरी के रूप में। मुझे लगता है कि यह सीधे एसक्यूएल का उपयोग करना संभव है, लेकिन मुझे नहीं पता कि कहां से शुरू करें। इस पर कोई सलाह बहुत सराहना की जाएगी ...

Answer 1

एक और विकल्प, फ़ंक्शन की आवश्यकता के बिना

update points set country = t1.country from 
(
    select points.oid, countries.name as country from
    countries INNER JOIN points on st_contains(countries.wkb_geometry,points.wkb_geometry)
) t1 
where t1.oid = points.oid

मुझे संदेह है (हालांकि मैंने परीक्षण नहीं किया है) कि यह आपके उदाहरण में नेस्टेड फ़ंक्शन का उपयोग करने से अधिक तेज़ होगा।

रनिंग एक्सपट्र्स से मेरा आउटपुट (उम्मीद है कि आपका लुक ऐसा ही है)। यदि आपको अधिक Seq स्कैन परिणाम मिले हैं, तो यह देखने के लिए कुछ है, शायद अनुक्रमणिका ठीक से सेट नहीं हैं।

Update on points  (cost=1.18..29.40 rows=121 width=129)"
  ->  Nested Loop  (cost=1.18..29.40 rows=121 width=129)"
        Join Filter: _st_contains(countries.geometry, public.points.geometry)"
        ->  Hash Join  (cost=1.18..2.37 rows=28 width=220)"
              Hash Cond: (public.points.oid = public.points.oid)"
              ->  Seq Scan on points  (cost=0.00..1.08 rows=28 width=114)"
              ->  Hash  (cost=1.08..1.08 rows=28 width=110)"
                    ->  Seq Scan on points  (cost=0.00..1.08 rows=28 width=110)"
        ->  Index Scan using "countries_Idx" on countries  (cost=0.00..0.91 rows=1 width=414)"
              Index Cond: (geometry && public.points.geometry)"

Comments

बहुत बढ़िया! यह एक बहुत तेज प्रतीत होता है। धन्यवाद!

— एडमएस्ट्राडा

Answer 2

ठीक है ... मैंने थोड़ी हैकिंग की और पाया कि एक SQL FUNCTION मुझे वहां सबसे ज्यादा मिलता है। किसी को पुल पर ले जाने पर कोई विचार है?

 CREATE OR REPLACE FUNCTION getcountry (
       country_geom geometry
    ) RETURNS TABLE(country text) AS $$
        SELECT b.name as country FROM  
                    geonames d, world_borders b WHERE
                    $1 && b.wkb_geometry 
                    AND intersects($1, b.wkb_geometry) ;
  $$ LANGUAGE SQL;

UPDATE geonames 
    SET country = val
    FROM (SELECT getcountry(point_geom) FROM geonames) AS val