Erro no Qgis do aplicativo R
Jose Antonio Maior Bono
Carregando camada resultante
As seguintes camadas não foram geradas corretamente.
• C:/Temp/processing_gkgjtJ/1b7dcbd05f234ff78971c33b9c2c01bb/kriging_prediction.tif
• C:/Temp/processing_gkgjtJ/2233b1ba6fe84179960b381305a9f33f/kriging_variance.tif
Você
pode verificar o 'Painel de Registro de Mensagens' na janela principal
do QGIS para encontrar mais informações sobre a realização do algoritmo.
abaixo todo resultado do processamento.
Sera que alguém sabe o que poderia estar acontecendo?
Algoritmo 'Ordinary Kriging' iniciando…
Entrada de parâmetros:
{ 'Estimate_range_and_psill_initial_values_from_sample_variogram' : True, 'Extent' : 0, 'Local_kriging' : False, 'Number_of_nearest_observations' : 25, 'RPLOTS' : 'TEMPORARY_OUTPUT', 'R_CONSOLE_OUTPUT' : 'TEMPORARY_OUTPUT', 'Resolution' : 0, 'Show_Sum_of_Square_Errors' : False, 'field' : 'altura', 'kriging_prediction' : 'TEMPORARY_OUTPUT', 'kriging_variance' : 'TEMPORARY_OUTPUT', 'layer' : 'D:/BONO/Teste QGis/pontoslageado.xlsx.pontos.vrt', 'model' : 1, 'nugget' : 0, 'psill' : 0, 'range' : 0 }
R execution commands
options("repos"="http://cran.at.r-project.org/")
.libPaths("C:/Users/José Antônio/AppData/Roaming/QGIS/QGIS3/profiles/default/processing/rlibs")
tryCatch(find.package("sf"), error = function(e) install.packages("sf", dependencies=TRUE))
library("sf")
tryCatch(find.package("raster"), error = function(e) install.packages("raster", dependencies=TRUE))
library("raster")
tryCatch(find.package("'gstat'"), error = function(e) install.packages("'gstat'", dependencies=TRUE))
library("'gstat'")
tryCatch(find.package("'sp'"), error = function(e) install.packages("'sp'", dependencies=TRUE))
library("'sp'")
layer <- st_read("C:/Temp/processing_gkgjtJ/8329b99c20b1415c9038943d2a209b8c/layer.shp", quiet = TRUE, stringsAsFactors = FALSE)
field <- "altura"
Estimate_range_and_psill_initial_values_from_sample_variogram <- TRUE
nugget <- 0.0
model <- 1
range <- 0.0
psill <- 0.0
Local_kriging <- FALSE
Number_of_nearest_observations <- 25.0
Show_Sum_of_Square_Errors <- FALSE
Extent <- 0
Resolution <- 0.0
RPLOTS <- "C:/Temp/processing_gkgjtJ/5572f1f377694b6691eff94abf0d2cb4/RPLOTS.html"
R_CONSOLE_OUTPUT <- "C:/Temp/processing_gkgjtJ/9a189c739ba247b9800ca185a382ccff/R_CONSOLE_OUTPUT.html"
png("C:/Temp/processing_gkgjtJ/62bb9fab47de42359f43650e07a558f1/RPLOTS.png")
library('gstat')
library('sp')
Models<-c("Exp","Sph","Gau","Mat")
model2<-Models[model+1]
create_new_data_ch <- function (layer)
{
convex_hull = chull(coordinates(layer)[, 1], coordinates(layer)[,2])
convex_hull = c(convex_hull, convex_hull[1])
d = Polygon(layer[convex_hull, ])
if(!is.projected(layer) | Resolution== 0){new_data = spsample(d, 5000,
type = "regular")}
if(is.projected(layer) & Resolution!= 0){
new_data = spsample(d, n= 1, cellsize=c(Resolution,Resolution),
type="regular")}
gridded(new_data) = TRUE
attr(new_data, "proj4string") <-layer@proj4string
return(new_data)
}
create_new_data_ext <- function (layer){
bottomright <- c(layer@bbox[1], layer@bbox[2])
topleft <- c(layer@bbox[3], layer@bbox[4])
d <- SpatialPolygons(
list(Polygons(list(Polygon(coords = matrix(
c(topleft[1],bottomright[1], bottomright[1],topleft[1],topleft[1],
topleft[2], topleft[2], bottomright[2],
bottomright[2],topleft[2]), ncol=2, nrow= 5))), ID=1)))
if(!is.projected(layer) | Resolution== 0){new_data = spsample(d, 5000,
type = "regular")}
if(is.projected(layer) & Resolution != 0){
new_data = spsample(d, n= 1, cellsize=c(Resolution,Resolution),
type="regular")}
gridded(new_data) = TRUE
attr(new_data, "proj4string") <-layer@proj4string
return(new_data)
}
if(Extent==0){mask<-create_new_data_ch(layer)}
if(Extent==1){mask<-create_new_data_ext(layer)}
field <- make.names(field)
names(layer)[names(layer)==field]="field"
layer$field <- as.numeric(as.character(layer$field))
str(layer)
layer <- remove.duplicates(layer)
layer <- layer[!is.na(layer$field),]
g = gstat(id = field, formula = field~1, data = layer)
vg = variogram(g)
if(Estimate_range_and_psill_initial_values_from_sample_variogram){range=NA}
if(Estimate_range_and_psill_initial_values_from_sample_variogram){psill=NA}
vgm = vgm(nugget=nugget, psill=psill, range=range, model=model2)
vgm = fit.variogram(vg, vgm)
vgm
plot(vg, vgm, plot.numbers = TRUE)
if(Local_kriging==FALSE){prediction = krige(field~1, layer, newdata = mask, vgm)}
if(Local_kriging==TRUE){prediction = krige(field~1, layer, newdata = mask, vgm, nmax=Number_of_nearest_observations)}
if(Show_Sum_of_Square_Errors==TRUE){paste("SSE:", attr(vgm, "SSErr"))}
#>if(!is.projected(layer)){warning(paste0("'layer' isn't projected.\n", "Resolution was not used. Interpolation was done over 5000 cells"))}
#>if(is.projected(layer) & Resolution == 0){warning("Resolution was set to 0. Final resolution estimated from data")}
kriging_prediction = raster(prediction)
kriging_variance = raster(prediction["var1.var"])
writeRaster(kriging_variance, "C:/Temp/processing_gkgjtJ/2233b1ba6fe84179960b381305a9f33f/kriging_variance.tif", overwrite = TRUE)
writeRaster(kriging_prediction, "C:/Temp/processing_gkgjtJ/1b7dcbd05f234ff78971c33b9c2c01bb/kriging_prediction.tif", overwrite = TRUE)
dev.off()
R execution console output
[1] "C:/Program Files/R/R-4.0.3/library/sf"
Erro: package or namespace load failed for 'sf' in loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck = vI[[j]]):
there is no package called 'DBI'
Execução interrompida
Tarefa executada em 0.51 segundos
Resultados:
{'RPLOTS': 'C:/Temp/processing_gkgjtJ/5ca1633e9af64532b5ee3339a5aac9be/RPLOTS.html',
'R_CONSOLE_OUTPUT': 'C:/Temp/processing_gkgjtJ/29ba42751c3841eeae0cb61dbf11558b/R_CONSOLE_OUTPUT.html',
'kriging_prediction': 'C:/Temp/processing_gkgjtJ/1b7dcbd05f234ff78971c33b9c2c01bb/kriging_prediction.tif',
'kriging_variance': 'C:/Temp/processing_gkgjtJ/2233b1ba6fe84179960b381305a9f33f/kriging_variance.tif'}
Carregando camada resultante
As seguintes camadas não foram geradas corretamente.
• C:/Temp/processing_gkgjtJ/1b7dcbd05f234ff78971c33b9c2c01bb/kriging_prediction.tif
• C:/Temp/processing_gkgjtJ/2233b1ba6fe84179960b381305a9f33f/kriging_variance.tif
Você pode verificar o 'Painel de Registro de Mensagens' na janela principal do QGIS para encontrar mais informações sobre a realização do algoritmo.
Hebert Azevedo Consultor Geo
[1] "C:/Program Files/R/R-4.0.3/library/sf"
Erro: package or namespace load failed for 'sf' in loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck = vI[[j]]):
there is no package called 'DBI'
Execução interrompida
no console do R instalar o DBI
install.packages("DBI")e no script incluir a biblioteca
library("DBI")
Não sei se vai resolver, mas vale testar.
Att--
Comunidade QGIS Brasil: http://qgisbrasil.org/
---
You received this message because you are subscribed to the Google Groups "qgisbrasil" group.
To unsubscribe from this group and stop receiving emails from it, send an email to qgisbrasil+...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/qgisbrasil/ead4f9f3-02fb-4534-82d6-22defa1dfcb1o%40googlegroups.com.
|
Remetente notificado por Mailtrack 29/01/21 17:37:56 |
maio...@gmail.com
Obrigado Herbert!!
Vou tentar, se der certo te aviso!!
Att
José Antonio
To view this discussion on the web visit https://groups.google.com/d/msgid/qgisbrasil/CACuh6VrYxMCzxxAa9uiFWxtQ%3DedVnHbq1Nx9avBsm-Wt6moUZg%40mail.gmail.com.