Как извлечь значения отражательной способности из оптического стека SAR в R
Я поместил изображения часового-1 и стража-2 в QGIS. Я выполнил предварительную обработку в SNAP и использовал инструмент "Выровнять растры" в QGIS, чтобы повторно сэмплировать изображения до 20 м и зарегистрировать их в общей проекции. Я хочу выполнить классификацию в R, но для извлечения значений отражательной способности / обратного рассеяния требуется вечность.
Прочитав мой шейп-файл и мой растр в 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)
Я ожидаю, что последняя строка будет извлекать значения отражательной способности, используя шейп-файл для создания кадра данных