Le portail francophone de la géomatique

SangJun · Wed 23 May 2007 16:38

Bonjour,

J'espère que je poste dans la bonne rubrique, si ce n'est pas le cas, je m'en excuse par avance.

J'essaye de calculer une zone autour d'une (multi)ligne sur un SIG.
Je pensais m'en sortir avec un buffer, mais cela ne correspond pas à ma demande.
Pour faire simple, ce que j'aimerai avoir c'est :

[img]http://www.hebergement-images.com/06/1179929327_exemple.jpg[/img]

La ligne centrale étant le chemin.

J'ai essayé (pour voir le résultat que j'aurai obtenu) :
CREATE TABLE Essai as select MemGeomUnion(buffer) from buffer2030;
alter table Essai add gid serial;
alter table Essai add primary key (gid);

Mais Quantum GIS se ferme lorsque je souhaite visualiser le résultat.

Quelqu'un pourrait-il m'aider svp?

Merci

Yves · Wed 23 May 2007 16:44

Bonjour,

QGIS ne permet pas de lancer une requête SQL sauf lors du chargement de données, dans ce cas là c'est un select (CLAUSE SELECT WHERE)
[Edit] je précise un peu : pour faire cela, vous pouvez passer par une console avec psql -d nomBase -U postgres -c "CREATE TABLE ..."

La prochaine version de QGIS améliorera la gestion de postgis.

Y.

Dernière modification par Yves (Wed 23 May 2007 16:45)

cbredel · Wed 23 May 2007 20:24

Bonjour,

Avec quelle version de Qgis travailles-tu ?
La géométrie de la table Essai que tu crée est de quel type ?

Au boulot, j'ai une version de qgis qui ferme direct, comme toi, si j'essai d'ouvrir une table dont la géométrie est de type geometry. Mais je n'ai pas la version en tête. Les versions que j'ai sous la main (0.8 et 0.81) affiche un message d'erreur donc je ne peux pas verifier que c'est bien cela.

A+

Yves · Wed 23 May 2007 21:44

Bonjour,

QGIS ne gère pas le type geometry

Y.

SangJun · Thu 24 May 2007 10:29

Je lance mes requêtes par pgAdmin III.

Qu'appeles-tu géométrie de table?
CREATE TABLE essai
(
memgeomunion geometry,
gid int4 NOT NULL DEFAULT nextval('essai_gid_seq'::regclass),
CONSTRAINT essai_pkey PRIMARY KEY (gid)
)
WITHOUT OIDS;
ALTER TABLE essai OWNER TO fsauvage;
Cela répond-t-il à ta question?

Ma version de qgis est la 0.7.4

SangJun · Thu 24 May 2007 10:34

@ Yves : Je ne comprend pas ta réponse. En effet, ma table Buffer s'affiche, alors que son type est :

CREATE TABLE buffer2030
(
gid int4 NOT NULL,
buffer geometry,
CONSTRAINT buffer2030_pkey PRIMARY KEY (gid)
)
WITHOUT OIDS;
ALTER TABLE buffer2030 OWNER TO fsauvage;

Yves · Thu 24 May 2007 11:04

Bonjour,

Essaye avec un WITH OIDS; à la place de WITHOUT OIDS;

Quand tu dis que ta table s'affiche, tu parles dans pgAdmin, pas dans QGIS ? Il me semble, mais je peut me tromper, que QGIS ne gérant pas ce type de géométrie (geometry, il gère seulement les type (MULTI)Polygon, point et line) il aura tendance à planter

Y.

Yves · Thu 24 May 2007 11:07

SangJun a écrit:

Qu'appeles-tu géométrie de table?

Il parle de memgeomunion geometry, tu peux définir un autre type de géometrie : MULTIPOLYGON par exemple (ce que je te conseil)

Y.

Dernière modification par Yves (Thu 24 May 2007 11:08)

SangJun · Thu 24 May 2007 11:52

ma table buffer2030 s'affiche dans qgis et un select gid, astext(buffer) from buffer2030;
m'affiche des résultats de type : 130546;"POLYGON((3.55758651896317 50.407053429159,3.55708285671878 50.4067498022559,3.55652963753288 50.406550269088,3.55594812130853 50.4064624976031,3.55536065538812 50.4064898608101,3.55478981575751 50.4066313071564,3.55425753946359 50.4068814009381,3.553784281 (..)"

tandis que select gid, astext(memgeomunion) from essaiu; m'affiche des résultats de type : 1;"POLYGON((3.55956587428577 50.4085222612187,3.5595877798921 50.408335560334,3.55954104353758 50.4077493174953,3.55938083490484 50.4071834569588,3.55911331072177 50.4066597244277,3.55874875179271 50.4061982466492,3.55830116791285 50.4058167579558,3.557787759 (..)" mais fait planter QGIS

De plus, j'ai effectué
drop table EssaiInter;
CREATE TABLE EssaiInter AS select intersection(B1.the_geom,B2.buffer) from buildup B1, buffer2030 B2 where isempty(intersection(B1.the_geom,B2.buffer))=false;
alter table EssaiInter add gid serial;
alter table EssaiInter add primary key (gid);
au résultat similaire (POLYGON) et j'obtiens le même plantage.

Yves · Thu 24 May 2007 14:34

Bonjour,

Vous pouvez essayer de voir si vos polygones sont valides avec la fonction isvalid(). Sinon je n'ai pas d'autres pistes pour l'instant.

Y.

SangJun · Thu 24 May 2007 15:00

select gid, isvalid(intersection) from essaiinter;

select gid, isvalid(memgeomunion) from essai;

Seconde colonne : t

:'(

Yves · Thu 24 May 2007 17:57

Bon je commence à peine à comprendre le problème :

Tu as une table buffer2030 qui s'affiche dans QGIS. Ta requête :

CREATE TABLE Essai AS SELECT MemGeomUnion(buffer) FROM buffer2030;
alter table Essai add gid serial;
alter table Essai add primary key (gid);

Tente de faire une Union sur toute des géométries pour en faire une seule. Elle fonctionne et renvoie des objets valides par contre QGIS plante quand tu essayes de la charger.

Tu utilises QGIS 0.7.4 ce qui n'est pas la version la plus récente (tente une mise à jour

La table buffer2030 est :

CREATE TABLE buffer2030
(
gid int4 NOT NULL,
buffer geometry,
CONSTRAINT buffer2030_pkey PRIMARY KEY (gid)
)
WITHOUT OIDS;
ALTER TABLE buffer2030 OWNER TO fsauvage;

Ai je bien résumé ton problème ?

Quelle est ta version de postgis ? Un geomUnion à la place de MemGeomUnion donne quoi ?

Y.

Dernière modification par Yves (Thu 24 May 2007 17:57)

cbredel · Thu 24 May 2007 18:37

Pour continuer sur la piste du champs de type geometry.
Quand tu crée ta table ... tu lui passe aussi les trois commandes suivantes :

select addgeometrycolumn ('Essai','geom',0,'MULTIPOLYGON',2);
update Essai set geom=geometry;
select dropgeometrycolumn('Essai','geometry');

Tu ajoute une colonne geometrique de type MULTIPOLYGON... tu rempli avec le contenu de ta colonne et tu supprimes l'ancienne colonne...
... J'aurai préféré une colonne de type alter geometry columns... mais j'ai pas trouvé.!!!

Normalement cela devrait être bon...

SangJun · Fri 25 May 2007 10:21

Yves a écrit:

Bon je commence à peine à comprendre le problème :

Tu as une table buffer2030 qui s'affiche dans QGIS. Ta requête :
CREATE TABLE Essai AS SELECT MemGeomUnion(buffer) FROM buffer2030;
alter table Essai add gid serial;
alter table Essai add primary key (gid);

Tente de faire une Union sur toute des géométries pour en faire une seule. Elle fonctionne et renvoie des objets valides par contre QGIS plante quand tu essayes de la charger.

Tu utilises QGIS 0.7.4 ce qui n'est pas la version la plus récente (tente une mise à jour

La table buffer2030 est :

CREATE TABLE buffer2030
(
gid int4 NOT NULL,
buffer geometry,
CONSTRAINT buffer2030_pkey PRIMARY KEY (gid)
)
WITHOUT OIDS;
ALTER TABLE buffer2030 OWNER TO fsauvage;

Ai je bien résumé ton problème ?

Quelle est ta version de postgis ? Un geomUnion à la place de MemGeomUnion donne quoi ?

Y.

Pour mettre à jour, je vais voir, je ne suis pas admin de ma machine, je ne peux rien installer.

GeomUnion et MemGeomUnion me renvoie le même problème.

SangJun · Fri 25 May 2007 10:30

cbredel a écrit:

Pour continuer sur la piste du champs de type geometry.
Quand tu crée ta table ... tu lui passe aussi les trois commandes suivantes :

select addgeometrycolumn ('Essai','geom',0,'MULTIPOLYGON',2);
update Essai set geom=geometry;
select dropgeometrycolumn('Essai','geometry');

Tu ajoute une colonne geometrique de type MULTIPOLYGON... tu rempli avec le contenu de ta colonne et tu supprimes l'ancienne colonne...
... J'aurai préféré une colonne de type alter geometry columns... mais j'ai pas trouvé.!!!

Normalement cela devrait être bon...

J'ai suivi tes recommandations, mais dès la première ligne, je me retrouve avec
ERROR: relation "public.Essai" does not exist
CONTEXT: SQL statement "ALTER TABLE public."Essai" ADD COLUMN geom geometry "
PL/pgSQL function "addgeometrycolumn" line 76 at execute statement
SQL statement "SELECT AddGeometryColumn('','', $1 , $2 , $3 , $4 , $5 )"
PL/pgSQL function "addgeometrycolumn" line 4 at select into variables
Etant plutôt novice, je me sens un peu perdu

Dernière modification par SangJun (Fri 25 May 2007 10:31)

SangJun · Fri 25 May 2007 10:46

C'est bon, je suis parvenu à mes fins avec :

DROP TABLE Essai;
CREATE TABLE Essai as select MemGeomUnion(buffer) from buffer2030;
alter table Essai add gid serial;
alter table Essai add primary key (gid);

drop table EssaiInter;
CREATE TABLE EssaiInter AS select B.gid, memgeomunion(intersection(B.the_geom,E.memgeomunion)) from buildup B, Essai E where isempty(intersection(B.the_geom,E.memgeomunion))=false GROUP BY B.gid;
alter table EssaiInter add primary key (gid);

J'ai donc l'intersection entre le buffer et buildup.

Bon, maintenant il ne me reste plus qu'à échanger la fonction buffer par une fonction qui me retournera l'image que j'ai mise dans mon premier post.

Yves · Fri 25 May 2007 10:49

Bonjour,

Le problème vient de la majuscule dans Essai :

select addgeometrycolumn ('essai','geom',0,'MULTIPOLYGON',2);
update essai set geom=geometry;
select dropgeometrycolumn('essai','geometry');

Le problème de l'absence d'entré dans la table geometry est bien vu ;-)
À ce sujet lire la FAQ de postgis (section 3.6) Pourquoi dois je utiliser la fonction AddGeometryColumn() et toutes les autres fonctions de l'OpenGIS ?
sur le site postgis.fr (onglet Manuel). http://www.postgis.fr/node/168

Y.

SangJun · Fri 25 May 2007 11:39

Yves a écrit:

Bonjour,

Le problème vient de la majuscule dans Essai :
select addgeometrycolumn ('essai','geom',0,'MULTIPOLYGON',2);
update essai set geom=geometry;
select dropgeometrycolumn('essai','geometry');

Le problème de l'absence d'entré dans la table geometry est bien vu ;-)
À ce sujet lire la FAQ de postgis (section 3.6) Pourquoi dois je utiliser la fonction AddGeometryColumn() et toutes les autres fonctions de l'OpenGIS ?
sur le site postgis.fr (onglet Manuel). http://www.postgis.fr/node/168

Y.

Vu que les tables que j'utilise sont le résultat de tables déjà existantes, je pense que je n'ai pas besoin d'AddGeometryColumn() (du moins pour le moment), non?

SangJun · Tue 29 May 2007 12:35

La ligne CREATE TABLE EssaiL as select MemGeomUnion(buffer) from bufferLong;
me renvoie comme erreur NOTICE: TopologyException: no outgoing dirEdge found
(3.84808,50.2911,50.2911) ERROR: GEOS union() threw an error!

Erreur que je ne comprend pas. Pourriez-vous m'expliquer ce que cela signifie?

SangJun · Thu 14 June 2007 10:15

Dans le même genre d'erreur, une de mes requêtes dans PostGIS me renvoie sans cesse :

Code:

NOTICE:  TopologyException: no outgoing dirEdge found (3.07069,50.2997,50.2997)
SELECT failed: ERROR:  GEOS union() threw an error!

maRequete= select GeomUnion(GeomFromText('POLYGON((3.064484 50.307842, 3.076605 50.291935, 3.076635 50.291958, 3.064514 50.307865, 3.064484 50.307842))'),'0103000000010000002500000095110D8268900840EB857D1B6226494065C5707500840840310BED9C6627494010ECF82F10840840E92ADD5D67274940D50500ADAC9008401B01558F5E2649407C596AA5A1910840009E184957264940946934B918A30840C576F700DD25494061C43E0114A308404E469561DC2549409835197DDF920840A415E5D84D264940EACE13CFD9A20840B743C362D4254940B7291E17D5A20840401361C3D3254940D33D81149A910840DFD35BBF562649401DED01CEC49008407EA0AE945C264940E5DD596BD490084014FACF4C5B2649400C40A374E99F084095F3C5DE8B254940BCE9961DE29F08405FB2F1608B254940BFECCD3EF39008404D568AC55826494003D29BEA16910840A34B88D85526494086FF7403059E0840274A42226D2549402672C119FC9D0840228D0A9C6C254940D784792179910840BE4344CA4D2649403EAE0D15E39C084065C746205E2549409487855AD39C0840ADA7565F5D2549406BC4EAB8CA900840067E390D5A2649408C04149994900840B84E57DB5D264940F4AFABFC92900840BBA67FF15D2649407C6308008E7D08400F2A711DE3264940A94BC631927D0840865AD3BCE3264940FC77D7715C900840394FFC2D60264940878BDCD3D57D0840DA39CD02ED264940BA30D28BDA7D0840516A2FA2ED26494055886B97599008400DD8ED0661264940E25817B7D1800840A9FB00A4362749402CF2EB87D8800840DE3CD521372749402D8F4D0160900840CBDCE48D6126494012C0CDE2C582084080272D5C56274940724D81CCCE82084086E464E25627494095110D8268900840EB857D1B62264940')

Les 2 éléments dans la fonction GeomUnion sont de type geometry (plus exactement des polygones) et ils sont valides. J'effectue plusieurs fois cette opération avec d'autres valeurs mais là ça bloque.

J'ai beau chercher je reste bloqué.

J'espère que quelqu'un pourra m'aider

SangJun · Tue 19 June 2007 09:54

Le problème venait de la version de GEOS

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 Wed 23 May 2007 16:38

QGIS et PostgreSQL

#2 Wed 23 May 2007 16:44

Re: QGIS et PostgreSQL

#3 Wed 23 May 2007 20:24

Re: QGIS et PostgreSQL

#4 Wed 23 May 2007 21:44

Re: QGIS et PostgreSQL

#5 Thu 24 May 2007 10:29

Re: QGIS et PostgreSQL

#6 Thu 24 May 2007 10:34

Re: QGIS et PostgreSQL

#7 Thu 24 May 2007 11:04

Re: QGIS et PostgreSQL

#8 Thu 24 May 2007 11:07

Re: QGIS et PostgreSQL

#9 Thu 24 May 2007 11:52

Re: QGIS et PostgreSQL

#10 Thu 24 May 2007 14:34

Re: QGIS et PostgreSQL

#11 Thu 24 May 2007 15:00

Re: QGIS et PostgreSQL

#12 Thu 24 May 2007 17:57

Re: QGIS et PostgreSQL

#13 Thu 24 May 2007 18:37

Re: QGIS et PostgreSQL

#14 Fri 25 May 2007 10:21

Re: QGIS et PostgreSQL

#15 Fri 25 May 2007 10:30

Re: QGIS et PostgreSQL

#16 Fri 25 May 2007 10:46

Re: QGIS et PostgreSQL

#17 Fri 25 May 2007 10:49

Re: QGIS et PostgreSQL

#18 Fri 25 May 2007 11:39

Re: QGIS et PostgreSQL

#19 Tue 29 May 2007 12:35

Re: QGIS et PostgreSQL

#20 Thu 14 June 2007 10:15

Re: QGIS et PostgreSQL

Code:

#21 Tue 19 June 2007 09:54

Re: QGIS et PostgreSQL

Pied de page des forums