[R-br] raster e shape file
Hélio Gallo Rocha
heliogallorocha em gmail.com
Quinta Abril 18 15:20:11 BRT 2013
Alexandre
Não sei se já tinha este resultado com valores de kc$predict, mas consegui
colocar numa tabela as coordenadas(var1 e var2) e o kc$predict.
Não tive muito tempo de trabalhar nisso, mas o ideal era exportar para
shape, mas não consegui.
Quem puder ajudar, se não teremos de lançar mão em soft pago para gerar um
shape.
ai vai o CRM:
require(geoR)
library(raster)
ml <- likfit(s100, ini = c(1,0.5), fix.nugget = T)
# defining the grid
pred.grid <- expand.grid(seq(0,1, l=100), seq(0,1, l=100))# ;pred.grid
#pred.grid
# kriging calculations
kc<- krige.conv(s100, loc = pred.grid, krige = krige.control(obj.m = ml))
#recorte=matrix(c(.2,.8,.6,0,.1,.2,.7,.4),ncol=2) #c(0.2, 0.6, 0.2, 1.1,
0.2, 0.3, 1.0, 1.1, ncol=2)
#recorte
#plot(recorte)
ima=image(kc, loc = pred.grid,col=gray(seq(1,0.1,l=10)))
##tutorial do paulo para incrementar o mapa sde krigagem
contour(kc)
image(kc)
contour(kc, add=TRUE, nlev=21)
persp(kc, theta=20, phi=20)
contour(kc, filled=TRUE)
contour(kc, filled=TRUE, color=terrain.colors)
contour(kc, filled=TRUE, col=gray(seq(1,0,l=21)))
coordinates(pred.grid) <- ~ Var1+Var2;plot(pred.grid)## Transformando as
coordenadas do grid
#fim
kc$predict<-as.data.frame(kc$predict)### Valores dos pixels preditos
kc_output<-SpatialPixelsDataFrame(points=pred.grid, data=kc$predict)###
Criando o objeto SpatialPixelsDataFrame
#kc_output
#plot(x,y,xlab="predição da krigagem",ylab="dados")
edit(data.frame(kc_output))
Seria bom conseguir apenas os valores de uma área com borda, também não
consegui
Abraço
Hélio
Em 11 de abril de 2013 10:37, ASANTOS [via R-br] <
ml-node+s2285057n4658942h70 em n4.nabble.com> escreveu:
> 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 <[hidden email]<http://user/SendEmail.jtp?type=node&node=4658942&i=0>
> > 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: <a moz-do-not-send="true" href="tel:%28612%29%20323%209898" value="
>> +16123239898" target="_blank">(612) 323 9898
>>
>>
>> ________________________________
>>
>> From: Hélio Gallo Rocha <[hidden email]<http://user/SendEmail.jtp?type=node&node=4658942&i=1>
>> >
>> To: r-br <[hidden email]<http://user/SendEmail.jtp?type=node&node=4658942&i=2>
>> >
>> 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
>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4658942&i=3>
>>
>> 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[hidden email] <http://user/SendEmail.jtp?type=node&node=4658942&i=4>
> 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)[hidden email] <http://user/SendEmail.jtp?type=node&node=4658942&i=5>
> [hidden email] <http://user/SendEmail.jtp?type=node&node=4658942&i=6>
> ======================================================================
>
>
> _______________________________________________
> R-br mailing list
> [hidden email] <http://user/SendEmail.jtp?type=node&node=4658942&i=7>
> 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.
>
> ------------------------------
> If you reply to this email, your message will be added to the discussion
> below:
>
> http://r-br.2285057.n4.nabble.com/R-br-raster-e-shape-file-tp4658937p4658942.html
> To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4>
> .
> NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml>
>
--
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130418/a2f3fd24/attachment.html>
Mais detalhes sobre a lista de discussão R-br