<div dir="ltr"><span style="font-family:arial,sans-serif;font-size:13px">Alexandro</span><div style="font-family:arial,sans-serif;font-size:13px">Vi sua msg sobre valores do predict, me interessa também, vou tentar.</div><div style="font-family:arial,sans-serif;font-size:13px">
Olhando seu CRM em:</div><div class="im" style="font-family:arial,sans-serif;font-size:13px"><div><span style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif">pixelsize <- (1.0-0.0/30)*(1.0-0.0/30)#Tamanho do pixel é dado por seq </span><br>
</div><div><span style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif"><br></span></div></div><div style="font-family:arial,sans-serif;font-size:13px"><span style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif">não teria de colocar um parentesis antes de dividir? (</span><span style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif">1.0-0.0)/30</span></div>
<div style="font-family:arial,sans-serif;font-size:13px"><span style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif"><br></span></div><div style="font-family:arial,sans-serif;font-size:13px"><span style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif">Hélio</span></div>
</div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 7 de abril de 2013 03:55, Thiago Veloso [via R-br] <span dir="ltr"><<a href="mailto:ml-node+s2285057n4658914h2@n4.nabble.com" target="_blank">ml-node+s2285057n4658914h2@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 style="font-size:12pt;font-family:times new roman,new york,times,serif"><div><div><div style="font-size:12pt;font-family:'times new roman','new york',times,serif"><div style="font-family:'times new roman','new york',times,serif;font-size:12pt">
<span>Alexandre,</span></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:'times new roman','new york',times,serif"><span><br></span></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:'times new roman','new york',times,serif">
Você já tem o polígono gridPol de onde você quer criar a média da krigagem. Não conheço bem as classes do pacote geoR, mas as coordenadas do kc$predict parecem estar no objeto pred.grid. <span style="font-size:12pt">A abordagem que eu usaria seria criar um objeto raster da krigagem e em seguida contar os valores usando a função extract (raster) usando como argumentos o raster da krigagem e o gridPol.</span></div>
<div style="font-style:normal;font-size:16px;background-color:transparent;font-family:'times new roman','new york',times,serif"><br></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:'times new roman','new york',times,serif">
Além da ajuda para
 a função extract no manual do raster, talvez o último exemplo dessa página possa te ajudar: <a href="http://stackoverflow.com/questions/14682606/extract-value-from-raster-stack-from-spatialpolygondataframe" style="font-size:12pt" rel="nofollow" link="external" target="_blank">http://stackoverflow.com/questions/14682606/extract-value-from-raster-stack-from-spatialpolygondataframe</a></div>
<div style="font-style:normal;font-size:16px;background-color:transparent;font-family:'times new roman','new york',times,serif"><br></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:'times new roman','new york',times,serif">
Mas não estou certo se era isso que você estava buscando. Acho que o Paulo vai te ajudar melhor.</div><div class="im"><div></div><div style="font-family:'times new roman','new york',times,serif;font-size:12pt">
 </div><div style="font-family:'times new roman','new york',times,serif;font-size:12pt">Saudações,<br>--<br>Thiago V. dos Santos<br>PhD student<br>Land and Atmospheric Science<br>University of Minnesota<br>
<a href="http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm" target="_blank">http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm</a><br>Phone: <a href="tel:%28612%29%20323%209898" value="+16123239898" target="_blank">(612) 323 9898</a><br>
</div>  </div><div style="font-family:'times new roman','new york',times,serif;font-size:12pt"> <div style="font-family:'times new roman','new york',times,serif;font-size:12pt"> <div dir="ltr">
 <font face="Arial"><div class="im"> <hr size="1">  <b><span style="font-weight:bold">From:</span></b> ASANTOS <<a href="http://user/SendEmail.jtp?type=node&node=4658914&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a>><br>
 </div><b><span style="font-weight:bold">To:</span></b> <a href="http://user/SendEmail.jtp?type=node&node=4658914&i=1" rel="nofollow" link="external" target="_blank">[hidden email]</a>; Thiago V. dos Santos <<a href="http://user/SendEmail.jtp?type=node&node=4658914&i=2" rel="nofollow" link="external" target="_blank">[hidden email]</a>> <br>
 <b><span style="font-weight:bold">Sent:</span></b> Saturday, April 6, 2013 10:26 PM<div><div class="h5"><br> <b><span style="font-weight:bold">Subject:</span></b> Re: [R-br] Obter valor médio dos pixels de uma imagem gerada por krigagem<br>
 </div></div></font> </div><div><div class="h5"> <br><div>
  

    
  
  <div>
    Thiago,<br>
    <br>
        Já pensei isto, mas para SpatialPoints(vp), preciso fazer um
    list() das coordenadas dos pixels (vp$x e vp$y) e isto não tenho e
    nem sei como recuperar o valor de um pixel à partir de uma
    coordenada conhecida.<br>
    <br>
    Abraço,<br>
    <br>
    Alexandre<br>
    <br>
    <br>
    <div>Em 06/04/2013 21:32, Thiago V. dos
      Santos escreveu:<br>
    </div>
    <blockquote style="border-left:2px solid #cccccc;padding:0 1em" type="cite">
      <div style="font-size:12pt;font-family:'times new roman','new york',times,serif">
        <div><span>Nao conheco bem todas as funcoes do sp, mas para pela
            mensagem de erro eu transformaria a tabela vp em
            SpatialPoints, que eh o argumento esperado pela funcao over.</span></div>
        <div> </div>
        <div>Saudações,<br>
          --<br>
          Thiago V. dos Santos<br>
          PhD student<br>
          Land and Atmospheric Science<br>
          University of Minnesota<br>
<a rel="nofollow" href="http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm" link="external" target="_blank">http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm</a><br>

          Phone: <a href="tel:%28612%29%20323%209898" value="+16123239898" target="_blank">(612) 323 9898</a><br>
        </div>
        <div style="font-family:'times new roman','new york',times,serif;font-size:12pt">
          <div style="font-family:'times new roman','new york',times,serif;font-size:12pt">
            <div dir="ltr"> <font face="Arial">
                <hr size="1"> <b><span style="font-weight:bold">From:</span></b>
                ASANTOS <a href="http://user/SendEmail.jtp?type=node&node=4658914&i=3" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>
                <b><span style="font-weight:bold">To:</span></b>
                <a href="http://user/SendEmail.jtp?type=node&node=4658914&i=4" rel="nofollow" link="external" target="_blank">[hidden email]</a> <br>
                <b><span style="font-weight:bold">Sent:</span></b>
                Saturday, April 6, 2013 8:12 PM<br>
                <b><span style="font-weight:bold">Subject:</span></b>
                Re: [R-br] Obter valor médio dos pixels de uma imagem
                gerada por krigagem<br>
              </font> </div>
            <br>
            Boa noite pessoal,<br>
            <br>
                  Acho evoluí um pouco na questão, mas agora a
            dificuldade esta em contar os pixels contidos no interior
            das malhas, quando tento um over(), ocorre o seguinte erro:<br>
            <br>
            Erro em function (classes, fdef, mtable)  :<br>
              unable to find an inherited method for function "over",
            for signature "table", "SpatialPolygons"<br>
            <br>
            <br>
            segue CRM:<br>
            #<br>
            require(geoR)<br>
            plot(s100$coords, xlim=c(0,1.2), ylim=c(0,1.2), xlab="Coord
            X", ylab="Coord Y")<br>
            loci <- matrix(c(0.2, 0.6, 0.2, 1.1, 0.2, 0.3, 1.0, 1.1),
            ncol=2)<br>
            text(loci, as.character(1:4), col="red")<br>
            polygon(x=c(0,1,1,0), y=c(0,0,1,1), lty=2)<br>
            <br>
            # Fitting models with nugget fixed to zero<br>
            ml <- likfit(s100, ini = c(1,0.5), fix.nugget = T)<br>
            <br>
            # defining the grid<br>
            pred.grid <-  expand.grid(seq(0,1, l=51), seq(0,1, l=51))<br>
            # kriging calculations<br>
            kc <- krige.conv(s100, loc = pred.grid, krige =
            krige.control(obj.m = ml))<br>
            # displaying predicted values<br>
            image(kc, loc = pred.grid, col=gray(seq(1,0.1,l=30)),
            xlab="Coord X", ylab="Coord Y")<br>
            #<br>
            <br>
            #Centro das 25 parcela 0.2 x 0.2<br>
            p0<-c(0.1,0.1)<br>
            #<br>
            c0 <- seq(p0[1],p0[1]+0.9,by=0.2)<br>
            d0 <- seq(p0[2],p0[2]+0.9,by=0.2)<br>
            d1 <- expand.grid(x=c0, y=d0)<br>
            points(d1[,1],d1[,2],col="red")<br>
            ver<-c(0.0,0.2,0.4,0.6,0.8,1.0)<br>
            hor<-c(0.0,0.2,0.4,0.6,0.8,1.0)<br>
            abline(v=ver, col="red")<br>
            abline(h=hor, col="red")<br>
            #<br>
            require(sp)<br>
            grid=GridTopology(c(0.1,0.1), c(0.2,0.2), c(5,5))### Cria as
            área<br>
            grid=SpatialGrid(grid=grid)<br>
            class(grid)  ## SpatialPixels<br>
            gridPol <- as(grid, "SpatialPolygons") ## Transfoma as
            areas em polígonos espaciais<br>
            options(digits=5)<br>
            pixelsize <- (1.0-0.0/30)*(1.0-0.0/30)#Tamanho do pixel é
            dado por seq Xmaior - Xmenor/l * Ymaior-Ymenor/l<br>
            vp <- table(area <- kc$predict)## Valor dos pixels<br>
            count.pix<-table(over(vp, gridPol)) #contar quantos
            pontos tem em cada area<br>
            <br>
            <br>
            Obrigado,<br>
            <br>
            Alexandre<br>
            <br>
            Em 06/04/2013 17:11, ASANTOS escreveu:<br>
            > Boa noite membros da lista,<br>
            > <br>
            >        Gostaria de receber ideias para calcular o valor
            médio dos pixel em uma imagem gerada por krigagem em
            sub-regiões de interesse na imagem, por exemplo usando um
            exemplo do Prof. Paulo, onde ajusta-se um semivariograma e
            faz a krigagem ordinária dos dados:<br>
            > <br>
            > #<br>
            > require(geoR)<br>
            > plot(s100$coords, xlim=c(0,1.2), ylim=c(0,1.2),
            xlab="Coord X", ylab="Coord Y")<br>
            > loci <- matrix(c(0.2, 0.6, 0.2, 1.1, 0.2, 0.3, 1.0,
            1.1), ncol=2)<br>
            > text(loci, as.character(1:4), col="red")<br>
            > polygon(x=c(0,1,1,0), y=c(0,0,1,1), lty=2)<br>
            > <br>
            > # Fitting models with nugget fixed to zero<br>
            > ml <- likfit(s100, ini = c(1,0.5), fix.nugget = T)<br>
            > <br>
            > # defining the grid<br>
            > pred.grid <-  expand.grid(seq(0,1, l=51), seq(0,1,
            l=51))<br>
            > # kriging calculations<br>
            > kc <- krige.conv(s100, loc = pred.grid, krige =
            krige.control(obj.m = ml))<br>
            > # displaying predicted values<br>
            > image(kc, loc = pred.grid, col=gray(seq(1,0.1,l=30)),
            xlab="Coord X", ylab="Coord Y")<br>
            > #<br>
            > <br>
            > #Agora vou dividir esta área em 25 partes e lançar os
            centroides destas subáreas:<br>
            > <br>
            > #Centro das 25 parcela 0.2 x 0.2<br>
            > p0<-c(0.1,0.1)<br>
            > #<br>
            > c0 <- seq(p0[1],p0[1]+0.9,by=0.2)<br>
            > d0 <- seq(p0[2],p0[2]+0.9,by=0.2)<br>
            > d1 <- expand.grid(x=c0, y=d0)<br>
            > points(d1[,1],d1[,2],col="red")<br>
            > ver<-c(0.0,0.2,0.4,0.6,0.8,1.0)<br>
            > hor<-c(0.0,0.2,0.4,0.6,0.8,1.0)<br>
            > abline(v=ver, col="red")<br>
            > abline(h=hor, col="red")<br>
            > #<br>
            > <br>
            >  E aqui vem minha pergunta, como calcular o valor médio
            dos pixels contidos em cada subárea de 0.2 x 0.2?<br>
            > <br>
            > Obrigado,<br>
            > <br>
            <br>
            --
            ======================================================================<br>
            Alexandre dos Santos<br>
            Proteção Florestal<br>
            Coordenador do curso Técnico em Florestas<br>
            Vice Coordenador do curso de Engenharia Florestal<br>
            IFMT - Instituto Federal de Educação, Ciência e Tecnologia
            de Mato Grosso<br>
            Campus Cáceres<br>
            Caixa Postal 244<br>
            Avenida dos Ramires, s/n<br>
            Bairro: Distrito Industrial<br>
            Cáceres - MT                      CEP: 78.200-000<br>
            Fone: <a href="tel:%28%2B55%29%2065%208132-8112" value="+556581328112" target="_blank">(+55) 65 8132-8112</a> (TIM)  <a href="tel:%28%2B55%29%2065%209686-6970" value="+556596866970" target="_blank">(+55) 65 9686-6970</a> (VIVO)<br>

            e-mails:<a href="http://user/SendEmail.jtp?type=node&node=4658914&i=5" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>
                    <a href="http://user/SendEmail.jtp?type=node&node=4658914&i=6" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>
======================================================================<br>
            <br>
            _______________________________________________<br>
            R-br mailing list<br>
            <a href="http://user/SendEmail.jtp?type=node&node=4658914&i=7" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>
            <a rel="nofollow" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
            Leia o guia de postagem (<a rel="nofollow" href="http://www.leg.ufpr.br/r-br-guia" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>)
            e forneça código mínimo reproduzível.<br>
            <br>
            <br>
          </div>
        </div>
      </div>
      <br>
      <fieldset></fieldset>
      <br>
      <pre>_______________________________________________
R-br mailing list
<a href="http://user/SendEmail.jtp?type=node&node=4658914&i=8" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<a rel="nofollow" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a rel="nofollow" href="http://www.leg.ufpr.br/r-br-guia" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
    </blockquote>
    <br>
    <pre>-- 
======================================================================
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: <a href="tel:%28%2B55%29%2065%208132-8112" value="+556581328112" target="_blank">(+55) 65 8132-8112</a> (TIM)   <a href="tel:%28%2B55%29%2065%209686-6970" value="+556596866970" target="_blank">(+55) 65 9686-6970</a> (VIVO)
<a href="http://user/SendEmail.jtp?type=node&node=4658914&i=9" rel="nofollow" link="external" target="_blank">[hidden email]</a> 
        <a href="http://user/SendEmail.jtp?type=node&node=4658914&i=10" rel="nofollow" link="external" target="_blank">[hidden email]</a>                    
======================================================================
</pre>
  </div>

</div><br><br> </div></div></div> </div>  </div></div></div></div><div><div class="h5"><br>_______________________________________________
<br>R-br mailing list
<br><a href="http://user/SendEmail.jtp?type=node&node=4658914&i=11" 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><div style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif"><div><div class="h5">
                <div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
                </div></div><a href="http://r-br.2285057.n4.nabble.com/R-br-Obter-valor-medio-dos-pixels-de-uma-imagem-gerada-por-krigagem-tp4658903p4658914.html" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-Obter-valor-medio-dos-pixels-de-uma-imagem-gerada-por-krigagem-tp4658903p4658914.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>