constraining Error while using lasclip

79 views
Skip to first unread message

Basem Nouh

unread,
Mar 28, 2021, 1:26:47 PM3/28/21
to LAStools - efficient tools for LiDAR processing
Good evening,

i am trying to remove all building from Lidar Data and i have shape files from the City of Bochum (see the photo)
i tried to clip the building but i became always an error.
the code is :
lasclip -cpu64 -lof file_list.17032.txt -merged -poly "H:\Masterarbeit\Masterarbeit\Masterarbeit\HsBo\Geobasis_laz\prozessdatei1_fullarea\clipped_area\2Drop_class\3Unclass\3Unclass\shp\WallSurfaceLod2.shp" -interior -o "merged1.laz"
 
some times i became different constrains number like 2,3,7
i wish any one could give an advice cuz this step is the most important one to get the best DTM

best regards
Basem
3dm_32_377_5700_1_nw_1.txt
shapefiles_Bochum.JPG
Error_lasclip_shp.JPG

Martin Isenburg

unread,
Mar 29, 2021, 9:24:49 AM3/29/21
to LAStools - efficient command line tools for LIDAR processing
Hello,

two possibilities:

(1) Likely: The polygons in your SHP file are overlapping or their polygonal lines are crossing. Run some tool that cleans up 2D topologies.

(2) Less likely: After quantizing in a local (!) coordinate system from double to float (as done by lasclip) some line segments become overlapping/crossing. Run some tool that assures the some distances between all polygonal lines.

Regards,

Martin


--
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
---
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/dd525aea-d9d9-4ff0-8884-48f4d11dd92fn%40googlegroups.com.

Terje Mathisen

unread,
Mar 29, 2021, 12:43:51 PM3/29/21
to last...@googlegroups.com
Martin Isenburg wrote:
Hello,

two possibilities:

(1) Likely: The polygons in your SHP file are overlapping or their polygonal lines are crossing. Run some tool that cleans up 2D topologies.

That seems likely:

After merging several shape files I run the following command to clean up the geometry:

sprintf(q(c:\OSGeo4W64\bin\ogr2ogr "%s" "%s" -dialect sqlite -sql "select st_union(geometry) from '%s'"), $target, $temp_shape, $layer);

Here $temp_layer is the temporary file that results from several update/append operations, $target is the final name, and $layer is the name I've used for the shape file layer.

The key is the sqlite SQL/GIS operation that creates a minimum union of all the source polygons.

Terje

(2) Less likely: After quantizing in a local (!) coordinate system from double to float (as done by lasclip) some line segments become overlapping/crossing. Run some tool that assures the some distances between all polygonal lines.

Regards,

Martin


On Sun, Mar 28, 2021 at 7:26 PM Basem Nouh <basem...@gmail.com> wrote:
Good evening,

i am trying to remove all building from Lidar Data and i have shape files from the City of Bochum (see the photo)
i tried to clip the building but i became always an error.
the code is :
lasclip -cpu64 -lof file_list.17032.txt -merged -poly "H:\Masterarbeit\Masterarbeit\Masterarbeit\HsBo\Geobasis_laz\prozessdatei1_fullarea\clipped_area\2Drop_class\3Unclass\3Unclass\shp\WallSurfaceLod2.shp" -interior -o "merged1.laz"
 
To view this discussion on the web visit https://groups.google.com/d/msgid/lastools/CABSWR-F9UvFH%2ByKa%3DdJNJz2terak1A%3DE3m39YxPzMMBPh71nvg%40mail.gmail.com.


-- 
- <Terje.M...@tmsw.no>
"almost all programming can be viewed as an exercise in caching"

Basem Nouh

unread,
Mar 30, 2021, 7:20:18 AM3/30/21
to last...@googlegroups.com
Hallo,

which kind of tools should i run to clean up those 2D topologies. is those in Lastools!

or it is better to ask the Office in City of Bochum (Stadtmodelle der stadt Bochum, Amt für Geoinformation) to send me another shape file. Some shape files were not supported according to an error 'unsupported fileshape 31' 

regards,
Basem

Knoor, Claudia

unread,
Mar 30, 2021, 8:35:47 AM3/30/21
to last...@googlegroups.com

Hello,

 

I suspect that a citygml format (3D buildings) was converted to shape format.

But that doesn't work in your case, because the redundant citygml information was written multiple times in the shape files, so there are many overlays. Furthermore these shapefiles are probably 3D. I guess LAStools software requires 2D-shapes (not sure whether LAStools is able to work with 3D-shapes).

So it will be best to take the „BuildingLoD1Solid-shape“ and convert it to 2D, the other shapefiles are not necessary.

By doing that, clipping should work. But even then there could last some sparse overlays, e.g. in the case of patio roofs.

 

Regards

Claudia

Basem Nouh

unread,
Apr 28, 2021, 3:50:32 PM4/28/21
to last...@googlegroups.com
Hallo,

i have tried several solution to clean the Shapefile but always i have the same problem 
i tried through Arcgis to repair the Geometry, Multipatch footprint and in Qgis vector clean but Lasclip cannot read the file

which kind of tools i may use !

thanks
Basem

Reply all
Reply to author
Forward
0 new messages