Max konig: I came across a bug in the Sentinel Toolbox, which has been present in Nest for a long time, not allowing to use many EPSG-projections. But I also see where the problem lies: For many projections I get this error:
cid:image001.jpg@01CF8ACE.9C6710D0
I compared a map projection that works with one where “Axis too complex” error is reported. Map projections that work contain these lines:
AXIS["Easting", EAST], AXIS["Northing", NORTH],
While the non-working ones contain these lines, which apparently are wrong or unprocessable by NEST:
AXIS["Easting", "North along 90 deg East"], AXIS["Northing", "North along 0 deg"],
If I manually replace these lines in the graph file (xml file) for the projection definition, then it works!
So the task for the next NEST version would be to correct all map projections in NEST for this issue.
Here one working and the one not: EPSG3033 works: PROJCS["WGS 84 / Australian Antarctic Lambert", GEOGCS["WGS 84", DATUM["World Geodetic System 1984", SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], UNIT["degree", 0.017453292519943295], AXIS["Geodetic longitude", EAST], AXIS["Geodetic latitude", NORTH], AUTHORITY["EPSG","4326"]], PROJECTION["Lambert_Conformal_Conic_2SP", AUTHORITY["EPSG","9802"]], PARAMETER["central_meridian", 70.0], PARAMETER["latitude_of_origin", -50.0], PARAMETER["standard_parallel_1", -68.5], PARAMETER["false_easting", 6000000.0], PARAMETER["false_northing", 6000000.0], PARAMETER["scale_factor", 1.0], PARAMETER["standard_parallel_2", -74.5], UNIT["m", 1.0], AXIS["Easting", EAST], AXIS["Northing", NORTH], AUTHORITY["EPSG","3033"]]
Now, if the selected CRS is an instance of AbstractReferenceSystem then the none strict mode is used to format the WKT.
Luis Veci 21 January 2015 at 15:57
The problem seems to be in formatting the Axis AXIS["Easting", "North along 90 deg East"], AXIS["Northing", "North along 0 deg"],
with toWkt(), the strict parameter is used which will throw if unknown text exists If it falls back to use toString() then the resulting projection appears fine.
This will need to be fixed in geotools. While we are still using an older version of geotools, we could fallback to toString() when needed
Max konig: I came across a bug in the Sentinel Toolbox, which has been present in Nest for a long time, not allowing to use many EPSG-projections. But I also see where the problem lies:
For many projections I get this error:
cid:image001.jpg@01CF8ACE.9C6710D0
I compared a map projection that works with one where “Axis too complex” error is reported. Map projections that work contain these lines:
AXIS["Easting", EAST],
AXIS["Northing", NORTH],
While the non-working ones contain these lines, which apparently are wrong or unprocessable by NEST:
AXIS["Easting", "North along 90 deg East"],
AXIS["Northing", "North along 0 deg"],
If I manually replace these lines in the graph file (xml file) for the projection definition, then it works!
So the task for the next NEST version would be to correct all map projections in NEST for this issue.
Here one working and the one not:
EPSG3033 works:
PROJCS["WGS 84 / Australian Antarctic Lambert",
GEOGCS["WGS 84",
DATUM["World Geodetic System 1984",
SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
UNIT["degree", 0.017453292519943295],
AXIS["Geodetic longitude", EAST],
AXIS["Geodetic latitude", NORTH],
AUTHORITY["EPSG","4326"]],
PROJECTION["Lambert_Conformal_Conic_2SP", AUTHORITY["EPSG","9802"]],
PARAMETER["central_meridian", 70.0],
PARAMETER["latitude_of_origin", -50.0],
PARAMETER["standard_parallel_1", -68.5],
PARAMETER["false_easting", 6000000.0],
PARAMETER["false_northing", 6000000.0],
PARAMETER["scale_factor", 1.0],
PARAMETER["standard_parallel_2", -74.5],
UNIT["m", 1.0],
AXIS["Easting", EAST],
AXIS["Northing", NORTH],
AUTHORITY["EPSG","3033"]]
EPSG3031 works not:
PROJCS["WGS 84 / Antarctic Polar Stereographic",
GEOGCS["WGS 84",
DATUM["World Geodetic System 1984",
SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
UNIT["degree", 0.017453292519943295],
AXIS["Geodetic longitude", EAST],
AXIS["Geodetic latitude", NORTH],
AUTHORITY["EPSG","4326"]],
PROJECTION["Polar Stereographic (variant B)", AUTHORITY["EPSG","9829"]],
PARAMETER["central_meridian", 0.0],
PARAMETER["Standard_Parallel_1", -71.0],
PARAMETER["false_easting", 0.0],
PARAMETER["false_northing", 0.0],
UNIT["m", 1.0],
AXIS["Easting", "North along 90 deg East"],
AXIS["Northing", "North along 0 deg"],
AUTHORITY["EPSG","3031"]]