Postgis, question basique d'unités

Importer dans OSM, contrôler, suivre et surveiller (osmose, keepright, ...)
JBosm
Messages : 69
Inscription : lun. avr. 22, 2013 12:43 pm

Postgis, question basique d'unités

Message par JBosm » dim. févr. 10, 2019 4:30 pm

Bonjour,
J'ai une question qui paraitra probablement évidente, mais je n'ai pas manipulé beaucoup les fonctions de Postgis.
Si un projet est défini avec :

Code : Tout sélectionner

<Layer name="rotated_buildings" srs="+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0.0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs +over" status="on">
puis que la fonction suivante est appelée :

Code : Tout sélectionner

ST_SimplifyPreserveTopology(l.way,100)
l'unité 100 est-elle en mètres ou dans une autre unité ? 
(l'origine des fonctions est le projet OpenArdenneMap, https://github.com/nobohan/OpenArdenneM ... .xml#L6613)
JB.

Avatar de l’utilisateur
cquest
Messages : 1872
Inscription : ven. avr. 16, 2010 12:22 am
Localisation : Val de Marne
Contact :

Re: Postgis, question basique d'unités

Message par cquest » lun. févr. 11, 2019 10:06 am

Le paramètre tolerance est exprimé dans la même unité que la géométrie concernée.

Ici il s'agit de "metre" mercator... qui ne sont pas vraiment des metre (variable en fonction de la latitude).

Pour faire simple, un noeud est éliminé si la distance projetée perpendiculairement au segment sans ce noeud est inférieure à cette tolérance.

Voir ce pseudo code: https://en.wikipedia.org/wiki/Ramer%E2% ... Pseudocode

JBosm
Messages : 69
Inscription : lun. avr. 22, 2013 12:43 pm

Re: Postgis, question basique d'unités

Message par JBosm » lun. févr. 11, 2019 8:56 pm

Il s'agit donc bien de 100"mètres mercators", ça me semblait assez grand comparé à la taille d'une église. Est-ce qu'il existe un moyen rapide de retrouver la distance réelle "au sol", ou est-ce qu'il faut repasser par lat/lon et la calculer (pas bien compliqué non plus) ?
En fait, la simplification ne fait pas dans la dentelle avec ces paramètres-là. En passant par Shapely, avec le même paramètre et dans la même zone géographique, j'obtiens ça :
Image Image
Pour la méthode, je n'avais pas de soucis, il me semble que c'est la même utilisée dans JOSM.
Merci Christian,
JB.

Avatar de l’utilisateur
cquest
Messages : 1872
Inscription : ven. avr. 16, 2010 12:22 am
Localisation : Val de Marne
Contact :

Re: Postgis, question basique d'unités

Message par cquest » lun. févr. 11, 2019 10:50 pm

Le ratio entre "metre mercator" et réels, est d'environ 1.5 sous nos latitudes.
Donc si tu veux simplifier à 1m près ça donne une tolérance de 1.5.

Répondre

Qui est en ligne ?

Utilisateurs parcourant ce forum : Aucun utilisateur inscrit et 2 invités