R语言学习记录：用矢量裁剪(掩膜)栅格

raster包和terra包都有crop和mask函数，前者生成shp文件范围的矩形，后者生成shp范围的栅格，范围之外都为NA。

用法为：crop(栅格，矢量shp）；mask(栅格，矢量shp）

今天发现了一个关于裁剪的问题：

我要用某个大的shp文件提取中间的某个县，然后用这个县的shp去裁剪栅格。

首先我想的是用vect函数读取大shp文件，然后再提取里面的县，比如这样

看起来挺好，但是，当我用这个函数去掩膜栅格的时候

结果并没有铺满整个屏幕，而是在原来的位置上。我猜测新掩膜后的栅格仍然保留着原来的行列数和范围。

那要如何解决这个问题呢？

很简单，先裁剪再掩膜。