[OSGeo-Discuss] Parsing well-known text

[Pierre Abbat]

I just added code to PerfectTIN to read the variable-length records of a LAS 
file. I found a WKT record in a square of West Virginia terrain, with no line 
feeds; I added line feeds and indentation and attached it. The other point 
clouds do not have variable-length records. I'd like to heed the WKT units 
when loading a point cloud, if there is a WKT.

I'm also working on a program called Wolkenbase (not public yet) which will 
separate ground from non-ground in a point cloud. If you load a point cloud 
with no unit information, I'd like to add a WKT that indicates only the unit 
and nothing else. What would that look like?

Do you have code to parse and manipulate WKT? The WKT is in a std::string.

Pierre
-- 
The Black Garden on the Mountain is not on the Black Mountain.
-------------- next part --------------
COMPD_CS["NAD83 / UTM zone 17N + NAVD88 height",
  PROJCS["NAD83 / UTM zone 17N",
    GEOGCS["NAD83",
      DATUM["North_American_Datum_1983",
	SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],
	TOWGS84[0,0,0,0,0,0,0],
	AUTHORITY["EPSG","6269"]],
      PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],
      UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],
      AUTHORITY["EPSG","4269"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",-81],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["meter",1,AUTHORITY["EPSG","9001"]],
    AXIS["X",EAST],
    AXIS["Y",NORTH],
    AUTHORITY["EPSG","26917"]],
  VERT_CS["NAVD88 height",
    VERT_DATUM["North American Vertical Datum 1988",2005,AUTHORITY["EPSG","5103"]],
    UNIT["meter",1,AUTHORITY["EPSG","9001"]],
    AXIS["Up",UP],
    AUTHORITY["EPSG","5703"]]]