From: Jan J. <h....@ce...> - 2008-08-29 10:07:31
|
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]]"); |