Setting Geographic Coordinate Systems in las files

[Michael Perdue]

Hello,

What is the proper way to assign geographic coordinate systems to la[sz] files with lastools? I can assign epsg 4326 with the following command

las2las -epsg 4326 -i foo.laz -o bar.laz

But it fails when I try to assign epsg 4617.

las2las -epsg 4617 -i foo.laz -o bar.laz
EPSG code 4617 not found in 'pcs.csv' file
ERROR: unknown EPSG code in '-epsg 4617'

Neither 4326 or 4617 are in the pcs.csv file. Both are in the gcs.csv file.

$ grep 4326 /usr/local/LAStools/bin/serf/geo/gcs.csv
4326,WGS 84,6326,World Geodetic System 1984,6326,9122,7030,8901,1,0,6422,,0,,,,,,,,

$ grep 4617 /usr/local/LAStools/bin/serf/geo/gcs.csv
4617,"NAD83(CSRS)",6140,NAD83 Canadian Spatial Reference System,6140,9122,7019,8901,1,0,6422,1842,1,9603,0,0,0,,,,

Cheers,

Mike

[Martin Isenburg]

Hello,

for historic reasons only projected coordinate systems use the EPSG codes ... except for geographic WGS 84.

Try this instead:

E:\software\LAStools\bin>las2las -i fusa.laz -longlat -nad83_csrs -o mist.laz

E:\software\LAStools\bin>lasinfo mist.laz
lasinfo (210720) report for 'mist.laz'
reporting all LAS header entries:
  file signature:             'LASF'
  file source ID:             0
  global_encoding:            0
  project ID GUID data 1-4:   00000000-0000-0000-0000-000000000000
  version major.minor:        1.1
  system identifier:          'LAStools (c) by rapidlasso GmbH'
  generating software:        'las2las (version 210729)'
  file creation day/year:     40/2010
  header size:                227
  offset to point data:       367
  number var. length records: 2
  point data format:          1
  point data record length:   28
  number of point records:    277573
  number of points by return: 263413 13879 281 0 0
  scale factor x y z:         0.01 0.01 0.01
  offset x y z:               0 0 0
  min x y z:                  277750.00 6122250.00 42.21
  max x y z:                  277999.99 6122499.99 64.35
variable length header record 1 of 2:
  reserved             0
  user ID              'LASF_Projection'
  record ID            34735
  length after header  32
  description          'by LAStools of rapidlasso GmbH'
    GeoKeyDirectoryTag version 1.1.0 number of keys 3
      key 1024 tiff_tag_location 0 count 1 value_offset 2 - GTModelTypeGeoKey: ModelTypeGeographic
      key 2048 tiff_tag_location 0 count 1 value_offset 4617 - GeographicTypeGeoKey: GCS_NAD83_CSRS
      key 4099 tiff_tag_location 0 count 1 value_offset 9001 - VerticalUnitsGeoKey: Linear_Meter
variable length header record 2 of 2:

Regards,

Martin

--
Download LAStools at
http://lastools.org
http://rapidlasso.com
Be social with LAStools at
http://facebook.com/LAStools.Guy
http://twitter.com/LAStools
http://linkedin.com/groups/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/CAP-P_3LRYsQ2w977oVEk_4J3hGR2AgDyvyEmabMK674kL6pw-g%40mail.gmail.com.