When the raster package was still active, I would usually have to do the following to crop the rasterlayer to an irregular polygon (with some functions updated for terra). I am wondering there is a faster/more efficient way to do this same thing?:

kml <- vect("pathtokml.kml")
r <- rast("pathtoraster.tif")

crp <- mask(r, project(kml, r)) #fix for different CRS's
crp2 <- crop(crp, ext(kml))
plot(crp2)

When the raster package was still active, I would usually have to do the following to crop the rasterlayer to an irregular polygon (with some functions updated for terra). I am wondering there is a faster/more efficient way to do this same thing?:

kml <- vect("pathtokml.kml")
r <- rast("pathtoraster.tif")

crp <- mask(r, project(kml, r)) #fix for different CRS's
crp2 <- crop(crp, ext(kml))
plot(crp2)

• r
• geospatial
• r-raster
• terra

Share Improve this question Follow asked Nov 19, 2024 at 4:07 AndrewAndrew 17122 silver badges66 bronze badges 1

• 1 For some performance gain try cropping first and and then mask. Or jus set mask = TRUE in crop() - stackoverflow/questions/77998525/… – margusl Commented Nov 19, 2024 at 8:02

Add a comment  | 

1 Answer 1

Sorted by: Reset to default 2

It is more efficient to first crop and then mask, and with terra you can do this in one step.

x <- crop(r, project(kml, r), mask=TRUE)