Dear Researchers,

I have two "SpatialPointsDataFrame" layer "p" and "onePLOT"  and I wish to
do an overlay, but I have an error message

> str(p)
Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots
  ..@ data       :'data.frame': 980000 obs. of  1 variable:
  .. ..$ V3: num [1:980000] 1 1 1 1 1 1 1 1 1 1 ...
  ..@ coords.nrs : int [1:2] 1 2
  ..@ coords     : num [1:980000, 1:2] 0.5 1.5 2.5 3.5 4.5 5.5 6.5 7.5 8.5
9.5 ...
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : NULL
  .. .. ..$ : chr [1:2] "x" "y"
  ..@ bbox       : num [1:2, 1:2] 0.5 0.5 1399.5 699.5
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : chr [1:2] "x" "y"
  .. .. ..$ : chr [1:2] "min" "max"
  ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots
  .. .. ..@ projargs: chr NA


> str(onePLOT)
Formal class 'SpatialPolygonsDataFrame' [package "sp"] with 5 slots
  ..@ data       :'data.frame': 1 obs. of  2 variables:
  .. ..$ Id      : Factor w/ 98 levels "LIDAR_Sicily_Area1_Plot1",..: 1
  .. ..$ firecode: num 1
  ..@ polygons   :List of 1
  .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots
  .. .. .. ..@ Polygons :List of 1
  .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots
  .. .. .. .. .. .. ..@ labpt  : num [1:2] 59 648
  .. .. .. .. .. .. ..@ area   : num 900
  .. .. .. .. .. .. ..@ hole   : logi FALSE
  .. .. .. .. .. .. ..@ ringDir: int 1
  .. .. .. .. .. .. ..@ coords : num [1:5, 1:2] 39.9 50 78.2 68.1 39.9 ...
  .. .. .. .. .. .. .. ..- attr(*, "dimnames")=List of 2
  .. .. .. .. .. .. .. .. ..$ : NULL
  .. .. .. .. .. .. .. .. ..$ : chr [1:2] "x" "y"
  .. .. .. ..@ plotOrder: int 1
  .. .. .. ..@ labpt    : num [1:2] 59 648
  .. .. .. ..@ ID       : chr "1"
  .. .. .. ..@ area     : num 900
  ..@ plotOrder  : int 1
  ..@ bbox       : num [1:2, 1:2] 39.9 628.7 78.2 667.1
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : chr [1:2] "x" "y"
  .. .. ..$ : chr [1:2] "min" "max"
  ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots
  .. .. ..@ projargs: chr NA


when I try overlay I have this error message

> p.ov <- overlay[onePLOT["Id"],p]
Error in overlay[onePLOT["Id"], p] :
  object of type 'closure' is not subsettable

thanks for help or suggestions
Gianni

	[[alternative HTML version deleted]]

