仅在多边形区域内剪切栅格

Question

我正在尝试使用leaflet根据大量多边形覆盖剪切的栅格。但是，我使用的光栅非常大并且不会映射，更不用说重新投影了。我有一种感觉，我实际上并没有“剪裁”＃39;完全是光栅。任何有关减少栅格大小或根据多边形实际剪切栅格的帮助都将不胜感激。

project_shp <- gUnion(spgeom1 = area1,spgeom2 = area2)

cr<- crop(raster, extent(project_shp), snap = "out")
fr<- rasterize(project_shp, cr)
lr<-mask(x=cr,mask=fr)

lr<-projectRaster(lr,crs=("+init=EPSG:4326"))

  

错误：无法分配大小为1.1 Gb的矢量

leaflet() %>% addTiles() %>%
    addPolygons(data= provincial,weight = 2, opacity = 0.8)%>%
    addRasterImage(data= lr, colors = lr$FBP_Fuel_Type1)

  

错误：无法分配大小为1.1 Gb的矢量