Я вычисление тренда временного ряда растровых данных в отношении времени, но из-NA и нулевых значений я не получаю правильные результаты ..Расчета склона в R после удаления ноля или НС значения
В примере ниже значений дается пиксель. Из-за нулевых значений я получаю значение наклона, как -23.275 вместо 3.897
time Value
1 289
2 289
3 353
4 305
5 0
6 0
7 385
8 0
9 0
10 305
11 0
12 0
13 0
14 0
15 0
correct form should be like this
Time Value
1 289
2 289
3 353
4 305
7 385
10 305
мой пример кода
library(raster)
setwd("F:\\PHD_work\\Paper_work\\EOS\\New folder")
eos = stack(list.files(pattern='*.img'))
time = 1:nlayers(eos)
f_slope=function(x) { if (is.na(x[1])){ NA } else { m = lm(x ~ time); summary(m)$coefficients[2] }}
z = calc(eos, f_slope)
y = writeRaster(z, filename = "slope", format = "HFA")
так что вы хотели бы иметь вектор с ненулевыми значениями? – mabdrabo
, пожалуйста, поделитесь выходом 'str (eos)' – mabdrabo
, который не смог опубликовать результаты всей str (eos), потому что он очень длинный. –