library(sf) library(sp) library(raster) shp <- read_sf("E:/waterconstraint/LoessPlateauRegion/LoessPlateauRegion.shp") r<-raster("E:/waterconstraint/GIMMS_NDVI/MonthMax/1982-01-01.tif") crs(shp) crs(r) shp_proj <- st_transform(shp,crs = (r) #r_crop <- crop(r,shp_proj) plot(r) plot(r_crop,main = "crop")

这段代码的作用是引入R语言中的sf、sp和raster三个包，然后读取一个shp文件和一个.tif文件（分别位于E:/waterconstraint/LoessPlateauRegion/ 和 E:/waterconstraint/GIMMS_NDVI/MonthMax/ 这两个文件夹下），并获取它们对应的坐标系信息。最后，将读取的shp文件转换为与.tif文件相同的投影坐标系（即广义选择不仅仅是地理坐标系）。