2

我正在使用卫星跟踪的动物并拥有大量的搬迁数据。

所以我有我的地图数据和重定位作为 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)
4

1 回答 1

2

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

鉴于评论中的proj4string(map), proj4string(locs),proj4string(SpatialPoints(locs))特别proj4string(SpatialPoints(locs))是NA, 我会说当它更改为对象时会count.points删除它。我认为在保留投影字符串的同时强制 a 的方法是通过.... 尝试使用将其插入?proj4stringSpatialPointsSPDFSPas(x,"SpatialPoints")tracecount.points

于 2013-05-13T11:01:44.193 回答