Hello Richard,
When I look at the resulting MIF file in a text editor, the coordinates are indeed correctly transformed but the header is CoordSys Earth Projection 3, 104, "m", 4.36748666666667, 90, 51.1666672333333, 49.8333339, 150000.013, 5400088.438 instead of CoordSys Earth Projection 3, 1019, "m", 4.3674866667, 90, 49.8333339, 51.1666672333, 150000.01300000001, 5400088.4380000001 Bounds (0, 0) (300000, 300000) Thanks for your help ! Jonathan Beliën GEO-6 De : Richard Greenwood [mailto:[email protected]] Envoyé : mardi 9 mai 2017 01:32 À : Jonathan Beliën <[email protected]> Cc : [email protected]; GDAL List <[email protected]>; Hellemans, Philippe <[email protected]> Objet : Re: [gdal-dev] ogr2ogr mapinfo reporjection When you say that the coordsys header is wrong, are you looking at it in MapInfo or in a text editor? If you look at the MIF file coordinates in a text editor do they look like they have been correctly transformed, and it is just the header that is wrong? On Mon, May 8, 2017 at 7:39 AM, Jonathan Beliën <[email protected] <mailto:[email protected]> > wrote: Hello everyone, I struggle to reproject GeoJSON file to TAB (or MIF) file. The source projection is EPSG:4326 and the destination projection is EPSG:31370 (Belge 1972 / Belgian Lambert 72 – Belgium : http://epsg.io/31370) I searched forums and mailing lists and saw that there was some issues with reprojection using Mitab but that it is supposed to be fixed. When I try the following command, the coordinates seems to be transformed but the « coordsys » header in the TAB (or MIF) file is wrong (WGS84 – Non Earth) ogr2ogr test.mif test.json -f "MapInfo File" -dsco FORMAT=MIF -s_srs EPSG:4326 -t_srs EPSG:31370 If I understood correctly what I read on the Internet, that’s because Mitab does not recognise « EPSG:31370 » ; I tried to replace « -t_srs EPSG:31370 » by « -t_srs "31370.txt" » where « 31370.txt » contained the WKT definition of EPSG:31370 ; I saw in the ogr2ogr documentation that there is a « MITAB_BOUNDS_FILE » parameter but did not find out how to use it (http://www.gdal.org/drv_mitab.html) ; I also searched in « /ogr/ogrsf_frmts/mitab/mitab_spatialref.cpp » to see if « Belgian Lambert 1972 » existed and it does : { 0, 1019,"Belgian 1972 7 Parameter", 4, -99.059, 53.322, -112.486, -0.419, 0.83, -1.885, 0.999999, 0}, So I suppose it should be possible to create a MapInfo file projected in « Belgiam Lambert 1972 » but I do not understand how ?! Any help or idea is welcome ! Thanks a lot ! Jonathan Beliën GEO-6 _______________________________________________ gdal-dev mailing list [email protected] <mailto:[email protected]> https://lists.osgeo.org/mailman/listinfo/gdal-dev -- Richard W. Greenwood, PLS www.greenwoodmap.com <http://www.greenwoodmap.com>
_______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
