banniere

Le portail francophone de la géomatique


Toujours pas inscrit ? Mot de passe oublié ?
Nom d'utilisateur    Mot de passe              Toujours pas inscrit ?   Mot de passe oublié ?

Annonce

Pour sécuriser votre compte sur les forums du GeoRezo, nous demandons de changer votre mot de passe.

Vous allez recevoir un message pour effectuer ce changement de mot de passe.

Merci de bien respecter les règles préconisées.

#1 Wed 19 August 2009 09:31

Maxime31
Participant assidu
Date d'inscription: 3 Sep 2008
Messages: 175

[MapServer / PostGIS] Problème de projection

Bonjour à tous,

J'ai un petit souci au niveau de PostGIS apparemment du à un problème de projection.

J'avais intégré dans PostGIS des shapefiles avec l'utilitaire shp2pgsql en ligne de commande. Cependant, je n'avais pas précisé le système de projection (avec l'option -S). Donc, en regardant dans la table geometry_columns de ma base PostGIS, j'avais tout les SRID à -1.

Pour plus de rigueur, j'ai remis tous les SRID à 27572 avec une fonction setSrid() sur toutes mes tables.

Mais, maintenant je veux réaliser cette requête :

Code:

SELECT * FROM matable WHERE Intersects(the_geom,'POINT("......,......")')

Et PgAdmin me renvoie cette erreur : ERROR:  Operation on mixed SRID geometries.

Alors j'ai essayé de réintégré toutes mes tables avec l'utilitaire shp2pgsql, cette fois-ci en précisant bien le système de projection que je désirais, mais j'ai toujours la même erreur.

Je précise que lorsque j'avais mes SRID à -1, cette requête marchait très bien.

Est-ce normal ?

Hors ligne

 

#2 Wed 19 August 2009 09:41

Yves
Membre du bureau
Lieu: Aix-les-Bains
Date d'inscription: 22 Mar 2006
Messages: 9828
Site web

Re: [MapServer / PostGIS] Problème de projection

Oui c'est normal il faut préciser la projection pour la géométrie point que tu créés :

Code:

SELECT * FROM matable WHERE Intersects(the_geom,GeomFromText('POINT("......,......")',27572))

Y.


Yves Jacolin, bénévole de l'association GeoRezo.net, agit au nom et pour le compte de l'association - Partageons ce qui nous départage !!  - GeoRezo vous aide ? Aidez GeoRezo !

Hors ligne

 

#3 Wed 19 August 2009 09:54

Maxime31
Participant assidu
Date d'inscription: 3 Sep 2008
Messages: 175

Re: [MapServer / PostGIS] Problème de projection

Bonjour Yves, merci de ta réponse rapide.

Alors je ne dois pas avoir la bonne syntaxe parce que je n'arrive pas à préciser le système de projections dasn ma requête.

A la base, ma requête était celle-ci :

Code:

SELECT * FROM communes WHERE Intersects(the_geom,'POINT(528251.35714286 1845315.7142857)')

Mais j'ai essayé :

Code:

SELECT * FROM communes WHERE Intersects(the_geom,'POINT(528251.35714286 1845315.7142857)',27572)

ou :

Code:

SELECT * FROM communes WHERE Intersects(the_geom,'POINT(528251.35714286 1845315.7142857),27572')

Mais sans succès...

Hors ligne

 

#4 Wed 19 August 2009 09:58

Maxime31
Participant assidu
Date d'inscription: 3 Sep 2008
Messages: 175

Re: [MapServer / PostGIS] Problème de projection

Ah non pardon j'avais mal lu ta requête...désolé.

Donc je viens d'essayer cela :

Code:

SELECT * FROM communes WHERE Intersects(the_geom,GeomFromText('POINT("528251.35714286,1845315.7142857")',27572))

Mas j'ai une erreur : ERROR:  parse error - invalid geometry
CONTEXT:  SQL function "geomfromtext" statement 1

Hors ligne

 

#5 Wed 19 August 2009 10:02

Maxime31
Participant assidu
Date d'inscription: 3 Sep 2008
Messages: 175

Re: [MapServer / PostGIS] Problème de projection

Ah c'ets bon ca marche !

Il fallait juste enlever les guillemets et la virgule pour les coordonnées du point.

Code:

SELECT * FROM communes WHERE Intersects(the_geom,GeomFromText('POINT(528251.35714286 1845315.7142857)',27572))

Voila, désolé pour les posts à répétition.

Dernière modification par Maxime31 (Wed 19 August 2009 10:09)

Hors ligne

 

#6 Wed 19 August 2009 10:06

nicolas.leclerc
Participant actif
Date d'inscription: 8 Apr 2009
Messages: 78

Re: [MapServer / PostGIS] Problème de projection

Bonjour

Et si tu fais :

Code:

SELECT * FROM communes WHERE Intersects(the_geom,('POINT("528251.35714286,1845315.7142857")',27572))

ha croisement...

Dernière modification par nicolas.leclerc (Wed 19 August 2009 10:07)

Hors ligne

 

#7 Wed 09 September 2009 15:24

Maxime31
Participant assidu
Date d'inscription: 3 Sep 2008
Messages: 175

Re: [MapServer / PostGIS] Problème de projection

Bonjour à tous,

Je poste ce message sur l'une de mes dernières discussions parce que je ne savais pas exactement où poster d'autre.
Donc juste pour dire que mon projet Webmapping est terminé et que je remercie toutes les personnes qui m'ont aidé sur ce forum, et notamment Yves, par leurs conseils précieux.

Encore une fois merci à tous et à bientôt.

Hors ligne

 

Pied de page des forums

Powered by FluxBB