[geos-devel] Does geos handle 3d points

[David Blasby]

The output of the included program is "geom = POINT (0.0000000000000000 
1.0000000000000000)".  Is the Z ordinate just not being printed or is it 
lost?  

If its not being printed, then we should fix this.  Also, the 
WKTReader() should handle 3d points as well.


// g++ test2.cpp -I/usr/local/include/geos -lgeos

#include "stdio.h"
#include "io.h"


int main(int argC, char* argV[])
{

    Coordinate *c;
    GeometryFactory *geomFactory;
    Geometry *g ;

    geomFactory= new GeometryFactory( new PrecisionModel(), -1);
    c = new Coordinate(0,1,2);
    g= geomFactory->createPoint(*c);

    string s= g->toString();
    cout << "geom = "<<s;
}