Skip to main content
4 of 6
added 2802 characters in body

lasclipCircle with catalog

I'm trying to clip a catalog using in package lidR (v2.0.1), I've tried with lasclipCircle and lasclip. It gives me the error: no point in ROI, but when I'm viewing the shapefile using mapview (I've exported the X and Y coordinates from this shapefile) and the catalog they are at the same place.

Here's my code.

page1 <- catalog("O:/LiDAR/NB_LiDAR_16/page1")
Buffer_16 <- readOGR(dsn = "O:/2018_Regen/Derived GIS", layer = "Buffer_16")
X_16 <- Buffer_16@data[["DDLONG"]]
Y_16 <- Buffer_16@data[["DDLAT"]]
page1_points <- lasclipCircle(page1, xcenter = X_16,ycenter = Y_16, radius = 10)

or

page1_points <- lasclip(page1, Buffer_16, radius = 10)

I get the following error

Error : No point found for within region of interest 2.

Summary of Buffer_16 Object of class SpatialPolygonsDataFrame Coordinates: min max x 2457220 2613228 y 7490508 7662516 Is projected: TRUE proj4string : [+init=epsg:2953 +proj=sterea +lat_0=46.5 +lon_0=-66.5 +k=0.999912 +x_0=2500000 +y_0=7500000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs] Data attributes: OBJECTID MAPNO PLOTID X Y DDLAT DDLONG SCHEDULE Min. : 62.0 4721 : 11 23239 : 1 2541224: 11 7600512: 9 Min. :46.41 Min. :-67.07 1:179 1st Qu.:257.5 4921 : 11 26157 : 1 2539224: 10 7554512: 8 1st Qu.:47.09 1st Qu.:-66.46 Median :454.0 5230 : 10 26159 : 1 2543224: 9 7604512: 8 Median :47.35 Median :-66.14 Mean :415.1 4425 : 9 26160 : 1 2505224: 8 7646512: 8 Mean :47.34 Mean :-66.12 3rd Qu.:578.5 4714 : 9 26369 : 1 2501224: 7 7602512: 7 3rd Qu.:47.61 3rd Qu.:-65.88 Max. :790.0 5523 : 9 26370 : 1 2507224: 7 7648512: 7 Max. :47.96 Max. :-65.00 (Other):120 (Other):173 (Other):127 (Other):132 ISPSP OBJECTID_1 objectid_2 PHOTO_YEAR LABEL NUMBERING Shape_Leng Shape_Area 0:131 Min. :5 Min. :2 2016:179 2016:179 0:179 Min. :749732 Min. :1.876e+10 1: 48 1st Qu.:5 1st Qu.:2 1st Qu.:749732 1st Qu.:1.876e+10 Median :5 Median :2 Median :749732 Median :1.876e+10 Mean :5 Mean :2 Mean :749732 Mean :1.876e+10 3rd Qu.:5 3rd Qu.:2 3rd Qu.:749732 3rd Qu.:1.876e+10 Max. :5 Max. :2 Max. :749732 Max. :1.876e+10

Summary of page1 : `class : LAScatalog extent : 2457000 , 2590000 , 7542000 , 7647000 (xmin, xmax, ymin, ymax) coord. ref. : +init=epsg:2953 +proj=sterea +lat_0=46.5 +lon_0=-66.5 +k=0.999912 +x_0=2500000 +y_0=7500000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs area : 43 km² points : 844.56 million points density : 19.6 points/m² num. files : 43 Summary of the processing options:

  • Catalog will be processed by file respecting the original tiling pattern
  • Catalog will be processed using 1 core(s). Summary of the output options:
  • Outputs will be returned in R objects. Summary of the input options:
  • readLAS will be called with the following select option: *
  • readLAS will be called with the following filter option: `