Le portail francophone de la géomatique

thierry974 · Thu 17 May 2018 11:28

Bonjour,

J’effectue une requête spatiale complexe (intersection puis agrégation de polygone suivant la valeur d’un attribut) sous via le Db manager. (cf requête ci-dessous) + contrôle sur la géométrie.
La requête fonctionne (4053 objets) .
Mais je n’arrive pas à la charger systématiquement sous QGIS. J’ai l’impression que cela dépend du nombre d’objets créés.

En effet si :
Je change la valeur de parcelle.code_com, (=> nombre d’objets moindre) j’arrive à charger la couche
Je rajoute en fin de requête : limit 500, j’arrive à charger la couche.

Quelqu’un a-t-il une solution ?
Je précise, je suis sur Qgis 3.0.2, Windows 7 64bits avec 16 Go de mémoire vive.

-------------------------------------------------------------------------------------------------------------------------------
With parcelle_union as (
SELECT ST_Intersection(obs_friche."saup_cultivable_moins_parcelle_bati_de_1000m2_v2".geom, "2016novcadastre".parcelle.geom) as geom_intersect,
"2016novcadastre".parcelle.code_com as code_com,
'97' || substring("2016novcadastre".parcelle.idu from 1 for 3) || substring("2016novcadastre".parcelle.idu from 7 for 11) as num_parcelle,
st_area(ST_Intersection(obs_friche."saup_cultivable_moins_parcelle_bati_de_1000m2_v2".geom, "2016novcadastre".parcelle.geom)) as sau,
"2016novcadastre".parcelle.supf as supparc,
round ((st_area(ST_Intersection(obs_friche."saup_cultivable_moins_parcelle_bati_de_1000m2_v2".geom, "2016novcadastre".parcelle.geom)) /"2016novcadastre".parcelle.supf)*100) as pc_sau
from "2016novcadastre".parcelle, obs_friche."saup_cultivable_moins_parcelle_bati_de_1000m2_v2"
WHERE
ST_Intersects(obs_friche."saup_cultivable_moins_parcelle_bati_de_1000m2_v2".geom, "2016novcadastre".parcelle.geom)=TRUE
AND st_area(ST_Intersection(obs_friche."saup_cultivable_moins_parcelle_bati_de_1000m2_v2".geom, "2016novcadastre".parcelle.geom)) > 500
AND "2016novcadastre".parcelle.supf >0
AND "2016novcadastre".parcelle.code_com ='410')
SELECT
row_number() OVER(ORDER BY num_parcelle DESC) as oid,
num_parcelle,
st_union (geom_intersect) as geom_intersect,
ST_IsValidReason(st_union (geom_intersect)) as code_erreur,
sum(sau) as sau,
sum(pc_sau) as pc_sau,
min (supparc) as suppar
from parcelle_union
Where ST_IsValid(geom_intersect)= true
group by num_parcelle

SANTANNA · Thu 17 May 2018 14:23

Bonjour,
Je vois des intersections dans la requête, et aussi des vérifications de validité de géométrie mais les géométries renvoyées sont-elles toutes du même type?

thierry974 · Thu 17 May 2018 15:51

deux type de géométrie sont renvoyés (polygone et multi polygone)
requête complète que je n'arrive pas à chargé dans qgis ou celle que je limite à 500
le problème ne vient donc pas de là.
Merci quand même

SANTANNA · Thu 17 May 2018 16:22

Quid lorsque la requête est faite sous PGAdmin et génère soit une nouvelle table soit une vue? QGIS lit-il le résultat?

trovez · Fri 18 May 2018 09:12

Bonjour,

La requête comporte
st_union (geom_intersect) as geom_intersect,
qui peut produire des géométrie invalides.

Avez-vous inspecté les valeurs du champ ST_IsValidReason(st_union (geom_intersect)) as code_erreur ?

Cdlt

thierry974 · Fri 18 May 2018 10:39

Bonjour,

Toutes les géométries sont valides
ST_IsValidReason(st_union (geom_intersect))

j'ai rajouté une condition pour en être sûr : Where ST_IsValid(geom_intersect)= true

trovez · Fri 18 May 2018 11:08

Le where porte sur St_Intersection, mais pas sur St_union, St_Union pouvant aussi générer des erreurs.

thierry974 · Fri 18 May 2018 11:12

SANTANNA

Si je crée une vue via pgadmin et que je l’affiche sous qgis : ça fonctionne

merci.
je ne sais toujours pas quel est le problème. j'aurais souhaité ne pas créer de vue. une vue par commune, c'est beaucoup.

SANTANNA · Fri 18 May 2018 11:20

j'aurais souhaité ne pas créer de vue. une vue par commune, c'est beaucoup.

Je ne comprends pas. N'est-ce pas la même requête que dans DB Manager (ou celle-là portait déjà sur une commune)? Il est aussi possible de faire une vue agrégée dans PGAdmin (à voir si la structure de la BD permet de faire cela facilement) et après soit la filtrer dans QGIS directement ou dans DB Manager avant ajout dans le projet. De ce fait, la BD PostGIS ne contiendrait qu'une seule vue.

Ben22 · Fri 18 May 2018 12:06

thierry974 a écrit:

deux type de géométrie sont renvoyés (polygone et multi polygone)
requête complète que je n'arrive pas à chargé dans qgis ou celle que je limite à 500
le problème ne vient donc pas de là.
Merci quand même

Essayez de tout passer en multipolygone en forçant avec st_multi, il me semble que le DBManager n'aime pas ajouter une couche avec deux types de géométries différentes.

dhaulagiri · Tue 28 August 2018 09:05

Bonjour,

Je parcours ce post en diagonale et il me semble rencontrer un problème similaire pour lequel je n'ai pas trouvé de solution:

J'effectue une requête simple pour découper les objets d'une table par ceux d'une autre table, du type:
SELECT table1.champ, st_multi(st_intersection(table2.the_geom, table1.the_geom)) AS the_geom
FROM table1 JOIN table2 ON st_intersects(table1.the_geom, table2.the_geom)

Dans PgAdmin ou dans le DBManager (QGIS 3.0.2), elle fonctionne bien. Si j'ajoute un CREATE TABLE, elle fonctionne toujours bien. Par contre, et c'est malheureusement mon objectif, elle ne fonctionne plus si je veux créer une vue avec un CREATE OR REPLACE VIEW. Concrètement, la vue est créée mais je ne peux pas l'afficher dans QGIS: elle est dédoublée dans la liste des couches (deux géométries?) mais aucune des deux couches ne peut être affichée.

J'ai aussi pensé à la géométrie mais je constate qu'un SELECT distinct ST_GeometryType(the_geom) renvoie parfois deux types différents pour certaines tables sans que cela m'empêche de les afficher dans QGIS. Manifestement, la création d'une table initialise la géométrie des objets d'une manière différente de la vue mais j'en suis réduit à cette hypothèse. Auriez-vous une idée de solution?

Je vous remercie,
Nicolas

Dernière modification par dhaulagiri (Tue 28 August 2018 09:22)

gglafouine · Wed 29 August 2018 08:34

dhaulagiri a écrit:

Auriez-vous une idée de solution?

Bonjour,
la pemière choseà faire c'est de vérifier les messages. Si il y a une erreur elle sera visible lors de l'exécution dans le panneau "Journal des messages"

Il est possible que le client postgresql + pgadmin utilisé dans la compilation de QGIS diffère de la version de votre base de données et de sa cartouche spatiale

Donc 4 choses à vérifier:

- messages du journal
- version des bases et cartouches utilisé dans QGIS
- version des bases et des cartouches du serveur SGBD
- version de GDAL OGR car il sert à lire et écrire les données dans la base car cela peut résulté d'un problème aussi à ce niveau.

Après si tu veux forcer une géométrie en particulier tu peux aussi filtrer sur des types et dire que le résultat attendu sera de type polygon

dhaulagiri · Thu 13 September 2018 12:56

Merci pour ton retour,

Les messages n'étaient pas très explicites mais bon.

Le problème était double:
- la table source avait deux entrées distinctes dans la table geometry_columns du schéma. J'en ai supprimé un
- j'ai généré un champ gid lors de la création de la vue: CAST OVER () AS VARCHAR) AS gid

Communauté	Ressources	L'emploi	Le marché	L'association	Fils RSS
Forum Blogs Wiki Planet Sites amis	Agenda Biblio Fiches AMO-CNIG Paris PCGIS Téléchargements	Offres d'emploi Banque de CV Formations Démarche Métiers	Geo-Entreprises Geo-Communiqué Appels d'offres	L'association Qui sommes-nous ? Nous aider Bienfaiteurs Partenaires Historique Statuts	Tout le site Job Géomatique Services GeoBlogs ...	Faire un don Mentions légales

Le portail francophone de la géomatique

Annonce

#1 Thu 17 May 2018 11:28

QGIS 3.02/DB Manager: requete postgis/Impossible de charger la couche

#2 Thu 17 May 2018 14:23

Re: QGIS 3.02/DB Manager: requete postgis/Impossible de charger la couche

#3 Thu 17 May 2018 15:51

Re: QGIS 3.02/DB Manager: requete postgis/Impossible de charger la couche

#4 Thu 17 May 2018 16:22

Re: QGIS 3.02/DB Manager: requete postgis/Impossible de charger la couche

#5 Fri 18 May 2018 09:12

Re: QGIS 3.02/DB Manager: requete postgis/Impossible de charger la couche

#6 Fri 18 May 2018 10:39

Re: QGIS 3.02/DB Manager: requete postgis/Impossible de charger la couche

#7 Fri 18 May 2018 11:08

Re: QGIS 3.02/DB Manager: requete postgis/Impossible de charger la couche

#8 Fri 18 May 2018 11:12

Re: QGIS 3.02/DB Manager: requete postgis/Impossible de charger la couche

#9 Fri 18 May 2018 11:20

Re: QGIS 3.02/DB Manager: requete postgis/Impossible de charger la couche

#10 Fri 18 May 2018 12:06

Re: QGIS 3.02/DB Manager: requete postgis/Impossible de charger la couche

#11 Tue 28 August 2018 09:05

Re: QGIS 3.02/DB Manager: requete postgis/Impossible de charger la couche

#12 Wed 29 August 2018 08:34

Re: QGIS 3.02/DB Manager: requete postgis/Impossible de charger la couche

#13 Thu 13 September 2018 12:56

Re: QGIS 3.02/DB Manager: requete postgis/Impossible de charger la couche

Pied de page des forums