Я скопировал изображения часового-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)
Я ожидаю, что последняя строка извлечет значения отражательной способности, используя шейп-файл для создания кадра данных