我的目标是使用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
处的材料创建的