[Proj] Programmatically obtaining the domain of a projection
pelson.pub at gmail.com
Mon Apr 22 10:23:00 EST 2013
I've got a case where I would like to identify the domain of a projection
programmatically. Take the geos projection for example:
$> invproj +proj=geos +h=35785831 +units=m
$> invproj +proj=geos +h=50000000 +units=m
Are the mathematical limits of Proj4 projections generally available (the C
interface is fine), or is using a numerical method the only way to
calculate the domain?
* There are some domains/bounds available on spatialreference.org (e.g.,
http://spatialreference.org/ref/epsg/4277/) but they are not exhaustive,
are transformed to WGS84 and have presumably been calculated manually.
* I've had a quick look at PJ_geos to see if this particular example has a
simple solution (and indeed grepped for "limits", "bounds", and "range" in
the whole proj.4 source) without much luck.
* I've already got a programmatic solution to this, I'm really looking to
find out if there is a more elegant mathematical answer to my problem.
Many thanks in advance,
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Proj