Hillmuth, Eric a écrit :
CoordinateReferenceSystem geoCRS = CRS.decode("EPSG:4269"); //geographic NAD83
CoordinateReferenceSystem vspCRS = CRS.decode("EPSG:32145"); //Vermont State
Plane NAD83
I'm using 2.1.0. Interestingly, the transform is a bit different than what
ArcView/ArcMap thinks it should be by a few decimal places, I'm not sure why
this would be the case but its not a show stopper.
The projection for "EPSG:32145" is "Transverse Mercator", which is a
somewhat complicated one. Many algorithms exist for this projection with
different rounding error / speed / area of validity tradeoff. There is
at least 2 or 3 different published algorithms if my memory serve me
right (Rueben would know better).
Some other projections like "Cylindrical Mercator" have less variations
in their algorithms.
I think that Geotools implements two "Transverse Mercator" algorithms,
but only one is the default one. Maybe we should provide an API in some
future Geotools version for allowing users to select the default
"Transverse Mercator" algorithm they want to use for instructions like
CRS.decode("EPSG:32145")...
Martin.
-------------------------------------------------------
This SF.net email is sponsored by: Splunk Inc. Do you grep through log files
for problems? Stop! Download the new AJAX search engine that makes
searching your log files as easy as surfing the web. DOWNLOAD SPLUNK!
http://sel.as-us.falkag.net/sel?cmd=lnk&kid3432&bid#0486&dat1642
_______________________________________________
Geotools-gt2-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/geotools-gt2-users
