Interpolating and writing results into model input files
Usage
interpolate(
  meteo_lst,
  catchment_boundary_path,
  dem_data_path,
  grid_spacing,
  p_vector = c("PCP", "SLR", "RELHUM", "WNDSPD", "TMP_MAX", "TMP_MIN"),
  idw_exponent = 2
)Arguments
- meteo_lst
 nested list of lists with dataframes. Nested structure meteo_lst -> data -> Station ID -> Parameter -> Dataframe (DATE, PARAMETER).
- catchment_boundary_path
 path to basin boundary shape file.
- dem_data_path
 path to DEM raster data in same projection as weather station.
- grid_spacing
 numeric value for distance in grid. Units of coordinate system should be used.
- p_vector
 character vector representing weather variables to interpolate (optional, default all variables selected c("PCP", "SLR", "RELHUM", "WNDSPD", "TMP_MAX", "TMP_MIN" ).
- idw_exponent
 numeric value for exponent parameter to be used in interpolation (optional, default value is 2).
Value
nested list of lists with dataframes for interpolation results. Nested structure lst -> data -> Station ID -> Parameter -> Dataframe (DATE, PARAMETER). Function also writes all SWAT weather text input files from the interpolation results.
Examples
if (FALSE) {
temp_path <- system.file("extdata", "weather_data.xlsx", package = "svatools")
DEM_path <- system.file("extdata", "GIS/DEM.tif", package = "svatools")
basin_path <- system.file("extdata", "GIS/basin.shp", package = "svatools")
met_lst <- load_template(temp_path, 3035)
interpolate(met_lst, basin_path, DEM_path, 2000) 
}
