Krigagem: Obter valores pontuais a partir da sup. interpolada
Bom dia pessoal, Para código reproduzível desta dúvida pode-se utilizar o exemplo disponível na página dos tutoriais do geoR: http://leg.ufpr.br/geoR/tutorials/ctc.R Ao final deste tutorial chegamos ao mapa da distribuição espacial obtido por krigagem. [image: Imagem inline 1] Gostaria de saber se é possível, obter o valor da variável para coordenadas específicas por exemplo na figura acima: x=5; y=8 x=25; y=10 x=15; y=15 Se é possível, como posso implementar no R. Desde já agradeço toda colaboração, *Hélder Gramacho * Recife-PE / *agrohelder@gmail.com <agrohelder@hotmail.com>*
defina estas coordenadas como pontos a serem interpolados e roda a função novamente: myLOCS <- cbind(c(5, 25, 15), c(8, 10, 15)) krige.conv(..., loc=myLOCS, ...) On Mon, 20 Oct 2014, Helder Gramacho wrote:
Bom dia pessoal, Para código reproduzível desta dúvida pode-se utilizar o exemplo disponível na página dos tutoriais do geoR: http://leg.ufpr.br/geoR/tutorials/ctc.R Ao final deste tutorial chegamos ao mapa da distribuição espacial obtido por krigagem. Imagem inline 1 Gostaria de saber se é possível, obter o valor da variável para coordenadas específicas por exemplo na figura acima: x=5; y=8 x=25; y=10 x=15; y=15 Se é possível, como posso implementar no R.
Desde já agradeço toda colaboração,
Hélder Gramacho Recife-PE / agrohelder@gmail.com
Olá Prof. Paulo Justiniano, Funcionou por aqui, muito obrigado pela ajuda! Outra dúvida: É possível exportar o mapa da krigagem para abrir num SIG por exemplo? Um dos formatos que o SIG que utilizo abre é o formato ASCII grided XYZ. Obrigado, *Hélder Gramacho * Recife-PE / *agrohelder@gmail.com <agrohelder@hotmail.com>* Em 20 de outubro de 2014 10:06, Paulo Justiniano <paulojus@leg.ufpr.br> escreveu:
defina estas coordenadas como pontos a serem interpolados e roda a função novamente:
myLOCS <- cbind(c(5, 25, 15), c(8, 10, 15))
krige.conv(..., loc=myLOCS, ...)
On Mon, 20 Oct 2014, Helder Gramacho wrote:
Bom dia pessoal,
Para código reproduzível desta dúvida pode-se utilizar o exemplo disponível na página dos tutoriais do geoR: http://leg.ufpr.br/geoR/ tutorials/ctc.R Ao final deste tutorial chegamos ao mapa da distribuição espacial obtido por krigagem. Imagem inline 1 Gostaria de saber se é possível, obter o valor da variável para coordenadas específicas por exemplo na figura acima: x=5; y=8 x=25; y=10 x=15; y=15 Se é possível, como posso implementar no R.
Desde já agradeço toda colaboração,
Hélder Gramacho Recife-PE / agrohelder@gmail.com
_______________________________________________ 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.
Sim neste caso basta voce exportar o objeto no formato desejado por exemplo suponha que voce tenha um grid de predição em um objeto predLOCS que voce passou para krige.conv() em: kr <- krige.conv(... , loc=predLOCS, ...) basta montar o data-frame(ou marix) e exportar como voce faria com qualquer dado do R nestes formatos cbind(predLOCS, kr$predict) write.table(...) alternativamente voce pode exportar como raster, shapefile ou outros formatos On Mon, 20 Oct 2014, Helder Gramacho wrote:
Olá Prof. Paulo Justiniano, Funcionou por aqui, muito obrigado pela ajuda! Outra dúvida: É possível exportar o mapa da krigagem para abrir num SIG por exemplo? Um dos formatos que o SIG que utilizo abre é o formato ASCII grided XYZ.
Obrigado,
Hélder Gramacho Recife-PE / agrohelder@gmail.com
Em 20 de outubro de 2014 10:06, Paulo Justiniano <paulojus@leg.ufpr.br> escreveu: defina estas coordenadas como pontos a serem interpolados e roda a função novamente:
myLOCS <- cbind(c(5, 25, 15), c(8, 10, 15))
krige.conv(..., loc=myLOCS, ...)
On Mon, 20 Oct 2014, Helder Gramacho wrote:
Bom dia pessoal, Para código reproduzível desta dúvida pode-se utilizar o exemplo disponível na página dos tutoriais do geoR: http://leg.ufpr.br/geoR/tutorials/ctc.R Ao final deste tutorial chegamos ao mapa da distribuição espacial obtido por krigagem. Imagem inline 1 Gostaria de saber se é possível, obter o valor da variável para coordenadas específicas por exemplo na figura acima: x=5; y=8 x=25; y=10 x=15; y=15 Se é possível, como posso implementar no R.
Desde já agradeço toda colaboração,
Hélder Gramacho Recife-PE / agrohelder@gmail.com
_______________________________________________ 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.
Olá Prof. Paulo, Consegui, primeiro converti a saída do krige.conv() para "SpatialPixelsDataFrame" e em seguida exportei para formato tif com writeGDAL() Obrigado pela ajuda mais uma vez, *Hélder Gramacho * Recife-PE / *agrohelder@gmail.com <agrohelder@hotmail.com>* Em 20 de outubro de 2014 14:27, Paulo Justiniano <paulojus@leg.ufpr.br> escreveu:
Sim neste caso basta voce exportar o objeto no formato desejado
por exemplo suponha que voce tenha um grid de predição em um objeto predLOCS que voce passou para krige.conv() em:
kr <- krige.conv(... , loc=predLOCS, ...)
basta montar o data-frame(ou marix) e exportar como voce faria com qualquer dado do R nestes formatos
cbind(predLOCS, kr$predict) write.table(...)
alternativamente voce pode exportar como raster, shapefile ou outros formatos
On Mon, 20 Oct 2014, Helder Gramacho wrote:
Olá Prof. Paulo Justiniano,
Funcionou por aqui, muito obrigado pela ajuda! Outra dúvida: É possível exportar o mapa da krigagem para abrir num SIG por exemplo? Um dos formatos que o SIG que utilizo abre é o formato ASCII grided XYZ.
Obrigado,
Hélder Gramacho Recife-PE / agrohelder@gmail.com
Em 20 de outubro de 2014 10:06, Paulo Justiniano <paulojus@leg.ufpr.br> escreveu: defina estas coordenadas como pontos a serem interpolados e roda a função novamente:
myLOCS <- cbind(c(5, 25, 15), c(8, 10, 15))
krige.conv(..., loc=myLOCS, ...)
On Mon, 20 Oct 2014, Helder Gramacho wrote:
Bom dia pessoal, Para código reproduzível desta dúvida pode-se utilizar o exemplo disponível na página dos tutoriais do geoR: http://leg.ufpr.br/geoR/tutorials/ctc.R Ao final deste tutorial chegamos ao mapa da distribuição espacial obtido por krigagem. Imagem inline 1 Gostaria de saber se é possível, obter o valor da variável para coordenadas específicas por exemplo na figura acima: x=5; y=8 x=25; y=10 x=15; y=15 Se é possível, como posso implementar no R.
Desde já agradeço toda colaboração,
Hélder Gramacho Recife-PE / agrohelder@gmail.com
_______________________________________________ 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.
_______________________________________________ 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)
-
Helder Gramacho -
Paulo Justiniano