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