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
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.,
R-br mailing list
R-br@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.