Le portail francophone de la géomatique

Bonjour Vincent,

Merci pour l'aide. Le résultat est malheureusement toujours le même. J'ai testé avec un autre jeu de données similaires (seule le nombre de données change) et une autre grille de même taille : le résultat est mieux mais toujours pas comme je le souhaiterais.

Au lieu d'avoir un unique nombre par espèce pour toutes les mailles, j'ai l'impression que seules les mailles où il n'y a pas de données de l'espèce ont ce nombre (qui devraient normalement être zéro mais qui là semble être le minimum de données comptées dans une maille pour l'espèce en question).

Je ne comprends pas pourquoi il attribut un nombre de données aux mailles qui n'en n'ont pas (si c'est bien ça) ?

Je ne comprends pas non plus pourquoi la requête marche mieux avec ce jeu de données ?

Moi j'écrirais quelque chose comme  ça :

SELECT    d.nom_latin,
        m.FID,
        count(d.n_bd_dep) as nombre_donnees,
        m.geometry
FROM    "données 2021" d
JOIN    "maillles 5k Loire" m ON ST_Intersects(d.geometry, m.geometry)
GROUP BY d.nom_latin,
        m.FID,
        m.geometry

N'hésite pas à mettre des alias aux tables, c'est tout de suite plus lisible

Dernière modification par Sylvain M. (Thu 18 November 2021 15:30)

Loup, tu n'aurais pas des entités multiparties dans ton jeu de données  ? Ça fausserait les résultats.

J'ai vais tester ma solution dans QGIS dès que possible.

Je vais regarder le fichier, mais comme je viens de le télécharger, je te donne un petit conseil que normalement donne tout formateur SIG à ses élèves : jamais d'espaces, d'accents ou de caractères spéciaux dans le nom (et le chemin) de ses fichiers SIG, shapefile en particulier !

C'est la raison pour laquelle tu dois mettre des guillemets partout dans tes requêtes, alors qu'ils seraient inutiles si tu avais respecté cette consigne de base en SIG
Le code de ta requête serait alors beaucoup plus lisible

Merci Sylvain de passer du temps pour tenter de m'aider ! :-)

J'ai testé ça (après avoir renommé mes couches "D" et "M") :

SELECT    D.nom_latin,
        M.FID,
        count(d.n_bd_dep) as nombre_donnees,
        M.geometry
FROM    D
JOIN    M ON ST_Intersects(D.geometry, M.geometry)
GROUP BY D.nom_latin,
        M.FID,
        M.geometry

Le résultat est presque bon, le nombre de données par espèce et par maille est bien affiché mais toujours le même problème avec les mailles vides :
Le minimum de chaque espèce est attribué au mailles vides au lieu de zéro. Par exemple : toutes les mailles sans donnée de A. affinis ont "8" en "nombre_donnees".

En effet : bizarre un SHP avec des géométries nulles : ce n'est pas très propre...
Cette requête corrige le problème :

SELECT    d.nom_latin,
        m.FID,
        count(d.n_bd_dep) as nombre_donnees,
        m.geometry
FROM    "données actualisées 2021" d
JOIN    "maillles 5k Loire" m ON ST_Intersects(d.geometry, m.geometry)
WHERE d.geometry is not null
GROUP BY d.nom_latin,
        m.FID,
        m.geometry;

WHERE d.geometry is not null

Par contre, vous devriez retirer le lien vers le fichier de données : ce sont des données nominatives des observateurs, et aussi avec potentiellement des espèces sensibles : ce fichier ne devrait pas être partagé tel quel.