NAD83(CSRS) for las2las
113 views
Skip to first unread message
Lukas Fraser
Aug 17, 2016, 2:38:27 PM8/17/16
to LAStools - efficient tools for LiDAR processing
Hi All,
I am trying to change the las header information to UTM 20 NAD83 CSRS (Canadian Spatial Reference System) using las2las. The EPSG code is 2961, which is currently not an option for las2las. Here is a link for that code: http://spatialreference.org/ref/epsg/2961/
I attempted to use -nad83 -utm 20N, but I believe that this is for NAD83(2011) which is not the same thing. Is there currently any way to do this?
Lukas
Evon Silvia
Aug 17, 2016, 2:47:27 PM8/17/16
to last...@googlegroups.com
LAStools by default uses the vanilla version of NAD83 (unspecified realization). You should just be able to use las2las and the argument -epsg 2691 as in the following:
las2las -i my.las -o csrs.las -epsg 3691
If that doesn't work, you need to update to the latest version of LAStools.
Evon
--
Download LAStools at
http://lastools.org
http://rapidlasso.com
Be social with LAStools at
http://facebook.com/LAStools
http://twitter.com/LAStools
http://linkedin.com/groups/LAStools-4408378
Manage your settings at
http://groups.google.com/group/lastools/subscribe
Jeffrey Osborne
Aug 26, 2016, 1:30:15 PM8/26/16
to LAStools - efficient tools for LiDAR processing
There is also support for the Canadian vertical datums and OGC WKT instead of the GeoTiff keys. https://groups.google.com/d/msg/lastools/2mIhkedcp9I/jxrGfGIxAQAJ
las2las -i in.las -epsg 2961 -vertical_cgvd28 -o out.las gives:
reporting all LAS header entries:
file signature: 'LASF'
file source ID: 0
global_encoding: 0
project ID GUID data 1-4: 00000000-0000-0000-F500-DD41E34B89C2
version major.minor: 1.1
system identifier: 'LAStools (c) by rapidlasso GmbH'
generating software: 'las2las (version 160606)'
file creation day/year: 247/2014
header size: 227
offset to point data: 329
number var. length records: 1
point data format: 1
point data record length: 28
number of point records: 8095202
number of points by return: 7001538 808595 233731 45175 5731
scale factor x y z: 0.01 0.01 0.01
offset x y z: 500000 5000000 0
min x y z: 579514.93 5051218.26 0.40
max x y z: 581098.40 5052359.28 66.93
variable length header record 1 of 1:
reserved 43707
user ID 'LASF_Projection'
record ID 34735
length after header 48
description 'by LAStools of rapidlasso GmbH'
GeoKeyDirectoryTag version 1.1.0 number of keys 5
key 1024 tiff_tag_location 0 count 1 value_offset 1 - GTModelTypeGeoKey: ModelTypeProjected
key 3072 tiff_tag_location 0 count 1 value_offset 2961 - ProjectedCSTypeGeoKey: NAD83(CSRS) / UTM zone 20N
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
key 4096 tiff_tag_location 0 count 1 value_offset 5114 - VerticalCSTypeGeoKey: VertCS_Canadian_Geodetic_Vertical_Datum_1928
reporting minimum and maximum for all LAS point record entries ...
X 7951493 8109840
Y 5121826 5235928
Z 40 6693
intensity 8 65520
return_number 1 7
number_of_returns 1 7
edge_of_flight_line 0 0
scan_direction_flag 0 0
classification 1 31
scan_angle_rank -23 25
user_data 205 205
point_source_ID 0 3
gps_time 567369.656098 568774.437919
number of first returns: 7001538
number of intermediate returns: 287217
number of last returns: 6983027
number of single returns: 6176580
WARNING: there are 417 points with return number 6
WARNING: there are 15 points with return number 7
overview over number of returns of given pulse: 6176580 1155340 571515 161219 27695 2688 165
histogram of classification of points:
2127 unclassified (1)
5865989 ground (2)
497044 low vegetation (3)
386198 medium vegetation (4)
1341057 high vegetation (5)
2787 Reserved for ASPRS Definition (31)
And las2las -i in.las -epsg 2961 -vertical_cgvd2013 -set_ogc_wkt -o out.las
reporting all LAS header entries:
file signature: 'LASF'
file source ID: 0
global_encoding: 0
project ID GUID data 1-4: 00000000-0000-0000-F500-DD41E34B89C2
version major.minor: 1.1
system identifier: 'LAStools (c) by rapidlasso GmbH'
generating software: 'las2las (version 160606)'
file creation day/year: 247/2014
header size: 227
offset to point data: 1230
number var. length records: 2
point data format: 1
point data record length: 28
number of point records: 8095202
number of points by return: 7001538 808595 233731 45175 5731
scale factor x y z: 0.01 0.01 0.01
offset x y z: 500000 5000000 0
min x y z: 579514.93 5051218.26 0.40
max x y z: 581098.40 5052359.28 66.93
variable length header record 1 of 2:
reserved 43707
user ID 'LASF_Projection'
record ID 34735
length after header 48
description 'by LAStools of rapidlasso GmbH'
GeoKeyDirectoryTag version 1.1.0 number of keys 5
key 1024 tiff_tag_location 0 count 1 value_offset 1 - GTModelTypeGeoKey: ModelTypeProjected
key 3072 tiff_tag_location 0 count 1 value_offset 2961 - ProjectedCSTypeGeoKey: NAD83(CSRS) / UTM zone 20N
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
key 4096 tiff_tag_location 0 count 1 value_offset 1127 - VerticalCSTypeGeoKey: VertCS_Canadian_Geodetic_Vertical_Datum_2013
variable length header record 2 of 2:
reserved 43707
user ID 'LASF_Projection'
record ID 2112
length after header 847
description 'by LAStools of rapidlasso GmbH'
WKT OGC COORDINATE SYSTEM:
PROJCS["NAD83(CSRS) / UTM zone 20N",GEOGCS["NAD83(CSRS)",DATUM["NAD83_Canadian_Spatial_Reference_System",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6140"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4617"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-63],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","2961"]]VERT_CS["CGVD2013",VERT_DATUM["Canadian Geodetic Vertical Datum of 2013",2005,AUTHORITY["EPSG","1127"]],UNIT["metre",1.0,AUTHORITY["EPSG","9001"]]AXIS["Gravity-related height",UP],AUTHORITY["EPSG","6647"]]
reporting minimum and maximum for all LAS point record entries ...
X 7951493 8109840
Y 5121826 5235928
Z 40 6693
intensity 8 65520
return_number 1 7
number_of_returns 1 7
edge_of_flight_line 0 0
scan_direction_flag 0 0
classification 1 31
scan_angle_rank -23 25
user_data 205 205
point_source_ID 0 3
gps_time 567369.656098 568774.437919
number of first returns: 7001538
number of intermediate returns: 287217
number of last returns: 6983027
number of single returns: 6176580
WARNING: there are 417 points with return number 6
WARNING: there are 15 points with return number 7
overview over number of returns of given pulse: 6176580 1155340 571515 161219 27695 2688 165
histogram of classification of points:
2127 unclassified (1)
5865989 ground (2)
497044 low vegetation (3)
386198 medium vegetation (4)
1341057 high vegetation (5)
2787 Reserved for ASPRS Definition (31)
Reply all
Reply to author
Forward
0 new messages