HortonMachine Spatial Toolbox - Lesto

7 views
Skip to first unread message

Afshin Ghafoori

unread,
Dec 12, 2022, 4:14:47 AM12/12/22
to jgrasstools
Hello,

New to GIS and the functions of gvSIG. Is there anyone who could help och share a guide on how till initiate the PointCloudMaximaFinder in Lesto?

Have filtered the pointcloud in order to get the vegetation seperated in FME. DEM surface (as TIFF) is exported and imported into a layer (so far so good).

When i start the The HortonMachine Spatial Toolbox and choose Lesto - vegetation - PointCloudMaximaFinder, the prarameters are set to:

The input LAS: Las (of the vegetation (original points from DEM))
A dtm raster is set to the imported DEM surface earlier
A set of polygons to use as region of interest is set to epmty
An optional dsm-dtm difference raster to use to check on the extracted tops is set to epmty
Radius for which a point can be local maxima is set to 4
Use an adaptive radius based on height is uncheked
Max permitted elevation difference around maxima is set to 3
The elevation threshold to apply to the chm is set to 0
The comma seperated list of classes to filter is set to 0
The output local maxima is set to Test.shp

Load experimental is checked
The heap is set to 24000

The error message from the log:

Warning, not writing an empty vector to file: 1.shp
No points processed in tile: POLYGON ((189377.6488098003 6375176.06587686, 189377.6488098003 6375778.06587686, 190150.6488098003 6375778.06587686, 190150.6488098003 6375176.06587686, 189377.6488098003 6375176.06587686))

Any help would be nice!

/Afshin

andrea antonello

unread,
Dec 12, 2022, 5:10:12 AM12/12/22
to jgras...@googlegroups.com
Hello Afshin,

> New to GIS and the functions of gvSIG. Is there anyone who could help och share a guide on how till initiate the PointCloudMaximaFinder in Lesto?

sure

> Have filtered the pointcloud in order to get the vegetation seperated in FME. DEM surface (as TIFF) is exported and imported into a layer (so far so good).
>
> When i start the The HortonMachine Spatial Toolbox and choose Lesto - vegetation - PointCloudMaximaFinder, the prarameters are set to:
>
> The input LAS: Las (of the vegetation (original points from DEM))
> A dtm raster is set to the imported DEM surface earlier
> A set of polygons to use as region of interest is set to epmty

ok

> An optional dsm-dtm difference raster to use to check on the extracted tops is set to epmty

ok

> Radius for which a point can be local maxima is set to 4

ok, this would be in the units of the las projection.

> Use an adaptive radius based on height is unchecked

ok

> Max permitted elevation difference around maxima is set to 3

ok

> The elevation threshold to apply to the chm is set to 0

ok

> The comma seperated list of classes to filter is set to 0

this means that only the points with class 0 are taken. If you want to
use all the points, leave it empty.

> The output local maxima is set to Test.shp

ok

> Load experimental is checked
> The heap is set to 24000
>
> The error message from the log:
>
> Warning, not writing an empty vector to file: 1.shp
> No points processed in tile: POLYGON ((189377.6488098003 6375176.06587686, 189377.6488098003 6375778.06587686, 190150.6488098003 6375778.06587686, 190150.6488098003 6375176.06587686, 189377.6488098003 6375176.06587686))

Ok, not knowing the data I have a few guesses:

- the dtm and las are in different projections: in that case the
search areas would not overlap
- there are no point with class 0

What is currently happening, is that for some reason no point is taken
from the input las to calculate the maxima.
We just need to understand why.

> Any help would be nice!

Let us know if the above guesses apply. Otherwise we need to investigate deeper.

Cheers,
Andrea


>
> /Afshin
>
> --
> You received this message because you are subscribed to the Google Groups "jgrasstools" group.
> To unsubscribe from this group and stop receiving emails from it, send an email to jgrasstools...@googlegroups.com.
> To view this discussion on the web visit https://groups.google.com/d/msgid/jgrasstools/99fac343-8d2d-43fe-bcf7-827430b4ac2en%40googlegroups.com.

Afshin Ghafoori

unread,
Dec 12, 2022, 6:25:21 AM12/12/22
to jgrasstools
Checked the coordinates of the las-file which seemed a bit of. Changed the  "The comma seperated list of classes to filter" to empty and set  "The elevation threshold to apply to the chm to 3. Works very good!

Thank you for quick help and good work!

/Afshin

andrea antonello

unread,
Dec 12, 2022, 8:22:28 AM12/12/22
to jgras...@googlegroups.com
> Checked the coordinates of the las-file which seemed a bit of. Changed the "The comma seperated list of classes to filter" to empty and set "The elevation threshold to apply to the chm to 3. Works very good!
>
> Thank you for quick help and good work!

Great to hear!

Cheers,
Andrea
> To view this discussion on the web visit https://groups.google.com/d/msgid/jgrasstools/55fbaf78-a2c2-44c2-ab07-bdbd367c00adn%40googlegroups.com.
Reply all
Reply to author
Forward
0 new messages