Yesterday I successfully requested a time series from the WeatherSOS service! The code below shows how that works: First requesting the data (including using some help-functions to create time interval and access the features of the SOS), second plotting it using R’s
weathersos = SOS("http://v-swe.uni-muenster.de:8080/WeatherSOS/sos") go.offering = sosOfferings(weathersos)[] # ATHMOSPHERIC_TEMPERATURE go.observedProperty = sosObservedProperties(weathersos)[] # temperature urn go.eventTime3 = sosCreateEventTime(sosCreateTimePeriod(sos = weathersos, begin = as.POSIXct("2010-09-16 18:00"), end = as.POSIXct("2010-09-20 18:00"))) obs3 <- getObservation(sos = weathersos,observedProperty = list(go.observedProperty), procedure = list(sosProcedures(weathersos)[]), eventTime = go.eventTime3, offering = go.offering@id) obs3@result # heureka! summary(obs3@result) # finally! plot(x = obs3@result[["Time"]], y = obs3@result[["urn:ogc:def:property:OGC::Temperature"]], type = "l", main = "Temperature in MÃ¼nster", xlab = "Time", ylab = "Temperature (Â°C)")
This is quite a milestone for the project, and a good start in this week. For today, spatial queries and testing (also other SOS instances) are on the agenda.