[Liblas-devel] can't detect projection

[Peter Tittmann]

All,

Working on Linux <> 4.9.0-2-amd64 #1 SMP Debian 4.9.10-1 (2017-02-17) x86_64 GNU/Linux

with 

liblas-config --version
1.8.1

lasinfo can’t detect the spatial reference of one of my .las tiles

> ---------------------------------------------------------
>  Header Summary
> ---------------------------------------------------------
> 
>  Version:                     1.2
>  Source ID:                   0
>  Reserved:                    1
>  Project ID/GUID:             '00000000-0000-0000-4b65-6e6169000000'
>  System ID:                   'Quantum Spatial'
>  Generating Software:         'LasMonkey 2.2.3'
>  File Creation Day/Year:      351/2016
>  Header Byte Size             227
>  Data Offset:                 329
>  Header Padding:              0
>  Number Var. Length Records:  1
>  Point Data Format:           1
>  Number of Point Records:     1222629
>  Compressed:                  False
>  Number of Points by Return:  1200340 20588 1548 141 11 
>  Scale Factor X Y Z:          0.01000000000000 0.01000000000000 0.01000000000000
>  Offset X Y Z:                -0.00 -0.00 -0.00
>  Min X Y Z:                   1368090.03 1867637.48 -0.69
>  Max X Y Z:                   1370722.43 1868999.99 425.33
>  Spatial Reference:           
> LOCAL_CS["unnamed",
>    GEOGCS["unknown",
>        DATUM["unknown",
>            SPHEROID["unretrievable - using WGS84",6378137,298.257223563]],
>        PRIMEM["Greenwich",0],
>        UNIT[,0.0174532925199433]],
>    AUTHORITY["EPSG","6397"],
>    UNIT["US survey foot",0.3048006096012192]]

However when I look at the same file on OS X : Darwin <> 15.6.0 Darwin Kernel Version 15.6.0: Mon Jan  9 23:07:29 PST 2017; :xnu-3248.60.11.2.1~1/RELEASE_X86_64 x86_64

using

liblas-config --version
1.8.0

and voila:

> ---------------------------------------------------------
>  Header Summary
> ---------------------------------------------------------
> 
>  Version:                     1.2
>  Source ID:                   0
>  Reserved:                    1
>  Project ID/GUID:             '00000000-0000-0000-4b65-6e6169000000'
>  System ID:                   'Quantum Spatial'
>  Generating Software:         'LasMonkey 2.2.3'
>  File Creation Day/Year:      351/2016
>  Header Byte Size             227
>  Data Offset:                 329
>  Header Padding:              0
>  Number Var. Length Records:  1
>  Point Data Format:           1
>  Number of Point Records:     2862909
>  Compressed:                  False
>  Number of Points by Return:  1419410 962642 389876 81545 8863 
>  Scale Factor X Y Z:          0.01000000000000 0.01000000000000 0.01000000000000
>  Offset X Y Z:                -0.00 -0.00 -0.00
>  Min X Y Z:                   1330298.86 1870307.91 9.78
>  Max X Y Z:                   1331999.99 1871999.99 657.02
>  Spatial Reference:           
> PROJCS["NAD83(2011) / Alaska zone 4",
>    GEOGCS["NAD83(2011)",
>        DATUM["NAD83_National_Spatial_Reference_System_2011",
>            SPHEROID["GRS 1980",6378137,298.2572221010002,
>                AUTHORITY["EPSG","7019"]],
>            AUTHORITY["EPSG","1116"]],
>        PRIMEM["Greenwich",0],
>        UNIT["degree",0.0174532925199433],
>        AUTHORITY["EPSG","6318"]],
>    PROJECTION["Transverse_Mercator"],
>    PARAMETER["latitude_of_origin",54],
>    PARAMETER["central_meridian",-150],
>    PARAMETER["scale_factor",0.9999],
>    PARAMETER["false_easting",1640416.666666667],
>    PARAMETER["false_northing",0],
>    UNIT["US survey foot",0.3048006096012192,
>        AUTHORITY["EPSG","9003"]],
>    AUTHORITY["EPSG","6397"]]
> 


Of note, the EPSG:6397 reference is actually wrong as EPSG:6397 is in metes not feet, I’m not sure if that screwing things up. 

Any idea what I can do to get libels on Debian to read the projection? Would changing/removing the erroneous EPSG authority make a difference?

Thanks!

Peter