Discussion :

[Bettylolz]

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

[ced]

Bonjour,

Première explication possible : le système de projection de chacune des deux couches. Est-ce que c'est le même système de projection ?

[Bettylolz]

Bonjour,

Première explication possible : le système de projection de chacune des deux couches. Est-ce que c'est le même système de projection ?

Bonjour,
Merci de votre réponse.
Normalement c'est bon.
D'ailleurs quand je transforme la géométrie de mon vecteur ca me renvoie ca:



Par contre pour le raster j'avoue que je me perds dans les moyens de vérifier la résolution et la projection.
Est-ce que vous avez quelques astuces pour que je vérifie si ca vient de là?

Merci

[SQLpro]

Votre topologie, vraisemblablement une ligne brisée, se recoupe. Une telle topologie n'existe pas dans la nature. Par exemple, les limites d'un champs de patates (ligne définissant le périmètre, aussi appelée frontière) ne peuvent se croiser. Corrigez votre topologie, sinon certaines opérations sont impossibles.

Lisez mon cours sur le sujet...
https://blog.developpez.com/sqlpro/p...on_geographiqu

A +

[Bettylolz]

Merci, j'explore ca et je reviens vers vous.

Bonne journée