Hello,
indeed. Thanks to Jeff for connecting me to the latest *.csv files. They will be part of the next release. Here they are until then
https://www.dropbox.com/s/magapow18hakmlr/LAStools_bin_serf_geo.zipSimply replace the three files in the ./LAStools/bin/serf/geo folder with these ones. Then the EPSG codes are recognized,
E:\LAStools\bin>las2las -i 476_5026_201901.laz -epsg 2961 -target_epsg 8083 -odix _8083 -olaz
SERIOUS WARNING: horizontal datum of source and target incompatible.
WARNING: reprojecting file with lastiling VLR invalidates this VLR
But there are two issues.
(1) The horizontal datum of source and target are incompatible. Do you need a more complex transformation? Do you need to use a grid of shift vectors? Your use case and your national geospatial agency (or other professionals) should have an answer for you here.
(2) The pre-existing LAStiling VLR in this LAZ file (added by lastile) no longer makes sense. It should be removed with option '-remove_tiling_vlr'.
Below the resulting lasinfo reports with and without option '-remove_tiling_vlr'.
E:\LAStools\bin>lasinfo -i 476_5026_201901_8083.laz
lasinfo (201111) report for '476_5026_201901_8083.laz'
reporting all LAS header entries:
file signature: 'LASF'
file source ID: 0
global_encoding: 17
project ID GUID data 1-4: 00000000-0000-0000-0000-000000000000
version major.minor: 1.4
system identifier: 'LAStools (c) by rapidlasso GmbH'
generating software: 'las2las (version 201030)'
file creation day/year: 15/2019
header size: 375
offset to point data: 1205
number var. length records: 2
point data format: 6
point data record length: 30
number of point records: 0
number of points by return: 0 0 0 0 0
scale factor x y z: 0.01 0.01 0.01
offset x y z: 25500000 5000000 0
min x y z: 25593433.53 5028155.58 -125.22
max x y z: 25594452.33 5029174.38 113.36
start of waveform data packet record: 0
start of first extended variable length record: 0
number of extended_variable length records: 0
extended number of point records: 17514732
extended number of points by return: 14366216 2429785 563140 127123 23629 4178 661 0 0 0 0 0 0 0 0
variable length header record 1 of 2:
reserved 0
user ID 'LASF_Projection'
record ID 34735
length after header 40
description 'by LAStools of rapidlasso GmbH'
GeoKeyDirectoryTag version 1.1.0 number of keys 4
key 1024 tiff_tag_location 0 count 1 value_offset 1 - GTModelTypeGeoKey: ModelTypeProjected
key 3072 tiff_tag_location 0 count 1 value_offset 8083 - ProjectedCSTypeGeoKey: NAD83(CSRS)v6 / MTM Nova Scotia zone 5
key 3076 tiff_tag_location 0 count 1 value_offset 9001 - ProjLinearUnitsGeoKey: Linear_Meter
key 4099 tiff_tag_location 0 count 1 value_offset 9001 - VerticalUnitsGeoKey: Linear_Meter
variable length header record 2 of 2:
reserved 0
user ID 'LASF_Projection'
record ID 2112
length after header 682
description 'by LAStools of rapidlasso GmbH'
WKT OGC COORDINATE SYSTEM:
PROJCS["NAD83(CSRS)v6 / NAD83(CSRS)v6 / MTM Nova Scotia zone 5",GEOGCS["NAD83(CSRS)v6",DATUM["North American Datum of 1983 (CSRS) version 6",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORIT
SG","7019"]],AUTHORITY["EPSG","1197"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","8252"]],PROJECTION["Transverse_
tor"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-64.5],PARAMETER["scale_factor",0.9999],PARAMETER["false_easting",25500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHO
"EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","8083"]]
LASzip compression (version 3.4r3 c3 50000): POINT14 3
LAStiling (idx 205963, lvl 9, sub 0, bbox 187000 4695000 699000 5207000) (size 1000 x 1000, buffer 0)
reporting minimum and maximum for all LAS point record entries ...
X 9343353 9445233
Y 2815558 2917438
Z -12522 11336
intensity 6 21432
return_number 1 7
number_of_returns 1 7
edge_of_flight_line 0 1
scan_direction_flag 0 1
classification 2 18
scan_angle_rank -33 32
user_data 0 255
point_source_ID 10011 30145
gps_time 246362877.637832 248873384.067001
extended_return_number 1 7
extended_number_of_returns 1 7
extended_classification 2 18
extended_scan_angle -5500 5333
extended_scanner_channel 1 3
number of first returns: 14366216
number of intermediate returns: 718552
number of last returns: 14366903
number of single returns: 11936939
overview over extended number of returns of given pulse: 11936939 3733305 1307530 413793 97960 20415 4790 0 0 0 0 0 0 0 0
histogram of classification of points:
2080670 ground (2)
6408505 low vegetation (3)
926788 medium vegetation (4)
8093469 high vegetation (5)
2712 noise (7)
2525 bridge deck (17)
63 Reserved for ASPRS Definition (18)
+-> flagged as withheld: 1186536
+---> 461025 of those are low vegetation (3)
+---> 57107 of those are medium vegetation (4)
+---> 665609 of those are high vegetation (5)
+---> 2712 of those are noise (7)
+---> 20 of those are bridge deck (17)
+---> 63 of those are Reserved for ASPRS Definition (18)
+-> flagged as extended overlap: 5071114
+---> 446228 of those are ground (2)
+---> 1995469 of those are low vegetation (3)
+---> 288884 of those are medium vegetation (4)
+---> 2339986 of those are high vegetation (5)
+---> 542 of those are noise (7)
+---> 5 of those are Reserved for ASPRS Definition (18)
E:\LAStools\bin>las2las -i 476_5026_201901.laz -epsg 2961 -remove_tiling_vlr -target_epsg 8083 -odix _8083 -olaz
SERIOUS WARNING: horizontal datum of source and target incompatible.
E:\LAStools\bin>lasinfo -i 476_5026_201901_8083.laz
lasinfo (201111) report for '476_5026_201901_8083.laz'
reporting all LAS header entries:
file signature: 'LASF'
file source ID: 0
global_encoding: 17
project ID GUID data 1-4: 00000000-0000-0000-0000-000000000000
version major.minor: 1.4
system identifier: 'LAStools (c) by rapidlasso GmbH'
generating software: 'las2las (version 201030)'
file creation day/year: 15/2019
header size: 375
offset to point data: 1205
number var. length records: 2
point data format: 6
point data record length: 30
number of point records: 0
number of points by return: 0 0 0 0 0
scale factor x y z: 0.01 0.01 0.01
offset x y z: 25500000 5000000 0
min x y z: 25593433.53 5028155.58 -125.22
max x y z: 25594452.33 5029174.38 113.36
start of waveform data packet record: 0
start of first extended variable length record: 0
number of extended_variable length records: 0
extended number of point records: 17514732
extended number of points by return: 14366216 2429785 563140 127123 23629 4178 661 0 0 0 0 0 0 0 0
variable length header record 1 of 2:
reserved 0
user ID 'LASF_Projection'
record ID 34735
length after header 40
description 'by LAStools of rapidlasso GmbH'
GeoKeyDirectoryTag version 1.1.0 number of keys 4
key 1024 tiff_tag_location 0 count 1 value_offset 1 - GTModelTypeGeoKey: ModelTypeProjected
key 3072 tiff_tag_location 0 count 1 value_offset 8083 - ProjectedCSTypeGeoKey: NAD83(CSRS)v6 / MTM Nova Scotia zone 5
key 3076 tiff_tag_location 0 count 1 value_offset 9001 - ProjLinearUnitsGeoKey: Linear_Meter
key 4099 tiff_tag_location 0 count 1 value_offset 9001 - VerticalUnitsGeoKey: Linear_Meter
variable length header record 2 of 2:
reserved 0
user ID 'LASF_Projection'
record ID 2112
length after header 682
description 'by LAStools of rapidlasso GmbH'
WKT OGC COORDINATE SYSTEM:
PROJCS["NAD83(CSRS)v6 / NAD83(CSRS)v6 / MTM Nova Scotia zone 5",GEOGCS["NAD83(CSRS)v6",DATUM["North American Datum of 1983 (CSRS) version 6",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY
SG","7019"]],AUTHORITY["EPSG","1197"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","8252"]],PROJECTION["Transverse_M
tor"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-64.5],PARAMETER["scale_factor",0.9999],PARAMETER["false_easting",25500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHOR
"EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","8083"]]
LASzip compression (version 3.4r3 c3 50000): POINT14 3
reporting minimum and maximum for all LAS point record entries ...
X 9343353 9445233
Y 2815558 2917438
Z -12522 11336
intensity 6 21432
return_number 1 7
number_of_returns 1 7
edge_of_flight_line 0 1
scan_direction_flag 0 1
classification 2 18
scan_angle_rank -33 32
user_data 0 255
point_source_ID 10011 30145
gps_time 246362877.637832 248873384.067001
extended_return_number 1 7
extended_number_of_returns 1 7
extended_classification 2 18
extended_scan_angle -5500 5333
extended_scanner_channel 1 3
number of first returns: 14366216
number of intermediate returns: 718552
number of last returns: 14366903
number of single returns: 11936939
overview over extended number of returns of given pulse: 11936939 3733305 1307530 413793 97960 20415 4790 0 0 0 0 0 0 0 0
histogram of classification of points:
2080670 ground (2)
6408505 low vegetation (3)
926788 medium vegetation (4)
8093469 high vegetation (5)
2712 noise (7)
2525 bridge deck (17)
63 Reserved for ASPRS Definition (18)
+-> flagged as withheld: 1186536
+---> 461025 of those are low vegetation (3)
+---> 57107 of those are medium vegetation (4)
+---> 665609 of those are high vegetation (5)
+---> 2712 of those are noise (7)
+---> 20 of those are bridge deck (17)
+---> 63 of those are Reserved for ASPRS Definition (18)
+-> flagged as extended overlap: 5071114
+---> 446228 of those are ground (2)
+---> 1995469 of those are low vegetation (3)
+---> 288884 of those are medium vegetation (4)
+---> 2339986 of those are high vegetation (5)
+---> 542 of those are noise (7)
+---> 5 of those are Reserved for ASPRS Definition (18)