raster包和terra包都有crop和mask函數，前者生成shp文件范圍的矩形，后者生成shp范圍的柵格，范圍之外都為NA。

用法為：crop(柵格，矢量shp）；mask(柵格，矢量shp）

今天發(fā)現(xiàn)了一個關于裁剪的問題：

我要用某個大的shp文件提取中間的某個縣，然后用這個縣的shp去裁剪柵格。

首先我想的是用vect函數讀取大shp文件，然后再提取里面的縣，比如這樣

看起來挺好，但是，當我用這個函數去掩膜柵格的時候

結果并沒有鋪滿整個屏幕，而是在原來的位置上。我猜測新掩膜后的柵格仍然保留著原來的行列數和范圍。

那要如何解決這個問題呢？

很簡單，先裁剪再掩膜。