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] < [email protected]> 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
_______________________________________________ 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.
