Simple plotting of a meteorological series for a given point.

meteoplot(object, index=1, var="MeanTemperature", 
          fun=NULL, freq=NULL, dates = NULL, months = NULL, add = FALSE,...)



A data frame with daily meteorological data (in this case index is not used) or an object of class SpatialPointsMeteorology. Alternatively, an object of class SpatialPointsDataFrame containing the meta data (columns dir, filename and possibly format) of meteorological files.


An integer to indicate the point in the SpatialPointsMeteorology object (or the SpatialPointsDataFrame object).


The meteorological variable to be plotted.


The name of a function to be calculated for summaries (only valid if freq is specified).


A string giving an interval specification for summaries (e.g., "week", "month", "quarter" or "year").


An object of class Date to define the period to be plotted. If dates = NULL then all dates in object are processed.


A numeric vector to indicate the subset of months for which plotting is desired (e.g. c(7,8) for July and August). When combined with fun and freq, this parameter allows plotting summaries for particular seasons. For example fun = "sum" freq = "years" and months = 6:8 leads to plotting the sum over summer months of each year.


A flag to indicate wether drawing should be done on the current plot (using function lines).


Additional parameters for functions plot or lines.


Daily precipitation is plotted using bars (i.e. type = "h" when calling plot). Otherwise the function draws lines (i.e. type = "l" when calling plot). If object is of class SpatialPointsDataFrame-class then the function reads the meteorological data to be plotted from the disk.


Miquel De Cáceres Ainsa, CREAF

See also



#Creates spatial topography points from the grid
p = 1:2
spt = as(examplegridtopography, "SpatialPointsTopography")[p]

#Interpolation of two points for the whole time period (2000-2003)
mp = interpolationpoints(exampleinterpolationdata, spt)
#> Processing point '1' (1/2) - done.
#> Processing point '2' (2/2) - done.

#Plot interpolated meteorological series
meteoplot(mp,1, ylab="Daily mean temperature")

meteoplot(mp,1, ylab="Monthly mean temperature", fun=mean, freq="months")