Le portail francophone de la géomatique

Syphys · Mon 10 May 2010 00:39

Bonjour,

J'ai une base de donnée postgis contenant les villes de France mais je n'arrive pas à comprendre en quelle unité sont retournés les résultats

J'ai utilisé le srid 2154 pour mettre mes données que j'ai pris à cette adresse : http://download.geofabrik.de/osm/europe/france/

quand je fais une requéte style

select distinct distance((select the_geom from villes v where v.name = 'Strasbourg' and type= 'city'),(select the_geom from villes v where v.name = 'Brest' and v.type= 'city')) from villes;

j'obtiens :
12.2371737413677

Si je fais des essaies avec plein de distance inter-ville je vois une correspondance dans l'ordre de grandeur mais ce n'est pas proportionnel.

Je vous remercie ne-serait-que pour l'espoir que l'attente d'une réponse me procure... Je sais je suis mélodramatique quoique... :'s

Guillaume Sueur · Mon 10 May 2010 01:13

Bonsoir,

Les données doivent être en WGS84, EPSG:4326, car cette distance Brest - Strasbourg ressemble à des degrés décimaux ! Essayez :
UPDATE villes SET the_geom = ST_Transform(ST_SetSRID(the_geom,4326),2154);
et refaites vos calculs de distance.

Guillaume

Sylvie · Mon 10 May 2010 07:50

Bonjour,

J'ai déplacé ce sujet sur Postgis sur le forum Géo'BD

Sylvie, pour la modération

Syphys · Mon 10 May 2010 09:29

C' est super ça marche!!!!!

Merci beaucoup!

Nicolas Ribot · Mon 10 May 2010 09:52

Syphys a écrit:

Bonjour,

J'ai une base de donnée postgis contenant les villes de France mais je n'arrive pas à comprendre en quelle unité sont retournés les résultats

J'ai utilisé le srid 2154 pour mettre mes données que j'ai pris à cette adresse : http://download.geofabrik.de/osm/europe/france/

quand je fais une requéte style

select distinct distance((select the_geom from villes v where v.name = 'Strasbourg' and type= 'city'),(select the_geom from villes v where v.name = 'Brest' and v.type= 'city')) from villes;

j'obtiens :
12.2371737413677

Si je fais des essaies avec plein de distance inter-ville je vois une correspondance dans l'ordre de grandeur mais ce n'est pas proportionnel.

Je vous remercie ne-serait-que pour l'espoir que l'attente d'une réponse me procure... Je sais je suis mélodramatique quoique... :'s

Bonjour,

Postgis fait les calculs dans les unités des données.

Nico

Yves · Mon 10 May 2010 10:53

Nicolas Ribot a écrit:

Syphys a écrit:
Bonjour,

J'ai une base de donnée postgis contenant les villes de France mais je n'arrive pas à comprendre en quelle unité sont retournés les résultats

J'ai utilisé le srid 2154 pour mettre mes données que j'ai pris à cette adresse : http://download.geofabrik.de/osm/europe/france/

quand je fais une requéte style

select distinct distance((select the_geom from villes v where v.name = 'Strasbourg' and type= 'city'),(select the_geom from villes v where v.name = 'Brest' and v.type= 'city')) from villes;

j'obtiens :
12.2371737413677

Si je fais des essaies avec plein de distance inter-ville je vois une correspondance dans l'ordre de grandeur mais ce n'est pas proportionnel.

Je vous remercie ne-serait-que pour l'espoir que l'attente d'une réponse me procure... Je sais je suis mélodramatique quoique... :'s

Bonjour,

Postgis fait les calculs dans les unités des données.

Nico

Je dirais même plus : dans les unités de la projection

Y.

Nicolas Ribot · Mon 10 May 2010 11:23

Yves a écrit:

Nicolas Ribot a écrit:
Syphys a écrit:
Bonjour,

J'ai une base de donnée postgis contenant les villes de France mais je n'arrive pas à comprendre en quelle unité sont retournés les résultats

J'ai utilisé le srid 2154 pour mettre mes données que j'ai pris à cette adresse : http://download.geofabrik.de/osm/europe/france/

quand je fais une requéte style

select distinct distance((select the_geom from villes v where v.name = 'Strasbourg' and type= 'city'),(select the_geom from villes v where v.name = 'Brest' and v.type= 'city')) from villes;

j'obtiens :
12.2371737413677

Si je fais des essaies avec plein de distance inter-ville je vois une correspondance dans l'ordre de grandeur mais ce n'est pas proportionnel.

Je vous remercie ne-serait-que pour l'espoir que l'attente d'une réponse me procure... Je sais je suis mélodramatique quoique... :'s

Bonjour,

Postgis fait les calculs dans les unités des données.

Nico

Je dirais même plus : dans les unités de la projection

Y.

Non, pas tout a fait: si une colonne geometrique contient des données metriques, (disons lambert 2), mais que lors de l'import un mauvais SRID a ete affecté a ses données (disons WGS84), le calcul se fera bien en metres. Autrement dit, postgis ne tient pas compte du SRID de la geometries pour connaitres les unités de calcul.
Au detail pres du nouveau type GEOGRAPHY, qui "sait" que les données sont dans un systeme spherique, avec un modele particulier de calcul, a condition qu'un SRID geodétique soit associé.

C'est a l'utilisateur de garantir que le SRID associé a des données soit valide et corresponde aux données.
Il faut donc y faire attention au moment de l'import de données dans postgis.

Nico

Syphys · Mon 10 May 2010 16:10

Toutes les questions que je me posais

Par contre je me demande comment on fait pour savoir le srid avec lequel on doit importer les données?

Dernière modification par Syphys (Mon 10 May 2010 16:12)

Yves · Mon 10 May 2010 16:20

Merci Nico, ma précision a pour origine la doc de postgis qui effectivement doit supposer que les données sont dans le srid correcte.

Y.

Nicolas Ribot · Mon 10 May 2010 16:33

Syphys a écrit:

Toutes les questions que je me posais

Par contre je me demande comment on fait pour savoir le srid avec lequel on doit importer les données?

Le fournisseur du shape DOIT avoir livré des metadonnées décrivant ledit shape (projection entr'autre)
Si un fichier .prj est present, il contient le WKT de la projection du shapefile.

Si aucune info n'est presente, associée au shape, alors ce n'est pas facile: données d'un pays => essayer les projections du pays, eventuellement en superposant des données dont on connait la projection.

Les coordonnées peuvent aider (metrique, lat-long ?)

Mais bon, faut vraiment gronder celui qui produit le shape s'il n'a pas les info de projections associées.

Nico

Syphys · Mon 10 May 2010 16:52

:$ :$ :$

C'est moi qui savait pas où regarder je voyais bien les autres fichiers autour mais je savais pas à quoi ils servaient :s

Syphys · Mon 10 May 2010 17:01

Comme vous avez l'air bien au courant de tout ça, je me permet de vous poser une question sans avoir vraiment regardé (parce que je sais pas trop où chercher...)

Savez-vous comment je peux entrer dans ma base des points dont j'ai la longitude et la latitude?

J'en aurais besoin plus tard mais ça me turlupine déjà...

Bon avec un tout petit peu de recherche il semble que j'ai trouvé ça ici:

http://www.postgis.fr/node/351

Merci encore de m'avoir débloqué et pour les explications

Dernière modification par Syphys (Mon 10 May 2010 17:10)

Nicolas Ribot · Mon 10 May 2010 17:43

Syphys a écrit:

Comme vous avez l'air bien au courant de tout ça, je me permet de vous poser une question sans avoir vraiment regardé (parce que je sais pas trop où chercher...)

Savez-vous comment je peux entrer dans ma base des points dont j'ai la longitude et la latitude?

J'en aurais besoin plus tard mais ça me turlupine déjà...

Bon avec un tout petit peu de recherche il semble que j'ai trouvé ça ici:

http://www.postgis.fr/node/351

Merci encore de m'avoir débloqué et pour les explications

De rien

L'ideal pour la fabrication de POINT postgis a partir de colonnes lat long est d'utiliser makePoint(x, y):

update matable set macolonnegeo = setSrid(st_makepoint(colLongitude, colLatitude), <le srid de ma colonne>);

Plus rapide et plus sur que la concatenation de string.

Nicolas

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 10 May 2010 00:39

Unité distance postgis

#2 Mon 10 May 2010 01:13

Re: Unité distance postgis

#3 Mon 10 May 2010 07:50

Re: Unité distance postgis

#4 Mon 10 May 2010 09:29

Re: Unité distance postgis

#5 Mon 10 May 2010 09:52

Re: Unité distance postgis

#6 Mon 10 May 2010 10:53

Re: Unité distance postgis

#7 Mon 10 May 2010 11:23

Re: Unité distance postgis

#8 Mon 10 May 2010 16:10

Re: Unité distance postgis

#9 Mon 10 May 2010 16:20

Re: Unité distance postgis

#10 Mon 10 May 2010 16:33

Re: Unité distance postgis

#11 Mon 10 May 2010 16:52

Re: Unité distance postgis

#12 Mon 10 May 2010 17:01

Re: Unité distance postgis

#13 Mon 10 May 2010 17:43

Re: Unité distance postgis

Pied de page des forums