multiple_interpolate_points: Interpolates DWD values

In whamer/papros: PAthogen PROgnosis System

Description

This functions interpolates points based on a dataframe containing coordinates, an aim variable, an aim dataset and a parameter such as a date, by which the interpolations are divided

Usage

multiple_interpolate_points(
  dataframe,
  coords,
  epsg,
  splitter,
  aim_variable,
  outputfile,
  trans_epsg = FALSE,
  co_variables = FALSE,
  procedure = c("ked", "ok", "idw"),
  progressbar = TRUE
)

Arguments

dataframe	dataframe containing the aim variable and if "ked" should be applied the covariables
coords	vector containing names of the columns containing x and y coordinate values
epsg	number of the EPSG code of the "coords" information
splitter	name of the splitter column
aim_variable	Character string with the name of the aim variable
outputfile	SpatialPointsDataframe, SpatialGridDataFrame or raster which should be filled with predictions; requires covariables for "ked"
trans_epsg	default = FALSE, number of the EPSG code the "coords" information should be transformed to
co_variables	default = FALSE, vector of covariables if needed
procedure	default = c("ked","ok","idw"); vector containing the interpolation technic to be used; the first method is used and if this does not work out, the second, and so on
progressbar	default = TRUE; should a progressbar be generated?

Value

a dataframe or a SpatialPointsDataFrame containing information about DWD locations in Germany

Author(s)

Wolfgang Hamer

Examples

# Download example data
shdat <- download_statewide_hourly_station_data(state = "Schleswig-Holstein", coord = TRUE)

shdat2 <- shdat %>%
    filter(DateTime < sort(unique(shdat$DateTime))[5])

example <- multiple_interpolate_points(dataframe = shdat2,
                                       coords = c("lon","lat"),
                                       splitter = "DateTime",
                                       aim_variable ="Temperature",
                                       outputfile=c(1000,1000),
                                       co_variables = FALSE,
                                       procedure = c("ked","ok","idw"),
                                       epsg = 4326,
                                       trans_epsg = 25832)
plot(stack(example))

whamer/papros documentation built on Feb. 6, 2021, 8:54 a.m.