使用自己的 crs 进行投影时，projectRaster 会产生不同的输出

Question

对于一个函数，我想使用用户集 CRS 重新投影光栅输入 - 本身是使用裁剪和蒙版后的输出。我认为用现有的 crs 投影不会做任何事情，只是 return 输入栅格。令我惊讶的是，情况并非如此。下面是一个可重现的例子

创建虚拟栅格：

library(raster)

# Download country polygin, in this case Malawi
mwi <- raster::getData("GADM", country = "MWI", level = 0)

# Create dummy raster
grid <- raster::raster() # 1 degree raster
grid <- raster::disaggregate(grid, fact = 12) # 5 arcmin
grid <- raster::crop(grid, mwi)
values(grid) <- rep(1, ncell(grid)) # Set a value

# The input raster with dimensions 94,39,3666
grid <- raster::mask(grid, mwi)
plot(grid)
grid

# Reproject the raster using its own crs. I use ngb as it is a categorical variable.
# This raster has dimensions 102, 47, 4794 so it seems a lot of white space (NA) is added.
own_crs <- crs(grid)
grid_reproj <- raster::projectRaster(grid, crs = own_crs, method = "ngb")
plot(grid_reproj)
grid_reproj

# To remove the white space I use trim
# This results in a waster with dimensions 93, 39, 3627
grid_trim <- raster::trim(grid_reproj)
plot(grid_trim)
grid_trim

# I also decided to compare the maps visually with mapview
library(mapview)

# There seems to be a trim function in mapview which I set to FALSE
# Also use the browser for easy viewing
options(viewer = NULL)
mapviewOptions(trim = FALSE)
mapviewGetOption("trim")
mapview(grid, col.regions = "green", legend = F) +
  mapview(grid_reproj, col.regions = "red", legend = F) +
  mapview(grid_trim, col.regions = "blue", legend = F)

比较地图我观察到两件事：

(1) 网格和 grid_trim 除了单个网格单元外几乎完全相同在上面。也许这是由于四舍五入，

(2) grid_reproj 具有更大的维度和不同的范围。它也是与其他地图相比，地图似乎略有偏移。这个由 trim 校正，所以我假设这些实际上是 NA 细胞，差异可能是与 mapview 如何显示地图有关。

因此，我的主要问题是，当 rasterProject 投影时发生了什么一样的程度？为什么结果会有所不同，甚至在 trim 之后？

Answer 1

要以正确方式投影光栅数据，您必须提供另一个光栅* 对象，而不仅仅是一个 crs。

如果您提供 crs，则会计算出一个新范围，它比避免丢失数据所需的预期范围大一点。当然，在这种 crs 实际上相同的奇怪情况下，return 输入不变是有意义的。

使用自己的 crs 进行投影时，projectRaster 会产生不同的输出

projectRaster results in different output when projecting with its own crs

gis

r

raster

map-projections