Realiza a krigagem em uma área delimitada

A função realiza, a partir de um objeto da classe GeoModel, a krigagem em uma área delimitada. O efeito pepita é incluído por padrão. Mas pode ser removido se nugget = FALSE.

# S4 method for class 'GeoModel'
krig(
  object,
  cellsize = 100,
  newarea,
  nugget,
  krig.type = c("ordinary", "universal"),
  fast = TRUE,
  crs_grid,
  buffer = 5,
  ...
)

Arguments

object

Objeto da classe GeoModel.

cellsize

Tamanho do pixel utilizado para criação da malha, em metros.

newarea

Polígono da classe sf ou sfc contendo a região delimitada para krigagem. Veja detalhes.

nugget

Lógico opcional informando se o efeito pepita deve ser incluído. Por padrão é FALSE.

krig.type

Caractere informando o tipo de krigagem entre "ordinary" e "universal". Padrão é "ordinary".

fast

Lógico informando se a predição é usando C++ ou não. Padrão é TRUE.

crs_grid

Numérico contendo o código epsg do sistema projetado de coordenadas para criação do grid de predição.

buffer

Inteiro informando o buffer em metros aplicado à envoltória convexa. Só utilizado quando envoltória convexa é calculada. Veja detalhes.

...

Demais argumentos a serem passados para predict.

Value

Retorna um objeto da classe GeoModelKrig.

Details

Por padrão, a krigagem é realizada em uma área delimitada por GeoModelObj. Veja que tal objeto já vem de um outro objeto Inv ou suas herdeiras, que devem obrigatoriamente possuir objetos espaciais associados para realizar o ajuste do modelo. Um desses objetos é a bordadura, incluída via função talhao() aplicada ao objeto Inv. Veja Inv, talhao, GeoModel para mais informações.

Caso não existe bordadura pré-especificada, então a função criará uma envoltória convexa a partir dos pontos observados, com adição de um buffer com padrão de 5 metros. Essa envoltória convexa será utilizada para criar a malha de predição.

Alternativamente é possível realizar a krigagem em outra região delimitada. Se for o caso, a nova área deve ser passada no argumento newarea. Veja que newarea deve ser um polígono regular da classe sfc ou sf. O sistema de coordenadas será automaticamente transformado para o mesmo dos pontos usados no ajuste de GeoModelObj, mas nenhuma conferência será feita acerca da proximidade dos pontos. Se newarea for informado, então esse polígono será utilizado na krigagem.

O argumento cellsize define a resolução/densidade da malha de predição. Se os objetos espaciais estiverem com sistemas de coordenadas geográficas, então a função encontra internamente a zona UTM e cria o grid a partir das coordenadas transformadas em WGS 84 UTM. Se a área de predição contemplar mais de uma zona, então uma projeção policônica/Albers baseada em SIRGAS 2000 é utilizada. O usuário pode incluir outra projeção para cálculo via argumento crs_grid.

Author

Gabriel Agostini Orso