[Liblas-devel] can't detect projection

Previous Topic Next Topic
 
classic Classic list List threaded Threaded
4 messages Options
Reply | Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

[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
_______________________________________________
Liblas-devel mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/liblas-devel
Reply | Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

Re: [Liblas-devel] can't detect projection

Mateusz Loskot
On 27 March 2017 at 20:27, Peter Tittmann <[hidden email]> wrote:

> 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
> [...]
>>  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]]

Does setting GDAL_DATA env make a difference?

https://trac.osgeo.org/gdal/wiki/FAQInstallationAndBuilding#WhatisGDAL_DATAenvironmentvariable

Best regards,
--
Mateusz Loskot, http://mateusz.loskot.net
_______________________________________________
Liblas-devel mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/liblas-devel
Reply | Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

Re: [Liblas-devel] can't detect projection

Peter Tittmann
Thanks Mateusz,

The following describes the current GDAL_DATA environment on the Debian machine (which can’t detect the CRS). Would there be another source for the data directory I should set?

$ gdal-config --datadir
/usr/share/gdal/2.1

and

$ ls -lrt /usr/share/gdal/2.1/
total 2704
-rw-r--r-- 1 root root   1081 Feb  2 14:15 vertcs.override.csv
-rw-r--r-- 1 root root  11931 Feb  2 14:15 vertcs.csv
-rw-r--r-- 1 root root   2854 Feb  2 14:15 vdv452.xsd
-rw-r--r-- 1 root root  25816 Feb  2 14:15 vdv452.xml
-rw-r--r-- 1 root root  21480 Feb  2 14:15 unit_of_measure.csv
-rw-r--r-- 1 root root  12071 Feb  2 14:15 trailer.dxf
-rw-r--r-- 1 root root  10360 Feb  2 14:15 stateplane.csv
-rw-r--r-- 1 root root   2048 Feb  2 14:15 seed_3d.dgn
-rw-r--r-- 1 root root   9216 Feb  2 14:15 seed_2d.dgn
-rw-r--r-- 1 root root  37030 Feb  2 14:15 s57objectclasses_iw.csv
-rw-r--r-- 1 root root  31226 Feb  2 14:15 s57objectclasses.csv
-rw-r--r-- 1 root root  63984 Feb  2 14:15 s57objectclasses_aml.csv
-rw-r--r-- 1 root root  20885 Feb  2 14:15 s57expectedinput.csv
-rw-r--r-- 1 root root   8745 Feb  2 14:15 s57attributes_iw.csv
-rw-r--r-- 1 root root   7254 Feb  2 14:15 s57attributes.csv
-rw-r--r-- 1 root root  18115 Feb  2 14:15 s57attributes_aml.csv
-rw-r--r-- 1 root root  13304 Feb  2 14:15 s57agencies.csv
-rw-r--r-- 1 root root  64880 Feb  2 14:15 ruian_vf_v1.gfs
-rw-r--r-- 1 root root  45972 Feb  2 14:15 ruian_vf_st_v1.gfs
-rw-r--r-- 1 root root   2600 Feb  2 14:15 ruian_vf_st_uvoh_v1.gfs
-rw-r--r-- 1 root root  46775 Feb  2 14:15 ruian_vf_ob_v1.gfs
-rw-r--r-- 1 root root 272735 Feb  2 14:15 projop_wparm.csv
-rw-r--r-- 1 root root   1685 Feb  2 14:15 prime_meridian.csv
-rw-r--r-- 1 root root   1525 Feb  2 14:15 pcs.override.csv
-rw-r--r-- 1 root root 719271 Feb  2 14:15 pcs.csv
-rw-r--r-- 1 root root   3466 Feb  2 14:15 pci_ellips.txt
-rw-r--r-- 1 root root  35105 Feb  2 14:15 pci_datum.txt
-rw-r--r-- 1 root root   1349 Feb  2 14:15 ozi_ellips.csv
-rw-r--r-- 1 root root   8482 Feb  2 14:15 ozi_datum.csv
-rw-r--r-- 1 root root   4659 Feb  2 14:15 osmconf.ini
-rw-r--r-- 1 root root  25079 Feb  2 14:15 ogrvrt.xsd
-rw-r--r-- 1 root root   6144 Feb  2 14:15 nitf_spec.xsd
-rw-r--r-- 1 root root  91275 Feb  2 14:15 nitf_spec.xml
-rw-r--r-- 1 root root   7490 Feb  2 14:15 netcdf_config.xsd
-rw-r--r-- 1 root root   4812 Feb  2 14:15 inspire_cp_CadastralZoning.gfs
-rw-r--r-- 1 root root   2450 Feb  2 14:15 inspire_cp_CadastralParcel.gfs
-rw-r--r-- 1 root root   1650 Feb  2 14:15 inspire_cp_CadastralBoundary.gfs
-rw-r--r-- 1 root root   1740 Feb  2 14:15 inspire_cp_BasicPropertyUnit.gfs
-rw-r--r-- 1 root root   9120 Feb  2 14:15 header.dxf
-rw-r--r-- 1 root root   1719 Feb  2 14:15 gt_ellips.csv
-rw-r--r-- 1 root root  15804 Feb  2 14:15 gt_datum.csv
-rw-r--r-- 1 root root   3670 Feb  2 14:15 gml_registry.xml
-rw-r--r-- 1 root root  12761 Feb  2 14:15 geoccs.csv
-rw-r--r-- 1 root root  19428 Feb  2 14:15 gdalvrt.xsd
-rw-r--r-- 1 root root  12305 Feb  2 14:15 GDALLogoGS.svg
-rw-r--r-- 1 root root  12305 Feb  2 14:15 GDALLogoColor.svg
-rw-r--r-- 1 root root  13022 Feb  2 14:15 GDALLogoBW.svg
-rw-r--r-- 1 root root   2021 Feb  2 14:15 gdalicon.png
-rw-r--r-- 1 root root 216876 Feb  2 14:15 gdal_datum.csv
-rw-r--r-- 1 root root    513 Feb  2 14:15 gcs.override.csv
-rw-r--r-- 1 root root  52729 Feb  2 14:15 gcs.csv
-rw-r--r-- 1 root root  75845 Feb  2 14:15 esri_Wisconsin_extra.wkt
-rw-r--r-- 1 root root 332546 Feb  2 14:15 esri_StatePlane_extra.wkt
-rw-r--r-- 1 root root     50 Feb  2 14:15 epsg.wkt
-rw-r--r-- 1 root root  13924 Feb  2 14:15 ellipsoid.csv
-rw-r--r-- 1 root root 177325 Feb  2 14:15 datum_shift.csv
-rw-r--r-- 1 root root   7978 Feb  2 14:15 coordinate_axis.csv
-rw-r--r-- 1 root root   8918 Feb  2 14:15 compdcs.csv


> On Mar 27, 2017, at 11:41 AM, Mateusz Loskot <[hidden email]> wrote:
>
> On 27 March 2017 at 20:27, Peter Tittmann <[hidden email]> wrote:
>> 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
>> [...]
>>> 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]]
>
> Does setting GDAL_DATA env make a difference?
>
> https://trac.osgeo.org/gdal/wiki/FAQInstallationAndBuilding#WhatisGDAL_DATAenvironmentvariable
>
> Best regards,
> --
> Mateusz Loskot, http://mateusz.loskot.net

_______________________________________________
Liblas-devel mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/liblas-devel
Reply | Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

Re: [Liblas-devel] can't detect projection

Mateusz Loskot
On 27 March 2017 at 20:59, Peter Tittmann <[hidden email]> wrote:
> Thanks Mateusz,
>
> The following describes the current GDAL_DATA environment on the Debian machine (which can’t detect the CRS). Would there be another source for the data directory I should set?

GDAL_DATA points to single location, unless something has changed lately.

The output you are getting is generated by these lines, so it's either
GeoTIFF or GDAL:
https://github.com/libLAS/libLAS/blob/master/src/header.cpp#L886-L895

No idea how debian packages are generated.
You may try to copy over GDAL_DATA files to Debian and see if it helps.
Perhaps the shared SRS files changed somehow.
Having no runtime env handy, I can only keep brainstorming.

Best regards,
--
Mateusz Loskot, http://mateusz.loskot.net
_______________________________________________
Liblas-devel mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/liblas-devel
Loading...