我是一名小学生，试图传授我的论文。
我有一个非常大的土地覆盖栅格（所有纽约市，分辨率为3x3英尺）。
为了测试我的脚本，我想使用mask函数来提取栅格的一部分，例如，Bronx或其中的一部分。
为此，我有一个具有低级别管理区域的shapefile，例如社区或人口普查区域。导入它们之后，我可以看到两个图层在空间上是如何重叠的（使用bbox和plot）。但是，当我尝试做面具时，需要很长时间才能看到任何结果。
在ArcMap上尝试此操作时，相同的操作需要很短的时间。
所执行的操作是：
#importing the different layers
land_cover <- raster(land_cover.img)
Bronx <- readOGR(filename)
#to make sure they have the same CRS, i change the projection of Bronx
Bronx <- spTransform(Bronx, crs(land_cover)
land_cover_bronx <- mask (land_cover, Bronx)
这是最后一行需要很长时间才能处理，并让我怀疑某些东西运行不好。每层的投影会有问题吗？ [compareCRS返回TRUE，图表显示它们是verlap]。