extract WKT (Well known text string) from laz vlrs, if available
store it in the header as WKT string
change encoding of the metadata json to newtonsoft to better format the string.
Note:
With the previous serializing of the metatdata json, the WKT string would end up with many "\u0022" characters written in it.
I seem to have corrected that out by importing and using newtonsoft json for serializing but perhaps there is a better way?
Here is an example of a WKT that may or may not be in a laz file, and may or may not be stored in the vlr.RecordID 2112
I can email a sample laz file containing the data for testing.
COMPD_CS["NAD83(2011) / UTM zone 15N / Geoid2018",PROJCS["NAD83(2011) / UTM zone 15N / Geoid2018",GEOGCS["NAD83(2011) / UTM zone 15N / Geoid2018",DATUM["NAD83(2011)",SPHEROID["GRS 1980",6378137.000,298.257222101,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","1116"]],PRIMEM["Greenwich",0.0000000000000000,AUTHORITY["EPSG","8901"]],UNIT["Degree",0.01745329251994329547,AUTHORITY["EPSG","9102"]],AUTHORITY["EPSG","6344"]],PROJECTION["Transverse_Mercator",AUTHORITY["EPSG","9807"]],PARAMETER["latitude_of_origin",0.0000000000000000],PARAMETER["central_meridian",-93.0000000000000142],PARAMETER["scale_factor",0.9996000000000000],PARAMETER["false_easting",500000.000],PARAMETER["false_northing",0.000],UNIT["Meter",1.00000000000000000000,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","6344"]],VERT_CS["Orthometric Heights",VERT_DATUM["Geoid2018",2005,AUTHORITY["EPSG","0"]],UNIT["Meter",1.00000000000000000000,AUTHORITY["EPSG","9001"]],AXIS["Height",UP]]]
unitycoder
commented
5 months ago
Note: With the previous serializing of the metatdata json, the WKT string would end up with many "\u0022" characters written in it. I seem to have corrected that out by importing and using newtonsoft json for serializing but perhaps there is a better way?
Here is an example of a WKT that may or may not be in a laz file, and may or may not be stored in the vlr.RecordID 2112 I can email a sample laz file containing the data for testing.
COMPD_CS["NAD83(2011) / UTM zone 15N / Geoid2018",PROJCS["NAD83(2011) / UTM zone 15N / Geoid2018",GEOGCS["NAD83(2011) / UTM zone 15N / Geoid2018",DATUM["NAD83(2011)",SPHEROID["GRS 1980",6378137.000,298.257222101,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","1116"]],PRIMEM["Greenwich",0.0000000000000000,AUTHORITY["EPSG","8901"]],UNIT["Degree",0.01745329251994329547,AUTHORITY["EPSG","9102"]],AUTHORITY["EPSG","6344"]],PROJECTION["Transverse_Mercator",AUTHORITY["EPSG","9807"]],PARAMETER["latitude_of_origin",0.0000000000000000],PARAMETER["central_meridian",-93.0000000000000142],PARAMETER["scale_factor",0.9996000000000000],PARAMETER["false_easting",500000.000],PARAMETER["false_northing",0.000],UNIT["Meter",1.00000000000000000000,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","6344"]],VERT_CS["Orthometric Heights",VERT_DATUM["Geoid2018",2005,AUTHORITY["EPSG","0"]],UNIT["Meter",1.00000000000000000000,AUTHORITY["EPSG","9001"]],AXIS["Height",UP]]]