[OSGeo-Discuss] Parsing well-known text
Pierre Abbat
phma at bezitopo.org
Wed Sep 2 01:14:00 PDT 2020
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"]]]
More information about the Discuss
mailing list