<div dir="ltr">Caros Prof. Paulo e Éder.<div><br></div><div>O CMR funcionou muito bem, resolvido</div><div><br></div><div>Muito obrigado aos dois</div><div><br></div><div>Hélio<br><div><br></div><div><br></div></div></div>
<div class="gmail_extra"><br><br><div class="gmail_quote">Em 6 de novembro de 2013 18:33, Eder Comunello [via R-br] <span dir="ltr"><<a href="mailto:ml-node+s2285057n4660867h45@n4.nabble.com" target="_blank">ml-node+s2285057n4660867h45@n4.nabble.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">Pessoal, boa tarde!<div><br></div><div>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.</div>
<div><br></div><div>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.</div>
<div><br></div><div>Os dados do CMR não tem projeção definida, mas a ideia de aplicação é a mesma.</div><div><br></div><div><div><font face="courier new, monospace">### <BEGIN></font></div><div><font face="courier new, monospace">### Areas de representação matricial (raster)</font></div>
<div><font face="courier new, monospace">### Usando o dataset s100 da geoR</font></div><div><font face="courier new, monospace">require(geoR); require(sp); require(raster) </font></div><div><font face="courier new, monospace">data(s100)</font></div>
<div><font face="courier new, monospace">vModel <- likfit(s100, ini=c(1,0.5), fix.nugget=T) ### ajuste de modelo (não avaliado)</font></div><div><font face="courier new, monospace">pGrid <- expand.grid(seq(0,1, l=30), seq(0,1, l=30)) ### grid de predição</font></div>
<div><font face="courier new, monospace">krig1 <- krige.conv(s100, loc=pGrid, krige=krige.control(obj.m=vModel)) </font></div><div><font face="courier new, monospace">image(krig1, col=2:5, asp=1)</font></div><div><font face="courier new, monospace"><br>
</font></div><div><font face="courier new, monospace">### Definição de classes</font></div><div><font face="courier new, monospace">range(krig1$pred) ### observa intervalo das classes</font></div><div><font face="courier new, monospace">classes <- -1:3*1; classes ### 5 classes definidas</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### Alternativa 1</font></div><div><font face="courier new, monospace">predFat <- cut(krig1$pred, breaks=classes)</font></div>
<div><font face="courier new, monospace">predFatTab <- rbind(pixels=table(predFat), perc=round(prop.table(table(predFat))*100,2))</font></div><div><font face="courier new, monospace">predFatTab ### quantificação</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### Alternativa 2</font></div><div><font face="courier new, monospace">pred <- cbind(pGrid,krig1$pred) ### data.frame com dados da predição</font></div>
<div><font face="courier new, monospace">coordinates(pred) <- ~Var1+Var2 ### transforma em SPointsDF</font></div><div><font face="courier new, monospace">gridded(pred) = TRUE ### transforma em SPixelsDF</font></div>
<div><font face="courier new, monospace">rPred <- raster(pred) ### transforma em raster</font></div><div><font face="courier new, monospace">plot(rPred)</font></div><div><font face="courier new, monospace"><br>
</font></div><div><font face="courier new, monospace">rPredFat <- cut(rPred, breaks=classes) ### fatia o objeto raster</font></div><div><font face="courier new, monospace">plot(rPredFat) ### visualiza o fatiamento</font></div>
<div><font face="courier new, monospace">freq(rPredFat, useNA='no') ### quantificação</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### Comparando as alternativas</font></div>
<div><font face="courier new, monospace">tmp <- t(as.table(freq(rPredFat, useNA='no')))[2,]; tmp</font></div><div><font face="courier new, monospace">round(prop.table(tmp)*100,2)</font></div><div><font face="courier new, monospace">rPredFatTab <- rbind(pixels=tmp, perc=round(prop.table(tmp)*100,2))</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">predFatTab; rPredFatTab</font></div><div><font face="courier new, monospace">### <END></font></div></div><div><br></div>
<div class="gmail_extra"><br clear="all"><div><div dir="ltr">Éder Comunello <<a href="http://user/SendEmail.jtp?type=node&node=4660867&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a><a href="http://user/SendEmail.jtp?type=node&node=4660867&i=1" rel="nofollow" link="external" target="_blank">[hidden email]</a>> <br>
Dourados, MS - [22 16.5'S, 54 49'W]<br></div></div>
<br><br></div></div><div class="im">
<br>_______________________________________________
<br>R-br mailing list
<br><a href="http://user/SendEmail.jtp?type=node&node=4660867&i=2" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.
<br>
<br>
<hr noshade size="1" color="#cccccc">
</div><div style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif"><div class="im">
<div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
</div><a href="http://r-br.2285057.n4.nabble.com/R-br-de-cada-area-na-krigagem-tp4660853p4660867.html" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-de-cada-area-na-krigagem-tp4660853p4660867.html</a>
</div><div class="HOEnZb"><div class="h5">
<div style="color:#666;font:11px tahoma,geneva,helvetica,arial,sans-serif;margin-top:.4em;line-height:1.5em">
To unsubscribe from R-br, <a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4" target="_blank">click here</a>.<br>
<a href="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" rel="nofollow" style="font:9px serif" target="_blank">NAML</a>
</div></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br>
</div>