Coordinate System Transformation Error?
---------------------------------------
Key: GEOT-1425
URL: http://jira.codehaus.org/browse/GEOT-1425
Project: GeoTools
Issue Type: Bug
Components: core referencing
Affects Versions: 2.3.3
Environment: Windows XP Pro, JAI 1.1.3, JDK 1.5.12
Reporter: Rafal Glowinski
Assignee: Martin Desruisseaux
Attachments: test.rar
Hi,
I have a problem with transformation between two coordinate systems. Both are
national coordinate systems in Poland. They are defined by PRJ files provided
by ESRI. Their definitions are as follows (WKT):
PROJCS["UWPP_1992",
GEOGCS["GCS_ETRF_1989",
DATUM["D_ETRF_1989",
SPHEROID["WGS_1984",6378137.0,298.257223563]],
PRIMEM["Greenwich",0.0],
UNIT["Degree",0.0174532925199433],
AXIS["Longitude",EAST],
AXIS["Latitude",NORTH]],
PROJECTION["Gauss-Kruger"],
PARAMETER["False_Easting",500000.0],
PARAMETER["False_Northing",-5300000.0],
PARAMETER["Central_Meridian",19.0],
PARAMETER["Scale_Factor",0.9993],
PARAMETER["Latitude_Of_Origin",0.0],
UNIT["Metre",1.0],
AXIS["x",EAST],
AXIS["y",NORTH]]
PROJCS["UWPP_2000_pas_5",
GEOGCS["GCS_ETRF_1989",
DATUM["D_ETRF_1989",
SPHEROID["WGS_1984",6378137.0,298.257223563]],
PRIMEM["Greenwich",0.0],
UNIT["Degree",0.0174532925199433],
AXIS["Longitude",EAST],
AXIS["Latitude",NORTH]],
PROJECTION["Gauss-Kruger"],
PARAMETER["False_Easting",5500000.0],
PARAMETER["False_Northing",0.0],
PARAMETER["Central_Meridian",15.0],
PARAMETER["Scale_Factor",0.999923],
PARAMETER["Latitude_Of_Origin",0.0],
UNIT["Metre",1.0],
AXIS["x",EAST],
AXIS["y",NORTH]]
Now let us consider a point from EUREF-POL, which number is 304
(http://www.cbk.waw.pl/zgp/proj/eurefpol_pl.html#Tab). Its Coordinates in BLH
form are as follows:
B: 53° 42' 14.42353"
L: 16° 07' 46.50919"
H: 238.935 m
Its coordinates in UWPP_1992 system are (652629.45871, 310593.99113)
Its coordinates in UWPP_2000_pas_5 system are (5953105.55795, 5574590.75485)
However when I try to transform a point from UWPP_1992 to UWPP_2000_pas_5 I get
weird results: (5935488.945, 5629996.686)
I also think that even the first phase of transformation (Inverse Mercator) to
decimal degrees yields wrong results. Your code gives the following coordinates
(in decimal degrees):
B: 50.6423
L: 21.1593
while coordinates taken from the EUREF-POL website are:
B: 53.0703
L: 16.1294
Attached is a JAVA source I used to reproduce the behaviour.
--
This message is automatically generated by JIRA.
-
If you think it was sent incorrectly contact one of the administrators:
http://jira.codehaus.org/secure/Administrators.jspa
-
For more information on JIRA, see: http://www.atlassian.com/software/jira
-------------------------------------------------------------------------
This SF.net email is sponsored by: Splunk Inc.
Still grepping through log files to find problems? Stop.
Now Search log events and configuration files using AJAX and a browser.
Download your FREE copy of Splunk now >> http://get.splunk.com/
_______________________________________________
Geotools-devel mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/geotools-devel
