R 使用栅格裁剪栅格

## clip
 if(projection(thistif) != projection(mo.range)){
   thistif = projectRaster(from = thistif,to = mo.range)
 }
 
 r2 <- crop(thistif, extent(mo.range))
 #plot(r2)
 r3 <- mask(r2, mo.range)

 

r4 = setExtent(r3,ext = extent(mo.range))
r4

 

然后writeraster就可以了