Copy extra byte values into R,G,B Attributes

159 views
Skip to first unread message

Aicke Damrau

unread,
Jun 15, 2020, 7:34:30 AM6/15/20
to LAStools - efficient tools for LiDAR processing
Hi all,

My task is to copy the "Reflectance" and "Pulse width" from Extra bytes into the Red and Green Attribute.

  description          'RIEGL Extra Bytes'
    Extra Byte Descriptions
      data type: 3 (unsigned short), name "Amplitude", description: "Echo signal amplitude [dB]", min: 0, max: 10000, scale: 0.01, offset: 0 (not set)
      data type: 3 (unsigned short), name "Pulse width", description: "Full width at half maximum [ns]", min: 1, max: 10000, scale: 0.1, offset: 0 (not set)
      data type: 4 (short), name "Reflectance", description: "Echo signal reflectance [dB]", min: -5000, max: 15000, scale: 0.01, offset: 0 (not set)
the header is followed by 2 user-defined bytes
LAStiling (idx 1013, lvl 6, sub 0, bbox 590000 5.751e+006 654000 5.815e+006)
reporting minimum and maximum for all LAS point record entries ...
  X          -1516000000 -1512000001
  Y          -884000000 -880000001
  Z              367267     631862
  intensity           0       4494
  return_number       0          7
  number_of_returns   0          7
  edge_of_flight_line 0          1
  scan_direction_flag 0          0
  classification      1         24
  scan_angle_rank   -31         32
  user_data           0          1
  point_source_ID    14        999
  gps_time -40787201.000000 263655956.894434
  Color R 0 0
        G 0 0
        B 0 0

Has any body an idea, how to solve?

Thanks 
Aicke

Martin Isenburg

unread,
Jun 18, 2020, 9:03:59 PM6/18/20
to LAStools - efficient command line tools for LIDAR processing
Hello,

I have added a capability to copy default field or extra attributes of a point to so called "registers" that can then be scaled or translated and then be copied back into the same or some other field. Here some examples. I expect this functionality to grow in the coming months based on new user demands:

lasinfo -i 2621_5779_1_20200121.laz -copy_attribute_into_register 1 0 -copy_register_into_R 0  
lasinfo -i 2621_5779_1_20200121.laz -copy_attribute_into_register 1 0 -scale_register 0 10.0 -copy_register_into_R 0
lasinfo -i 2621_5779_1_20200121.laz -copy_attribute_into_register 2 1 -translate_register 1 20 -scale_register 1 10.0 -copy_register_into_G 1
lasinfo -i 2621_5779_1_20200121.laz -copy_attribute_into_register 2 1 -translate_register 1 20 -scale_register 1 100.0 -copy_register_into_G 1  

Here the full example:

E:\software\LAStools\bin>lasinfo -i 2621_5779_1_20200121.laz -copy_attribute_into_register 1 0 -copy_register_into_R 0
lasinfo (200619) report for '2621_5779_1_20200121.laz'
reporting all LAS header entries:
  file signature:             'LASF'
  file source ID:             0
  global_encoding:            1
  project ID GUID data 1-4:   00000000-0000-0000-0000-000000000000
  version major.minor:        1.2
  system identifier:          'LAStools (c) by rapidlasso GmbH'
  generating software:        'las2las (version 200509)'
  file creation day/year:     164/2020
  header size:                227
  offset to point data:       2385
  number var. length records: 5
  point data format:          3
  point data record length:   40
  number of point records:    13576466
  number of points by return: 10539428 1766996 909481 333208 25380
  scale factor x y z:         0.00025 0.00025 0.00025
  offset x y z:               1000000 6000000 0
  min x y z:                  621000.00000 5779000.00000 91.81675
  max x y z:                  621999.99975 5779999.99975 157.96550
variable length header record 1 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            34735
  length after header  208
  description          'GeoKeyDirectoryTag (mandatory)'
    GeoKeyDirectoryTag version 1.1.0 number of keys 25
      key 1024 tiff_tag_location 0 count 1 value_offset 1 - GTModelTypeGeoKey: ModelTypeProjected
      key 1025 tiff_tag_location 0 count 1 value_offset 2 - GTRasterTypeGeoKey: RasterPixelIsPoint
      key 1026 tiff_tag_location 34737 count 33 value_offset 0 - GTCitationGeoKey: ETRS89 / UTM zone 32N/ORTH_GCG16
      key 2048 tiff_tag_location 0 count 1 value_offset 32767 - GeographicTypeGeoKey: user-defined
      key 2049 tiff_tag_location 34737 count 7 value_offset 66 - GeogCitationGeoKey: ETRS89
      key 2050 tiff_tag_location 0 count 1 value_offset 32767 - GeogGeodeticDatumGeoKey: user-defined
      key 2051 tiff_tag_location 0 count 1 value_offset 32767 - GeogPrimeMeridianGeoKey: user-defined
      key 2052 tiff_tag_location 0 count 1 value_offset 9001 - GeogLinearUnitsGeoKey: Linear_Meter
      key 2054 tiff_tag_location 0 count 1 value_offset 9102 - GeogAngularUnitsGeoKey: Angular_Degree
      key 2056 tiff_tag_location 0 count 1 value_offset 32767 - GeogEllipsoidGeoKey: user-defined
      key 2057 tiff_tag_location 34736 count 1 value_offset 0 - GeogSemiMajorAxisGeoKey: 6378137
      key 2059 tiff_tag_location 34736 count 1 value_offset 1 - GeogInvFlatteningGeoKey: 298.2572221
      key 2060 tiff_tag_location 0 count 1 value_offset 9102 - GeogAzimuthUnitsGeoKey: Angular_Degree
      key 2061 tiff_tag_location 34736 count 1 value_offset 2 - GeogPrimeMeridianLongGeoKey: 0
      key 3072 tiff_tag_location 0 count 1 value_offset 32767 - ProjectedCSTypeGeoKey: user-defined
      key 3073 tiff_tag_location 34737 count 33 value_offset 33 - PCSCitationGeoKey: ETRS89 / UTM zone 32N/ORTH_GCG16
      key 3074 tiff_tag_location 0 count 1 value_offset 32767 - ProjectionGeoKey: user-defined
      key 3075 tiff_tag_location 0 count 1 value_offset 1 - ProjCoordTransGeoKey: CT_TransverseMercator
      key 3076 tiff_tag_location 0 count 1 value_offset 9001 - ProjLinearUnitsGeoKey: Linear_Meter
      key 3080 tiff_tag_location 34736 count 1 value_offset 4 - ProjNatOriginLongGeoKey: 9
      key 3081 tiff_tag_location 34736 count 1 value_offset 3 - ProjNatOriginLatGeoKey: 0
      key 3082 tiff_tag_location 34736 count 1 value_offset 5 - ProjFalseEastingGeoKey: 500000
      key 3083 tiff_tag_location 34736 count 1 value_offset 6 - ProjFalseNorthingGeoKey: 0
      key 3092 tiff_tag_location 34736 count 1 value_offset 7 - ProjScaleAtNatOriginGeoKey: 0.9996
      key 4099 tiff_tag_location 0 count 1 value_offset 9001 - VerticalUnitsGeoKey: Linear_Meter
variable length header record 2 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            34736
  length after header  64
  description          'GeoDoubleParamsTag (optional)'
    GeoDoubleParamsTag (number of doubles 8)
      6.37814e+006 298.257 0 0 9 500000 0 0.9996
variable length header record 3 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            34737
  length after header  73
  description          'GeoASCIIParamsTag (optional)'
    GeoAsciiParamsTag (number of characters 73)
      ETRS89 / UTM zone 32N/ORTH_GCG16|ETRS89 / UTM zone 32N/ORTH_GCG16|ETRS89|
variable length header record 4 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            2112
  length after header  965
  description          'OGC COORDINATE SYSTEM WKT'
    WKT OGC COORDINATE SYSTEM:
    COMPD_CS["ETRS89 / UTM zone 32N/ORTH_GCG16",PROJCS["ETRS89 / UTM zone 32N/ORTH_GCG16",GEOGCS["ETRS89 / UTM zone 32N/ORTH_GCG16",DATUM["ETRS89",SPHEROID["GRS 1980",6378137.000,298.257222101,AUTHORI
TY["EPSG","7019"]],AUTHORITY["EPSG","6258"]],PRIMEM["Greenwich",0.0000000000000000,AUTHORITY["EPSG","8901"]],UNIT["Degree",0.01745329251994329547,AUTHORITY["EPSG","9102"]],AUTHORITY["EPSG","25832"]],P
ROJECTION["Transverse_Mercator",AUTHORITY["EPSG","9807"]],PARAMETER["latitude_of_origin",0.0000000000000000],PARAMETER["central_meridian",9.0000000000000000],PARAMETER["scale_factor",0.999600000000000
0],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","2
5832"]],VERT_CS["Orthometric Heights",VERT_DATUM["GCG2016",2005,AUTHORITY["EPSG","0"]],UNIT["Meter",1.00000000000000000000,AUTHORITY["EPSG","9001"]],AXIS["Height",UP]]]
variable length header record 5 of 5:
  reserved             0
  user ID              'LASF_Spec'
  record ID            4
  length after header  576

  description          'RIEGL Extra Bytes'
    Extra Byte Descriptions
      data type: 3 (unsigned short), name "Amplitude", description: "Echo signal amplitude [dB]", min: 0, max: 10000, scale: 0.01, offset: 0 (not set)
      data type: 3 (unsigned short), name "Pulse width", description: "Full width at half maximum [ns]", min: 1, max: 10000, scale: 0.1, offset: 0 (not set)
      data type: 4 (short), name "Reflectance", description: "Echo signal reflectance [dB]", min: -5000, max: 15000, scale: 0.01, offset: 0 (not set)
the header is followed by 2 user-defined bytes
LASzip compression (version 3.4r3 c2 50000): POINT10 2 GPSTIME11 2 RGB12 2 BYTE 2
LAStiling (idx 1013, lvl 6, sub 0, bbox 590000 5751000 654000 5815000) (size 1000 x 1000, buffer 0)

reporting minimum and maximum for all LAS point record entries ...
  X          -1516000000 -1512000001
  Y          -884000000 -880000001
  Z              367267     631862
  intensity           0       4494
  return_number       0          7
  number_of_returns   0          7
  edge_of_flight_line 0          1
  scan_direction_flag 0          0
  classification      1         24
  scan_angle_rank   -31         32
  user_data           0          1
  point_source_ID    14        999
  gps_time -40787201.000000 263655956.894434
  Color R 0 34

        G 0 0
        B 0 0
  attribute0          0      28.08  ('Amplitude')
  attribute1          0       34.5  ('Pulse width')
  attribute2     -18.54      26.14  ('Reflectance')
number of first returns:        10540061
number of intermediate returns: 1268951
number of last returns:         10542131
number of single returns:       8774677
WARNING: there are 633 points with return number 0
WARNING: there are 1272 points with return number 6
WARNING: there are 68 points with return number 7
overview over number of returns of given pulse: 8774044 1714674 1727941 1230537 120921 7246 470
WARNING: there are 633 points with a number of returns of given pulse of 0
histogram of classification of points:
           21376  unclassified (1)
         6062893  ground (2)
            2190  noise (7)
             632  keypoint (8)
            8005  water (9)
         4214965  overlap (12)
         3256456  Reserved for ASPRS Definition (20)
            9949  Reserved for ASPRS Definition (24)

E:\software\LAStools\bin>lasinfo -i 2621_5779_1_20200121.laz -copy_attribute_into_register 1 0 -scale_register 0 10.0 -copy_register_into_R 0
lasinfo (200619) report for '2621_5779_1_20200121.laz'
reporting all LAS header entries:
  file signature:             'LASF'
  file source ID:             0
  global_encoding:            1
  project ID GUID data 1-4:   00000000-0000-0000-0000-000000000000
  version major.minor:        1.2
  system identifier:          'LAStools (c) by rapidlasso GmbH'
  generating software:        'las2las (version 200509)'
  file creation day/year:     164/2020
  header size:                227
  offset to point data:       2385
  number var. length records: 5
  point data format:          3
  point data record length:   40
  number of point records:    13576466
  number of points by return: 10539428 1766996 909481 333208 25380
  scale factor x y z:         0.00025 0.00025 0.00025
  offset x y z:               1000000 6000000 0
  min x y z:                  621000.00000 5779000.00000 91.81675
  max x y z:                  621999.99975 5779999.99975 157.96550
variable length header record 1 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            34735
  length after header  208
  description          'GeoKeyDirectoryTag (mandatory)'
    GeoKeyDirectoryTag version 1.1.0 number of keys 25
      key 1024 tiff_tag_location 0 count 1 value_offset 1 - GTModelTypeGeoKey: ModelTypeProjected
      key 1025 tiff_tag_location 0 count 1 value_offset 2 - GTRasterTypeGeoKey: RasterPixelIsPoint
      key 1026 tiff_tag_location 34737 count 33 value_offset 0 - GTCitationGeoKey: ETRS89 / UTM zone 32N/ORTH_GCG16
      key 2048 tiff_tag_location 0 count 1 value_offset 32767 - GeographicTypeGeoKey: user-defined
      key 2049 tiff_tag_location 34737 count 7 value_offset 66 - GeogCitationGeoKey: ETRS89
      key 2050 tiff_tag_location 0 count 1 value_offset 32767 - GeogGeodeticDatumGeoKey: user-defined
      key 2051 tiff_tag_location 0 count 1 value_offset 32767 - GeogPrimeMeridianGeoKey: user-defined
      key 2052 tiff_tag_location 0 count 1 value_offset 9001 - GeogLinearUnitsGeoKey: Linear_Meter
      key 2054 tiff_tag_location 0 count 1 value_offset 9102 - GeogAngularUnitsGeoKey: Angular_Degree
      key 2056 tiff_tag_location 0 count 1 value_offset 32767 - GeogEllipsoidGeoKey: user-defined
      key 2057 tiff_tag_location 34736 count 1 value_offset 0 - GeogSemiMajorAxisGeoKey: 6378137
      key 2059 tiff_tag_location 34736 count 1 value_offset 1 - GeogInvFlatteningGeoKey: 298.2572221
      key 2060 tiff_tag_location 0 count 1 value_offset 9102 - GeogAzimuthUnitsGeoKey: Angular_Degree
      key 2061 tiff_tag_location 34736 count 1 value_offset 2 - GeogPrimeMeridianLongGeoKey: 0
      key 3072 tiff_tag_location 0 count 1 value_offset 32767 - ProjectedCSTypeGeoKey: user-defined
      key 3073 tiff_tag_location 34737 count 33 value_offset 33 - PCSCitationGeoKey: ETRS89 / UTM zone 32N/ORTH_GCG16
      key 3074 tiff_tag_location 0 count 1 value_offset 32767 - ProjectionGeoKey: user-defined
      key 3075 tiff_tag_location 0 count 1 value_offset 1 - ProjCoordTransGeoKey: CT_TransverseMercator
      key 3076 tiff_tag_location 0 count 1 value_offset 9001 - ProjLinearUnitsGeoKey: Linear_Meter
      key 3080 tiff_tag_location 34736 count 1 value_offset 4 - ProjNatOriginLongGeoKey: 9
      key 3081 tiff_tag_location 34736 count 1 value_offset 3 - ProjNatOriginLatGeoKey: 0
      key 3082 tiff_tag_location 34736 count 1 value_offset 5 - ProjFalseEastingGeoKey: 500000
      key 3083 tiff_tag_location 34736 count 1 value_offset 6 - ProjFalseNorthingGeoKey: 0
      key 3092 tiff_tag_location 34736 count 1 value_offset 7 - ProjScaleAtNatOriginGeoKey: 0.9996
      key 4099 tiff_tag_location 0 count 1 value_offset 9001 - VerticalUnitsGeoKey: Linear_Meter
variable length header record 2 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            34736
  length after header  64
  description          'GeoDoubleParamsTag (optional)'
    GeoDoubleParamsTag (number of doubles 8)
      6.37814e+006 298.257 0 0 9 500000 0 0.9996
variable length header record 3 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            34737
  length after header  73
  description          'GeoASCIIParamsTag (optional)'
    GeoAsciiParamsTag (number of characters 73)
      ETRS89 / UTM zone 32N/ORTH_GCG16|ETRS89 / UTM zone 32N/ORTH_GCG16|ETRS89|
variable length header record 4 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            2112
  length after header  965
  description          'OGC COORDINATE SYSTEM WKT'
    WKT OGC COORDINATE SYSTEM:
    COMPD_CS["ETRS89 / UTM zone 32N/ORTH_GCG16",PROJCS["ETRS89 / UTM zone 32N/ORTH_GCG16",GEOGCS["ETRS89 / UTM zone 32N/ORTH_GCG16",DATUM["ETRS89",SPHEROID["GRS 1980",6378137.000,298.257222101,AUTHORI
TY["EPSG","7019"]],AUTHORITY["EPSG","6258"]],PRIMEM["Greenwich",0.0000000000000000,AUTHORITY["EPSG","8901"]],UNIT["Degree",0.01745329251994329547,AUTHORITY["EPSG","9102"]],AUTHORITY["EPSG","25832"]],P
ROJECTION["Transverse_Mercator",AUTHORITY["EPSG","9807"]],PARAMETER["latitude_of_origin",0.0000000000000000],PARAMETER["central_meridian",9.0000000000000000],PARAMETER["scale_factor",0.999600000000000
0],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","2
5832"]],VERT_CS["Orthometric Heights",VERT_DATUM["GCG2016",2005,AUTHORITY["EPSG","0"]],UNIT["Meter",1.00000000000000000000,AUTHORITY["EPSG","9001"]],AXIS["Height",UP]]]
variable length header record 5 of 5:
  reserved             0
  user ID              'LASF_Spec'
  record ID            4
  length after header  576

  description          'RIEGL Extra Bytes'
    Extra Byte Descriptions
      data type: 3 (unsigned short), name "Amplitude", description: "Echo signal amplitude [dB]", min: 0, max: 10000, scale: 0.01, offset: 0 (not set)
      data type: 3 (unsigned short), name "Pulse width", description: "Full width at half maximum [ns]", min: 1, max: 10000, scale: 0.1, offset: 0 (not set)
      data type: 4 (short), name "Reflectance", description: "Echo signal reflectance [dB]", min: -5000, max: 15000, scale: 0.01, offset: 0 (not set)
the header is followed by 2 user-defined bytes
LASzip compression (version 3.4r3 c2 50000): POINT10 2 GPSTIME11 2 RGB12 2 BYTE 2
LAStiling (idx 1013, lvl 6, sub 0, bbox 590000 5751000 654000 5815000) (size 1000 x 1000, buffer 0)

reporting minimum and maximum for all LAS point record entries ...
  X          -1516000000 -1512000001
  Y          -884000000 -880000001
  Z              367267     631862
  intensity           0       4494
  return_number       0          7
  number_of_returns   0          7
  edge_of_flight_line 0          1
  scan_direction_flag 0          0
  classification      1         24
  scan_angle_rank   -31         32
  user_data           0          1
  point_source_ID    14        999
  gps_time -40787201.000000 263655956.894434
  Color R 0 345

        G 0 0
        B 0 0
  attribute0          0      28.08  ('Amplitude')
  attribute1          0       34.5  ('Pulse width')
  attribute2     -18.54      26.14  ('Reflectance')
number of first returns:        10540061
number of intermediate returns: 1268951
number of last returns:         10542131
number of single returns:       8774677
WARNING: there are 633 points with return number 0
WARNING: there are 1272 points with return number 6
WARNING: there are 68 points with return number 7
overview over number of returns of given pulse: 8774044 1714674 1727941 1230537 120921 7246 470
WARNING: there are 633 points with a number of returns of given pulse of 0
histogram of classification of points:
           21376  unclassified (1)
         6062893  ground (2)
            2190  noise (7)
             632  keypoint (8)
            8005  water (9)
         4214965  overlap (12)
         3256456  Reserved for ASPRS Definition (20)
            9949  Reserved for ASPRS Definition (24)

E:\software\LAStools\bin>lasinfo -i 2621_5779_1_20200121.laz -copy_attribute_into_register 2 1 -translate_register 1 20 -scale_register 1 10.0 -copy_register_into_G 1
lasinfo (200619) report for '2621_5779_1_20200121.laz'
reporting all LAS header entries:
  file signature:             'LASF'
  file source ID:             0
  global_encoding:            1
  project ID GUID data 1-4:   00000000-0000-0000-0000-000000000000
  version major.minor:        1.2
  system identifier:          'LAStools (c) by rapidlasso GmbH'
  generating software:        'las2las (version 200509)'
  file creation day/year:     164/2020
  header size:                227
  offset to point data:       2385
  number var. length records: 5
  point data format:          3
  point data record length:   40
  number of point records:    13576466
  number of points by return: 10539428 1766996 909481 333208 25380
  scale factor x y z:         0.00025 0.00025 0.00025
  offset x y z:               1000000 6000000 0
  min x y z:                  621000.00000 5779000.00000 91.81675
  max x y z:                  621999.99975 5779999.99975 157.96550
variable length header record 1 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            34735
  length after header  208
  description          'GeoKeyDirectoryTag (mandatory)'
    GeoKeyDirectoryTag version 1.1.0 number of keys 25
      key 1024 tiff_tag_location 0 count 1 value_offset 1 - GTModelTypeGeoKey: ModelTypeProjected
      key 1025 tiff_tag_location 0 count 1 value_offset 2 - GTRasterTypeGeoKey: RasterPixelIsPoint
      key 1026 tiff_tag_location 34737 count 33 value_offset 0 - GTCitationGeoKey: ETRS89 / UTM zone 32N/ORTH_GCG16
      key 2048 tiff_tag_location 0 count 1 value_offset 32767 - GeographicTypeGeoKey: user-defined
      key 2049 tiff_tag_location 34737 count 7 value_offset 66 - GeogCitationGeoKey: ETRS89
      key 2050 tiff_tag_location 0 count 1 value_offset 32767 - GeogGeodeticDatumGeoKey: user-defined
      key 2051 tiff_tag_location 0 count 1 value_offset 32767 - GeogPrimeMeridianGeoKey: user-defined
      key 2052 tiff_tag_location 0 count 1 value_offset 9001 - GeogLinearUnitsGeoKey: Linear_Meter
      key 2054 tiff_tag_location 0 count 1 value_offset 9102 - GeogAngularUnitsGeoKey: Angular_Degree
      key 2056 tiff_tag_location 0 count 1 value_offset 32767 - GeogEllipsoidGeoKey: user-defined
      key 2057 tiff_tag_location 34736 count 1 value_offset 0 - GeogSemiMajorAxisGeoKey: 6378137
      key 2059 tiff_tag_location 34736 count 1 value_offset 1 - GeogInvFlatteningGeoKey: 298.2572221
      key 2060 tiff_tag_location 0 count 1 value_offset 9102 - GeogAzimuthUnitsGeoKey: Angular_Degree
      key 2061 tiff_tag_location 34736 count 1 value_offset 2 - GeogPrimeMeridianLongGeoKey: 0
      key 3072 tiff_tag_location 0 count 1 value_offset 32767 - ProjectedCSTypeGeoKey: user-defined
      key 3073 tiff_tag_location 34737 count 33 value_offset 33 - PCSCitationGeoKey: ETRS89 / UTM zone 32N/ORTH_GCG16
      key 3074 tiff_tag_location 0 count 1 value_offset 32767 - ProjectionGeoKey: user-defined
      key 3075 tiff_tag_location 0 count 1 value_offset 1 - ProjCoordTransGeoKey: CT_TransverseMercator
      key 3076 tiff_tag_location 0 count 1 value_offset 9001 - ProjLinearUnitsGeoKey: Linear_Meter
      key 3080 tiff_tag_location 34736 count 1 value_offset 4 - ProjNatOriginLongGeoKey: 9
      key 3081 tiff_tag_location 34736 count 1 value_offset 3 - ProjNatOriginLatGeoKey: 0
      key 3082 tiff_tag_location 34736 count 1 value_offset 5 - ProjFalseEastingGeoKey: 500000
      key 3083 tiff_tag_location 34736 count 1 value_offset 6 - ProjFalseNorthingGeoKey: 0
      key 3092 tiff_tag_location 34736 count 1 value_offset 7 - ProjScaleAtNatOriginGeoKey: 0.9996
      key 4099 tiff_tag_location 0 count 1 value_offset 9001 - VerticalUnitsGeoKey: Linear_Meter
variable length header record 2 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            34736
  length after header  64
  description          'GeoDoubleParamsTag (optional)'
    GeoDoubleParamsTag (number of doubles 8)
      6.37814e+006 298.257 0 0 9 500000 0 0.9996
variable length header record 3 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            34737
  length after header  73
  description          'GeoASCIIParamsTag (optional)'
    GeoAsciiParamsTag (number of characters 73)
      ETRS89 / UTM zone 32N/ORTH_GCG16|ETRS89 / UTM zone 32N/ORTH_GCG16|ETRS89|
variable length header record 4 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            2112
  length after header  965
  description          'OGC COORDINATE SYSTEM WKT'
    WKT OGC COORDINATE SYSTEM:
    COMPD_CS["ETRS89 / UTM zone 32N/ORTH_GCG16",PROJCS["ETRS89 / UTM zone 32N/ORTH_GCG16",GEOGCS["ETRS89 / UTM zone 32N/ORTH_GCG16",DATUM["ETRS89",SPHEROID["GRS 1980",6378137.000,298.257222101,AUTHORI
TY["EPSG","7019"]],AUTHORITY["EPSG","6258"]],PRIMEM["Greenwich",0.0000000000000000,AUTHORITY["EPSG","8901"]],UNIT["Degree",0.01745329251994329547,AUTHORITY["EPSG","9102"]],AUTHORITY["EPSG","25832"]],P
ROJECTION["Transverse_Mercator",AUTHORITY["EPSG","9807"]],PARAMETER["latitude_of_origin",0.0000000000000000],PARAMETER["central_meridian",9.0000000000000000],PARAMETER["scale_factor",0.999600000000000
0],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","2
5832"]],VERT_CS["Orthometric Heights",VERT_DATUM["GCG2016",2005,AUTHORITY["EPSG","0"]],UNIT["Meter",1.00000000000000000000,AUTHORITY["EPSG","9001"]],AXIS["Height",UP]]]
variable length header record 5 of 5:
  reserved             0
  user ID              'LASF_Spec'
  record ID            4
  length after header  576

  description          'RIEGL Extra Bytes'
    Extra Byte Descriptions
      data type: 3 (unsigned short), name "Amplitude", description: "Echo signal amplitude [dB]", min: 0, max: 10000, scale: 0.01, offset: 0 (not set)
      data type: 3 (unsigned short), name "Pulse width", description: "Full width at half maximum [ns]", min: 1, max: 10000, scale: 0.1, offset: 0 (not set)
      data type: 4 (short), name "Reflectance", description: "Echo signal reflectance [dB]", min: -5000, max: 15000, scale: 0.01, offset: 0 (not set)
the header is followed by 2 user-defined bytes
LASzip compression (version 3.4r3 c2 50000): POINT10 2 GPSTIME11 2 RGB12 2 BYTE 2
LAStiling (idx 1013, lvl 6, sub 0, bbox 590000 5751000 654000 5815000) (size 1000 x 1000, buffer 0)

reporting minimum and maximum for all LAS point record entries ...
  X          -1516000000 -1512000001
  Y          -884000000 -880000001
  Z              367267     631862
  intensity           0       4494
  return_number       0          7
  number_of_returns   0          7
  edge_of_flight_line 0          1
  scan_direction_flag 0          0
  classification      1         24
  scan_angle_rank   -31         32
  user_data           0          1
  point_source_ID    14        999
  gps_time -40787201.000000 263655956.894434
  Color R 0 0
        G 14 461
        B 0 0
  attribute0          0      28.08  ('Amplitude')
  attribute1          0       34.5  ('Pulse width')
  attribute2     -18.54      26.14  ('Reflectance')
number of first returns:        10540061
number of intermediate returns: 1268951
number of last returns:         10542131
number of single returns:       8774677
WARNING: there are 633 points with return number 0
WARNING: there are 1272 points with return number 6
WARNING: there are 68 points with return number 7
overview over number of returns of given pulse: 8774044 1714674 1727941 1230537 120921 7246 470
WARNING: there are 633 points with a number of returns of given pulse of 0
histogram of classification of points:
           21376  unclassified (1)
         6062893  ground (2)
            2190  noise (7)
             632  keypoint (8)
            8005  water (9)
         4214965  overlap (12)
         3256456  Reserved for ASPRS Definition (20)
            9949  Reserved for ASPRS Definition (24)

E:\software\LAStools\bin>lasinfo -i 2621_5779_1_20200121.laz -copy_attribute_into_register 2 1 -translate_register 1 20 -scale_register 1 100.0 -copy_register_into_G 1
lasinfo (200619) report for '2621_5779_1_20200121.laz'
reporting all LAS header entries:
  file signature:             'LASF'
  file source ID:             0
  global_encoding:            1
  project ID GUID data 1-4:   00000000-0000-0000-0000-000000000000
  version major.minor:        1.2
  system identifier:          'LAStools (c) by rapidlasso GmbH'
  generating software:        'las2las (version 200509)'
  file creation day/year:     164/2020
  header size:                227
  offset to point data:       2385
  number var. length records: 5
  point data format:          3
  point data record length:   40
  number of point records:    13576466
  number of points by return: 10539428 1766996 909481 333208 25380
  scale factor x y z:         0.00025 0.00025 0.00025
  offset x y z:               1000000 6000000 0
  min x y z:                  621000.00000 5779000.00000 91.81675
  max x y z:                  621999.99975 5779999.99975 157.96550
variable length header record 1 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            34735
  length after header  208
  description          'GeoKeyDirectoryTag (mandatory)'
    GeoKeyDirectoryTag version 1.1.0 number of keys 25
      key 1024 tiff_tag_location 0 count 1 value_offset 1 - GTModelTypeGeoKey: ModelTypeProjected
      key 1025 tiff_tag_location 0 count 1 value_offset 2 - GTRasterTypeGeoKey: RasterPixelIsPoint
      key 1026 tiff_tag_location 34737 count 33 value_offset 0 - GTCitationGeoKey: ETRS89 / UTM zone 32N/ORTH_GCG16
      key 2048 tiff_tag_location 0 count 1 value_offset 32767 - GeographicTypeGeoKey: user-defined
      key 2049 tiff_tag_location 34737 count 7 value_offset 66 - GeogCitationGeoKey: ETRS89
      key 2050 tiff_tag_location 0 count 1 value_offset 32767 - GeogGeodeticDatumGeoKey: user-defined
      key 2051 tiff_tag_location 0 count 1 value_offset 32767 - GeogPrimeMeridianGeoKey: user-defined
      key 2052 tiff_tag_location 0 count 1 value_offset 9001 - GeogLinearUnitsGeoKey: Linear_Meter
      key 2054 tiff_tag_location 0 count 1 value_offset 9102 - GeogAngularUnitsGeoKey: Angular_Degree
      key 2056 tiff_tag_location 0 count 1 value_offset 32767 - GeogEllipsoidGeoKey: user-defined
      key 2057 tiff_tag_location 34736 count 1 value_offset 0 - GeogSemiMajorAxisGeoKey: 6378137
      key 2059 tiff_tag_location 34736 count 1 value_offset 1 - GeogInvFlatteningGeoKey: 298.2572221
      key 2060 tiff_tag_location 0 count 1 value_offset 9102 - GeogAzimuthUnitsGeoKey: Angular_Degree
      key 2061 tiff_tag_location 34736 count 1 value_offset 2 - GeogPrimeMeridianLongGeoKey: 0
      key 3072 tiff_tag_location 0 count 1 value_offset 32767 - ProjectedCSTypeGeoKey: user-defined
      key 3073 tiff_tag_location 34737 count 33 value_offset 33 - PCSCitationGeoKey: ETRS89 / UTM zone 32N/ORTH_GCG16
      key 3074 tiff_tag_location 0 count 1 value_offset 32767 - ProjectionGeoKey: user-defined
      key 3075 tiff_tag_location 0 count 1 value_offset 1 - ProjCoordTransGeoKey: CT_TransverseMercator
      key 3076 tiff_tag_location 0 count 1 value_offset 9001 - ProjLinearUnitsGeoKey: Linear_Meter
      key 3080 tiff_tag_location 34736 count 1 value_offset 4 - ProjNatOriginLongGeoKey: 9
      key 3081 tiff_tag_location 34736 count 1 value_offset 3 - ProjNatOriginLatGeoKey: 0
      key 3082 tiff_tag_location 34736 count 1 value_offset 5 - ProjFalseEastingGeoKey: 500000
      key 3083 tiff_tag_location 34736 count 1 value_offset 6 - ProjFalseNorthingGeoKey: 0
      key 3092 tiff_tag_location 34736 count 1 value_offset 7 - ProjScaleAtNatOriginGeoKey: 0.9996
      key 4099 tiff_tag_location 0 count 1 value_offset 9001 - VerticalUnitsGeoKey: Linear_Meter
variable length header record 2 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            34736
  length after header  64
  description          'GeoDoubleParamsTag (optional)'
    GeoDoubleParamsTag (number of doubles 8)
      6.37814e+006 298.257 0 0 9 500000 0 0.9996
variable length header record 3 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            34737
  length after header  73
  description          'GeoASCIIParamsTag (optional)'
    GeoAsciiParamsTag (number of characters 73)
      ETRS89 / UTM zone 32N/ORTH_GCG16|ETRS89 / UTM zone 32N/ORTH_GCG16|ETRS89|
variable length header record 4 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            2112
  length after header  965
  description          'OGC COORDINATE SYSTEM WKT'
    WKT OGC COORDINATE SYSTEM:
    COMPD_CS["ETRS89 / UTM zone 32N/ORTH_GCG16",PROJCS["ETRS89 / UTM zone 32N/ORTH_GCG16",GEOGCS["ETRS89 / UTM zone 32N/ORTH_GCG16",DATUM["ETRS89",SPHEROID["GRS 1980",6378137.000,298.257222101,AUTHORI
TY["EPSG","7019"]],AUTHORITY["EPSG","6258"]],PRIMEM["Greenwich",0.0000000000000000,AUTHORITY["EPSG","8901"]],UNIT["Degree",0.01745329251994329547,AUTHORITY["EPSG","9102"]],AUTHORITY["EPSG","25832"]],P
ROJECTION["Transverse_Mercator",AUTHORITY["EPSG","9807"]],PARAMETER["latitude_of_origin",0.0000000000000000],PARAMETER["central_meridian",9.0000000000000000],PARAMETER["scale_factor",0.999600000000000
0],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","2
5832"]],VERT_CS["Orthometric Heights",VERT_DATUM["GCG2016",2005,AUTHORITY["EPSG","0"]],UNIT["Meter",1.00000000000000000000,AUTHORITY["EPSG","9001"]],AXIS["Height",UP]]]
variable length header record 5 of 5:
  reserved             0
  user ID              'LASF_Spec'
  record ID            4
  length after header  576

  description          'RIEGL Extra Bytes'
    Extra Byte Descriptions
      data type: 3 (unsigned short), name "Amplitude", description: "Echo signal amplitude [dB]", min: 0, max: 10000, scale: 0.01, offset: 0 (not set)
      data type: 3 (unsigned short), name "Pulse width", description: "Full width at half maximum [ns]", min: 1, max: 10000, scale: 0.1, offset: 0 (not set)
      data type: 4 (short), name "Reflectance", description: "Echo signal reflectance [dB]", min: -5000, max: 15000, scale: 0.01, offset: 0 (not set)
the header is followed by 2 user-defined bytes
LASzip compression (version 3.4r3 c2 50000): POINT10 2 GPSTIME11 2 RGB12 2 BYTE 2
LAStiling (idx 1013, lvl 6, sub 0, bbox 590000 5751000 654000 5815000) (size 1000 x 1000, buffer 0)

reporting minimum and maximum for all LAS point record entries ...
  X          -1516000000 -1512000001
  Y          -884000000 -880000001
  Z              367267     631862
  intensity           0       4494
  return_number       0          7
  number_of_returns   0          7
  edge_of_flight_line 0          1
  scan_direction_flag 0          0
  classification      1         24
  scan_angle_rank   -31         32
  user_data           0          1
  point_source_ID    14        999
  gps_time -40787201.000000 263655956.894434
  Color R 0 0
        G 146 4614
        B 0 0
  attribute0          0      28.08  ('Amplitude')
  attribute1          0       34.5  ('Pulse width')
  attribute2     -18.54      26.14  ('Reflectance')
number of first returns:        10540061
number of intermediate returns: 1268951
number of last returns:         10542131
number of single returns:       8774677
WARNING: there are 633 points with return number 0
WARNING: there are 1272 points with return number 6
WARNING: there are 68 points with return number 7
overview over number of returns of given pulse: 8774044 1714674 1727941 1230537 120921 7246 470
WARNING: there are 633 points with a number of returns of given pulse of 0
histogram of classification of points:
           21376  unclassified (1)
         6062893  ground (2)
            2190  noise (7)
             632  keypoint (8)
            8005  water (9)
         4214965  overlap (12)
         3256456  Reserved for ASPRS Definition (20)
            9949  Reserved for ASPRS Definition (24)

E:\software\LAStools\bin>lasinfo -i 2621_5779_1_20200121.laz -copy_attribute_into_register 1 0 -scale_register 0 10.0 -copy_register_into_R 0 -copy_attribute_into_register 2
1 -translate_register 1 20 -scale_register 1 100.0 -copy_register_into_G 1
lasinfo (200619) report for '2621_5779_1_20200121.laz'
reporting all LAS header entries:
  file signature:             'LASF'
  file source ID:             0
  global_encoding:            1
  project ID GUID data 1-4:   00000000-0000-0000-0000-000000000000
  version major.minor:        1.2
  system identifier:          'LAStools (c) by rapidlasso GmbH'
  generating software:        'las2las (version 200509)'
  file creation day/year:     164/2020
  header size:                227
  offset to point data:       2385
  number var. length records: 5
  point data format:          3
  point data record length:   40
  number of point records:    13576466
  number of points by return: 10539428 1766996 909481 333208 25380
  scale factor x y z:         0.00025 0.00025 0.00025
  offset x y z:               1000000 6000000 0
  min x y z:                  621000.00000 5779000.00000 91.81675
  max x y z:                  621999.99975 5779999.99975 157.96550
variable length header record 1 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            34735
  length after header  208
  description          'GeoKeyDirectoryTag (mandatory)'
    GeoKeyDirectoryTag version 1.1.0 number of keys 25
      key 1024 tiff_tag_location 0 count 1 value_offset 1 - GTModelTypeGeoKey: ModelTypeProjected
      key 1025 tiff_tag_location 0 count 1 value_offset 2 - GTRasterTypeGeoKey: RasterPixelIsPoint
      key 1026 tiff_tag_location 34737 count 33 value_offset 0 - GTCitationGeoKey: ETRS89 / UTM zone 32N/ORTH_GCG16
      key 2048 tiff_tag_location 0 count 1 value_offset 32767 - GeographicTypeGeoKey: user-defined
      key 2049 tiff_tag_location 34737 count 7 value_offset 66 - GeogCitationGeoKey: ETRS89
      key 2050 tiff_tag_location 0 count 1 value_offset 32767 - GeogGeodeticDatumGeoKey: user-defined
      key 2051 tiff_tag_location 0 count 1 value_offset 32767 - GeogPrimeMeridianGeoKey: user-defined
      key 2052 tiff_tag_location 0 count 1 value_offset 9001 - GeogLinearUnitsGeoKey: Linear_Meter
      key 2054 tiff_tag_location 0 count 1 value_offset 9102 - GeogAngularUnitsGeoKey: Angular_Degree
      key 2056 tiff_tag_location 0 count 1 value_offset 32767 - GeogEllipsoidGeoKey: user-defined
      key 2057 tiff_tag_location 34736 count 1 value_offset 0 - GeogSemiMajorAxisGeoKey: 6378137
      key 2059 tiff_tag_location 34736 count 1 value_offset 1 - GeogInvFlatteningGeoKey: 298.2572221
      key 2060 tiff_tag_location 0 count 1 value_offset 9102 - GeogAzimuthUnitsGeoKey: Angular_Degree
      key 2061 tiff_tag_location 34736 count 1 value_offset 2 - GeogPrimeMeridianLongGeoKey: 0
      key 3072 tiff_tag_location 0 count 1 value_offset 32767 - ProjectedCSTypeGeoKey: user-defined
      key 3073 tiff_tag_location 34737 count 33 value_offset 33 - PCSCitationGeoKey: ETRS89 / UTM zone 32N/ORTH_GCG16
      key 3074 tiff_tag_location 0 count 1 value_offset 32767 - ProjectionGeoKey: user-defined
      key 3075 tiff_tag_location 0 count 1 value_offset 1 - ProjCoordTransGeoKey: CT_TransverseMercator
      key 3076 tiff_tag_location 0 count 1 value_offset 9001 - ProjLinearUnitsGeoKey: Linear_Meter
      key 3080 tiff_tag_location 34736 count 1 value_offset 4 - ProjNatOriginLongGeoKey: 9
      key 3081 tiff_tag_location 34736 count 1 value_offset 3 - ProjNatOriginLatGeoKey: 0
      key 3082 tiff_tag_location 34736 count 1 value_offset 5 - ProjFalseEastingGeoKey: 500000
      key 3083 tiff_tag_location 34736 count 1 value_offset 6 - ProjFalseNorthingGeoKey: 0
      key 3092 tiff_tag_location 34736 count 1 value_offset 7 - ProjScaleAtNatOriginGeoKey: 0.9996
      key 4099 tiff_tag_location 0 count 1 value_offset 9001 - VerticalUnitsGeoKey: Linear_Meter
variable length header record 2 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            34736
  length after header  64
  description          'GeoDoubleParamsTag (optional)'
    GeoDoubleParamsTag (number of doubles 8)
      6.37814e+006 298.257 0 0 9 500000 0 0.9996
variable length header record 3 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            34737
  length after header  73
  description          'GeoASCIIParamsTag (optional)'
    GeoAsciiParamsTag (number of characters 73)
      ETRS89 / UTM zone 32N/ORTH_GCG16|ETRS89 / UTM zone 32N/ORTH_GCG16|ETRS89|
variable length header record 4 of 5:
  reserved             0
  user ID              'LASF_Projection'
  record ID            2112
  length after header  965
  description          'OGC COORDINATE SYSTEM WKT'
    WKT OGC COORDINATE SYSTEM:
    COMPD_CS["ETRS89 / UTM zone 32N/ORTH_GCG16",PROJCS["ETRS89 / UTM zone 32N/ORTH_GCG16",GEOGCS["ETRS89 / UTM zone 32N/ORTH_GCG16",DATUM["ETRS89",SPHEROID["GRS 1980",6378137.000,298.257222101,AUTHORI
TY["EPSG","7019"]],AUTHORITY["EPSG","6258"]],PRIMEM["Greenwich",0.0000000000000000,AUTHORITY["EPSG","8901"]],UNIT["Degree",0.01745329251994329547,AUTHORITY["EPSG","9102"]],AUTHORITY["EPSG","25832"]],P
ROJECTION["Transverse_Mercator",AUTHORITY["EPSG","9807"]],PARAMETER["latitude_of_origin",0.0000000000000000],PARAMETER["central_meridian",9.0000000000000000],PARAMETER["scale_factor",0.999600000000000
0],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","2
5832"]],VERT_CS["Orthometric Heights",VERT_DATUM["GCG2016",2005,AUTHORITY["EPSG","0"]],UNIT["Meter",1.00000000000000000000,AUTHORITY["EPSG","9001"]],AXIS["Height",UP]]]
variable length header record 5 of 5:
  reserved             0
  user ID              'LASF_Spec'
  record ID            4
  length after header  576

  description          'RIEGL Extra Bytes'
    Extra Byte Descriptions
      data type: 3 (unsigned short), name "Amplitude", description: "Echo signal amplitude [dB]", min: 0, max: 10000, scale: 0.01, offset: 0 (not set)
      data type: 3 (unsigned short), name "Pulse width", description: "Full width at half maximum [ns]", min: 1, max: 10000, scale: 0.1, offset: 0 (not set)
      data type: 4 (short), name "Reflectance", description: "Echo signal reflectance [dB]", min: -5000, max: 15000, scale: 0.01, offset: 0 (not set)
the header is followed by 2 user-defined bytes
LASzip compression (version 3.4r3 c2 50000): POINT10 2 GPSTIME11 2 RGB12 2 BYTE 2
LAStiling (idx 1013, lvl 6, sub 0, bbox 590000 5751000 654000 5815000) (size 1000 x 1000, buffer 0)

reporting minimum and maximum for all LAS point record entries ...
  X          -1516000000 -1512000001
  Y          -884000000 -880000001
  Z              367267     631862
  intensity           0       4494
  return_number       0          7
  number_of_returns   0          7
  edge_of_flight_line 0          1
  scan_direction_flag 0          0
  classification      1         24
  scan_angle_rank   -31         32
  user_data           0          1
  point_source_ID    14        999
  gps_time -40787201.000000 263655956.894434
  Color R 0 345
        G 146 4614
        B 0 0
  attribute0          0      28.08  ('Amplitude')
  attribute1          0       34.5  ('Pulse width')
  attribute2     -18.54      26.14  ('Reflectance')
number of first returns:        10540061
number of intermediate returns: 1268951
number of last returns:         10542131
number of single returns:       8774677
WARNING: there are 633 points with return number 0
WARNING: there are 1272 points with return number 6
WARNING: there are 68 points with return number 7
overview over number of returns of given pulse: 8774044 1714674 1727941 1230537 120921 7246 470
WARNING: there are 633 points with a number of returns of given pulse of 0
histogram of classification of points:
           21376  unclassified (1)
         6062893  ground (2)
            2190  noise (7)
             632  keypoint (8)
            8005  water (9)
         4214965  overlap (12)
         3256456  Reserved for ASPRS Definition (20)
            9949  Reserved for ASPRS Definition (24)

--
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
---
You received this message because you are subscribed to the Google Groups "LAStools - efficient tools for LiDAR processing" group.
To unsubscribe from this group and stop receiving emails from it, send an email to lastools+u...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/lastools/829f00b2-289c-448d-a732-33256aa5be55o%40googlegroups.com.

Susana Gonzalez

unread,
Jun 18, 2020, 9:09:19 PM6/18/20
to last...@googlegroups.com

Excellent!! 😊

Thanks

 

From: last...@googlegroups.com <last...@googlegroups.com> On Behalf Of Martin Isenburg
Sent: Friday, 19 June 2020 1:03 PM
To: LAStools - efficient command line tools for LIDAR processing <last...@googlegroups.com>
Subject: Re: [LAStools] Copy extra byte values into R,G,B Attributes

 

[EXTERNAL]

Aicke Damrau

unread,
Jul 2, 2020, 8:21:12 AM7/2/20
to LAStools - efficient tools for LiDAR processing
HI Martin,

thank's a lot.

Scale Function is clear, but could you please explain the transform function?

Aicke
To unsubscribe from this group and stop receiving emails from it, send an email to last...@googlegroups.com.

Martin Isenburg

unread,
Jul 8, 2020, 12:34:21 AM7/8/20
to LAStools - efficient command line tools for LIDAR processing
Hello Aike,

sorry for the late reply ... but I do not understand the question. Please be more specific. 

Martin

Evon Silvia

unread,
Jul 10, 2020, 2:02:35 PM7/10/20
to last...@googlegroups.com
This looks like a Google Translate fail to me.

lasinfo -i 2621_5779_1_20200121.laz -copy_attribute_into_register 1 0 -copy_register_into_R 0  
lasinfo -i 2621_5779_1_20200121.laz -copy_attribute_into_register 1 0 -scale_register 0 10.0 -copy_register_into_R 0
lasinfo -i 2621_5779_1_20200121.laz -copy_attribute_into_register 2 1 -translate_register 1 20 -scale_register 1 10.0 -copy_register_into_G 1
lasinfo -i 2621_5779_1_20200121.laz -copy_attribute_into_register 2 1 -translate_register 1 20 -scale_register 1 100.0 -copy_register_into_G 1  

-scale_register will multiply the value by the specified amount.
-translate_register will add or offset the value by the specified amount - I'm guessing this got mistranslated/misunderstood as "transform."

This is similar to the LAS header's provision of a scale (multiply) and offset (add) value to the XYZ values, which are stored as integers.

You have to do this to the ExtraBytes because the RGB fields are unsigned 16-bit integers, while ExtraBytes can be anything at all. Storing an ExtraByte value like 3.1415926535 in an unsigned 16-bit integer would result in its truncation to the value of 3, which represents considerable data loss. Providing a scale factor of 100 will store the value of 314, and providing an offset lets you store negative numbers (for example) in a sensible way.

Hope that helps.

Evon
--

Evon Silvia PLS

QSI Solutions Architect

ASPRS LAS Working Group Chair

1100 NE Circle Blvd #126, Corvallis, OR 97330

D: 541-249-5818 | E: esi...@quantumspatial.com

     


--
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
---
You received this message because you are subscribed to the Google Groups "LAStools - efficient tools for LiDAR processing" group.
To unsubscribe from this group and stop receiving emails from it, send an email to lastools+u...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/lastools/CABSWR-FV1cFPeBBZArp59sCBdnSvYiE_v0S4sVU%2B%2Bgm7o6RDOA%40mail.gmail.com.

Martin Isenburg

unread,
Jul 10, 2020, 4:21:03 PM7/10/20
to LAStools - efficient command line tools for LIDAR processing
Some examples:

:: copy attribute 1 into R channel of RGB via register 0

las2las -i in.laz ^
            -copy_attribute_into_register 1 0 ^
            -copy_register_into_R 0  

:: copy attribute 1 into register 2, scale by factor 10.0, copy result into G channel of RGB

las2las -i in.laz ^
            -copy_attribute_into_register 1 2 ^
            -scale_register 2 10.0 ^
            -copy_register_into_G 2    

:: copy attribute 0 into register 3, translate by 20.0 (aka add 20), scale by factor 11.0, copy result into B channel of RGB

las2las -i in.laz ^
            -copy_attribute_into_register 0 3 ^
            -translate_register 3 20.0 ^
            -scale_register 3 11.0 ^
            -copy_register_into_B 3   

:: copy attribute 0 into register 4, scale by factor 5.0, translate by 2.5 (aka add 2.5), copy result into R channel of RGB

las2las -i in.laz ^
            -copy_attribute_into_register 0 4 ^
            -scale_register 4 5.0 ^
            -translate_register 4 2.5 ^
            -copy_register_into_R 4

More options to modify the registers will be implemented as needed. Let me know what you want.

Regards

Reply all
Reply to author
Forward
0 new messages