Zähle Punkte im Polygon mit Postgis

Ich habe ein einfaches Problem: Ich möchte die Anzahl der Punkte innerhalb eines Satzes von Polygonen zählen.

Ich habe bereits eine SQL, aber sie gibt nur die GID des Polygons zurück, das tatsächlich Punkte enthält.

Meine Tabellen: eine Polygonebene mit 19.000 Zeilen und eine Punktebene mit 450 Zeilen.

Die folgende SQL

select grid.gid, count(*) AS totale FROM grid, kioskdhd3 WHERE 
st_contains(grid.geom,kioskdhd3.geom) GROUP BY grid.gid;

Es werden nur 320 Polygone zurückgegeben, die tatsächlich Punkte enthalten. Aber ich möchte, dass alle Polygone zurückgegeben werden, auch wenn die Anzahl der Punkte 0 ist.

Das hat natürlich mit meiner WHERE-Klausel zu tun. Wo muss ich meine st_contains ablegen?

Vielen Dank Carsten

SELECT grid.gid, count(kioskdhd3.geom) AS totale 
FROM grid LEFT JOIN kioskdhd3 
ON st_contains(grid.geom,kioskdhd3.geom) 
GROUP BY grid.gid;

Nach diesem Link werden Joins in GiST-Indizes nicht unterstützt.

Darf ich empfehlen:

SELECT grid.gid, 
       SUM(CASE WHEN st_contains(grid.geom,kioskdhd3.geom) THEN 1 ELSE 0 END) AS total
FROM grid, kioskdhd3 
GROUP BY grid.gid;