Le portail francophone de la géomatique

Bonjour,
QGIS n'autorise l'affichage de couches PG que s'il détecte un champ clé primaire (pour lui, il faut que ce soit un integer non dupliqué).
Cette condition semble être respectée dans la première situation (le gid). Mais l'est-elle dans la seconde? En récupérant des valeurs de la seconde table (cours_eau_permanents_nommes), n'y a-t-il pas le risque de se retrouver, pour un même gid de captage, avec deux valeurs cours_eau_nom? Auquel cas, le gid n'est plus unique... (ce qui expliquerait le 55ème individu)

Rebonjour,
J'ai fait une seconde capture dans la foulée, qui vous montre
- la requête visant à récupérer les noms des cours d'eau associés aux captages
- une partie des données remontées
- les messages d'échec d'affichage des données graphiques
- la gestion, comme une couche Qgis entière de la première requête portant sur les seuls attributs des captages
En espérant que l'exposé de la situation soit suffisamment clair, même si cela ne vaut pas l'examen des données en direct
Cordialement

Pour le order by n_prelev_aep_p_068_echange.gid, j'ai bien vérifié, j'ai 55 gid différents. Sur une série réduite comme celle-ci, que j'ai pu examiner lentement, un doublon m'aurait vite interpellé.

J'ai essayé la requête, que je n'ai pas très bien comprise. En tout cas le message d'erreur de Echec_requete se rapporte à un élément du code qui fonctionnait sans problème avec la requête initiale (celle des 55 enregistrements)

"SELECT row_number() over ()::integer as num_unique, liste.* FROM" a-t-il pour objet de forcer, si elle ne se fait pas, la création d'un identifiant unique dans la requête ? A quel ensemble d'objets (captages de la table postgis, ou cours d'eau de la requête enregistrée dans le projet qgs) se rapporterait num_unique ?

"(SELECT DISTINCT n_prelev_aep_p_068_echange.gid, n_prelev_aep_p_068_echange.geom, n_prelev_aep_p_068_echange.indice_nat, n_prelev_aep_p_068_echange.date_dup_i, n_prelev_aep_p_068_echange.date_derni, n_prelev_aep_p_068_echange.uge_nom, n_prelev_aep_p_068_echange.identperi_, n_prelev_aep_p_068_echange.nom as captage_nom, cours_eau_permanents_nommes.nom as cours_eau_nom
FROM n_prelev_aep_p_068_echange, (SELECT* from n_troncon_cours_eau_bdt_068 where not nom='NR' and not regime='Intermittent') as cours_eau_permanents_nommes)" : j'ai du mal à comprendre où est l'erreur. Faudrait-il mettre un alias à l'extérieur de cet ensemble ? J'ai fait un essai, non concluant, cf la capture alias_inoperant.jpg

Et où cette requête serait-elle stockée ?

Dernière modification par Bertrand DORNER (Fri 07 November 2014 16:56)

Bonjour, et bonne année à tous,

J'ai finalement, étant revenu sur cette recherche, trouvé comment répondre à la série de questions :
trouver et cartographier les captages d'eau potable situés à moins de 40m de distance des tronçons de cours d'eau permanents et nommés de la BD Topo, en affichant les identifiants des captages, les noms des gestionnaires des captages, et les noms des tronçons de cours d'eau.

Voici déjà la requête, brute de commentaires
SELECT DISTINCT row_number() over() as myid, n_prelev_aep_p_068_echange.indice_nat, n_prelev_aep_p_068_echange.gid, n_prelev_aep_p_068_echange.uge_nom, n_prelev_aep_p_068_echange.geom, cours_eau_permanents_nommes.nom
FROM n_prelev_aep_p_068_echange, (SELECT* from n_troncon_cours_eau_bdt_068
where not nom='NR' and not regime='Intermittent'
) as cours_eau_permanents_nommes
WHERE ST_DWithin(n_prelev_aep_p_068_echange.geom, cours_eau_permanents_nommes.geom, 40)
ORDER BY n_prelev_aep_p_068_echange.indice_nat

Essayons de procéder par étapes successives
- Sélection des tronçons de cours d'eau permanents et nommés :
SELECT* (sélection de la clé primaire gid et de la géométrie pour un affichage cartographique)
FROM n_troncon_cours_eau_bdt_068
WHERE not nom='NR' and not regime='Intermittent' (absence de valeur 'NULL' dans le champ nom)

- Sélection des captages à l'intérieur d'une distance de 40m autour des cours d'eau recherchés au départ :
SELECT DISTINCT n_prelev_aep_p_068_echange.indice_nat, n_prelev_aep_p_068_echange.gid, n_prelev_aep_p_068_echange.uge_nom, n_prelev_aep_p_068_echange.geom (sélection de la clé primaire des captages et de la géométrie pour un affichage cartographique)
FROM n_prelev_aep_p_068_echange, (SELECT* from n_troncon_cours_eau_bdt_068
where not nom='NR' and not regime='Intermittent'
) as cours_eau_permanents_nommes (renommer en alias la première requête permet qu'elle soit considérée comme une table)
WHERE ST_DWithin(n_prelev_aep_p_068_echange.geom, cours_eau_permanents_nommes.geom, 40)
On obtient 54 enregistrements, cartographiables, sans relation attributaire avec les cours d'eau

- Question supplémentaire, les points de captages répondant à ces conditions doivent afficher le nom du tronçon de cours d'eau
Dans la clause SELECT DISTINCT, ajout de "row_number() over() as myid" (histoire d'avoir un identifiant unique pour les enregistrements de captages présents en doublons, car à moins de 40m de plusieurs cours d'eau, ou tronçons différents de mêmes cours d'eau) et "cours_eau_permanents_nommes.nom" (noms des tronçons de cours d'eau de la BD Topo)
La clause ORDER BY n'apparaît pas comme indispensable à la cartographie
En rajoutant les noms de cours d'eau j'obtiens 55 enregistrements, non cartographiables (un captage en doublon car à moins de 40m de deux cours d'eau nommés)
En rajoutant le champ myid, j'obtiens 60 enregistrements, cartographiables. Les enregistrements supplémentaires correspondent à deux tronçons d'un cours d'eau de part et d'autre d'un affluent lui aussi nommé et permanent, et à un/des captages situés à cheval de deux tronçons d'un même cours d'eau nommé et permanent, peut-être au confluent de cours d'eau temporaires et sans nom.

En espérant que cette démarche soit utile

Cordialement