retrieve_predictor: read predictor values at the point locations, returns a...
In mengluchu/APMtools: statistical air pollution mapping tools

read predictor values at the point locations, returns a dataframe or save into a csv

retrieve_predictor(
  predictorstack,
  pointdatafile,
  coordslonlat,
  proj,
  csvname = "predictors"
)

`predictorstack`	raster stack to read point values from
`pointdatafile`	the file contains point locations to be read
`coordslonlat`	a vector specify the names of the longitude and latitude variable names.
`proj`	the projection of the point file
`csvname`	the name to save the csv file

lus = raster("/data/lu01/NWA/predictor/NLstack.grd")
lf_lo = list.files("/data/lu01/NWA/Bakfietsdata", pattern = "^.*morning.*.csv$", full.names = T)
bakfile1 = read.csv(lf_lo[1])
proj = "+proj=longlat +datum=WGS84"
retrieve_predictor(lus, bakfile1, c("Lon", "Lat"), proj)