亲爱的人群

问题

我试图用包nfc,pgirmess,SpatialPack和spdep来计算空间相关图.但是,我很难定义距离的起点和终点.我只对较小距离的空间自相关感兴趣,但有较小的箱子.此外,由于光栅非常大(1.8百万像素),我遇到了这些软件包的内存问题,但是SpatialPack.

所以我尝试使用包栅格中的Moran函数生成我自己的代码.但我必须有一些错误,因为完整数据集的结果与其他包中的结果略有不同.如果我的代码中没有错误,它至少可以帮助其他类似问题的人.

题

我不确定,我的焦点矩阵是否错误.你能否告诉我中心像素是否需要合并？使用testdata我无法显示方法之间的差异,但在我的完整数据集中,可见差异,如下图所示.但是,箱子并不完全相同(50米对69米),所以这可以解释部分差异.然而,在第一个箱子里,这个解释似乎对我来说似乎不合理.或者我的光栅的不规则形状,以及处理NA的不同方式会导致差异吗？

自己的方法与SpatialPack的比较

可运行的例子

测试数据

计算测试数据的代码取自http://www.petrkeil.com/?p=1050#comment-416317

# packages used for the data generation
library(raster)
library(vegan) # will be used for PCNM

# empty matrix and spatial coordinates of its cells
side=30
my.mat <- matrix(NA, nrow=side, ncol=side)
x.coord <- rep(1:side, each=side)*5
y.coord <- rep(1:side, times=side)*5
xy <- data.frame(x.coord, y.coord)

# all paiwise euclidean distances between the cells
xy.dist <- dist(xy)

# PCNM axes of the dist. matrix (from 'vegan' package)
pcnm.axes <- pcnm(xy.dist)$vectors

# using 8th …