Question

我正在使用卫星跟踪的动物并且有大量的重新定位数据。

所以我将地图数据和重定位作为SpatialPointsDataFrames，当我问

时

if proj4string(map)==proj4string(locs)我知道了。

但是当我尝试count.points函数时如下

cp <- count.points(locs, map)

我收到以下错误

Error in count.points(SpatialPoints(x), w) : 
  different proj4string in w and xy

有没有人对为什么会出现这种情况有任何想法？

修改代码：

load("mydata") 
map = mydata$map 
map 
mimage(map) 
locs= mydata$relocs 
locs 
image(map) 
points(locs, col=as.numeric(slot(locs, "data")[,1]), pch=16) 
cp <- count.points(locs, map)

Answer 1

可重复的例子在这里会有很长的路要走。但一般而言，R对投影字符串的比较大致是逐字的。因此，如果有额外的空间等，它就会失败。

考虑到评论中的proj4string(map)，proj4string(locs)，proj4string(SpatialPoints(locs))，特别是proj4string(SpatialPoints(locs))是NA，我会说{{1}当count.points更改为proj4string对象时，它会删除它。我认为在保留投影字符串的同时强制SpatialPoints到SPDF的方法是通过SP ....尝试使用as(x,"SpatialPoints")将其插入trace ？

为什么R认为投影数据不同？

1 个答案: