where句の交差面積の割合を計算する

postgresにポリゴン（国勢調査ブロックグループ）のテーブルがあります。各ブロックグループに、主に存在する町（別のポリゴンテーブル）をタグ付けします。これは可能ですか？私は本質的に次のようなものを作成する必要があると考えています：

select b.*,t.name  
from blockgroups b, towns t  
where (st_area(st_intersection(b.wkb_geometry, t.wkb_geometry))/st_area(b.wkb_geometry)) > .5  

しかし、このクエリには永遠に時間がかかります（約5,000のブロックグループと375の町があります...）。このクエリが正しくない場合はまったく機能しない、または正しい場合は高速になる方法についての提案はありますか？

あなたがそれをしている方法は機能しますが、postgisは触れていなくても、すべての「ブロックグループと町」の組み合わせの交差点のジオメトリを作成しようとしているため、時間がかかりすぎます。

WHERE句に別の条件チェックを追加して、2つのジオメトリがインターセプトするかどうかを確認し、既存のジオメトリの前に配置します。

select b.*,t.name
from blockgroups b, towns t
where st_intersects(b.wkb_geometry, t.wkb_geometry) and    
    (st_area(st_intersection(b.wkb_geometry, t.wkb_geometry))/st_area(b.wkb_geometry)) > .5

SQLでは、WHERE句に条件のリストがある場合、書き込まれた順序でテストされます。初期操作の1つでFALSEが返された場合、クエリは結果が常にFALSEになるため、他の条件チェックをスキップします。

また、blockgroups.wkb_geometryおよびtowns.wkb_geometryに空間インデックスがあることを確認してください。

アレクサンドルの非常に有用な答えに加えて、国勢調査ユニットのいくつかがあなたの町の3つにまたがる可能性がある場合（したがって、どの町でも50％以上の落下が保証できない場合）、これを行うことができます：

select distinct on (b.id)
b.*,t.name,
(st_area(st_intersection(b.wkb_geometry, t.wkb_geometry))/st_area(b.wkb_geometry)) as proportion
from blockgroups b, towns t
where st_intersects(b.wkb_geometry, t.wkb_geometry) 
order by b.id, proportion desc;

これは、基本的に次の状況から保護します-青の領域が消えます：

st_intersectsおよび&&演算子では、これを使用できます：Calcolo_perc_intersez_postgis