Skip to main content
Geographic Information Systems

Return to Question

Tweeted twitter.com/StackGIS/status/1103445295172587520
deleted 1 character in body; edited tags; edited title
Source Link
Andre Silva
Andre Silva
  • 10.3k
  • 12
  • 55
  • 109

lasclipCircle with Clipping a catalog with lasclipCircle returns "no point found within region of interest"

I'm trying to clip a catalog using in package lidR (v2.0.1), package. 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.

I get the following error:

I got the same message for ROI 1 and all the other ROI as well.

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.

I get the following error

Clipping a catalog with lasclipCircle returns "no point found within region of interest"

I'm trying to clip a catalog using lidR (v2.0.1) package. 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.

I get the following error:

I got the same message for ROI 1 and all the other ROI as well.

Use code block
Source Link
edit approved Mar 1, 2019 at 17:37
JRR
edit approved Mar 1, 2019 at 17:37
JRR
  • 9.7k
  • 1
  • 14
  • 29
page1_points <- lasclip(page1, Buffer_16, radius = 10)

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+10Buffer_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 optionspage1:

  • 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: `
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:    `  

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

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: `
page1_points <- lasclip(page1, Buffer_16)

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:    `
added 2802 characters in body
Source Link
Stéphanie Landry
Stéphanie Landry
  • 33
  • 3

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: `

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: `
Format code, fix code, format Q
Source Link
edit approved Feb 28, 2019 at 8:43
JRR
edit approved Feb 28, 2019 at 8:43
JRR
  • 9.7k
  • 1
  • 14
  • 29
added 28 characters in body
Source Link
Bera
Bera
  • 77.8k
  • 14
  • 78
  • 188
Source Link
Stéphanie Landry
Stéphanie Landry
  • 33
  • 3