<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="mailto:comunello.eder@gmail.com" target="_blank">c</a><a href="mailto:omunello.eder@gmail.com" target="_blank">omunello.eder@gmail.com</a>> <br>
Dourados, MS - [22 16.5'S, 54 49'W]<br></div></div>
<br><br></div></div>