Daugiau

Ogr2ogr - prarasta projekcija konvertuojant geojsoną į „MapInfo“ skirtuką


Turiu šį geojsoną, kurį norėčiau konvertuoti į „MapInfo“ skirtuką su ogr2ogr:ogr2ogr -f „MapInfo skirtukas“ output.tab input.geojson

geojsonas:

{"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}}]}

Tačiau išvestyje „MapInfo Tab“ nėra geojsone apibrėžtos projekcijos.

Naudojant ogrinfo output.tab:ogrinfo -ro -al output.tab

Sluoksnio pavadinimas: išvesties geometrija: eilutės eilutės funkcijų skaičius: 1 apimtis: (-7560915.990000, 6260488.995000) - (-7551132.000000, 6293815.995000) Sluoksnis SRS WKT: PROJCS ["nenurodytas", GEOGCS ["neįvardytas", DATUM [[WGS_ERO] ["WGS 84", 6378137,298.257223563], TOWGS84 [0,0,0,0,0,0,0]], PRIMEM ["Greenwich", 0], UNIT ["laipsnis", 0.0174532925199433]], PROJEKCIJA [ "Mercator_1SP"], PARAMETER ["central_meridian", 0], PARAMETER ["scale_factor", 1], PARAMETER ["false_easting", 0], PARAMETER ["false_northing", 0], UNIT ["Meter", 1.0]] id: sveikasis skaičius (0.0) OGRFeature (išvestis2): 1 id (sveikas skaičius) = 36 stilius = PEN (w: 1px, c: # 000000, id: "mapinfo-pen-2, ogr-pen-0") LINESTRING (- 7551132.0 6293815.995, -7560915.99 6260488.995)

Atrodo, kad tai yra žinoma „MapInfo Tab“ (mitab) formato problema, kaip matome šiame sename „osgeo“ biliete: https://trac.osgeo.org/gdal/ticket/481

„MapInfo Tab“ formatu gali kilti problemų atpažinti projekciją pagal pavadinimą.

REDAGUOTI: Geojsone vietoj pavadintų crs naudoti susietus crs taip pat neveikia.

2 REDAGAVIMAS: ši problema buvo ištaisyta „gdal“ bagažinės versijoje nuo 2015-09-18. Šis pleistras bus pateiktas kitame oficialiame leidime po šios datos.


Apie tai rašyčiau laiškus į „gdal-dev“ adresų sąrašą, nes tai gali būti „MapInfo“ tvarkyklės klaida. Palyginkite su „OGC WKT“ išvestimi iš „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 + vienetai = m + nadgrids = @ null + 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 "]], AUTHORITY [" EPSG "," 6326 "]], PRIMEM [" Greenwich ", 0, AUTHORITY [" EPSG "," 8901 "]], UNIT [" laipsnis ", 0.0174532925199433, AUTHORITY [" EPSG " , "9122"]], AUTHORITY ["EPSG", "4326"]], PROJECTION ["Mercator_1SP"], PARAMETER ["central_meridian", 0], PARAMETER ["scale_factor", 1], PARAMETER ["false_easting", 0], PARAMETRAS ["false_northing", 0], UNIT ["skaitiklis", 1, AUTHORITY ["EPSG", "9001"]], AXIS ["X", RYTAS], AXIS ["Y", NORTH], IŠPLĖTIMAS ["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 + vienetai = m + nadgrids = @ null + wktext + no_defs "], AUTHORITY [" EPSG "," 3857 "]]


Žiūrėti video įrašą: 2012-04-14 TV laida Gimtoji žemė pirma dalis (Spalio Mėn 2021).