Martin,
>
>Given that if ESRI just ignore the AXIS elements, then it would be wrong for
>any
>legal WKT having (WEST,SOUTH) orientation without the ESRI-specific extra
>parameters (in which case GeoTools already produces different result than
>ESRI),
>then maybe it is okay to consider that we would get different result in any
>case
>where the AXIS["x", WEST], etc. element are present, no matter if
>PARAMETER["X_Scale", -1]... are present or not?
>
You are almost right - but when using ESRI Krovak without the parameters
like X_SCALE (etc..) it produces CRS with westing southing (!!) directions for
projected axis.
(Sorry - I did not express 100% right in my previous post about this).
In other worlds - this ESRI definition:
PROJCS["S-JTSK_Ferro_Krovak_East_North",
GEOGCS["GCS_S_JTSK_Ferro",
DATUM["D_S_JTSK",
SPHEROID["Bessel_1841",6377397.155,299.1528128]],
PRIMEM["Ferro",-17.66666666666667],
UNIT["Degree",0.0174532925199432955]],
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",42.5],
PARAMETER["Latitude_Of_Center",49.5],
PARAMETER["X_Scale",-1],
PARAMETER["Y_Scale",1],
PARAMETER["XY_Plane_Rotation",90],
UNIT["Meter",1]]
equals to this in GeoTools:
PROJCS["S-JTSK (Ferro) / Krovak",
GEOGCS["S-JTSK (Ferro)",
DATUM["S_JTSK_Ferro",
SPHEROID["Bessel 1841", 6377397.155, 299.1528128,
AUTHORITY["EPSG","7004"]],
AUTHORITY["EPSG","6818"]],
PRIMEM["Ferro", -17.66666666666667, AUTHORITY["EPSG","8909"]],
UNIT["degree", 0.017453292519943295],
AXIS["Longitude", EAST],
AXIS["Latitude", NORTH],
AUTHORITY["EPSG","4818"]],
PROJECTION["Krovak"],
PARAMETER["latitude_of_center", 49.5],
PARAMETER["longitude_of_center", 42.5],
PARAMETER["azimuth", 30.288139722222223],
PARAMETER["pseudo_standard_parallel_1", 78.5],
PARAMETER["scale_factor", 0.9999],
PARAMETER["false_easting", 0.0],
PARAMETER["false_northing", 0.0],
UNIT["m", 1.0],
AXIS["x", EAST],
AXIS["y", NORTH],
AUTHORITY["EPSG","2065"]]
Or eventually this (using X_scale etc...):
PROJCS["S-JTSK (Ferro) / Krovak",
GEOGCS["S-JTSK (Ferro)",
DATUM["S_JTSK_Ferro",
SPHEROID["Bessel 1841", 6377397.155, 299.1528128,
AUTHORITY["EPSG","7004"]],
AUTHORITY["EPSG","6818"]],
PRIMEM["Ferro", -17.66666666666667, AUTHORITY["EPSG","8909"]],
UNIT["degree", 0.017453292519943295],
AXIS["Longitude", EAST],
AXIS["Latitude", NORTH],
AUTHORITY["EPSG","4818"]],
PROJECTION["Krovak"],
PARAMETER["latitude_of_center", 49.5],
PARAMETER["longitude_of_center", 42.5],
PARAMETER["azimuth", 30.288139722222223],
PARAMETER["pseudo_standard_parallel_1", 78.5],
PARAMETER["scale_factor", 0.9999],
PARAMETER["false_easting", 0.0],
PARAMETER["false_northing", 0.0],
PARAMETER["X_Scale",-1.0],
PARAMETER["Y_Scale",1.0],
PARAMETER["XY_Plane_Rotation",90.0]
UNIT["m", 1.0],
AXIS["x", SOUTH],
AXIS["y", WEST],
AUTHORITY["EPSG","2065"]]
All these definitions project
source pt = (12 , 48 )
to
target pt = (-953172.26, -1245573.32)
So the last WKT should work for both systems, but is not the on that is created
by these systems automatically and has to be created by hand. Actually I can't
see other working solution (except of forcing ESRI to use AXIS parameter).
I hope I did not cause you headache :-)
Bests,
Jan.
------------------------------------------------------------------------------
SF.Net email is Sponsored by MIX09, March 18-20, 2009 in Las Vegas, Nevada.
The future of the web can't happen without you. Join us at MIX09 to help
pave the way to the Next Web now. Learn more and register at
http://ad.doubleclick.net/clk;208669438;13503038;i?http://2009.visitmix.com/
_______________________________________________
Geotools-devel mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/geotools-devel
