Hélio,
Postei o problema no R-sig-geo internacional (Extract a pixel
values in regular grid on a krigging
object<https://stat.ethz.ch/pipermail/r-sig-geo/2013-April/018003.html>), acompanhe
os desdobramentos em:
https://stat.ethz.ch/pipermail/r-sig-geo/2013-April/subject.html
Abraço,
Alexandre
Em 11/04/2013 07:12, Hélio Gallo Rocha escreveu:
Thiago,
Obrigado pelo retorno.
Consegui plotar, mas não consegui extrair os dados kc$predict com
coordenadas para fazer uma recomendação de calagem
Mas dá mais uma luz.
--
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho
Em 10 de abril de 2013 22:40, Thiago V. dos Santos
<[email protected] <mailto:[email protected]>> escreveu:
Hélio,
Algumas sugestões foram postadas por mim e pelo Paulo, que
envolvem converter o resultado da krigagem em um raster usando as
funções do sp ou do raster.
Certas perguntas podem ajudar na solução do problema: qual o nome
do objeto produzido pela krigagem em que você está interessado?
qual a classe desse objeto? onde estão as coordenadas dele? como
converter um objeto dessa classe para um objeto da classe raster?
Seguindo o código que eu Alexandre postou recentemente, eu faria
uma matriz ou dataframe a partir de lat lon e krig (baseado no
output da krigagem) e montaria um raster usando a função
rasterFromXYZ. Uma vez em raster, a conversão para polígono
(shapefile) é direta. Veja código abaixo, que pode ser acoplado ao
CMR da última mensagem do Alexandre:
#--------------------------------------
# First, add kriging values to the coordinates list
pred.grid$krig<-kc$predict
# Then, create a raster
r <- rasterFromXYZ(pred.grid)
# Last, check the result
plot(r)
# Converting to shapefile is straightforward
s <- rasterToPolygons(r)
# Plot shapefile (BEWARE - pretty slow)
plot(s)
#--------------------------------------
Saudações,
--
Thiago V. dos Santos
PhD student
Land and Atmospheric Science
University of Minnesota
http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm
Phone: (612) 323 9898 <tel:%28612%29%20323%209898>
________________________________
From: Hélio Gallo Rocha <[email protected]
<mailto:[email protected]>>
To: r-br <[email protected] <mailto:[email protected]>>
Sent: Wednesday, April 10, 2013 2:27 PM
Subject: [R-br] raster e shape file
Caros,
Também estou as voltas com resultados de krigagen.
Estou acompanhando a msg do Alexandre Santos.
Já pesquisei bastante, mas não achei solução para transformar o
resultado de uma krigagem em um raster ou num shape file.
Agradeço desde já
--
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho
_______________________________________________
R-br mailing list
[email protected] <mailto:[email protected]>
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.
--
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho
_______________________________________________
R-br mailing list
[email protected]
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.
--
======================================================================
Alexandre dos Santos
Proteção Florestal
Coordenador do curso Técnico em Florestas
Vice Coordenador do curso de Engenharia Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)
e-mails:[email protected]
[email protected]
======================================================================
_______________________________________________
R-br mailing list
[email protected]
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.
