[R-br] % de cada area na krigagem

[Hélio Gallo Rocha]

Caros Prof. Paulo e Éder.

O CMR funcionou muito bem, resolvido

Muito obrigado aos dois

Hélio




Em 6 de novembro de 2013 18:33, Eder Comunello [via R-br] <
ml-node+s2285057n4660867h45 em n4.nabble.com> escreveu:

> Pessoal, boa tarde!
>
> Também precisei desse procedimento anteriormente e já tinha pesquisado
> alguma coisa. Embora a dúvida já tenha sido respondida, pode ser que o
> trecho de código que tenho aqui possa ser de interesse.
>
> A alternativa 1 é a solução colocada pelo Prof. Paulo. A alternativa 2 é
> mais trabalhosa, mas possibilita visualizar mais facilmente os fatiamentos
> efetuados. Além disso, haveria a possibilidade de vetorizar o objeto
> 'raster' para então obter as áreas ou armazenar em um shapefile.
>
> Os dados do CMR não tem projeção definida, mas a ideia de aplicação é a
> mesma.
>
> ### <BEGIN>
> ### Areas de representação matricial (raster)
> ### Usando o dataset s100 da geoR
> require(geoR); require(sp); require(raster)
> data(s100)
> vModel <- likfit(s100, ini=c(1,0.5), fix.nugget=T)    ### ajuste de modelo
> (não avaliado)
> pGrid  <- expand.grid(seq(0,1, l=30), seq(0,1, l=30)) ### grid de predição
> krig1  <- krige.conv(s100, loc=pGrid, krige=krige.control(obj.m=vModel))
> image(krig1, col=2:5, asp=1)
>
> ### Definição de classes
> range(krig1$pred)          ### observa intervalo das classes
> classes <- -1:3*1; classes ### 5 classes definidas
>
> ### Alternativa 1
> predFat    <- cut(krig1$pred, breaks=classes)
> predFatTab <- rbind(pixels=table(predFat),
> perc=round(prop.table(table(predFat))*100,2))
> predFatTab ### quantificação
>
> ### Alternativa 2
> pred <- cbind(pGrid,krig1$pred) ### data.frame com dados da predição
> coordinates(pred) <- ~Var1+Var2 ### transforma em SPointsDF
> gridded(pred) = TRUE            ### transforma em SPixelsDF
> rPred <- raster(pred)           ### transforma em raster
> plot(rPred)
>
> rPredFat <- cut(rPred, breaks=classes) ### fatia o objeto raster
> plot(rPredFat) ### visualiza o fatiamento
> freq(rPredFat, useNA='no') ### quantificação
>
> ### Comparando as alternativas
> tmp <- t(as.table(freq(rPredFat, useNA='no')))[2,]; tmp
> round(prop.table(tmp)*100,2)
> rPredFatTab <- rbind(pixels=tmp, perc=round(prop.table(tmp)*100,2))
>
> predFatTab; rPredFatTab
> ### <END>
>
>
> Éder Comunello <[hidden email]<http://user/SendEmail.jtp?type=node&node=4660867&i=0>[hidden
> email] <http://user/SendEmail.jtp?type=node&node=4660867&i=1>>
> Dourados, MS - [22 16.5'S, 54 49'W]
>
>
>
> _______________________________________________
> R-br mailing list
> [hidden email] <http://user/SendEmail.jtp?type=node&node=4660867&i=2>
> 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-de-cada-area-na-krigagem-tp4660853p4660867.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/20131106/f64f5cad/attachment.html>