[Liblas-devel] Vertical Coordinate Systems in LAS format

[Frank Warmerdam]

Guys,

Similar to Mateusz' suggestion I have implemented GetWKT() as:

     enum WKTModeFlag
     {
         eHorizontalOnly = 1,
         eCompoundOK = 2
     };

     /// Returns the OGC WKT describing Spatial Reference System.
     /// If GDAL is linked, it uses GDAL's operations and methods to determine
     /// the WKT.  If GDAL is not linked, no WKT is returned.
     /// \param mode_flag May be eHorizontalOnly indicating the WKT will not
     /// include vertical coordinate system info (the default), or
     /// eCompoundOK indicating the the returned WKT may be a compound
     /// coordinate system if there is vertical coordinate system info
     /// available.
     std::string GetWKT( WKTModeFlag mode_flag = eHorizontalOnly ) const;

We seem to have consensus so I have committed the code implementing
the proposal at:

   http://liblas.org/wiki/VerticalCS

and pushed it up to the master repository.  I have also implemented
python bindings for the new methods, and added some very minimal tests
in the python bindings test suite (SRS-GDAL.txt).

For the time being, I consider the "vertical datum activity" complete,
though I'm still open to feedback.

Howard, I'd appreciate it if you could add an appropriate link to
the VerticalCS wiki topic from whereever you consider appropriate.

Best regards,
-- 
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Programmer for Rent