<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.Shkpostityyli17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 2.0cm 70.85pt 2.0cm;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="FI" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal">Hi,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">This command<o:p></o:p></p>
<p class="MsoNormal">ogr2ogr -f "MapInfo file" -s_srs epsg:3067 -t_srs epsg:2393 2393.mif tm35.jml<o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US">writes into MIF<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">CoordSys Earth Projection 8, 104, "m", 27, 0, 1, 3500000, 0<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">ogrinfo shows<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">ogrinfo 2393.mif -al -so<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">INFO: Open of `2393.mif'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> using driver `MapInfo File' successful.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Layer name: 2393<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Geometry: Unknown (any)<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Feature Count: 3<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Extent: (3284260.027390, 6962722.914217) - (3284610.369420, 6962941.953889)<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Layer SRS WKT:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">PROJCRS["unnamed",<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> BASEGEOGCRS["unnamed",<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> DATUM["World Geodetic System 1984",<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ELLIPSOID["WGS 84",6378137,298.257223563,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> LENGTHUNIT["metre",1,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ID["EPSG",9001]]]],<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> PRIMEM["Greenwich",0,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ANGLEUNIT["degree",0.0174532925199433,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ID["EPSG",9122]]]],<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> CONVERSION["unnamed",<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> METHOD["Transverse Mercator",<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ID["EPSG",9807]],<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> PARAMETER["Latitude of natural origin",0,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ANGLEUNIT["degree",0.0174532925199433],<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ID["EPSG",8801]],<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> PARAMETER["Longitude of natural origin",27,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ANGLEUNIT["degree",0.0174532925199433],<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ID["EPSG",8802]],<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> PARAMETER["Scale factor at natural origin",1,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> SCALEUNIT["unity",1],<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ID["EPSG",8805]],<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> PARAMETER["False easting",3500000,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> LENGTHUNIT["metre",1],<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ID["EPSG",8806]],<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> PARAMETER["False northing",0,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> LENGTHUNIT["metre",1],<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ID["EPSG",8807]]],<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> CS[Cartesian,2],<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> AXIS["easting",east,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ORDER[1],<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> LENGTHUNIT["metre",1,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ID["EPSG",9001]]],<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> AXIS["northing",north,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ORDER[2],<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> LENGTHUNIT["metre",1,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ID["EPSG",9001]]]]<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Data axis to CRS axis mapping: 1,2<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">When I convert this data back into EPSG:3067 the coordinates are shifted by a hundred meters probably because what was called towgs84 parameters are not applied. Parameters appear in the mitab source code
</span><a href="https://github.com/OSGeo/gdal/blob/master/gdal/ogr/ogrsf_frmts/mitab/mitab_spatialref.cpp#L251"><span lang="EN-US">https://github.com/OSGeo/gdal/blob/master/gdal/ogr/ogrsf_frmts/mitab/mitab_spatialref.cpp#L251</span></a><span lang="EN-US"><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">{ 0, 1016,"Finnish_KKJ", 4, -96.062, -82.428, -121.754, -4.801, -0.345, 1.376, 1.496, 0},<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">I tried to fix the MIF by changing the CoordSys line into<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">CoordSys Earth Projection 24, 1016, "m", 27, 0, 1, 3500000, 0<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">that is to my knowledge what MapInfo itself is using but that does not work at all with GDAL<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">ogr2ogr -f "jml" -t_srs epsg:3067 2393_alt_from_mif.jml 2393_alt.mif<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">ERROR 1: PROJ: proj_as_proj_string: Unsupported conversion method: Transverse_Mercator_MapInfo_24<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">ERROR 6: Cannot find coordinate operations from `PROJCRS["unnamed",BASEGEOGCRS["unnamed",DATUM["Finnish_KKJ",ELLIPSOID["International 1924",6378388,297,LENGTHUNIT["metre",1,ID["EPSG",9001]]]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]]],CONVERSION["unnamed",METHOD["Transverse_Mercator_MapInfo_24"],PARAMETER["latitude_of_origin",0,ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]],PARAMETER["central_meridian",27,ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]],PARAMETER["scale_factor",1,SCALEUNIT["unity",1,ID["EPSG",9201]]],PARAMETER["false_easting",3500000,LENGTHUNIT["metre",1,ID["EPSG",9001]]],PARAMETER["false_northing",0,LENGTHUNIT["metre",1,ID["EPSG",9001]]]],CS[Cartesian,2],AXIS["easting",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["northing",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]'
to `EPSG:3067'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">ERROR 1: Failed to create coordinate transformation between the<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">following coordinate systems. This may be because they<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">are not transformable.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Is there something that I could do for making the conversion to round-trip?<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">-Jukka Rahkonen-<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
</div>
</body>
</html>