<html><head><meta http-equiv="Content-Type" content="text/html charset=windows-1252"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">Bom dia Éder,<div><br></div><div>Obrigado pelo retorno.</div><div><br></div><div>Estou estudando as metodologias de interpolação não havia pensado na alternativa que mencionaste a respeito do TIN. Somente me pareceu lógico criar uma linha de borda com uma informação que sei que é zero para a batimetria.</div><div><br></div><div>A partir disto interpolar por kriging os valores. No entanto, confesso que preciso de um caminho para construir um mapa e que também consiga avaliar a qualidade da interpolação.</div><div><br></div><div>Estou estudando pelos exemplos apresentados pelo Bivand (2008) - Applied Spatial Data Analysis with R.</div><div><br></div><div>Grande abraço e obrigado mais uma vez.</div><div><br></div><div>Thiago</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br><div><div>On Jun 20, 2014, at 12:39 PM, Éder Comunello <<a href="mailto:comunello.eder@gmail.com">comunello.eder@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div dir="ltr"><div>Thiago, bom dia!</div><div><br></div><div>Olhei seus dados e pude notar que na verdade você tem 236 pontos cotados. Os 1263 restantes são na verdade da linha de borda, transformados em pontos e que receberam o valor zero. Na figura anexa, eu isolei os pontos cotados e usei os demais pra restaurar a linha de borda.</div>
<div><br></div><div>O uso de linhas de borda (um tipo de breaklines) é uma prática comum quando você utiliza grades triangulares (TIN) pra modelar tridimensionalmente uma área. A partir da grade triangular você interpola por outro método pra obter uma grade regular. Talvez seja o procedimento mais adequado pro seu caso.</div>
<div><br></div><div>Ao optar por outro método que não o uso de TIN, eu particularmente, nunca tinha visto adicionar os pontos da linha de borda aos dados para interpolar e, num primeiro momento, acredito que seria melhor interpolar os 236 pontos e utilizar a linha de borda apenas pra recortar a superfície gerada. Mas precisaria estudar melhor sobre a questão.</div>
<div><br></div><div>Ainda pretendo seguir a análise do seu código, mas queria adiantar esses pontos.</div><div><br></div><div>Atte.,</div><div><br></div><div><div><img width="653" height="653" apple-inline="yes" id="89FC2547-D45E-4E6D-B715-95512F23A05D" apple-width="yes" apple-height="yes" src="cid:ii_hwnepncq1_146b9022328d706f"><br>
</div><br></div><div><br></div><div><br></div></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">Em 19 de junho de 2014 04:39, Thiago Cesar Lima Silveira <span dir="ltr"><<a href="mailto:thiagoclsilveira@yahoo.com.br" target="_blank">thiagoclsilveira@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; position: static; z-index: auto;"><div style="word-wrap:break-word">Olá Éder, <div><br></div><div>Obrigado pela disponibilidade!</div><div><br></div><div>
Quando mencionei 1200 pontos, me referi aos pontos reais de batimetria amostrados. </div><div><br></div><div>O grid4 que me refiro é isso mesmo: </div><div>Features : 1630464 <div class=""><br>extent : 468333.3, 499983.3, 6634933, 6681223 (xmin, xmax, ymin, ymax)</div>
</div><div>Resolution: 30 x 30</div><div><br></div><div>Não cheguei a usar o IDW para batimetria. </div><div><br></div><div>É a primeira vez que tento construir mapas interpolados no R e a minha intenção é desenvolver uma linha de trabalho para gerar os mapas e avaliar a qualidade da interpolação.</div>
<div><br></div><div>Segue abaixo como construí o gri4, os dados da batimetria te mando para o teu e-mail pessoal.</div><div><br></div><div>#############</div><div># Grid4 </div><div>> grid <- raster( )</div><div>
> grid <- raster(ncol=1544, nrow=1056, xmn=468318.3, xmx=499998.3, ymn=6634918, ymx=6681238)</div><div>> projection (grid) <- '+proj=utm +zone=22 +south +ellps=aust_SA +units=m +no_defs ‘</div><div>> res(grid) <- 30</div>
<div>> grid4 <- as(grid, ‘SpatialPoints’)</div><div><br></div><div>Grande abraço,</div><div><br></div><div>Thiago</div><div class=""><div><br></div><div><br></div><div><br><div><div style="word-wrap:break-word"><div style="word-wrap:break-word">
<div>_________________________________<br>Thiago Cesar Lima Silveira<br>PhD candidate - Zoology PUCRS - Brazil<br>Visitor PhD. student - School of Ocean Sciences - Bangor University</div><div><span style="font-size:10pt;text-align:justify">Building Westbury Mouth</span></div>
<div><span style="font-size:13px;text-align:justify">Room Nautilus 327</span></div><div><span style="font-size:10pt;text-align:justify">Menai Bridge - Isle of Anglesey - UK</span></div><div>Telephone.: 07843 115244<br>e-mail: <a href="mailto:thiago.cesar@acad.pucrs.br" target="_blank">thiago.cesar@acad.pucrs.br</a><br>
CV: <a href="http://lattes.cnpq.br/5960267776845701" target="_blank">http://lattes.cnpq.br/5960267776845701</a></div></div></div></div></div><div><br></div><div><br></div><div><br></div></div><div><div>
</div>
<br><div><div><div class="h5"><div>On Jun 18, 2014, at 7:39 PM, Éder Comunello <<a href="mailto:comunello.eder@gmail.com" target="_blank">comunello.eder@gmail.com</a>> wrote:</div><br></div></div><blockquote type="cite">
<div><div class="h5"><div dir="ltr">Thiago, boa tarde!<div><br></div><div>Me disponho a dar uma olhada, mas vou precisar dos objetos <span style="font-size:13px;font-family:arial,sans-serif">interpol_pcsand (ou data_</span><span style="font-size:13px;font-family:arial,sans-serif">interpol) e grid4. Me chamou a atenção o fato de aparecer a menção de </span><font face="arial, sans-serif">1.630.464 features no seu objeto grid4. Se for isso mesmo são mais de 1200x1200pontos ao invés de 1200 pontos como mencionado.</font></div>
<div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">Além disso, talvez você possa utilizar IDW para batimetria. A não ser que você esteja avaliando os parâmetros do modelo ajustado.</font></div>
<div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">Atte.,</font></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></div></div><div class="">
_______________________________________________<br>R-br mailing list<br><a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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.</div></blockquote></div><br></div></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>
_______________________________________________<br>R-br mailing list<br><a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br<br>Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.</blockquote></div><br></div></body></html>