Discussion: Karttapullautin help

I need some help with Karttapullautin.

I have a set of lidar files that cause the error below. One of these files can be downloaded here: https://giovanniberlanda.name/5h680051100.laz

Is there any suggestion on what I could look at?

C:\Temp1\temp_pullauta>pullauta.exe in\5h680051100.laz
Karttapullautin v. 20190203 (c) Jarkko Ryyppo 2012-19  All rights reserved.
This exe is free for non commercial use or if used for for navsport
mapping (orienteering, rogaining, adventure racing mapping).
There is no warranty. Use it at your own risk!
Preparing input file
.......... done.
Knoll detection part 1
.. done.
Knoll detection part 2
..Illegal modulus zero at script/pullautin.pl line 6163.

Contour generation part 1
Modification of non-creatable array value attempted, subscript -660477 at script/pullautin.pl line 6985.

Contour generation part 2
.................... done.
Contour generation part 3
.................... done.
Contour generation part 4
.................... done.
Vegetation generation
......Modification of non-creatable array value attempted, subscript -680000 at script/pullautin.pl line 3713,  line 1.

Cliff generation
.....GD Warning: one parameter to a memory allocation multiplication is negative or zero, failing operation gracefully
gdImageCreate error at GD/Image.pm line 86.

Rendering png map with depressions
...Modification of non-creatable array value attempted, subscript -660500 at script/pullautin.pl line 8400.

Rendering png map without depressions
...Modification of non-creatable array value attempted, subscript -660500 at script/pullautin.pl line 8400.


All done!

D:\testfolder>lasinfo.exe 5h680051100.laz
lasinfo (191127) report for '5h680051100.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.0
  system identifier:          ''
  generating software:        'TerraScan'
  file creation day/year:     203/2011
  header size:                227
  offset to point data:       291
  number var. length records: 1
  point data format:          1
  point data record length:   28
  number of point records:    3500049
  number of points by return: 2138557 965123 323797 72572 0
  scale factor x y z:         0.01 0.01 0.01
  offset x y z:               0 0 0
  min x y z:                  680000.00 5110000.01 1039.29
  max x y z:                  680500.00 5110500.00 1247.20
variable length header record 1 of 1:
  reserved             0
  user ID              'ESRI_Export'
  record ID            3
  length after header  8
  description          'Local Time: 2020 04 30 11 44 13'
the header is followed by 2 user-defined bytes
LASzip compression (version 3.4r1 c2 50000): POINT10 2 GPSTIME11 2
reporting minimum and maximum for all LAS point record entries ...
  X            68000000   68050000
  Y           511000001  511050000
  Z              103929     124720
  intensity           1       4820
  return_number       1          4
  number_of_returns   1          4
  edge_of_flight_line 0          0
  scan_direction_flag 0          1
  classification      1          1
  scan_angle_rank   -17          6
  user_data          32         32
  point_source_ID     1          2
  gps_time 292203.833043 473345.298273
number of first returns:        2138557
number of intermediate returns: 397985
number of last returns:         2132735
number of single returns:       1169228
overview over number of returns of given pulse: 1169228 1282873 755937 292011 0 0 0
histogram of classification of points:
         3500049  unclassified (1)

Point could needs to be ground classified, these files are not, all of the points are classified as "unclassififed".

I took a different route with the LiDAR file and came to the same conclusion.
The 3.5 million points are 'unclassified'. Usually they are identified as ground points and others. The ground points give you the contours, the hill shading, etc.
So, working the LiDAR in OCAD I got the ground contours all right but also contour lines for each of the trees and the buildings.
So I think we have identified the problem. I'm sorry I do not know what is the solution.

I have not tried it but there is:

https://lastools.github.io/download/lasground_new_...

Yup, try lasground. Let us know how it works out!

to make it seamless between tiles you will need to re-tile dataset with something like 50 m buffer. You can do it with lastile's option -buffer 50. If I remember right lasground will automatically drop the buffer.

https://rapidlasso.com/2015/08/07/use-buffers-when...

(two first commands there are about what you need, first making tiles with a buffer and then ground classification)

I use lasground_new, even when the data is already classified. I choose the -wilderness option to reduce the suppression of upright objects the program thinks are buildings. More mounds and boulders appear.

Thanks everybody for the support, lasground_new seems to do what I need.
I'm trying to complete this project ( https://mapanttrentino.org/ ) and the files that cover one of the two big holes in the middle are different that everything else.

That is one huge project! I'm glad it will work out for you.
So far how many LiDAR tiles you have used in this project?

There are about 35000 tiles. We are processing them in groups of 12x12 and cropping to 10x10 after extracting contours, vegetation, and adding OSM data. Then we use QGIS to merge everything and generate the tiles for the web.

Images seem to have only about 32 colors or so (I think you are using Nearest Neighbour), but image format is 32bit png. You could drop file sizes and (bandwidth usage) down to 30% by converting them to 8bit png images. But I guess thats something you can do when all tiles are done.

You missed a gate.

Thanks for the suggestion about png file format!
I'm not sure if the QGIS web-tile-generator plugin uses the same format as its input from Karttapullautin or it has a relevant option on its own. Or I'll just add an extra step and convert them after generation.

> You missed a gate.
Can you please fix it in OpenStreetMap ?

A post step with imagemagick would be better, lets you iterate compression settings if you overwrite tiles.

I can remeber 5 days there especially the event I ran back loop wrong direction and and ended up drinking at a water station three times with always about 20 min in between. I wonder what they thougt, a crazy hopeless Finn having drink breaks while searchng a control I guess.