В настоящее время я работаю над проектом по сравнению эвапотранспирации MODIS между классами земного покрова с течением времени. У меня есть набор изображений MODIS ET в виде 16-дневных композитов за 11-летний период и набор классификаций земного покрова. Я хочу создать график коробочных диаграмм, суммирующих ET для пикселей в каждом классе растительного покрова за год, который представляет классификация растительного покрова, с датой (с шагом 16 дней) по оси x, значением ET по оси y,и коробочные участки, сгруппированные / окрашенные по классу земельного покрова.
Я попытался переформатировать данные в фрейм данных и использовать ggplot, но создание графика занимает много времени, и я получаю ошибку. Теперь я пытаюсь использовать boxplot или bwplot с растровым стеком, но полученные графики либо отображают только значения ET для всего изображения (не сгруппированные по классу земного покрова) с датой на оси x, либо отображают сгруппированные значения ETпо классу земного покрова только для 1 16-дневного изображения / даты.
Ниже приведен код, над которым я сейчас работаю, с образцом из 5 16-дневных составных изображений MODIS и 1 классификации земного покрова с 3 классами. Все изображения имеют одинаковую проекцию, разрешение и размер, поэтому выравнивание не является проблемой. Я относительно новичок в R, поэтому любые советы о том, как лучше подойти к этому, были бы полезны. Мне нужно будет сделать несколько итераций этого графика, поэтому мне нужен сравнительно быстрый подход / код.
#import MODIS images
ET_files <- list.files(ET_path, full.names=TRUE, pattern=".tif$")
ET <- stack(ET_files)
LCclass <- raster(lc_path)
#define NA values
NAvalue(ET) <- 0
NAvalue(LCclass) <- 0
#plot boxplot
boxplot(ET, LCclass)
, который создает это изображение с именем каждого изображения MODIS (которое содержит дату) наось х, но не сгруппирована по классу земного покрова. Я не уверен, почему класс земного покрова игнорируется. изображение результирующего графика коробчатого графика