我将shapefile裁剪到较小的范围（我的AOI - 感兴趣的区域） - 使用缩小的工作目录。在我的工作流程中，我光栅化shapefile。
这是我的问题：我保存了我的两个shapefile（越来越小）以比较栅格化结果（它应该与底层栅格具有相同范围（AOI）相同的较小shapefile（显然也是如此） AOI））。 但遗憾的是他们不是。 CRS和数字o单元是相同的 - 但是例如NA的数量不是。
我使用合成数据执行了相同的过程和工作流程，并且它工作得很好 - 所以问题必须是我的数据。这是dropboxlink，您可以在其中下载shapefile和栅格。 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)
你有什么建议可能出错吗？ 非常感谢！