我将 shapefile 裁剪到较小的程度（我的 AOI - 感兴趣的区域） - 以使用减少的工作目录。在我的工作流程中，我栅格化了 shapefile。
这是我的问题：我保存了我的两个 shapefile（更小和更大）以比较栅格化结果（这应该与底层栅格具有与较小 shapefile 相同的范围（AOI）相同（显然也是我的 AOI）） . 但不幸的是，他们不是。CRS 和 o 个单元格是相同的 - 但例如 NA 的数量不是。
我对合成数据做了同样的程序和工作流程，而且效果很好——所以问题可能出在我的数据上。这是您可以下载 shapefile 和栅格的 dropboxlink。https://www.dropbox.com/sh/btgt2rc7uzawtx5/AADJ2YrKOnPh8gM-PPF7rmIQa?dl=0
我把我的代码留在这里：
#load shp files
setwd("C:/Users/.../R")
TESTshp<-readOGR(dsn="test_crop_dropbox", layer="boden_ebod_reproj")
extent(TESTshp)
setwd("C:/Users/.../R/test_crop_dropbox")
TESTraster<-raster("testraster.tif")
extent(TESTraster)
TESTshp_small <- crop(TESTshp, extent(TESTraster))
TESTrasterize<- rasterize(TESTshp, TESTraster, field="BodTyp_gen")
TESTrasterize_small<- rasterize(TESTshp_small, TESTraster, field="BodTyp_gen")
identical(TESTrasterize, TESTrasterize_small)
你有什么建议可能是错的吗？非常感谢！