Olá pessoal, quero gerar um mapa em que a variável em questão é obtida por meio de uma função, na qual as variáveis de entrada são dados espacializados em um grid. Entretanto estou me deparando com este erro: Error in image.default(x, y, z) : increasing 'x' and 'y' values expected Vocês teriam uma solução? Segue o código exemplo Qp <- function(mu,sigma,area,p){ Qp<-(exp(mu+(sigma*qnorm(1-p))))*area } grid <- pred_grid(range(Longitude),range(Latitude), by=0.5) gr <- locations.inside(grid, limite) x <- gr$Var1 y <- gr$Var2 z <- Qp(kc.mu.ANO$predict,kc.sigma.ANO$predict,1,0.9) image(x,y,z) contour(z, add=T, nlev=5)
require(geoR) grid <- pred_grid(c(-55,-50),c(-25,-22), by=0.5) plot(grid) x <- unique(grid[,1]); nx <- length(x) y <- unique(grid[,2]); ny <- length(y) z <- matrix(runif(nx*ny, 10,50), nx, ny) image(x,y,z) Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Wagner, queira desculpar! Enviei o email anterior antes de finalizar. Basicamente faltou acrescentar que o z deve ser uma matriz e x e y devem corresponder às colunas e linhas da matriz. O erro mencionado é decorrente de haver valores fora de ordem no seu x (na verdade a sequência se repete). Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W] 2015-05-05 16:32 GMT-04:00 Éder Comunello <comunello.eder@gmail.com>:
require(geoR) grid <- pred_grid(c(-55,-50),c(-25,-22), by=0.5) plot(grid) x <- unique(grid[,1]); nx <- length(x) y <- unique(grid[,2]); ny <- length(y) z <- matrix(runif(nx*ny, 10,50), nx, ny) image(x,y,z)
Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Olá Éder obrigado novamente pela ajuda, entretanto, não consegui resolver o problema com essa solução, pois os parâmetros de entrada na função Qp são resultados do krig.conv() ou seja, são da classe "kriging" (...$predict), e me fornecem um vetor de tamanho diferente de nx*ny, sendo assim fiz um improviso: kc.sigma.ANO$z <- z image(kc.sigma.ANO,val=kc.sigma.ANO$z), contour(kc.sigma.ANO,val=kc.sigma.ANO$z, add=T, nlev=10) talvez não seja a forma mais adequada, mas resolvi assim. Em 5 de maio de 2015 17:58, Éder Comunello <comunello.eder@gmail.com> escreveu:
Wagner, queira desculpar! Enviei o email anterior antes de finalizar.
Basicamente faltou acrescentar que o z deve ser uma matriz e x e y devem corresponder às colunas e linhas da matriz. O erro mencionado é decorrente de haver valores fora de ordem no seu x (na verdade a sequência se repete).
Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
2015-05-05 16:32 GMT-04:00 Éder Comunello <comunello.eder@gmail.com>:
require(geoR) grid <- pred_grid(c(-55,-50),c(-25,-22), by=0.5) plot(grid) x <- unique(grid[,1]); nx <- length(x) y <- unique(grid[,2]); ny <- length(y) z <- matrix(runif(nx*ny, 10,50), nx, ny) image(x,y,z)
Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
participantes (2)
-
Wagner Wolff -
Éder Comunello