Pages: 1
- Sujet précédent - OGR : Une conversion de ESRI à MapInfo génère des recouvrements - Sujet suivant
#1 Wed 05 March 2014 16:48
OGR : Une conversion de ESRI à MapInfo génère des recouvrements
Bonjour,
Il y a une couche ESRI Shapefile dont j'ai traité automatiquement les recouvrements avec l'excellent outil pprepair
Je détecte les recouvrements en intégrant ma couche dans GRASS (La couche des recouvrements porte 2 comme identifiant)
La couche résultant de pprepair au format ESRI ne présente que 3 recouvrements contre plusieurs milliers auparavant
En la convertissant au format MapInfo avec ogr2ogr, la couche MapInfo obtenue a 3353 recouvrements, ce qui sous-entend que les coordonnées des noeuds ne seraient pas fidèlement reproduits suite à la conversion depuis ESRI par ogr2ogr ?
Les géométries sont-elles bien conservées suite au changement de format ?
Dernière modification par MathieuR (Wed 05 March 2014 16:51)
geodata au cerema et petits billets en géomatique
Hors ligne
#2 Wed 05 March 2014 17:25
- dominique.lys
- Participant assidu
- Date d'inscription: 5 Oct 2006
- Messages: 473
- Site web
Re: OGR : Une conversion de ESRI à MapInfo génère des recouvrements
Bonjour
voir ces discusions:
http://georezo.net/forum/viewtopic.php?id=89081
http://georezo.net/forum/viewtopic.php?id=88741
visiblement une perte de précision dans les coordonnées lorsque l'on passe au format Mapinfo lié à une mauvaise définition des bornes dans les projections utilisées par Mapinfo.
Hors ligne
#3 Wed 05 March 2014 17:46
Re: OGR : Une conversion de ESRI à MapInfo génère des recouvrements
Merci pour les liens.
Oui, en effet, j'obtiens ceci pour les coordonnées des noeuds
Pour le shape d'origine :
Code:
POLYGON((933459.60252128774300218 1820336.15873935609124601 [...]
Pour le TAB converti :
Code:
POLYGON((933459.59999999997671694 1820336.15999999991618097 [...]
Je ne vois pas trop comment procéder pour avoir une couche correcte au format MapInfo du coup..
et je n'ai pas trouvé de réponse dans les deux fils cités..
geodata au cerema et petits billets en géomatique
Hors ligne
#4 Wed 05 March 2014 18:06
- dominique.lys
- Participant assidu
- Date d'inscription: 5 Oct 2006
- Messages: 473
- Site web
Re: OGR : Une conversion de ESRI à MapInfo génère des recouvrements
En effet pas de solution en vue
Peut être essayer avec un format intermédiaire selon ce que le traducteur universelle de mapinfo est capable d'importer
Hors ligne
#5 Thu 06 March 2014 08:39
- Nicolas Granier
- Participant assidu
- Date d'inscription: 19 Apr 2007
- Messages: 271
Re: OGR : Une conversion de ESRI à MapInfo génère des recouvrements
Bonjour,
A l'époque j'avais contacté un revendeur Mapinfo justement pour savoir comment étaient gérées les coordonnées dans Mapinfo.
Le résultat c'est que Mapinfo utilise en effet ses systèmes de coordonnées pour affecter un nombre de décimales dans les coordonnées des objets vecteurs.
je m'explique : on retrouve 3 "formats" de projection (si mes souvenirs sont bons), borne Monde, borne Europe et borne France.
En borne Monde, les cordonnées seront stockées sans décimales c'est à dire avec une précision métrique, en borne Europe avec 2 décimales par exemple et en borne France avec 8 décimales. Je détaille ici le principe, le nombre de décimales n'est pas réel, c'est juste pour me faire comprendre.
Que ce passe t-il si vous convertissez du shape vers du tab en forçant la projection par un -a_srs suivi du prj d'Arcgis? La projection ne sera pas reconnue dans Mapinfo, mais est-ce que les coordonnées seront bonnes ?
Nicolas GRANIER
Hors ligne
#6 Fri 07 March 2014 11:05
Re: OGR : Une conversion de ESRI à MapInfo génère des recouvrements
En utilisant le prj issu du site EPSG
http://spatialreference.org/ref/epsg/2154/
Code:
ogr2ogr -f "mapinfo file" -a_srs D:\2154.prj prj_esri\sample.TAB sample.shp
J'obtiens les même résultats qu'en haut.
Sachant qu'en haut, j'avais utilisé ogr2ogr avec un en-tête reconnu par mapinfo.
Code:
PROJCS["unnamed", GEOGCS["unnamed", DATUM["GRS_80", SPHEROID["GRS 80",6378137,298.257222101], TOWGS84[0,0,0,0,0,0,0]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433]], PROJECTION["Lambert_Conformal_Conic_2SP"], PARAMETER["standard_parallel_1",44], PARAMETER["standard_parallel_2",49], PARAMETER["latitude_of_origin",46.5], PARAMETER["central_meridian",3], PARAMETER["false_easting",700000], PARAMETER["false_northing",6600000], UNIT["Meter",1.0]]
Autrement dit, c'est le format MapInfo qui semble poser problème, plus particulièrement le driver MITAB..(?)
En utilisant le traducteur universel et la borne "monde", j'obtiens les coordonnées suivantes :
Code:
POLYGON((933459.63435469125397503 1820336.14230716833844781
Normalement, les coordonnées auraient dû être stockées sans décimales, ce qui, au passage, aurait permis d'éviter les recouvrements. Etrange qu'elles apparaissent ici.
Dernière modification par MathieuR (Fri 07 March 2014 11:05)
geodata au cerema et petits billets en géomatique
Hors ligne
Pages: 1
- Sujet précédent - OGR : Une conversion de ESRI à MapInfo génère des recouvrements - Sujet suivant