Daily Archives: March 30, 2018

    Improvement of clip by extent example from Wednesday…  The following code clips the LAS to the boundary of the shapefile.  (example data) library(lidR) library(rgdal) setwd("c:/temp/sf/") tdat<- readLAS("sp2_pcloud.las") plot(tdat) tpoly<- readOGR( ".","AnalysisUnit") clpdat <- lasclipPolygon(tdat, tpoly@polygons[[1]]@Polygons[[1]]@coords[,1], tpoly@polygons[[1]]@Polygons[[1]]@coords[,2]) plot(clpdat)     COARSE WOODY DEBRIS EXAMPLE: Task at hand is to determine […]