On 04/05/2011 05:32 PM, gianni lavaredo wrote: > 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] :
Maybe you mean overlay(onePlot["Id"], p) ? Are you looking for coinciding points? In that case, there is also the function zerodist. Hth, > object of type 'closure' is not subsettable > > thanks for help or suggestions > Gianni > > [[alternative HTML version deleted]] > > _______________________________________________ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251 8333081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de http://www.52north.org/geostatistics e.pebe...@wwu.de _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo