我在QGIS中堆叠了哨兵1和哨兵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提取反射率值以生成数据框