r - 如何从 R 中的 SAR 光学堆栈中提取反射率值

Question

我在 QGIS 中堆叠了 sentinel-1 和 sentinel-2 图像。我在 SNAP 中进行了预处理，并使用 QGIS 中的“对齐栅格”工具将图像重新采样到 20m 并将它们注册到一个公共投影。我想在 R 中执行分类，但提取反射率/反向散射值需要很长时间。

在将我的 shapefile 和我的光栅读入 R 之后，我陷入了反射率值提取步骤。

training <- readOGR('C:/march2018/TrainingData.shp')

s12 <- brick('C:/march2018/s1s213bstack.tif')

names(s12) <- paste0("B", c(1:13))

s12ref <- extract(s12, training, df=TRUE)

我希望最后一行使用 shapefile 提取反射率值以生成数据帧