<div dir="ltr">Even,<div><br></div><div>I'm not sure if this helps, but in dealing with an issue with WKID 5514, this WKT workers to properly convert data from WGS84. This was confirmed against a customer's source data.</div><div><span style="color:rgb(24,54,145);font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace;font-size:12px;line-height:16.8px;white-space:pre-wrap">'</span><span style="color:rgb(24,54,145);font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace;font-size:12px;line-height:16.8px;white-space:pre-wrap">PROJCS["S-JTSK_Krovak_East_North",GEOGCS["GCS_S_JTSK",DATUM["Jednotne_Trigonometricke_Site_Katastralni",SPHEROID["Bessel_1841",6377397.155,299.1528128]],TOWGS84[589,76,480,0,0,0,0],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Krovak"],PARAMETER["False_Easting",0],PARAMETER["False_Northing",0],PARAMETER["Pseudo_Standard_Parallel_1",78.5],PARAMETER["Scale_Factor",0.9999],PARAMETER["Azimuth",30.28813975277778],PARAMETER["Longitude_Of_Center",24.83333333333333],PARAMETER["Latitude_Of_Center",49.5],PARAMETER["X_Scale",-1],PARAMETER["Y_Scale",1],PARAMETER["XY_Plane_Rotation",90],UNIT["Meter",1],AUTHORITY["EPSG","102067"]]</span><span style="color:rgb(24,54,145);font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace;font-size:12px;line-height:16.8px;white-space:pre-wrap">'</span><br></div><div></div><div><br></div><div>Daniel</div><div><br></div><div><i>Daniel Fenton</i></div><div><i>Software Engineer | Esri R&D</i></div><div><i><a href="https://koopjs.github.io">https://koopjs.github.io</a></i><br><div class="gmail_quote"><div dir="ltr">On Thu, Jan 28, 2016 at 12:05 PM Even Rouault <<a href="mailto:even.rouault@spatialys.com" target="_blank">even.rouault@spatialys.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
I'm having a look at the following all related issues :<br>
<a href="https://trac.osgeo.org/postgis/ticket/1940" rel="noreferrer" target="_blank">https://trac.osgeo.org/postgis/ticket/1940</a><br>
<a href="https://trac.osgeo.org/proj/ticket/185" rel="noreferrer" target="_blank">https://trac.osgeo.org/proj/ticket/185</a><br>
<a href="https://trac.osgeo.org/gdal/ticket/4762" rel="noreferrer" target="_blank">https://trac.osgeo.org/gdal/ticket/4762</a><br>
<br>
related to EPSG:2065 support, ie "S-JTSK (Ferro) / Krovak". There's an issue<br>
about the datum shift used, but that's an easy one compared to a more serious<br>
one regarding axis order.<br>
<br>
The proj.4 krovak method expects the first value of a coordinate pair to be a<br>
easting, and the second value a northing. Whereas original Krovak (EPSG:2065)<br>
is supposed to be southing first, westing second.<br>
<br>
<a href="https://trac.osgeo.org/gdal/ticket/4762" rel="noreferrer" target="_blank">https://trac.osgeo.org/gdal/ticket/4762</a> suggests to add a "+to_meter=-1"<br>
parameter. Which means that the coordinates supplied are first westing, second<br>
southing (instead of "+to_meter=-1", a cleaner equivalent way is to use<br>
"+axis=wsu", meaning westing,southing,up )<br>
But EPSG would mandate southing,westing.<br>
<br>
But admitedly, for less exotic coordinate systems with official EPSG axis order<br>
latitude,longitude (EPSG:4326...) or northing,easting (some Finland projected<br>
CRS for example) we have never introduced in GDAL WKT or proj.4 strings the<br>
+axis option, and the application/user does the coordinate switch itself<br>
before passing that to GDAL/proj.<br>
<br>
So, does that sound OK that a user providing coordinates advertized to be in<br>
EPSG:2065 does that in westing, southing order ? Or was this something<br>
particular to this user ? (I'd note that the dataset mentionned in the ticket<br>
is a GML resulting from a WFS server, and we all know the axis ordering issues<br>
related to WFS)<br>
<br>
What a mess...<br>
(Using EPSG:5221 which is the easting, northing variant avoids those<br>
headaches)<br>
<br>
Even<br>
<br>
--<br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a></blockquote></div></div></div>