Which tool(s) can I use to convert a LAS file into a binary grid (0 1)?

199 views
Skip to first unread message

Rosann Aarsen

unread,
Jul 10, 2015, 3:23:52 AM7/10/15
to last...@googlegroups.com
Hi,

I want to make a binary grid (0 1) of a point cloud, LAS file (see attachment).

I already looked at LASgrid, but I do not think that I can make a binary grid with this tool, right? Or is this possible by adding certain arguments?
 
My LAS file consists of a pedestrian path (with a part of the street) in the Netherlands, called Nieuwe Looiersstraat in Amsterdam, and covers an area of 33 x 13 meters (429 m2).
I want to detect in my point cloud data: where the width of the pedestrian path is < 900 mm (f.e. by using ‘opening’ with a circle as kernel of 900 mm width).
Therefore, first I need to have a binary grid to use 'opening'.
In the binary grid the values: 0 = no points are detected (as pedestrian path) ; 1 = points are detected in grid-cell. 
 
Hopefully you can help me, and recommend me which tool(s) I can use.

Thank in advance.
I am looking forward to see your response(s).
2. LASclip.las

Martin Isenburg

unread,
Jul 12, 2015, 7:35:37 AM7/12/15
to LAStools - efficient command line tools for LIDAR processing
Hello,

first off a couple of requests (for you but also in general):

1) If you send sample data to the list make sure it is small (max 500 KB)
2) Never send LAS files. Always convert them to LAZ first.
3) Do not use spaces or dots or other funny characters in your file names

You can use the lasgrid tool with the '-occupancy' option.

lasgrid -i "2. LASclip.las" -occupancy -step 2 -nodata 0 -o occu_2m.asc
set_epsg_code: look-up for 28992 not implemented
set_ProjectedCSTypeGeoKey: look-up for 28992 not implemented

more occu_2m.asc
ncols 16
nrows 5
xllcorner 121498.000000
yllcorner 485940.000000
cellsize 2.000000
NODATA_value 0
1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0
1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 0
0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1

lasgrid -i "2. LASclip.las" -occupancy -step 1 -nodata 0 -o occu_1m.asc
set_epsg_code: look-up for 28992 not implemented
set_ProjectedCSTypeGeoKey: look-up for 28992 not implemented

more occu_1m.asc
ncols 32
nrows 10
xllcorner 121498.000000
yllcorner 485940.000000
cellsize 1.000000
NODATA_value 0
0 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0
0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0
0 0 0 0 0 0 0 0 0 1 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 0 0 1 1


Reply all
Reply to author
Forward
0 new messages