Re: [OSM-dev-fr] Shapé les admin_level=boundary

2013-12-17 Par sujet kimaidou
Bonjour,

Pour info, j'avais fait un test rapide de Mapshaper pour simplifier les
régions françaises, et cela a été très rapide à faire et efficace. CF
fichier lien ci-dessous
http://we.tl/gnmRd84GMY

Bon après, comme a priori le nombre d’éléments après simplification peut
avoir baissé, il ne reste que le shp et le shx, le dbf a sauté. On peut
toujours retrouver les attributs via intersection du centroide des régions
d'origine avec les polygones simplifiés, mais ça demande une étape de plus.

Pour info en python il existe un outil qui fonctionne bien aussi :
https://github.com/mdouchin/Bloch
Testé et approuvé : cf lien ici avec tolérance utilisée : 1000 (après
reprojection en Lambert 93)
http://we.tl/OQsNXaO9dI
Dans ce cas les attributs sont bien conservés.

Bonne journée
Michael


Le 17 décembre 2013 10:01, Christian Quest cqu...@openstreetmap.fr a
écrit :

 Et ici ? http://neocarto.hypotheses.org/244 (CC-by-NC-SA)


 Le 16 décembre 2013 20:26, Rodolphe Quiédeville rodol...@quiedeville.orga 
 écrit :


 Yves yve...@gmail.com writes:

  Juste pour les frontières, natural earth ne te suffit pas ?

 Dans le jeu Admin 0 boundary lines, il y a les frontières, il faudrait
 alors que je merge avec un autre jeu pour obtenir des polygones pas
 pays.
 En tout cas natural earth est une bonne piste que j'avais oubliée.

 Merci



 
  Rodolphe Quiédeville rodol...@quiedeville.org a écrit :
 Bonjour,
 
 Est-ce que quelqu'un à par devers lui les limites administratives
 extraites d'OSM au format shapefile ? Soit les données, soit le script
 pour le faire, voir une méthode à implémenter.
 
 Pour les départements j'ai utilisé Geofla, mais il me faudrait les
 limites nationales des pays européens désormais.
 
 http://www.data.gouv.fr/DataSet/30383060
 
 Merci
 
 --
 Rodolphe Quiédeville
 Expert Tsung - Consulting en performance des SI
 http://blog.rodolphe.quiedeville.org, Tel : 06 13 79 63 41
 Travaillons libres, http://fr.lolix.org/
 
 ___
 dev-fr mailing list
 dev-fr@openstreetmap.org
 https://lists.openstreetmap.org/listinfo/dev-fr

 --
 Rodolphe Quiédeville
 Expert Tsung - Consulting en performance des SI
 http://blog.rodolphe.quiedeville.org, Tel : 06 13 79 63 41
 Travaillons libres, http://fr.lolix.org/

 ___
 dev-fr mailing list
 dev-fr@openstreetmap.org
 https://lists.openstreetmap.org/listinfo/dev-fr




 --
 Christian Quest - OpenStreetMap France
 Un nouveau serveur pour OSM... http://donate.osm.org/server2013/

 ___
 dev-fr mailing list
 dev-fr@openstreetmap.org
 https://lists.openstreetmap.org/listinfo/dev-fr


___
dev-fr mailing list
dev-fr@openstreetmap.org
https://lists.openstreetmap.org/listinfo/dev-fr


Re: [OSM-dev-fr] Shapé les admin_level=boundary

2013-12-17 Par sujet Christian Quest
Les extensions topo de postgis sont bien lentes. J'ai l'impression que ma
série de requête est quand même nettement plus rapide.

L'autre problème, c'est que j'ai déjà du mal à créer la topologie initiale
car ça me sort des erreurs de topo peu claires.
J'explore donc avec des choses plus light, comme les départements. Ca
devrait faire moins de géométries pour tester le principe.

Pour l'instant c'est avec les multipolygones que j'essaye de créer la
topologie, mais je pense passer sur les linestring de limites, au moins ça
aidera à identifier plus rapidement les endroits où il y a un souci de topo
dans les données de départ.

Le 17 décembre 2013 10:19, kimaidou kimai...@gmail.com a écrit :

 Oh tiens, je viens de voir que les régions simplifiées via Bloch ont
 laissé des plumes, comme l'enclave des Papes par exemple.. A creuser donc.
 Sinon, je pense que la solution PostGIS, si elle peut sembler lourde, sera
 sûrement la plus souple et robuste pour ce travail, comme l'a commencé
 Christian (tester aussi l'utilisation de la topologie dans PostGIS)

 Michael


-- 
Christian Quest - OpenStreetMap France
Un nouveau serveur pour OSM... http://donate.osm.org/server2013/
___
dev-fr mailing list
dev-fr@openstreetmap.org
https://lists.openstreetmap.org/listinfo/dev-fr


Re: [OSM-dev-fr] Shapé les admin_level=boundary

2013-12-17 Par sujet Christian Quest
Je confirme, très très lent... 69 secondes pour transformer 8 linestring en
topo... à ce rythme, vu qu'il y en a 109902, ça va mettre plus de 10j  !

Donc si je veux avoir des limites simplifiées d'ici demain 16h30 je
retourne à mes requêtes ;)


Le 17 décembre 2013 11:27, Christian Quest cqu...@openstreetmap.fr a
écrit :

 Les extensions topo de postgis sont bien lentes. J'ai l'impression que ma
 série de requête est quand même nettement plus rapide.

 L'autre problème, c'est que j'ai déjà du mal à créer la topologie initiale
 car ça me sort des erreurs de topo peu claires.
 J'explore donc avec des choses plus light, comme les départements. Ca
 devrait faire moins de géométries pour tester le principe.

 Pour l'instant c'est avec les multipolygones que j'essaye de créer la
 topologie, mais je pense passer sur les linestring de limites, au moins ça
 aidera à identifier plus rapidement les endroits où il y a un souci de topo
 dans les données de départ.

 Le 17 décembre 2013 10:19, kimaidou kimai...@gmail.com a écrit :

 Oh tiens, je viens de voir que les régions simplifiées via Bloch ont
 laissé des plumes, comme l'enclave des Papes par exemple.. A creuser donc.
 Sinon, je pense que la solution PostGIS, si elle peut sembler lourde,
 sera sûrement la plus souple et robuste pour ce travail, comme l'a commencé
 Christian (tester aussi l'utilisation de la topologie dans PostGIS)

 Michael


 --
 Christian Quest - OpenStreetMap France
 Un nouveau serveur pour OSM... http://donate.osm.org/server2013/




-- 
Christian Quest - OpenStreetMap France
Un nouveau serveur pour OSM... http://donate.osm.org/server2013/
___
dev-fr mailing list
dev-fr@openstreetmap.org
https://lists.openstreetmap.org/listinfo/dev-fr


Re: [OSM-dev-fr] Shapé les admin_level=boundary

2013-12-17 Par sujet Christophe Merlet
Le 17/12/2013 11:43, Christian Quest a écrit :
 Je confirme, très très lent... 69 secondes pour transformer 8 linestring
 en topo... à ce rythme, vu qu'il y en a 109902, ça va mettre plus de 10j  !
 
 Donc si je veux avoir des limites simplifiées d'ici demain 16h30 je
 retourne à mes requêtes ;)

Je viens de simplifier manuellement les limites de régions une par une à
1m (quasi l'épaisseur du trait dans le cadastre) dans JOSM en moins de 3 mn

Statistiquement ça a supprimé 130 833 nœuds sur 619 858. sur 7759 chemins.

Il serait sans doute intéressant de pouvoir extraire de JOSM cet algo,
ou rendre JOSM scriptable sans avoir besoin d'afficher le fichier osm à
l'écran.


2ème point. Il serait bon de simplifier dans OSM ces limites de communes
abusivement et inutilement détaillées.  Les limites de 2 communes
limitrophes ne se superposent déjà pas dans le cadastre, alors ce n'est
pas en les dessinant avec énormément de points, quelles seront plus
précises.

Je vais faire une simplification à 10m et à 100m pour voir...


Librement,
-- 
Christophe Merlet (RedFox)

___
dev-fr mailing list
dev-fr@openstreetmap.org
https://lists.openstreetmap.org/listinfo/dev-fr


Re: [OSM-dev-fr] Shapé les admin_level=boundary

2013-12-17 Par sujet Christophe Merlet
Le 17/12/2013 12:01, Christophe Merlet a écrit :
 Le 17/12/2013 11:43, Christian Quest a écrit :
 Je confirme, très très lent... 69 secondes pour transformer 8 linestring
 en topo... à ce rythme, vu qu'il y en a 109902, ça va mettre plus de 10j  !

 Donc si je veux avoir des limites simplifiées d'ici demain 16h30 je
 retourne à mes requêtes ;)
 
 Je viens de simplifier manuellement les limites de régions une par une à
 1m (quasi l'épaisseur du trait dans le cadastre) dans JOSM en moins de 3 mn
 
 Statistiquement ça a supprimé 130 833 nœuds sur 619 858. sur 7759 chemins.
 
 Il serait sans doute intéressant de pouvoir extraire de JOSM cet algo,
 ou rendre JOSM scriptable sans avoir besoin d'afficher le fichier osm à
 l'écran.
 
 
 2ème point. Il serait bon de simplifier dans OSM ces limites de communes
 abusivement et inutilement détaillées.  Les limites de 2 communes
 limitrophes ne se superposent déjà pas dans le cadastre, alors ce n'est
 pas en les dessinant avec énormément de points, quelles seront plus
 précises.
 
 Je vais faire une simplification à 10m et à 100m pour voir...

à 10m je supprime 435 367 noeuds...
à 100m, j'en supprime 578 441 soit plus de 93%

Pour faire mieux dans JOSM, il faudrait refusionner les chemins qui
composent les limites de régions avant de lancer la simplification.


Librement,
-- 
Christophe Merlet (RedFox)

___
dev-fr mailing list
dev-fr@openstreetmap.org
https://lists.openstreetmap.org/listinfo/dev-fr


Re: [OSM-dev-fr] Shapé les admin_level=boundary

2013-12-17 Par sujet Pieren
2013/12/17 Christophe Merlet red...@redfoxcenter.org:

 2ème point. Il serait bon de simplifier dans OSM ces limites de communes
 abusivement et inutilement détaillées.

Finalement, c'est comme dans OSM ;-) Ils ont des contributeurs à
précision géométrique variable ;-)

 Les limites de 2 communes
 limitrophes ne se superposent déjà pas dans le cadastre, alors ce n'est
 pas en les dessinant avec énormément de points, quelles seront plus
 précises.

Pour quelqu'un qui cherche sa parcelle, le niveau de détail peut avoir
son importance (c'est toujours une question d'échelle). Et on peut
toujours simplifier alors que l'inverse n'est pas possible. Le fait
qu'elles ne se superposent pas d'une commune à l'autre n'est pas lié à
leur niveau de détails. En plus, ça serait directement faire mentir le
communiqué de presse qui vante le nombre de points 15 fois supérieur à
geofla ;-)

Pieren

___
dev-fr mailing list
dev-fr@openstreetmap.org
https://lists.openstreetmap.org/listinfo/dev-fr


Re: [OSM-dev-fr] Shapé les admin_level=boundary

2013-12-17 Par sujet Christophe Merlet
Le 17/12/2013 12:15, Pieren a écrit :
 2013/12/17 Christophe Merlet red...@redfoxcenter.org:
 
 2ème point. Il serait bon de simplifier dans OSM ces limites de communes
 abusivement et inutilement détaillées.
 
 Finalement, c'est comme dans OSM ;-) Ils ont des contributeurs à
 précision géométrique variable ;-)
 
 Les limites de 2 communes
 limitrophes ne se superposent déjà pas dans le cadastre, alors ce n'est
 pas en les dessinant avec énormément de points, quelles seront plus
 précises.
 
 Pour quelqu'un qui cherche sa parcelle, le niveau de détail peut avoir
 son importance (c'est toujours une question d'échelle). Et on peut
 toujours simplifier alors que l'inverse n'est pas possible. Le fait
 qu'elles ne se superposent pas d'une commune à l'autre n'est pas lié à
 leur niveau de détails. En plus, ça serait directement faire mentir le
 communiqué de presse qui vante le nombre de points 15 fois supérieur à
 geofla ;-)

On peut dessiner une ligne droite avec 30 points, elle ne sera pas plus
précise que celle qui se contente de 2 points.


Librement,
-- 
Christophe Merlet (RedFox)

___
dev-fr mailing list
dev-fr@openstreetmap.org
https://lists.openstreetmap.org/listinfo/dev-fr


Re: [OSM-dev-fr] Shapé les admin_level=boundary

2013-12-17 Par sujet Christian Quest
Ah si... on est sûr qu'elle est droite et pas simplifiée ;)


Le 17 décembre 2013 12:23, Christophe Merlet red...@redfoxcenter.org a
écrit :

 Le 17/12/2013 12:15, Pieren a écrit :
  2013/12/17 Christophe Merlet red...@redfoxcenter.org:
 
  2ème point. Il serait bon de simplifier dans OSM ces limites de communes
  abusivement et inutilement détaillées.
 
  Finalement, c'est comme dans OSM ;-) Ils ont des contributeurs à
  précision géométrique variable ;-)
 
  Les limites de 2 communes
  limitrophes ne se superposent déjà pas dans le cadastre, alors ce n'est
  pas en les dessinant avec énormément de points, quelles seront plus
  précises.
 
  Pour quelqu'un qui cherche sa parcelle, le niveau de détail peut avoir
  son importance (c'est toujours une question d'échelle). Et on peut
  toujours simplifier alors que l'inverse n'est pas possible. Le fait
  qu'elles ne se superposent pas d'une commune à l'autre n'est pas lié à
  leur niveau de détails. En plus, ça serait directement faire mentir le
  communiqué de presse qui vante le nombre de points 15 fois supérieur à
  geofla ;-)

 On peut dessiner une ligne droite avec 30 points, elle ne sera pas plus
 précise que celle qui se contente de 2 points.


 Librement,
 --
 Christophe Merlet (RedFox)

 ___
 dev-fr mailing list
 dev-fr@openstreetmap.org
 https://lists.openstreetmap.org/listinfo/dev-fr




-- 
Christian Quest - OpenStreetMap France
Un nouveau serveur pour OSM... http://donate.osm.org/server2013/
___
dev-fr mailing list
dev-fr@openstreetmap.org
https://lists.openstreetmap.org/listinfo/dev-fr


Re: [OSM-dev-fr] Shapé les admin_level=boundary

2013-12-17 Par sujet Christophe Merlet
Le 17/12/2013 12:24, Christian Quest a écrit :
 Ah si... on est sûr qu'elle est droite et pas simplifiée ;)

Toute la subtile différence entre résolution et précision :)

 Le 17 décembre 2013 12:23, Christophe Merlet red...@redfoxcenter.org
 mailto:red...@redfoxcenter.org a écrit :
 
 Le 17/12/2013 12:15, Pieren a écrit :
  2013/12/17 Christophe Merlet red...@redfoxcenter.org
 mailto:red...@redfoxcenter.org:
 
  2ème point. Il serait bon de simplifier dans OSM ces limites de
 communes
  abusivement et inutilement détaillées.
 
  Finalement, c'est comme dans OSM ;-) Ils ont des contributeurs à
  précision géométrique variable ;-)
 
  Les limites de 2 communes
  limitrophes ne se superposent déjà pas dans le cadastre, alors ce
 n'est
  pas en les dessinant avec énormément de points, quelles seront plus
  précises.
 
  Pour quelqu'un qui cherche sa parcelle, le niveau de détail peut avoir
  son importance (c'est toujours une question d'échelle). Et on peut
  toujours simplifier alors que l'inverse n'est pas possible. Le fait
  qu'elles ne se superposent pas d'une commune à l'autre n'est pas lié à
  leur niveau de détails. En plus, ça serait directement faire mentir le
  communiqué de presse qui vante le nombre de points 15 fois supérieur à
  geofla ;-)
 
 On peut dessiner une ligne droite avec 30 points, elle ne sera pas plus
 précise que celle qui se contente de 2 points.
 
 
 Librement,
 --
 Christophe Merlet (RedFox)
 
 ___
 dev-fr mailing list
 dev-fr@openstreetmap.org mailto:dev-fr@openstreetmap.org
 https://lists.openstreetmap.org/listinfo/dev-fr
 
 
 
 
 -- 
 Christian Quest - OpenStreetMap France
 Un nouveau serveur pour OSM... http://donate.osm.org/server2013/
 
 
 ___
 dev-fr mailing list
 dev-fr@openstreetmap.org
 https://lists.openstreetmap.org/listinfo/dev-fr
 


-- 
Christophe Merlet (RedFox)

___
dev-fr mailing list
dev-fr@openstreetmap.org
https://lists.openstreetmap.org/listinfo/dev-fr


Re: [OSM-dev-fr] Shapé les admin_level=boundary

2013-12-17 Par sujet kimaidou
Un autre test rapide avec Mapshaper :
* conversion préalable du Shapefile des communes du département 34 via QGIS
en GeoJSON
* Simplification avec Mapshaper avec une tolérance de 10%
La conversion est instantannée et est réalisée en ligne. Comme c'est du
GeoJson, Mapshaper a pu conserver les attributs.

Fichier source : 6,5 Mo
Fichier résultat : 607 ko

http://we.tl/hoS16168Ii


Le 17 décembre 2013 13:55, Ab_fab gamma@gmail.com a écrit :

 Plus d'infos sur les scripts dans JOSM, mais pour un toute autre sujet,
 via ce message sur la liste talk-fr
 https://lists.openstreetmap.org/pipermail/talk-fr/2011-November/037702.html


 Le 17 décembre 2013 13:34, Ab_fab gamma@gmail.com a écrit :

 JOSM est scriptable (python et java), mais je ne sais pas si tu peux
 éviter l'affichage à l'écran
  Le 17/12/2013 11:43, Christian Quest a écrit :
  Je confirme, très très lent... 69 secondes pour transformer 8 linestring
  en topo... à ce rythme, vu qu'il y en a 109902, ça va mettre plus de
 10j  !
 
  Donc si je veux avoir des limites simplifiées d'ici demain 16h30 je
  retourne à mes requêtes ;)

 Je viens de simplifier manuellement les limites de régions une par une à
 1m (quasi l'épaisseur du trait dans le cadastre) dans JOSM en moins de 3
 mn

 Statistiquement ça a supprimé 130 833 nœuds sur 619 858. sur 7759 chemins.

 Il serait sans doute intéressant de pouvoir extraire de JOSM cet algo,
 ou rendre JOSM scriptable sans avoir besoin d'afficher le fichier osm à
 l'écran.


 2ème point. Il serait bon de simplifier dans OSM ces limites de communes
 abusivement et inutilement détaillées.  Les limites de 2 communes
 limitrophes ne se superposent déjà pas dans le cadastre, alors ce n'est
 pas en les dessinant avec énormément de points, quelles seront plus
 précises.

 Je vais faire une simplification à 10m et à 100m pour voir...


 Librement,
 --
 Christophe Merlet (RedFox)

 ___
 dev-fr mailing list
 dev-fr@openstreetmap.org
 https://lists.openstreetmap.org/listinfo/dev-fr




 --
 ab_fab http://wiki.openstreetmap.org/wiki/User:Ab_fab
 Il n'y a pas de pas perdus, Nadja

 ___
 dev-fr mailing list
 dev-fr@openstreetmap.org
 https://lists.openstreetmap.org/listinfo/dev-fr


___
dev-fr mailing list
dev-fr@openstreetmap.org
https://lists.openstreetmap.org/listinfo/dev-fr


Re: [OSM-dev-fr] Simplifier nos limites admin...

2013-12-17 Par sujet Philippe Verdy
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 dépend jusqu'où on fait
remonter l'estuaire maritime pour les limites des communes littorales.

En revanche pour les limites de régions et départements (i.e. les limites
préfectorales, pas celles des collectivités locales de la région et du
département, qui n'ont pas compétence du domaine maritime et même pas la
compétence totale sur le littoral terrestre protégé) et les limites
nationales, on n'a pas de mal à définir leur domaine maritime en utilisant
les limites dans les eaux territoriales.

C'est bien pour ça qu'on devrait supporter en France aussi ce que font les
Allemands depuis longtemps (et d'autres pays aussi) : les land_areas
définis sur les lignes de côtes et non les limites des eaux territoriales
nationales.

Peut-être que ces land_areas ne sont pas documentés sur le wiki en anglais
mais ça l'est parfaitement en allemand et il y a assez d'explication et de
notes d'alertes dans les données OSM pour expliquer la différence avec les
boundary=*, et ils en ont aussi discuté depuis longtemps de la nécessité
d'avoir ces land_areas (qui ne sont pas seulement un besoin administratif,
mais aussi un clair besoin au plan géographique/cartographique, même si en
théorie on nos lignes de côtes devraient être celles de la ligne de base
légale, qui traverse aussi les estuaires et inclut les ports presque fermés
par une digue de protection et toute zone d'eau de moins de 50 mètres de
large environ, soit l'ordre de grandeur des marées le long des plages
atlantiques et de la Manche ou la Mer du Nord, cette distance étant plus
réduite en Méditerranée, mais variant dans les estuaires selon les régimes
des fleuves pour ce qui concerne les limites des estuaires)..

Mais c'est vrai qu'un land_area pour la France demanderait un multipolygone
très complexe avec de nombreux membres, tellement nos cotes sont découpées
(surtout en Bretagne).

Actuellement on ne sait toujours pas si en France les limites de nos
régions doivent être celles des collectivités territoriales (excluant les
eaux territoriales), ce qui semble être le cas, ou les limites
(sous-)préfectorales (en fait de la compétence exclusive de l'Etat, lequel
inclut le domaine maritime). Selon les cas on aurait besoin d'avoir les
deux aussi en France, les limites administratives (de l'Etat) incluant les
eaux territoriales, et une autre relation pour les land_areas (dans la
compétence réelle des régions, départements et communes en tant que
collectivités locales mais qui ne sont pas réellement un découpage complet
du pays).
___
dev-fr mailing list
dev-fr@openstreetmap.org
https://lists.openstreetmap.org/listinfo/dev-fr