From: Ben A. <ad...@in...> - 2014-01-10 16:58:34
|
Charles, brilliant! Thank you for these quick and helpful answers! I tried implementing them, here's my problems: Why do the GeographicLib::UTMUPS::Forward/Reverse method not use const for input parameters? That breaks a lot of my code. Of course I could make relevant members mutable or use const_cast, but that seems weird. Now that I have the EPSG code, I need to get the WKT. Grepping for wkt in the GeographicLib sources didn't yield any results, so I guess I need to search and extract the relevant line from EPSG2WKT.TXT (google-hit at geospatialpython) myself? Comparing that file with the WKT from http://spatialreference.org, the text file seems to be missing the AXIS-entries. Also, of course, the height information is missing. Maybe it's out of GeographicLib's scope, but I wonder how exactly the WKT for e.g. UTM 11N with height above the WGS84 ellipsoid in meters should look like. Why does GeographicLib::UTMUPS::EncodeEPSG(11, true) return 32610? Shouldn't it be 32611 (http://spatialreference.org/ref/epsg/32611/)? Thanks again! Ben |