Я делаю анализ тренда для растра.Мне нужен этот график.
На графике анализ Манн-Кендалтренда годовых осадков ... Цвета показывают единицы оценки наклона Сена (мм у − 1) для осадков для каждой ячейки сетки (1 км регулярной сетки).Значимые тренды (α = 0,1) показаны в виде заштрихованных областей. Пока у меня есть следующий код.
##libraries
library(raster)
library(grDevices)
library(spatialEco)
##download the data for tile with lat long
datasets<-getData('worldclim', var='prec', res=10)
##cropping to some random smaller extent
e <- as(extent(55,65,30,35), 'SpatialPolygons')
crs(e) <- "+proj=longlat +datum=WGS84 +no_defs"
datasets <- crop(datasets, e)
names(datasets)<-month.abb
##plot the data to check
plot(datasets[[1]])
result<-raster.kendall(datasets, tau = TRUE, intercept = TRUE, p.value = TRUE,
z.value = TRUE, confidence = F)
names(result) <- c("slope","tau", "intercept", "p.value","z.value")
##make the significant p value shapefile
p.value <- result[[4]]
p.value[p.value >= 0.05] <- NA ##95%significance
p.value<-rasterToPolygons(p.value)
p.value <- aggregate(p.value)
##colour scale
Uniques <- cellStats(result[["slope"]],stat=unique)
Uniques.max <- max(Uniques)
Uniques.min <- min(Uniques)
my.at <- round(seq(ceiling(Uniques.max), floor(Uniques.min), length.out= 10),0)
myColorkey <- list(at=my.at, labels=list(at=my.at))
levelplot(result[["slope"]],margin=F,main="slope",par.settings=RdBuTheme(),colorkey=myColorkey,xlab=NULL,ylab=NULL,scales = list(draw=F,tck = 0))+
layer(sp.polygons(p.value,density=5,angle=45))
Это дает это.
Мне нужно проверить следующее.
Это правильный способ вычисления тренда?
Могу ли я интерпретировать этот наклонкак тренд?Могу ли я сказать, что тренд только для полигонов в левом нижнем углу значительный (<= 0,05).Даже при том, что более высокие значения уклона находятся в верхнем правом углу. </p>
- Как заполнить многоугольник наклонными штриховками?Я использовал плотность = 5, угол = 45, но это не график.