Mais

Salvando resultados no pacote R do automap para dados de série temporal

Salvando resultados no pacote R do automap para dados de série temporal


Estou planejando interpolar os dados mensais de precipitação de 1980 e 2010. Estou usando o seguinte código:

require (gstat) require (sp) require (automap) require (rgdal) require (raster) para (i em 1: ncol (data @ data)) {kriging_new <- autoKrige (data [, i], data, grid) plot (kriging_new)}

Quero salvar os valores previstos (estimativas de krigagem), variância da previsão de krigagem, desvio padrão de krigagem para cadaeue salve-o como um arquivo img. Alguém pode me ajudar a configurar isso dentro do loop.


Você pode usar rgdal ou raster (com uma etapa adicional de coersão) para exportar a previsão, a variância da previsão ou o desvio padrão da previsão como rasters. As somas dos quadrados são armazenadas como um vetor onde os modelos de ajuste e variograma experimental são objetos data.frame. Você terá que atribuir uma coluna e gerá-los em um único data.frame ou escrever cada um independentemente usando algo como write.csv.

Aqui está um exemplo prático que grava os resultados no diretório de trabalho atual.

library (sp) library (automap) library (rgdal) library (raster) setwd ("D: / TMP") dados (meuse) coordenadas (meuse) = ~ x + y data (meuse.grid) em grade (meuse.grid) = ~ x + y vars <- c ("cádmio", "cobre", "chumbo", "zinco") sum.squares <- vector () var.model <- data.frame () para (i em vars) {kriging_new <- autoKrige (meuse @ data [, i] ~ 1, meuse, meuse.grid) sum.squares <- append (sum.squares, kriging_new $ sserr) kriging_new $ var_model <- data.frame (y = i, kriging_new $ var_model) var.model <- rbind (var.model, kriging_new $ var_model) writeGDAL (kriging_new $ krige_output ["var1.pred"], paste (paste (i, "pred", sep = "_"), " img ", sep =". ")) writeGDAL (kriging_new $ krige_output [" var1.var "], paste (paste (i," var ", sep =" _ ")," img ", sep =". ") )} nomes do var.model (sum.squares) <- vars print (sum.squares) r <- raster ("cadmium_pred.img") plot (r)

Assista o vídeo: General Automap Concepts Novation SL MkII Explained