Crop a Geospatial Raster in R
Today you will learn how to crop a raster into a vector. In this task, we will use the raster and rgdal packages. Let’s load them (read here). Then we load a raster (in the example we use Sentinel-2 satellite imagery) into the variable r, which we will crop (read here). We display the loaded raster (read here):
We load the polygon to be clipped into the variable border (read here) and display it on the raster:
We crop the raster with the polygon using the crop function from the raster package:
|
1 |
r2 <- crop(r,border) |
We display the result:
The raster has been cropped to the extent of the polygon only. To mask the pixels outside the polygon, you still need to use the mask function:
|
1 |
r3 <- mask(r2,border) |
Show the result:
As a result of our action, we have obtained a raster that is clipped exactly to our polygon.
