Conta i punti in poligono con Postgis

Ho un semplice problema: voglio contare il numero di punti all'interno di una serie di poligoni.

Ho già un SQL ma restituisce solo il gid del poligono che in realtà contiene punti.

Le mie tabelle: uno strato poligonale con 19.000 righe e uno strato punto con 450 righe.

Il seguente SQL

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

restituisce solo circa 320 poligoni che contengono effettivamente punti. Voglio che tutti i poligoni vengano restituiti, anche se il numero di punti è 0.

Naturalmente ha a che fare con la mia clausola WHERE. Dove devo inserire i miei st_contains?

Grazie 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;

Per questo i join di sinistra non sono supportati negli indici GiST.

Posso raccomandare:

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;