栅格数据包中提取功能的使用不一致

时间:2020-02-04 12:15:31

标签: r r-raster

使用提取函数对协变量值从1分钟分辨率重新采样到5分钟分辨率时,我一直遇到奇怪的行为。我想为一个大区域建模数据,该区域的响应范围为5分钟乘5分钟的网格系统(这些单元称为pentads),但感兴趣的气候协变量的分辨率更高,为1分钟乘1分钟。我一直在使用提取函数来获取均值,总和和计数,但是均值与计数和总和不匹配。

此外,我还遇到了一些其他不一致的情况,具体取决于返回的对象类型,即数据框对象与空间多边形数据框。下面的输出显示了这两种对象类型的前10个均值,计数和总和,很明显,有时它们是相同的,但有时是不同的,并且不是系统的:

covs.mean.df [1:10,2]

[1] 215.2727 231.4000 221.1000 222.8333 387.8750 381.7857 390.6667 346.6000 322.3684 [10] 319.8000

covs.mean.sp@data [1:10,12]

[1] 381.7857 231.4000 215.2727 221.1000 222.8333 387.8750 390.6667 346.6000 322.3684 [10] 319.8000

covs.sum.df [1:10,2]

[1] 1157 2368 2211 1337 3103 5345 8204 6932 6125 7995

covs.sum.sp@data [1:10,12]

[1] 1157 2368 2211 1337 3103 6932 5345 6125 7995 6204

covs.count.sp@data [1:10,12]

[1] 10 5 11 6 8 14 21 20 19 25

covs.count.df [1:10,2]

[1] 5 11 10 6 8 14 21 20 19 25

几天前,我重新运行了代码,唯一的不同是我显式地插入了参数“ cellnumbers = F”,现在看起来似乎更加一致,这当然很奇怪,因为cellnumbers参数的默认值设置为FALSE。因此,现在大多数行是一致的(因为均值=总和/计数),但是仍然有一些行的总和和计数与均值不一致。任何见解/建议将不胜感激。

非常感谢

格雷格

0 个答案:

没有答案