如何从陆地卫星图像中提取值并绘制R中不同物质的光谱轮廓

我有一个区域的形状文件，我想得到该区域不同物质的光谱轮廓，例如水、道路或植被。像这样的东西

最好的方法是什么

我已经在r中加载了shapefile，但不确定如何从Landsat光栅堆栈中提取值

```{r}
landsat<-stack(B1, B2, B3, B4, B5, B6, B7)
roi <- readOGR("roi.shp")
shapes
```

`{r}
陆地卫星使用以下代码

df = raster::extract(landsat,     # raster layer
    roi,                       #shapefile
    fun=mean,         # what value to extract
    df=TRUE)          # return a dataframe

光栅对象和roi应该具有相同的坐标系