Le 18 décembre 2013 00:34, Christian Quest cqu...@openstreetmap.fr a
écrit :
Il va falloir trouver une astuce... j'ai testé département par département
et là c'est jouable.
... sauf si les limites côtières entre deux départements traversent
justement par un estuaire (ex. 35 et 44), mais ça
Le 15 décembre 2013 15:07, Christian Quest cqu...@openstreetmap.fr a
écrit :
Comme indiqué plus haut, celui-ci fonctionne bien dans la très grande
majorité des cas, mais sur certaines communes aux contours tarabiscotés, le
fait de simplifier les limites une à une peut faire qu'elles se
2013/12/16 Matthias Dietrich eiger@gmail.com:
J'ai des doutes sur Ferney-Voltaire. La route douanière ne me semble pas
faire partie de la commune, ni même du territoire français, même si l'usage
est réservé à la France. En tout cas, le cadastre ne contient pas cette
route.
On a le même
Bon... ça avance... j'ai maintenant l'intégralité des limites.
Le principe est le suivant:
- extraction des limites entre 2 communes (ST_Touches + ST_Intersection)
- extraction des limites côtières et frontalières
- extraction des communes non limitrophes (les iles)
- première génération de
Salut,
En complément, je pensais plutôt que si on avait des limites
topologiquement justes, on pourrait passer par MapShaper
http://www.mapshaper.org (c'est l'interface en ligne mais on peut aussi
faire de la ligne de commande https://github.com/mbloch/mapshaper)
L'outil corrige aussi les limites
Si quelqu'un veux essayer mapshaper, j'aimerai bien voir ce que ça donne
sur 36680 polygones et 9.5 million de nœuds...
Objectif (quasi) atteint de mon côté; j'ai 2 niveaux de simplification (100
et 250) qui sont ok et un autre avec encore un petit bug (simplification la
plus légère à 10).
J'ai
Le dimanche 15 décembre 2013 21:03:56, Christian Quest a écrit :
La seule chose absolument pas gérées ce sont des polygones en 8.
J'en ai identifié 2 et en reprenant le cadastre, le 8 n'avait pas lieu
d'être.
Attention à ne pas céder à la tentation du je corrige les données pour
compenser un
C'est pour ça que j'ai regardé le cadastre avant de corriger... mais le
manque dans l'algo est bien là.
Le 15 décembre 2013 21:36, sly (sylvain letuffe) lis...@letuffe.org a
écrit :
Le dimanche 15 décembre 2013 21:03:56, Christian Quest a écrit :
La seule chose absolument pas gérées ce sont
Bon, je progresse...
Si vous voulez voir la tête des requêtes c'est par ici (WARNING, ça peut
piquer les yeux):
https://gist.github.com/cquest/66797473e5663bb4ba43
Et graphiquement au final ça peut donner ça: http://cl.ly/image/2e2Y0l0K1Q0O
:)
mais aussi ça: http://cl.ly/image/3g2B2O2w3x2O :(
Salut Christian,
Je prends la discussion en cours de route donc excuse moi par avance si
la question t'a déjà été posée.
Mais as-tu pensé à avoir recours à la topologie ?
Ex -
http://blog.mathieu-leplatre.info/use-postgis-topologies-to-clean-up-road-networks.html
A.
On 13-12-13 03:27 PM,
En fait je voulais résoudre 2 problème en même temps.
Le premier: le rendu des limites admin ne me plait pas.
En effet, le rendu se base sur les polygones et pas sur les ways qui les
constituent. Du coup, pour chaque limite, on a en général 2 tracé, voir 4,
6 ou 8 pour les limites de département,
Je ne comprends pas le code cité, mais les communes limitrophes de bordure
de département ont des boundaries avec admin_level=6.
( (région, 4), respectivement).
Serait-ce une piste ?
Art.
Le 11 déc. 2013 16:59, Christian Quest cqu...@openstreetmap.fr a écrit
:
Les limites admins sont super
Je m'appuie uniquement sur les polygones de communes, pas sur les way
d'origine qui peuvent avoir de tout comme tag admin_level (y compris aucun).
Le 11 décembre 2013 17:04, Art Penteur art.pent...@gmail.com a écrit :
Je ne comprends pas le code cité, mais les communes limitrophes de bordure
J'me disait bien que c'était trop bête pour être ça.
Art.
Le 11 déc. 2013 17:07, Christian Quest cqu...@openstreetmap.fr a écrit :
Je m'appuie uniquement sur les polygones de communes, pas sur les way
d'origine qui peuvent avoir de tout comme tag admin_level (y compris aucun).
Le 11
De : Art Penteur
Le 11 déc. 2013 17:07, Christian Quest a écrit :
Le 11 décembre 2013 17:04, Art Penteur a écrit :
Le 11 déc. 2013 16:59, Christian Quest a écrit :
Ca semble ok sauf que j'ai des manques.
voici la tête de la requête si quelqu'un a une idée:
SELECT
En plus lisible :
SELECT
ST_Simplify(ST_LineMerge (st_intersection(p.way,p2.way)),10) as limite,
concat(p.name ,' - ', p2.name) as nom
from
planet_osm_polygon p
join planet_osm_polygon p2
on
(st_touches(p.way,p2.way)
and p2.ref:INSEE is not null
and p2.admin_level='8'
Je en comprends pas tout à ce que fait cette requête, mais pourquoi tu as
p2.osm_id ** p.osm_id dans ta clause ON ? Ca ne devrait pas être
p2.osm_id *=* p.osm_id, plutôt ?
Francescu
Le 11 décembre 2013 16:58, Christian Quest cqu...@openstreetmap.fr a
écrit :
Les limites admins sont super
!bbox! c'est pour TileMill (mapnik) car je l'utilise pour visualiser le
résultat.
La tolérance correspond au nombre de points intermédiaires pour définir un
quart de cercle. Donc avec 1 ça te simplifie un cercle en 8 points. Ca ne
joue que sur la simplification et pour être sûr je l'ai passé à
De : Christian Quest
!bbox! c'est pour TileMill (mapnik) car je l'utilise pour visualiser le
résultat.
Ah oui, le truc qui remplace QGis ;-)
La tolérance correspond au nombre de points intermédiaires pour définir un
quart de
cercle. Donc avec 1 ça te simplifie un cercle en 8 points. Ca
En remplaçant le ST_Touches par un simple j'ai ça:
http://cl.ly/image/1O3r3Z0y2K1l
Très bizarre...
J'ai donc regardé ce que me sortait mon ST_Intersection des polygones des
communes 14053 et 14206.
J'ai une GEOMETRYCOLLECTION avec plein de LINESTRING et deux POLYGON...
alors qu'en chargeant
Et voilà:
UPDATE planet_osm_rels SET pending=true FROM planet_osm_polygon where
id=-osm_id and ref:INSEE is not null;
41936 polygones seront recalculés au prochain update... donc d'ici quelques
minutes :)
Le 11 décembre 2013 18:19, Christian Quest cqu...@openstreetmap.fr a
écrit :
En
21 matches
Mail list logo