Le portail francophone de la géomatique

fnixou · Mon 06 March 2017 12:51

Salut à tous je me retrouve face à un problème un peu gênant, en effet quand j'essaye d’exécuter une requête pour créer une couche virtuelle mon QGIS plante ( minidump )
et le message dans le rapport de bug est le suivant : "Le thread a tenté de lire à partir d'une adresse virtuelle ou d'écrire vers une adresse virtuelle pour laquelle il ne dispose pas d'un accès appropriés".

Cela provient surement de la requête que je vous poste également.
J'espère que quelqu'un aura une idée car pour le moment c'est sacrément vide de mon côté

SELECT t3.label ,
SUM(t3.NBHAB) AS TOTHAB
FROM (SELECT t1.label,
area(intersection(t1.geometry, t2.geometry))/area(t2.geometry)*t2.DONNEESP10 as NBHAB
FROM test_fusion AS t1,
test_iris AS t2
WHERE st_intersects(t2.geometry, t1.geometry)) t3
GROUP BY t3.label

Merci d'avance.

michel56 · Mon 06 March 2017 16:07

Bonjour,

pour la fonction st_intersects il faut mettre un alias

WHERE st_intersects(alias1, alias2)

fnixou · Mon 06 March 2017 16:14

a moins que je n ai pas bien compris votre réponse mais t1 et t2 sont bien des alias dans ma requête . Si vous pouviez préciser du coup ce que vous souhaitiez pointer dans la requête je vous en serai gré , car la je suis un poil perdu .

michel56 · Mon 06 March 2017 17:52

OUI OK

mets juste l'alias

Code:

WHERE st_intersects(t2, t1)

et tu es sur que intersection() est une fonction?

Le deuxième SELECT ne serait pas plutôt dans la clause where

Tu devrais utiliser le gestionnaire de base de données pour faire ta requête, la fenêtre sql et l'assistant.

Une requête sur deux tables c'est déjà compliqué mais sur 3 bon courage.

Une requête sur deux tables qui fonctionne si ça peut t'aider.

Code:

SELECT "CouchePoints"."geometry" AS G_CouchePoints,
"CouchePoints"."annee",
"CouchePolygones"."geometry" AS G_CouchePolygones,
"CouchePolygones"."NomCommunes",
Sum("CouchePoints"."Nb") AS SumCouchePoints,
Count("CouchePoints"."geometry") AS CountCouchePoints 
FROM "CouchePoints", "CouchePolygones"
WHERE ST_Within("G_CouchePoints","G_CouchePolygones")
GROUP BY "CouchePoints"."annee", "CouchePolygones"."NomCommunes"
ORDER BY "CouchePoints"."annee" ASC, "CouchePolygones"."NomCommunes" ASC

Dernière modification par michel56 (Mon 06 March 2017 18:58)

fnixou · Tue 07 March 2017 10:16

la fonction intersection existe belle et bien elle permet de créer un nouvelle objet a partir de l intersection des objets 1 et 2 si je ne me trompe pas.

A priori non le select est bien dans le from l'idée étant de créer une table 3 à partir de 1 et 2 et d'en extraire la donnée nbhab pour chaque label.

Pour faire simple j'ai 2 couches qui se superposent et j'aimerai calculer la quantité d'habitants dans la seconde couche à partir de mes données dans la première couche.

Merci en tout cas du coup de main.

Dernière modification par fnixou (Tue 07 March 2017 10:16)

michel56 · Tue 07 March 2017 10:55

Bonjour,

intersection() c'est une fonction QGIS qui fonctionne dans une expression.

Dans une requête SQL ?.

Pour tes deux couches, quelle est la relation entre les deux tables.

Dernière modification par michel56 (Tue 07 March 2017 10:56)

fnixou · Tue 07 March 2017 12:38

J'utilise l'ajout de couche virtuelle de QGIS mes 2 couches sont des shp donc non contenu dans une quelconque base ( je peux éventuellement les ajouter à un postgre+postgis si nécessaire mais je doute que ça change la donne )

Pour ce qui est de intersection l'ajout de couche virtuelle de Qgis utilise SQLite et spatialite comme moteur et cette fonction est incluse dans Spatialite.

Voilà pour le petit topo j'espère que ça t'aidera.

michel56 · Tue 07 March 2017 12:58

Bonjour

Tu as essayé avec le gestionnaire de base de données et l'assistant.
Dans la liste des fonctions il y bien ST_intersection()

fnixou · Tue 07 March 2017 15:07

C'est ce sur quoi je planche actuellement et j'ai également décomposé ma requête en 2 étapes.
Dans un premier temps je crée ma table T3 et je vais la requêter ensuite (le traitement était surement trop lourd pour QGIS )

fnixou · Tue 07 March 2017 18:06

Finalement je suis passé sur postgre j'ai créé mes tables et ma requête ne me fait plus planter qgis mais je me retrouve avec une erreur de topologie sur mes données j'avance mais j'ai vraiment l'impression d'y aller à pas d'escargots.
Enfin bon je finirai bien par réussir :p

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 Mon 06 March 2017 12:51

QGIS: Plantage pendant la création d'une couche virtuelle

#2 Mon 06 March 2017 16:07

Re: QGIS: Plantage pendant la création d'une couche virtuelle

#3 Mon 06 March 2017 16:14

Re: QGIS: Plantage pendant la création d'une couche virtuelle

#4 Mon 06 March 2017 17:52

Re: QGIS: Plantage pendant la création d'une couche virtuelle

Code:

Code:

#5 Tue 07 March 2017 10:16

Re: QGIS: Plantage pendant la création d'une couche virtuelle

#6 Tue 07 March 2017 10:55

Re: QGIS: Plantage pendant la création d'une couche virtuelle

#7 Tue 07 March 2017 12:38

Re: QGIS: Plantage pendant la création d'une couche virtuelle

#8 Tue 07 March 2017 12:58

Re: QGIS: Plantage pendant la création d'une couche virtuelle

#9 Tue 07 March 2017 15:07

Re: QGIS: Plantage pendant la création d'une couche virtuelle

#10 Tue 07 March 2017 18:06

Re: QGIS: Plantage pendant la création d'une couche virtuelle

Pied de page des forums