pointgridding: Grid point data to spatial grid

In mjlaine/fastgrid: Gridding point data to a regular grid

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.