Archives for posts with tag: gstat

Within the UncertWeb project the ifgi provides a SOS with parts of the AirBase dataset provided by the EEA. Release 0.1-15 contains a great demo on what you can do with that data: How to build up a spatial interpolation, in this case inverse distance weighting (IDW), on top of a GetObservation result. This comprises some data formatting and spatial projection.

Thanks go to Edzer Pebesma for contributing the interpolation examples using the packages gstat and spacetime and pointing out bugs in the package!

Some plots generated by the demo are:

  1. overview plot of the available offerings,
  2. a histogram of all downloaded NO2 values,
  3. a bubble plot of the data for Germany,
  4. the stations in Germany in an UTM projection,
  5. two IDW interpolations, one aggregated over time and one for a specific time stamp, and last but not least
  6. a time series plot for a whole year for one station including a polynomically fitted regression line.

Go ahead and run the analysis yourself:


library("sos4R")

demo("airquality")

I intended to do a nice demo with CSIRO’s South Esk Testbed. But I got stuck with a long planned feature – the plotting of offerings and SOS. So this is what this demo is about now, and additionally it shows how to combine measurements of different SOSs and convert these to sp objects for further analysis. Below is a map of the available sensors.

CSIRO South Esk Testbed Map

It must be noted that the focus of the demo is on data consolidation and plotting. The demo contains code for two Kriging examples using the gstat and automap packages, which sadly I did not get to work yet. These analyses are not necessarily sensible!

Try it out:


library("sos4R")

demo("southesk")

Some example plots:

The demo will be in a soon to be released new version of sos4R, but the documentation needs to be updated before uploading it to CRAN. Read the rest of this entry »