[R-br] Krigagem: Obter valores pontuais a partir da sup. interpolada
Helder Gramacho
agrohelder em gmail.com
Segunda Outubro 20 18:14:12 BRST 2014
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 em gmail.com <agrohelder em hotmail.com>*
Em 20 de outubro de 2014 14:27, Paulo Justiniano <paulojus em 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 em gmail.com
>>
>>
>> Em 20 de outubro de 2014 10:06, Paulo Justiniano <paulojus em 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 em gmail.com
>>
>>
>>
>> _______________________________________________
>> R-br mailing list
>> R-br em 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 em 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141020/e3205231/attachment.html>
Mais detalhes sobre a lista de discussão R-br