Berechnung der prozentualen Schnittfläche in der where-Klausel

Ich habe eine Tabelle mit Polygonen (Zensusblockgruppen) in Postgres. Ich möchte jede Blockgruppe mit der Stadt (einer anderen Polygontabelle) versehen, in der sie sich hauptsächlich befindet. Ist das möglich? Ich denke, ich müsste im Wesentlichen Folgendes erstellen:

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  

aber diese Abfrage dauert ewig (ich habe ungefähr 5.000 Blockgruppen und 375 Städte ...). Irgendwelche Vorschläge, wie diese Abfrage funktioniert, wenn sie falsch ist oder schneller, wenn sie korrekt ist?

Die Art und Weise, wie Sie es tun, wird funktionieren, aber es wird zu viel Zeit in Anspruch nehmen, da postgis versucht, die Geometrie der Schnittpunkte jeder Kombination aus "Blockgruppe und Stadt" zu erstellen, selbst wenn sie sich nicht einmal berühren.

Fügen Sie Ihrer WHERE-Klausel eine weitere Bedingungsüberprüfung hinzu, um zu überprüfen, ob die beiden Geometrien sich überschneiden, und setzen Sie sie vor die vorhandene:

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

Wenn Sie in SQL eine Liste von Bedingungen in der WHERE-Klausel haben, werden diese in der Reihenfolge getestet, in der sie geschrieben werden. Wenn in einem der frühen Vorgänge ein FALSE zurückgegeben wird, überspringt die Abfrage nur die Prüfung der anderen Bedingungen, da das Ergebnis immer FALSE ist.

Stellen Sie außerdem sicher, dass Sie über räumliche Indizes für blockgroups.wkb_geometry und towns.wkb_geometry verfügen.

Wenn sich einige Ihrer Zensuseinheiten über drei Ihrer Städte erstrecken (und Sie daher nicht mehr als 50% Stürze in einer Stadt garantieren können), können Sie dies zu Alexandres sehr nützlicher Antwort hinzufügen:

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;

Dies schützt grundsätzlich vor folgender Situation - in der die Bereiche in Blau verschwinden würden:

Mit st_intersects und && können Sie Folgendes verwenden: Calcolo_perc_intersez_postgis