ポリゴン内にあるポイントでデータベース行を更新する

2つのテーブルを含むPostGIS / Postgresql DBがあります。1つはポイントジオメトリを持ち、もう1つは国の境界をポリゴンとして表します。各ポイントが交差する国名を、ポイントテーブルの各行に追加したいと思います。多分1つの大きな更新クエリとして。これは単純なSQLを使用して実行できると思いますが、どこから始めればよいかわかりません。これについてのアドバイスは大歓迎です...

機能を必要としない別のオプション

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

これはテストしたことはありませんが、例のようにネストされた関数を使用するよりも高速になると思います。

Explainを実行したときの私の出力（うまくいけば、あなたは似ているはずです）Seq Scanの結果がそれ以上ある場合は、それを確認する必要があります。おそらく、インデックスが適切に設定されていない可能性があります。

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)"

OK ...私は少しハッキングを行ったところ、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