LASHeight Troubleshooting

LHall

unread,

Aug 26, 2016, 10:34:25 AM8/26/16

to LAStools - efficient tools for LiDAR processing

I am having trouble successfully running lasheight to convert ellipsoid elevations to geoid elevations. I have checked both input files to make sure that they are both in the same coordinate system.

Here is the error I keep getting: "ERROR: using zero ground points from 'C:\Users\........"

Here is the command I am running: lasheight -i "C:\Users\jimyd\Desktop\TEMP\8_25_16\tile_254150_3280300_1.las" -ground_points "C:\Users\jimyd\Desktop\TEMP\8_25_16\g2012au6_1_UTM_clip.las" -replace_z -olas

Anybody have any suggestions on how to correct this issue?

Martin Isenburg

unread,

Aug 26, 2016, 2:53:33 PM8/26/16

to LAStools - efficient command line tools for LIDAR processing

Hello,

it would help if you could include two lasinfo reports. You can generate them as follows:

lasinfo -cd -i C:\Users\jimyd\Desktop\TEMP\8_25_16\tile_254150_3280300_1.las -odix _info -otxt
lasinfo -cd -i C:\Users\jimyd\Desktop\TEMP\8_25_16\g2012au6_1_UTM_clip.las -odix _info -otxt

http://rapidlasso.com/lasinfo
http://lastools.org/download/lasinfo_README.txt

Regards,

Martin @rapidlasso

> --
> 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

LHall

unread,

Aug 26, 2016, 5:19:10 PM8/26/16

to LAStools - efficient tools for LiDAR processing

g2012au6_1_UTM_clip_info.txt:

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: ''

generating software: 'Global Mapper'

file creation day/year: 239/2016

header size: 227

offset to point data: 506

number var. length records: 3

point data format: 0

point data record length: 20

number of point records: 69540

number of points by return: 0 0 0 0 0

scale factor x y z: 0.0000001 0.0000001 0.001

offset x y z: -95.671838309564293 29.485984095198067 0

min x y z: -95.6718383 29.4859841 -27.390

max x y z: -95.3987505 29.6899017 -27.080

variable length header record 1 of 3:

reserved 43707

user ID 'LASF_Projection'

record ID 34735

length after header 80

description 'Georeferencing Information'

GeoKeyDirectoryTag version 1.1.0 number of keys 9

key 1024 tiff_tag_location 0 count 1 value_offset 2 - GTModelTypeGeoKey: ModelTypeGeographic

key 1026 tiff_tag_location 34737 count 13 value_offset 0 - GTCitationGeoKey: GCS_WGS_1984

key 2048 tiff_tag_location 0 count 1 value_offset 4326 - GeographicTypeGeoKey: GCS_WGS_84

key 2054 tiff_tag_location 0 count 1 value_offset 9102 - GeogAngularUnitsGeoKey: Angular_Degree

key 2056 tiff_tag_location 0 count 1 value_offset 7030 - GeogEllipsoidGeoKey: Ellipse_WGS_84

key 2057 tiff_tag_location 34736 count 1 value_offset 0 - GeogSemiMajorAxisGeoKey: 6378137

key 2058 tiff_tag_location 34736 count 1 value_offset 1 - GeogSemiMinorAxisGeoKey: 6356752.314

key 2059 tiff_tag_location 34736 count 1 value_offset 2 - GeogInvFlatteningGeoKey: 298.2572236

key 4099 tiff_tag_location 0 count 1 value_offset 9001 - VerticalUnitsGeoKey: Linear_Meter

variable length header record 2 of 3:

reserved 43707

user ID 'LASF_Projection'

record ID 34736

length after header 24

description 'Double Param Array'

GeoDoubleParamsTag (number of doubles 3)

6.37814e+006 6.35675e+006 298.257

variable length header record 3 of 3:

reserved 43707

user ID 'LASF_Projection'

record ID 34737

length after header 13

description 'GeoAsciiParamsTag'

GeoAsciiParamsTag (number of characters 13)

GCS_WGS_1984|

reporting minimum and maximum for all LAS point record entries ...

X 0 2730878

Y 0 2039176

Z -27390 -27080

intensity 0 0

return_number 1 1

number_of_returns 0 0

edge_of_flight_line 0 0

scan_direction_flag 0 0

classification 0 0

scan_angle_rank 0 0

user_data 0 0

point_source_ID 0 0

number of first returns: 69540

number of intermediate returns: 0

number of last returns: 69540

number of single returns: 69540

covered area in square units/kilounits: 4/0.00

point density: all returns 17385.00 last only 17385.00 (per square units)

spacing: all returns 0.01 last only 0.01 (in units)

WARNING: for return 1 real number of points by return is 69540 but header entry was not set.

WARNING: there are 69540 points with a number of returns of given pulse of 0

histogram of classification of points:

69540 never classified (0)

tile_254150_3280300_1_info.txt:

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.2

system identifier: 'LAStools (c) by rapidlasso GmbH'

generating software: 'las2las (version 160228)'

file creation day/year: 208/2016

header size: 227

offset to point data: 1082

number var. length records: 3

point data format: 1

point data record length: 28

number of point records: 9512530

number of points by return: 9432140 80390 0 0 0

scale factor x y z: 0.001 0.001 0.001

offset x y z: 254095.81465054984 3280392.2110245381 -0.585747686600944

min x y z: 254140.001 3280290.000 -7.275

max x y z: 254210.000 3280359.999 71.540

variable length header record 1 of 3:

reserved 43707

user ID 'LASF_Projection'

record ID 34735

length after header 64

description 'GeoTIFF GeoKeyDirectoryTag'

GeoKeyDirectoryTag version 1.1.0 number of keys 7

key 1024 tiff_tag_location 0 count 1 value_offset 1 - GTModelTypeGeoKey: ModelTypeProjected

key 1025 tiff_tag_location 0 count 1 value_offset 1 - GTRasterTypeGeoKey: RasterPixelIsArea

key 1026 tiff_tag_location 34737 count 22 value_offset 0 - GTCitationGeoKey: WGS 84 / UTM zone 15N

key 2049 tiff_tag_location 34737 count 7 value_offset 22 - GeogCitationGeoKey: WGS 84

key 2054 tiff_tag_location 0 count 1 value_offset 9102 - GeogAngularUnitsGeoKey: Angular_Degree

key 3072 tiff_tag_location 0 count 1 value_offset 32615 - ProjectedCSTypeGeoKey: WGS 84 / UTM 15N

key 3076 tiff_tag_location 0 count 1 value_offset 9001 - ProjLinearUnitsGeoKey: Linear_Meter

variable length header record 2 of 3:

reserved 43707

user ID 'LASF_Projection'

record ID 34737

length after header 30

description 'GeoTIFF GeoAsciiParamsTag'

GeoAsciiParamsTag (number of characters 30)

WGS 84 / UTM zone 15N|WGS 84|

variable length header record 3 of 3:

reserved 43707

user ID 'liblas'

record ID 2112

length after header 599

description 'OGR variant of OpenGIS WKT SRS'

LAStiling (idx 6, lvl 2, sub 0, bbox 254050 3.28025e+006 254250 3.28045e+006, buffer)

reporting minimum and maximum for all LAS point record entries ...

X 44186 114185

Y -102211 -32212

Z -6689 72126

intensity 0 25600

return_number 1 2

number_of_returns 1 2

edge_of_flight_line 0 0

scan_direction_flag 0 0

classification 0 0

scan_angle_rank 0 0

user_data 0 0

point_source_ID 1 1

gps_time 311576.078234 311896.193117

number of first returns: 9432140

number of intermediate returns: 0

number of last returns: 9429801

number of single returns: 9349411

covered area in square meters/kilometers: 4900/0.00

point density: all returns 1941.33 last only 1924.45 (per square meter)

spacing: all returns 0.02 last only 0.02 (in meters)

overview over number of returns of given pulse: 9349411 163119 0 0 0 0 0

histogram of classification of points:

9512530 never classified (0)

Martin Isenburg

unread,

Aug 27, 2016, 5:34:07 AM8/27/16

to LAStools - efficient command line tools for LIDAR processing

Hello,

well ... now it should be obvious. The files are not in the same CRS. The 'g2012au6_1_UTM_clip.laz' file is - despite its name - in geographic coordinates and not in projected UTM coordinates like the 'tile_254150_3280300_1.laz' file. Hence they have completely different bounding boxes and fail the initial overlap test already.