# NOT RUN {
## The following example takes some time and download heavy data.
## It will therefore not be run unless you type:
## example(prepcipitate, run.dontrun = TRUE)
# }
# NOT RUN {
## We fit the models for Germany:
GNIPDataDEagg <- prepdata(data = GNIPDataDE)
GermanFit <- isofit(iso.data = GNIPDataDEagg,
mean.model.fix = list(elev = TRUE, lat.abs = TRUE))
elevation.raster <- relevate(
elevation.raster = ElevRasterDE,
isofit = GermanFit,
aggregation.factor = 0)
getprecip(path = "~/Desktop/")
precipitation.brick <- prepcipitate(path = "~/Desktop/",
elevation.raster = ElevRaster
)
if (require(rasterVis)) {
levelplot(precipitation.brick)
}
# }
Run the code above in your browser using DataLab