光栅：：焦点返回不正确的值

Question

我正在使用光栅软件包v2.5-8的焦距函数来获得3x3窗口中的最大值。我期望两行/列的边以NA的形式返回，相反，返回的输出是9,9,9。这是正确的吗？

示例:

library(raster); require(rasterVis)
r <- raster(nrows=3, ncols=3)
r[] <- 1:ncell(r)
plot(r);text(r);
r.class <- focal(r, w=matrix(1,nrow=3,ncol=3), fun=max) 
plot(r.class); text(r.class);

输出：

     [,1] [,2] [,3]
[1,]   NA   NA   NA
[2,]    9    9    9
[3,]   NA   NA   NA

预期输出：

       [,1] [,2] [,3]
[1,]   NA   NA   NA
[2,]   NA    9   NA
[3,]   NA   NA   NA

Answer 1

得到这个结果是因为地球的“左”和“右”边(经度= -180或180)是相同的位置。

library(raster)
r <- raster(nrows=3, ncols=3)
r[] <- 1:ncell(r)
as.matrix(r)

##     [,1] [,2] [,3]
## [1,]    1    2    3
## [2,]    4    5    6
## [3,]    7    8    9

rf <- focal(r, w=matrix(1,nrow=3,ncol=3), fun=max) 
as.matrix(rf)

##     [,1] [,2] [,3]
## [1,]   NA   NA   NA
## [2,]    9    9    9
## [3,]   NA   NA   NA

默认的CRS是lonlat。

crs(r)

## CRS arguments:
## +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 

有了平面CRS，你就能得到你所期望的：

crs(r) <- "+proj=utm +zone=1 +datum=WGS84"
rf2 <- focal(r, w=matrix(1,nrow=3,ncol=3), fun=max) 
as.matrix(rf2)

##     [,1] [,2] [,3]
## [1,]   NA   NA   NA
## [2,]   NA    9   NA
## [3,]   NA   NA   NA