pointgridding: Grid point data to spatial grid

View source: R/gridutils.R

pointgriddingR Documentation

Grid point data to spatial grid

Description

This is an easier to use version of fastkriege

Usage

pointgridding(
  y = NULL,
  pointdata,
  modelgrid,
  variable = "y",
  trend_model = NULL,
  priorfield = FALSE,
  cov.pars
)

Arguments

y

data matrix of point data to be gridded, alternative to giving variable

pointdata

SpatialPointsDataFrame of the spatial points

modelgrid

SpatialGrid or SpatialGridDataFrame of the output grid

variable

name of the variable in pointdata to be gridded

cov.par

Kriging parameters, sigma^2, phi, tau^2.

Examples

 

lonmin <- 19.0; lonmax <- 33.0
latmin <- 59.0; latmax <- 71.5
newgrid <- makelonlatgrid(lonmin=lonmin,lonmax=lonmax,latmin=latmin,latmax=latmax,londx=0.1,latdx=0.1)
nobs <- 50
obs <- data.frame(longitude=runif(nobs,lonmin,lonmax),
                 latitude=runif(nobs,latmin,latmax),
                 temperature=rnorm(nobs,mean=10,sd=4))
sp::coordinates(obs) <- c('longitude','latitude')
out <- pointgridding(pointdata=obs, modelgrid = newgrid, variable="temperature",cov.pars = c(0.5^2,2.0,0.0))
MOSplotting::MOS_plot_field(out,stations=obs,cmin=0,cmax=20)


mjlaine/fastgrid documentation built on July 24, 2022, 2:43 a.m.