Suite

Ogr2ogr - Projection perdue lors de la conversion de geojson en onglet MapInfo


J'ai le geojson suivant que je voudrais convertir en onglet MapInfo avec ogr2ogr :ogr2ogr -f "Onglet MapInfo" output.tab input.geojson

geojson :

{ "type":"FeatureCollection", "crs":{ "type":"name", "properties":{ "name":"urn:ogc:def:crs:EPSG:3857" } }, "features" :[ { "type":"Feature", "geometry":{ "type":"LineString", "coordinates":[[-7551132, 6293816],[-7560916, 6260489]], }, "properties": { "id":36 } } ] }

Cependant, l'onglet MapInfo en sortie n'inclut pas la projection définie dans le geojson.

Utilisation d'ogrinfo sur output.tab :ogrinfo -ro -al sortie.tab

Nom de la couche : sortie Géométrie : chaîne de ligne Nombre d'entités : 1 Etendue : (-7560915.990000, 6260488.995000) - (-7551132.000000, 6293815.995000) Couche SRS WKT : PROJCS["unnamed", GEOGCS["unnamed", DATUM["WGS_1984", SPHEROID ["WGS 84",6378137,298.257223563], TOWGS84[0,0,0,0,0,0,0,0]], PRIMEM["Greenwich",0], UNIT["degré",0,0174532925199433]], PROJECTION[ "Mercator_1SP"], PARAMETRE["central_meridian",0], PARAMETER["scale_factor",1], PARAMETER["false_easting",0], PARAMETER["false_northing",0], UNIT["Meter",1.0]] id: Entier (0.0) OGRFeature(output2):1 id (Entier) = 36 Style = PEN(w:1px,c:#000000,id:"mapinfo-pen-2,ogr-pen-0") LINESTRING (- 7551132.0 6293815.995,-7560915.99 6260488.995)

Il semble que ce soit un problème connu avec le format MapInfo Tab (mitab) comme nous pouvons le voir dans cet ancien ticket osgeo ici : https://trac.osgeo.org/gdal/ticket/481

Le format de l'onglet MapInfo peut avoir des problèmes pour reconnaître la projection par son nom.

EDIT: L'utilisation de crs liés au lieu de crs nommés dans le geojson ne fonctionne pas non plus.

EDIT 2 : Ce problème a été corrigé dans la version tronc de gdal depuis le 18/09/2015. Ce patch sera présent dans la prochaine version officielle après cette date.


J'écrirais un courrier à ce sujet à la liste de diffusion gdal-dev car cela peut être un bogue dans le pilote MapInfo. Comparez avec la sortie OGC WKT de gdalsrsinfo :

gdalsrsinfo epsg:3857 PROJ.4 : '+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0= 0 +k=1.0 +units=m [email protected] +wktext +no_defs' OGC WKT : PROJCS["WGS 84 / Pseudo-Mercator", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG"," 7030"]], AUTORITE["EPSG","6326"]], PRIMEM["Greenwich",0, AUTORITE["EPSG","8901"]], UNITÉ["degré",0.0174532925199433, AUTORITE["EPSG" ,"9122"]], AUTORITE["EPSG","4326"]], PROJECTION["Mercator_1SP"], PARAMETRE["central_meridian",0], PARAMETRE["scale_factor",1], PARAMETRE["false_easting", 0], PARAMETER["false_northing",0], UNIT["metre",1, AUTHORITY["EPSG","9001"]], AXIS["X",EAST], AXIS["Y",NORTH], EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 + x_0=0.0 +y_0=0 +k=1.0 +units=m [email protected] +wktext +no_defs"], AUTORITÉ["EPSG","3857"]]


Voir la vidéo: Mapinfo-16: How to create a mapinfo tab file and BTS map from excel data (Octobre 2021).