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 <thi_veloso@yahoo.com.br>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/inde... Phone: (612) 323 9898

________________________________ From: Hélio Gallo Rocha <heliogallorocha@gmail.com> To: r-br <r-br@listas.c3sl.ufpr.br> 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 R-br@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.

-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

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@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/inde... 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... 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

Boa tarde,   Estou tentando exportar para o R o quadro 14 da POF 2008/2009 usando o pacote dicionariosIBGE elaborado por vcs.   Gostaria de saber    

library(dicionariosIBGE) data(dicpof2008.2009) caminho.microdados <- file.path("C:","DADOS","T_INVENTARIO_S.txt") dados <- le.pesquisa(dicionario=dic5pof2008.2009, rotulos=rot5pof2008.2009, pathname.in=caminho.microdados, tbloco=577450, nlines=577450) Erro em data.frame(cod = codigos, stringsAsFactors = FALSE) :   argumento "codigos" ausente, sem padrão   Sabendo q falta o argumento “códigos” e vendo na discrição do help da função como “Array of characters with codes of the wanted variables to be read, as specified in dictionary”. No entanto, não consigo encontrar os códigos dessa variável.   Vcs podem me ajudar?   Abs   Fernando Cunha