prepare.functions: Functions to prepare input for 'kriging'
In jobonaf/pesco: Data fusion for air quality data

Functions to prepare input for kriging

prepare.day(day, obs.daily, ctm.daily,
            pollutant,
            emis.winter, emis.summer,
            elev = NULL, verbose = FALSE)                        

prepare.ctm(ctm.daily, day,
            x.pnt, y.pnt, x.grd, y.grd,
            conc.min = 10^-6)
prepare.emis(emis.winter, emis.summer, day,
             x.pnt, y.pnt, x.grd = NULL, y.grd = NULL)
prepare.elev(elev, x.pnt, y.pnt, 
             z.pnt = rep(NA,length(x.pnt)),
             x.grd = NULL, y.grd = NULL)   
prepare.obs(obs.daily, day, pollutant,
             conc.min = 10^-6)

`day`	required day, in the format "YYYY-MM-DD"
`obs.daily`	data frame with daily observations, as returned by `dailyObs`
`ctm.daily`	CTM output, aggregated on a daily basis, as returned by `dailyCtm`
`emis.winter`	winter emissions, in a list of 2 elements: coords coordinates, a list of 2 numeric vectors x and y data numeric vector
`emis.summer`	summer emissions, in a list of 2 elements: coords coordinates, a list of 2 numeric vectors x and y data numeric vector
`elev`	elevation, in a list of 2 elements: coords coordinates, a list of 2 numeric vectors x and y data numeric vector
`pollutant`	name of the column with pollutant concentations
`verbose`	logical; if `TRUE` some messages are given
`x.pnt`	numeric vector of x coordinates of the stations
`y.pnt`	numeric vector of y coordinates of the stations
`z.pnt`	numeric vector of elevation of the stations
`x.grd`	numeric vector of x coordinates of the grid cells
`y.grd`	numeric vector of y coordinates of the grid cells
`conc.min`	minimum concentration; if a concentration is less than this value, then it is set equal to it

To prepare the required input for the function kriging, you need only the function prepare.day, which includes all the functions prepare.ctm, prepare.emis, prepare.elev and prepare.obs.