[R-br] Predição lenta - Interpolação Kriging

[Éder Comunello]

Thiago, bom dia!

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.

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.

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.

Ainda pretendo seguir a análise do seu código, mas queria adiantar esses
pontos.

Atte.,






Éder Comunello <c <comunello.eder em gmail.com>omunello.eder em gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]


Em 19 de junho de 2014 04:39, Thiago Cesar Lima Silveira <
thiagoclsilveira em yahoo.com.br> escreveu:

> Olá Éder,
>
> Obrigado pela disponibilidade!
>
> Quando mencionei 1200 pontos, me referi aos pontos reais de batimetria
> amostrados.
>
> O grid4 que me refiro é isso mesmo:
> Features    : 1630464
>
> extent      : 468333.3, 499983.3, 6634933, 6681223  (xmin, xmax, ymin,
> ymax)
> Resolution: 30 x 30
>
> Não cheguei a usar o IDW para batimetria.
>
> É 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.
>
> Segue abaixo como construí o gri4, os dados da batimetria te mando para o
> teu e-mail pessoal.
>
> #############
> # Grid4
> > grid <- raster( )
> > grid <- raster(ncol=1544, nrow=1056, xmn=468318.3, xmx=499998.3,
> ymn=6634918, ymx=6681238)
> > projection (grid) <- '+proj=utm +zone=22 +south +ellps=aust_SA +units=m
> +no_defs '
> > res(grid) <- 30
> > grid4 <- as(grid, 'SpatialPoints')
>
> Grande abraço,
>
> Thiago
>
>
>
> _________________________________
> Thiago Cesar Lima Silveira
> PhD candidate - Zoology PUCRS - Brazil
> Visitor PhD. student - School of Ocean Sciences -  Bangor University
> Building Westbury Mouth
> Room Nautilus 327
> Menai Bridge - Isle of Anglesey - UK
> Telephone.: 07843 115244
> e-mail: thiago.cesar em acad.pucrs.br
> CV: http://lattes.cnpq.br/5960267776845701
>
>
>
>
> On Jun 18, 2014, at 7:39 PM, Éder Comunello <comunello.eder em gmail.com>
> wrote:
>
> Thiago, boa tarde!
>
> Me disponho a dar uma olhada, mas vou precisar dos objetos interpol_pcsand
> (ou data_interpol) e grid4. Me chamou a atenção o fato de aparecer a
> menção de 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.
>
> Além disso, talvez você possa utilizar IDW para batimetria. A não ser que
> você esteja avaliando os parâmetros do modelo ajustado.
>
> Atte.,
>
> Éder Comunello <c <comunello.eder em gmail.com>omunello.eder em gmail.com>
> Dourados, MS - [22 16.5'S, 54 49'W]
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> 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.
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140620/098b54da/attachment.html>
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: bat1.png
Tipo: image/png
Tamanho: 9771 bytes
Descrição: não disponível
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140620/098b54da/attachment.png>