Le portail francophone de la géomatique

Tu devrais y parvenir avec le script "Atlas par valeurs uniques", proposé par Jérémie Cornet sur le Portail SIG.

Couplé avec les fonctions d’agrégations de la 2.16, c'est très puissant !

Avant cela, créé ta couche de données synthétisées "à la maille", via une jointure spatiale. Les requêtes SQL sur des "couches virtuelles", via le menu "Bases de données > Gestionnaire de bases de données" sont très bien pour ce travail.

Bon courage !

A+

Sylvain M.

Bonjour Mathieu,
Merci une nouvelle fois pour ton aide.
J'ai donc réussi à obtenir une table de ma grille avec 3 colones "espece" (nom de l'espèce); "code"(max_code); "code10km" (code de la grille)
Pour cela, j'ai exécuté le code suivant après avoir créé ma BD dans QSpatialite avec mes 2 tables et exporté ma table résultat (avec geometrie)

SELECT espece, max(CODE), code10km, grille.geometry
FROM grille, obs
WHERE st_intersects(grille.geometry, obs.geometry)
GROUP BY espece,  code10km

Nickel, ca correspond à ceux dont j'ai besoin pour générer mes cartes (matière première). Même si j'ai 4800 lignes alors que je n'ai que 91 carrés (les carrés "positifs" se superposent).

Par contre ensuite, pour la génération de l'atlas comme expliqué dans le lien, ça correspond pas tout à fait à mon besoin.
J'ai besoin de "colorier" les carrés de ma grille suivant 3 couleurs en fonction du max_code (= indice de nidification, possible, probable, certain) alors qu'avec "ta" méthode, on peut juste faire de la présence absence.
J'ai fouillé l'outil Atlas, cela ne semble pas possible à réaliser?
Je vais me débrouiller autrement avec un style catégorisé couplé à des sélections sucessives pour chaque espèce

Un grand MERCI
A bientôt
Bérenger

Merci Mathieu,
Ta conbine est la bonne
J'y suis presque, il y a un hic, l'atlas génére 499 cartes. Voir copie_ecran2 en PJ.
La table en sortie que j'utilise pour faire l'atlas, compte 4840 enregistrements (nb espece x mailles >0). Voir copie_ecran en PJ.
Du coup, en générant l'atlas , il me fait X cartes par espece (autant de cartes que le nombre de mailles >0 par espece).
Est-ce que tu aurais une dernière astuce pour résoudre ce problème?
Merci

Dernière modification par BerengerREMY (Sat 10 December 2016 11:47)

Merci Lsam !

"grille_sp" ne doit pas être ta couche de couverture. Relis mon message, il faut que tu crées une autre couche que celle de ta couche de carrés, sans information de géométrie, simplement une couche texte listant tes espèces (avec comme nom de champ "espece"). C'est elle qu'il faudra mettre en couverture.

Dernière modification par Lsam (Tue 13 December 2016 00:36)

Bonjour,

Je suis en charge de l'application du code demandé par Bérenger Remy et de la génération de l'atlas cité ici.
Je rencontre un problème avec le code trouvé pour générer les cartes :

SELECT espece, max(CODE), code10km, grille.geometry
FROM grille, obs
WHERE st_intersects(grille.geometry, obs.geometry)
GROUP BY espece,  code10km

En effet ce code fonctionne très bien lorsque les noms de mes colonnes et de mes tables attributaires sont les mêmes que dans le code, seulement lorsque que j'essaye de changer les appellations ou de rajouter des champs dans "FROM", car j'ai en réalité 29 tableurs de données à traiter, des messages d'erreur apparaissent.

Idéalement, je souhaiterais que mon code puisse ressembler à ça :

SELECT ESPECE, max(CODE), Code10km, grille.geometry
FROM grille, obs, obs1, obs2, obs3....., obs29
WHERE st_intersects(grille.geometry, obs.geometry, obs1.geometry,...., obs29.geometry)
GROUP BY espece, code10km

Il y a sûrement une erreur d'écriture, mais je ne parviens pas à la trouver !

Jeanne.

Bonsoir Jeanne,

une solution pourrait consister à faire des UNION (29 en tout) :

SELECT ESPECE, max(CODE), Code10km, grille.geometry
FROM grille, obs
WHERE st_intersects(grille.geometry, obs.geometry)
GROUP BY espece, code10km
UNION
SELECT ESPECE, max(CODE), Code10km, grille.geometry
FROM grille, obs1
WHERE st_intersects(grille.geometry, obs1.geometry)
GROUP BY espece, code10km
UNION
SELECT ESPECE, max(CODE), Code10km, grille.geometry
FROM grille, obs2
WHERE st_intersects(grille.geometry, obs2.geometry)
GROUP BY espece, code10km
UNION...

Bonjour Jeanne,

la symbologie est appliquée à l'objet courant de l'atlas. Ce qui suppose que tu génères l'atlas.

Je pense que là tu visualises le composeur, pas un aperçu de l'atlas.

Il est accessible avec l'icône entourée en pièce jointe.

C'est donc que la règle de symbologie ne s'applique pas.
Je pense qu'il faut renommer le champ max dans la requête SQL. Là dans ta table de donnée tu dois avoir une colonne qui s'appelle max et non CODE... Je me trompe ?

SELECT ESPECE, max(CODE) AS CODE, Code10km, grille.geometry
FROM grille, obs
WHERE st_intersects(grille.geometry, obs.geometry)
GROUP BY espece, code10km
UNION
SELECT ESPECE, max(CODE) AS CODE, Code10km, grille.geometry
FROM grille, obs1
WHERE st_intersects(grille.geometry, obs1.geometry)
GROUP BY espece, code10km
UNION
SELECT ESPECE, max(CODE) AS CODE, Code10km, grille.geometry
FROM grille, obs2
WHERE st_intersects(grille.geometry, obs2.geometry)
GROUP BY espece, code10km
UNION...

OK il faut donc renommer les champs dans la requête (cf code ce-dessus).

Dans la requête SQL nous avons demandé à spaptialite de calculer la valeur max de ta colonne code, sans préciser de nom particulier pour cette colonne. Par défaut la base de données l'a appelé max(CODE). Pour la renommer il aurait fallu le dire à la base de données avec la syntaxe en gras ci-dessous :

SELECT ESPECE, max(CODE) AS CODE, Code10km, grille.geometry
FROM grille, obs2
...