<div dir="ltr">Elias, bom dia!<div><br></div><div>Reconheço o mérito das suas colocações, só complementei! </div><div><br></div><div>Aproveitando a deixa... A exemplo do material do Jackson Aquino, mantenho também para consulta um material seu, "Introdução à análise de dados espacialmente referenciados". Tenho aqui uma versão de 2012. Tens alguma edição mais atual ou intenção de transformar esse material em livro?</div>
<div><br></div><div>Grato,</div><div><br><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 class="gmail_quote">Em 29 de maio de 2014 07:18, Elias T. Krainski <span dir="ltr"><<a href="mailto:eliaskrainski@yahoo.com.br" target="_blank">eliaskrainski@yahoo.com.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">

  
    
  
  <div bgcolor="#FFFFFF" text="#000000">
    humilhou... :)<div><div class="h5"><br>
    <br>
    <div>On 29/05/14 12:36, Éder Comunello
      wrote:<br>
    </div>
    </div></div><blockquote type="cite"><div><div class="h5">
      <div dir="ltr">
        <div>Hélio, bom dia!</div>
        <div><br>
        </div>
        <div>Faz sentido interpolar? Acredito que sim, partindo do
          princípio que a dependência espacial pode existir, mas seu
          método não permitiu detectá-la. Acredito que a causa mais
          comum seja o esquema de amostragem utilizado, sobretudo no que
          se refere à distância entre amostras. No caso particular, você
          conseguiu modelar em algumas áreas e não em outras, dando
          margem à essa interpretação.</div>
        <div><br>
        </div>
        <div>O colega Elias já postou uma solução, mas acrescento o
          código abaixo, caso ainda tenha interesse em utilizar o
          {gstat}.</div>
        <div><br>
        </div>
        <div>Em termos gerais, a ideia de operação no {gstat} é similar
          a do {geoR}, no sentido em que você vai precisar criar um grid
          pra receber o resultado da interpolação. O espaçamento da
          grade será a resolução da interpolação. O ponto principal é
          que tem que trabalhar com objetos da classe 'sp'.</div>
        <div><br>
        </div>
        <div>Verifique os parâmetros número de vizinhos (nmax) e
          peso/potência da  distância (idp). No {gstat} por default
          utilizam-se todos os pontos e o idp=2.</div>
        <div><br>
        </div>
        <div>Atte.,</div>
        <br>
        <div>
          <div><font face="courier new, monospace">### <code r></font></div>
          <div><font face="courier new, monospace">sapply(c("gstat",
              "sp", "geoR", "RColorBrewer"), require, character.only=T)</font></div>
          <div><font face="courier new, monospace"><br>
            </font></div>
          <div><font face="courier new, monospace">data(parana);
              names(parana)</font></div>
          <div><font face="courier new, monospace">hist(parana$data,
              col=3)</font></div>
          <div><font face="courier new, monospace">points(parana,
              pt.divide='quart')</font></div>
          <div><font face="courier new, monospace"><br>
            </font></div>
          <div><font face="courier new, monospace">### Criar objeto 'sp'</font></div>
          <div><font face="courier new, monospace">pr <-
              data.frame(x=parana$coords[,1], y=parana$coords[,2],
              chuva=parana$data)</font></div>
          <div><font face="courier new, monospace">names(pr);
              coordinates(pr) <- ~x+y</font></div>
          <div><font face="courier new, monospace">class(pr)</font></div>
          <div><font face="courier new, monospace">plot(pr, asp=1,
              axes=T, pch=20); polygon(parana$borders, border=2)</font></div>
          <div><font face="courier new, monospace"><br>
            </font></div>
          <div><font face="courier new, monospace">### Criar 'grid'</font></div>
          <div><font face="courier new, monospace">bbox(pr)</font></div>
          <div><font face="courier new, monospace">lim   <-
              round(bbox(pr)+c(-1,-1, 1, 1)*50); lim ### amplia a área
              do bbox</font></div>
          <div><font face="courier new, monospace">grid  <-
              expand.grid(x=seq(lim[1,1],lim[1,2], by=10),
              y=seq(lim[2,1],lim[2,2], by=10))</font></div>
          <div><font face="courier new, monospace">grid  -> <a href="http://grid.pt" target="_blank">grid.pt</a></font></div>
          <div><font face="courier new, monospace">gridded(grid) = ~x+y
                     ### SpatialPixels</font></div>
          <div><font face="courier new, monospace">coordinates(<a href="http://grid.pt" target="_blank">grid.pt</a>)
              = ~x+y ### SpatialPoints</font></div>
          <div><font face="courier new, monospace"><br>
            </font></div>
          <div><font face="courier new, monospace">### IDW Default</font></div>
          <div><font face="courier new, monospace">idw    <-
              idw(pr$chuva~1, pr, grid)</font></div>
          <div>
            <font face="courier new, monospace"><a href="http://idw.pt" target="_blank">idw.pt</a>
              <- idw(chuva~1, pr, <a href="http://grid.pt" target="_blank">grid.pt</a>)</font></div>
          <div><font face="courier new, monospace"><br>
            </font></div>
          <div><font face="courier new, monospace">### Visualização</font></div>
          <div><font face="courier new, monospace">spplot(idw,  
               "var1.pred", main = "IDW Default")</font></div>
          <div><font face="courier new, monospace">spplot(<a href="http://idw.pt" target="_blank">idw.pt</a>,
              "var1.pred", main = "IDW default - Mapa Pontuado")</font></div>
          <div><font face="courier new, monospace">image(idw);
              polygon(parana$borders); points(pr)</font></div>
          <div><font face="courier new, monospace"><br>
            </font></div>
          <div><font face="courier new, monospace">### Variando
              parâmetros</font></div>
          <div><font face="courier new, monospace">idw.data <-
              data.frame(</font></div>
          <div><font face="courier new, monospace">          default =
              idw(chuva ~ 1, pr, grid)$var1.pred, ### nmax=todos &
              idp=2</font></div>
          <div>
            <font face="courier new, monospace">          idw6    =
              idw(chuva ~ 1, pr, grid, nmax=6)$var1.pred,</font></div>
          <div><font face="courier new, monospace">          idw9    =
              idw(chuva ~ 1, pr, grid, nmax=9)$var1.pred,</font></div>
          <div><font face="courier new, monospace">          idp1    =
              idw(chuva ~ 1, pr, grid, idp = 1)$var1.pred,</font></div>
          <div><font face="courier new, monospace">          idp4    =
              idw(chuva ~ 1, pr, grid, idp = 4)$var1.pred,</font></div>
          <div><font face="courier new, monospace">          idp8    =
              idw(chuva ~ 1, pr, grid, idp = 8)$var1.pred)</font></div>
          <div><font face="courier new, monospace"><br>
            </font></div>
          <div><font face="courier new, monospace">grid.data <-
              SpatialPixelsDataFrame(grid, idw.data)</font></div>
          <div><font face="courier new, monospace">dput(names(grid.data))</font></div>
          <div><font face="courier new, monospace"><br>
            </font></div>
          <div><font face="courier new, monospace">spplot(grid.data,
              c("default", "idw6", "idw9", "idp1", "idp4", "idp8"), main
              = "IDW", col.regions=rainbow(16)</font></div>
          <div><font face="courier new, monospace"># idp:  numeric;
              specify the inverse distance weighting power</font></div>
          <div><font face="courier new, monospace"># nmax: the number of
              nearest observations that should be used </font></div>
          <div><font face="courier new, monospace">### </code></font></div>
        </div>
        <div class="gmail_extra"><br>
        </div>
        <div class="gmail_extra">
          <div><img style="margin-right: 0px;"><br>
          </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 class="gmail_quote"><br>
          </div>
        </div>
      </div>
      <br>
      <fieldset></fieldset>
      <br>
      </div></div><div class=""><pre>_______________________________________________
R-br mailing list
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
    </div></blockquote>
    <br>
  </div>

<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" 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" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div></div></div>