[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