Problem projecting LAS --> non-EPSG Albers to WA State Plane

[gmcw49...@gmail.com]

I have some LAS files that were delivered in an Albers Projection (NAD83-1983-USFS-R6-Abers, also SR-ORG:25, or US Forest Service Region 6 Albers)

I am Trying to re-project them to our standard of Washington State Plane south, NAD83.

 

I have not been able to find a standard EPSG code for this projection, and have not been able to figure out if I can specify a user-defined Albers projection parameter.  Trying to use las2las.exe.

Does Las2Las have this projection coded in anywhere?  Or can you clue me in to the syntax of specifying an Albers input projection parameters?

 

(I tried NOT specifying the input projection, but it ended up in the middle of the Pacific Ocean.)

Syntax tried --> .\las2las.exe -v -i 48117B3101.las -target_sp83 WA_S -target_feet -target_elevation_feet -olas -o las_files

Lasinfo reports for a sample input file, and for the resultant projected file are below.

Thanks!   Any help will be greatly appreciated!

Guy McWethy

 

Input File report from lasinfo ....

lasinfo report for N:\projects\public\colville_2015\original_data\All_Return_Point_Cloud\48117B3101.las
reporting all LAS header entries:
  file signature:             'LASF'
  file source ID:             0
  global_encoding:            1
  project ID GUID data 1-4:   084B9F75-F1A0-42C1-0000-000000000000
  version major.minor:        1.2
  system identifier:          'ALS70'
  generating software:        'GeoCue LAS Updater'
  file creation day/year:     250/2015
  header size:                227
  offset to point data:       1059
  number var. length records: 5
  point data format:          1
  point data record length:   28
  number of point records:    2182817
  number of points by return: 1512970 483528 147854 33531 4559
  scale factor x y z:         0.001 0.001 0.001
  offset x y z:               0 3000000 0
  min x y z:                  794981.808 1578195.843 793.423
  max x y z:                  795339.265 1578606.996 1776.882
variable length header record 1 of 5:
  reserved             43707
  user ID              'LASF_Projection'
  record ID            34735
  length after header  216
  description          'GeoTiff Projection Keys'
    GeoKeyDirectoryTag version 1.1.0 number of keys 26
      key 1024 tiff_tag_location 0 count 1 value_offset 1 - GTModelTypeGeoKey: ModelTypeProjected
      key 1026 tiff_tag_location 34737 count 35 value_offset 54 - GTCitationGeoKey: PCS Name = NAD_1983_USFS_R6_A
      key 2048 tiff_tag_location 0 count 1 value_offset 4269 - GeographicTypeGeoKey: GCS_NAD83
      key 2049 tiff_tag_location 34737 count 103 value_offset 116 - GeogCitationGeoKey: GCS Name = GCS_North_Amer
83|Datum = North_American_1983|Ellipsoid = GRS_1980|Primem = Greenwich
      key 2050 tiff_tag_location 0 count 1 value_offset 6269 - GeogGeodeticDatumGeoKey: Datum_North_American_Datu
      key 2051 tiff_tag_location 0 count 1 value_offset 8901 - GeogPrimeMeridianGeoKey: PM_Greenwich
      key 2054 tiff_tag_location 0 count 1 value_offset 9102 - GeogAngularUnitsGeoKey: Angular_Degree
      key 2055 tiff_tag_location 34736 count 1 value_offset 10 - GeogAngularUnitSizeGeoKey: 0.01745329252
      key 2056 tiff_tag_location 0 count 1 value_offset 7019 - GeogEllipsoidGeoKey: Ellipse_GRS_1980
      key 2057 tiff_tag_location 34736 count 1 value_offset 7 - GeogSemiMajorAxisGeoKey: 6378137
      key 2059 tiff_tag_location 34736 count 1 value_offset 8 - GeogInvFlatteningGeoKey: 298.2572221
      key 2061 tiff_tag_location 34736 count 1 value_offset 9 - GeogPrimeMeridianLongGeoKey: 0
      key 3072 tiff_tag_location 0 count 1 value_offset 32767 - ProjectedCSTypeGeoKey: user-defined
      key 3073 tiff_tag_location 34737 count 54 value_offset 0 - PCSCitationGeoKey: NAD_1983_USFS_R6_Albers|proje
Albers Equal Area
      key 3074 tiff_tag_location 0 count 1 value_offset 32767 - ProjectionGeoKey: user-defined
      key 3075 tiff_tag_location 0 count 1 value_offset 11 - ProjCoordTransGeoKey: CT_AlbersEqualArea
      key 3076 tiff_tag_location 0 count 1 value_offset 9001 - ProjLinearUnitsGeoKey: Linear_Meter
      key 3077 tiff_tag_location 34736 count 1 value_offset 6 - ProjLinearUnitSizeGeoKey: 1
      key 3078 tiff_tag_location 34736 count 1 value_offset 0 - ProjStdParallel1GeoKey: 43
      key 3079 tiff_tag_location 34736 count 1 value_offset 1 - ProjStdParallel2GeoKey: 48
      key 3081 tiff_tag_location 34736 count 1 value_offset 2 - ProjNatOriginLatGeoKey: 34
      key 3082 tiff_tag_location 34736 count 1 value_offset 3 - ProjFalseEastingGeoKey: 600000
      key 3083 tiff_tag_location 34736 count 1 value_offset 4 - ProjFalseNorthingGeoKey: 0
      key 3088 tiff_tag_location 34736 count 1 value_offset 5 - ProjCenterLongGeoKey: -120
      key 4097 tiff_tag_location 34737 count 27 value_offset 89 - VerticalCitationGeoKey: NAVD88 - Geoid12A (Mete
      key 4099 tiff_tag_location 0 count 1 value_offset 9001 - VerticalUnitsGeoKey: Linear_Meter
variable length header record 2 of 5:
  reserved             43707
  user ID              'LASF_Projection'
  record ID            34736
  length after header  88
  description          'GeoTiff double parameters'
    GeoDoubleParamsTag (number of doubles 11)
      43 48 34 600000 0 -120 1 6.37814e+006 298.257 0 0.0174533
variable length header record 3 of 5:
  reserved             43707
  user ID              'LASF_Projection'
  record ID            34737
  length after header  220
  description          'GeoTiff ASCII parameters'
    GeoAsciiParamsTag (number of characters 220)
      NAD_1983_USFS_R6_Albers|projection: Albers Equal Area|PCS Name = NAD_1983_USFS_R6_Albers|NAVD88 - Geoid12A
)|GCS Name = GCS_North_American_1983|Datum = North_American_1983|Ellipsoid = GRS_1980|Primem = Greenwich|
variable length header record 4 of 5:
  reserved             43707
  user ID              'NIIRS10'
  record ID            4
  length after header  10
  description          'NIIRS10 Timestamp'
variable length header record 5 of 5:
  reserved             43707
  user ID              'NIIRS10'
  record ID            1
  length after header  26
  description          'NIIRS10 Tile Index'
the header is followed by 2 user-defined bytes
reporting minimum and maximum for all LAS point record entries ...
  X           794981808  795339265
  Y          -1421804157 -1421393004
  Z              793423    1776882
  intensity           0        257
  return_number       1          7
  number_of_returns   1          7
  edge_of_flight_line 0          1
  scan_direction_flag 0          1
  classification      1         18
  scan_angle_rank   -21         22
  user_data           0          1
  point_source_ID    85         87
  gps_time 120405375.745355 120406049.706210
number of first returns:        1512970
number of intermediate returns: 187250
number of last returns:         1510088
number of single returns:       1027491
WARNING: there are 353 points with return number 6
WARNING: there are 22 points with return number 7
overview over number of returns of given pulse: 1027491 671347 343891 116695 21195 2039 159
histogram of classification of points:
         2038620  unclassified (1)
          144009  ground (2)
             183  noise (7)
               5  Reserved for ASPRS Definition (18)

PROJECTED FILE lasinfo Report ...

lasinfo report for N:\projects\public\colville_2015\las_files\48117B3101_rp.las
reporting all LAS header entries:
  file signature:             'LASF'
  file source ID:             0
  global_encoding:            1
  project ID GUID data 1-4:   084B9F75-F1A0-42C1-0000-000000000000
  version major.minor:        1.2
  system identifier:          'LAStools (c) by rapidlasso GmbH'
  generating software:        'las2las (version 160124)'
  file creation day/year:     250/2015
  header size:                227
  offset to point data:       467
  number var. length records: 3
  point data format:          1
  point data record length:   28
  number of point records:    2182817
  number of points by return: 1512970 483528 147854 33531 4559
  scale factor x y z:         0.001 0.001 0.001
  offset x y z:               0 3000000 0
  min x y z:                  794981.808 1578195.843 793.423
  max x y z:                  795339.265 1578606.996 1776.882
variable length header record 1 of 3:
  reserved             43707
  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 32149 - ProjectedCSTypeGeoKey: NAD83 / Washington South
      key 3076 tiff_tag_location 0 count 1 value_offset 9002 - ProjLinearUnitsGeoKey: Linear_Foot
      key 4099 tiff_tag_location 0 count 1 value_offset 9002 - VerticalUnitsGeoKey: Linear_Foot
variable length header record 2 of 3:
  reserved             43707
  user ID              'NIIRS10'
  record ID            1
  length after header  26
  description          'NIIRS10 Tile Index'
variable length header record 3 of 3:
  reserved             43707
  user ID              'NIIRS10'
  record ID            4
  length after header  10
  description          'NIIRS10 Timestamp'
the header is followed by 2 user-defined bytes
reporting minimum and maximum for all LAS point record entries ...
  X           794981808  795339265
  Y          -1421804157 -1421393004
  Z              793423    1776882
  intensity           0        257
  return_number       1          7
  number_of_returns   1          7
  edge_of_flight_line 0          1
  scan_direction_flag 0          1
  classification      1         18
  scan_angle_rank   -21         22
  user_data           0          1
  point_source_ID    85         87
  gps_time 120405375.745355 120406049.706210
number of first returns:        1512970
number of intermediate returns: 187250
number of last returns:         1510088
number of single returns:       1027491
WARNING: there are 353 points with return number 6
WARNING: there are 22 points with return number 7
overview over number of returns of given pulse: 1027491 671347 343891 116695 21195 2039 159
histogram of classification of points:
         2038620  unclassified (1)
          144009  ground (2)
             183  noise (7)
               5  Reserved for ASPRS Definition (18)