C++ coordinate class that automates projection?

Question

Can anyone point me to a C++ class for representing coordinates that automatically handles projection transforms? What I'm thinking of is something that would allow me to create vars that incorporate a projection and have appropriate geographic operators such that comparisons make sense even when the coordinates are in different projections. Using template notation (though I can't figure out how to template on a projection)

Coord < SRID:1234> point1 (x1, y1); Coord < SRID:5678> point2 (x2, y2);

if (point1.east(point2)) ... // true if point1 is east of point2

or

point2 = point1 // This operation assigns the coordinates of point1 to point2, 
using the coordinate system of point2 (in this case SRID::5678).

Answer 1

Try to look at Proj.4 library.

You can reproject you points to common coordinate system and then perform the comparison.

Try to create some wrapper class that will perform automate reprojection to common CS and then compare coordinates.

Answer 2

If you want to support comparing points on different datums you have a lot of work ahead; there are often multiple "paths" between datums and it is often the case that the correct path is context dependent.

That being said, CS-MAP (C API) will do its best to pick a path between datums - if you aren't picky it will probably do something acceptable.

And many CS-MAP coordinate systems are tagged with EPSG #s, so you could look up the name associated with a given number.