使用over(spatialPoints,spatialPolygons)时出错

时间:2017-02-21 21:47:40

标签: r sp point-in-polygon

我的目标是使用sp中的over()函数来确定点存在于哪个多边形中。我在使用时遇到此错误(错误:意外'>'in“>”):

over(ptsUTM, hex_grid)

其中class(ptsUTM)=“SpatialPoints”和class(hex_grid)=“SpatialPolygons”。两者都有相同的@ proj4string =

CRS arguments: +proj=utm +zone=14 +datum=WGS84 +units=km +no_defs +ellps=WGS84

ptsUTM变量是根据rDGAL, Tiff Files, and WorldFile中的材料创建的,并添加了以下内容:

rr2[rr2<minDbzValue]  <- NA
pol <- rasterToPolygons(rr2, fun=function(x){ x >= minDbzValue } )
pol2                = spTransform( pol, "+proj=utm +zone=14 +datum=WGS84 +units=km" )
ptsUTM              = SpatialPoints( pol2 )
ptsUTM              = hex_grid@proj4string

和hex_grid是根据Fishnets and Honeycomb: Square vs. Hexagonal Spatial Grids

处的材料创建的

0 个答案:

没有答案