Bonjour à tous,

Je suis statisticienne dans le domaine de la santé et je m'intéresse depuis peu à l'utilisation de données environnementales à proximité du lieu d'habitation des patients.
Après avoir constaté que QGis était un logiciel "clic/bouton" qui ne permettait pas, entre autre, d'automatiser mes programmes : je me lance sur PostGIS.
J'ai commencé par les vecteurs et je commence à comprendre la logique mm si je pense que j'aurai l'occasion de rencontrer très prochainement de nouvelles difficultés.
Je me penche maintenant sur les rasters et j'avoue que je suis un peu perdue.

Je veux faire une intersection entre mon raster de pollution lumineuse nocturne moyennée sur l'année et mon vecteur des contours des départements.
Le but étant d'obtenir une valeur moyenne par département.

Je bosse à partir du bouquin PostGIS Cookbook dont j'ai sorti le programme suivant, que j'ai adapté par la suite:

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
select dep.dep_name,(ST_SummaryStats(ST_Union(ST_Clip(A.rast,1,dep.geom,TRUE)),1)).mean
From public."nighttime_light" as A
Join public."dep_4326" as dep ON ST_Intersects(A.rast,dep.geom)
group by dep_name;
Ca tourne!
Mais ca me sort 0 ligne.

Pourriez-vous m'aider?

D'avance un grand Merci