R 在上使用时出错(空间点、空间多边形)
我的目标是使用sp中的over()函数来确定点存在于哪个多边形中。我在使用时出现此错误(错误:在“>”中出现意外“>”):R 在上使用时出错(空间点、空间多边形),r,sp,point-in-polygon,R,Sp,Point In Polygon,我的目标是使用sp中的over()函数来确定点存在于哪个多边形中。我在使用时出现此错误(错误:在“>”中出现意外“>”): over(ptsUTM, hex_grid) 其中类(ptsUTM)=“空间点”,类(十六进制网格)=“空间多边形”。两者都有相同的@proj4string= CRS arguments: +proj=utm +zone=14 +datum=WGS84 +units=km +no_defs +ellps=WGS84 ptsUTM变量是根据中的材料创建的,添加了以下内容:
over(ptsUTM, hex_grid)
其中类(ptsUTM)=“空间点”,类(十六进制网格)=“空间多边形”。两者都有相同的@proj4string=
CRS arguments: +proj=utm +zone=14 +datum=WGS84 +units=km +no_defs +ellps=WGS84
ptsUTM变量是根据中的材料创建的,添加了以下内容:
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
rr2[rr2Hi,你能澄清最后两条指令应该做什么吗?(ptsUTM=SpatialPoints(pol2)
和ptsUTM=hex_grid@proj4string
)似乎您首先尝试将多边形转换为点,然后用proj4string覆盖结果。特别是最后一条指令似乎用字符串替换ptsUTM:因此,over
无法工作(除非我缺少某些内容)当然。如果不提供所有的代码,这是很难的:rr2是一个光栅对象;rasterToPolygons将rr2转换为空间多边形数据帧对象pol。然后,spTransform将CRS投影信息添加到pol2(过渡)。我使用空间点将pol2转换为空间点对象(ptsUTM)并添加与六角网格(是空间多边形的对象)具有相同CRS投影的投影信息.好的,但我的问题是如果执行ptsUTM=hex_grid@proj4string
在我看来,正如所写的,您似乎在覆盖变量,而不是将crs分配给它。在调用之前,您是否检查了ptsUTM
的类,超过了Mi scuzi,Lorenzo。我的错误。我被信息淹没了。谢谢。不用担心,很高兴能帮上忙。