Hi Martin,
>
>Just to make sure that I get the big picture right, could you confirm me
>that in the current state, invoking
>
> ProjectedCRS.getBaseToCRS().getMathTransform().toWKT()
>
>on the CRS you created gives you a CONCAT_MT with 2 transforms, the first
>one containing a few elements like PARAMETER["elt_0_0", 0] and the second
>one containing the map projection parameters without the
>PARAMETER["X_Scale", ...] parameters?
>
yes -
System.out.println(((ProjectedCRS)CRS.parseWKT(prj)).getConversionFromBase().getMathTransform());
gives:
CONCAT_MT[PARAM_MT["Krovak",
PARAMETER["semi_major", 6377397.155],
PARAMETER["semi_minor", 6356078.962818189],
PARAMETER["latitude_of_center", 49.5],
PARAMETER["longitude_of_center", 24.833333333333332],
PARAMETER["azimuth", 30.28813975277778],
PARAMETER["pseudo_standard_parallel_1", 78.5],
PARAMETER["scale_factor", 0.9999],
PARAMETER["false_easting", 0.0],
PARAMETER["false_northing", 0.0]],
PARAM_MT["Affine",
PARAMETER["num_row", 3],
PARAMETER["num_col", 3],
PARAMETER["elt_0_0", -0.00000000000000006123],
PARAMETER["elt_0_1", 1.0],
PARAMETER["elt_1_0", -1.0],
PARAMETER["elt_1_1", -0.00000000000000006123]]]
But CRS.parseWKT(prj) doesn't give the affine transform but just this:
PROJCS["S-JTSK_Krovak_East_North",
GEOGCS["GCS_S_JTSK",
DATUM["D_S_JTSK",
SPHEROID["Bessel_1841", 6377397.155, 299.1528128]],
PRIMEM["Greenwich", 0.0],
UNIT["degree", 0.017453292519943295],
AXIS["Longitude", EAST],
AXIS["Latitude", NORTH]],
PROJECTION["Krovak"],
PARAMETER["latitude_of_center", 49.5],
PARAMETER["longitude_of_center", 24.833333333333332],
PARAMETER["azimuth", 30.28813975277778],
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]]
The original wkt is:
String prj = new String(
"PROJCS[\"S-JTSK_Krovak_East_North\", \n" +
"GEOGCS[\"GCS_S_JTSK\",DATUM[\"D_S_JTSK\",\n" +
"SPHEROID[\"Bessel_1841\",6377397.155,299.1528128]],\n" +
"PRIMEM[\"Greenwich\",0.0],\n"+
"UNIT[\"Degree\",0.0174532925199433]],\n"+
"PROJECTION[\"Krovak\"],\n" +
"PARAMETER[\"Pseudo_Standard_Parallel_1\",78.5],\n" +
"PARAMETER[\"Scale_Factor\",0.9999],\n" +
"PARAMETER[\"Azimuth\",30.28813975277778],\n" +
"PARAMETER[\"Longitude_Of_Center\",24.83333333333333],\n" +
"PARAMETER[\"Latitude_Of_Center\",49.5],\n"+
"PARAMETER[\"X_Scale\",-1.0],\n"+
"PARAMETER[\"Y_Scale\",-1.0],\n" +
"PARAMETER[\"XY_Plane_Rotation\",90.0],\n" +
"AXIS[\"x\", EAST], \n" +
"AXIS[\"y\", NORTH], \n" +
"UNIT[\"Meter\",1.0]]");
-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Geotools-devel mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/geotools-devel
