simulariatools

simulariatools is a free collection collection of functions and tools
useful to pre and post process data for air quality modelling and
assessment. The package is developed and maintained by the people at
Simularia and it is extensively used for
their daily job.

If you use this package in your work, please consider citing it. Refer
to its Zenodo DOI to cite it.

Table of Contents

• Installation
• Brief examples
• List of Functions
• Contact
• Contributors
• License
• Changelog

Installation

Install the released version of simulariatools from CRAN:

install.packages("simulariatools")

Or install the development version from GitHub with:

pak::pkg_install("Simularia/simulariatools")

Note: in order to use importADSOBIN() to import ADSO/BIN data files,
a working installation of Python 3 is required. For more information
about R and Python interoperability, please refer to
reticulate documentation.

Brief examples

Contour plot

First, import air quality data from NetCDF or ADSO/BIN files with
the appropriate convenience function:

library(simulariatools)
mydata <- importRaster(file = "./development/conc_avg.nc",
                       k = 1000,
                       destaggering = TRUE,
                       variable = "nox",
                       verbose = TRUE)
#> 
#> Raster statistics -----------------------------------------------
#>        X (min, max, dx)  :  496000.000   519250.000      250.000
#>        Y (min, max, dy)  : 4943000.000  4955250.000      250.000
#>      nox (min, max, mean):    0.00e+00     2.71e+00     1.52e-01
#> -----------------------------------------------------------------

A quick contour plot with default configuration can be easily obtained
by running contourPlot2() without any argument:

contourPlot2(mydata)

The plot is customisable by using contourPlot2() arguments and by
piping ggplot2 instructions:

library(ggplot2)
contourPlot2(mydata,
             domain = c(502000, 519000, 4943125, 4955125, 5, 5),
             levels = c(-Inf, 0.5, 1, 1.5, 2, Inf),
             legend = "NOx [ug/m3]") +
    labs(x = NULL, y = NULL) +
    theme_minimal()

Use ggsave() to save the last plot to file:

ggsave(filename = "~/path/to/myplot.png", width = 7, height = 6, dpi = 300)

Use tile optional argument to produce a plot without interpolation:

library(ggplot2)
contourPlot2(mydata,
             tile = TRUE,
             legend = "NOx [ug/m3]") +
    labs(x = NULL, y = NULL) +
    theme_minimal()

List of functions

Available functions are listed below:

• contpourPlot2()
• downloadBasemap()
• importRaster()
• importADSOBIN()
• importSurferGrd()
• plotAvgRad()
• plotAvgTemp()
• plotStabilityClass()
• removeOutliers()
• rollingMax()
• stabilityClass()
• vectorField()
• contourPlot() (deprecated)
• createBaseMap() (deprecated)

Deprecated functions will be removed in the near future.

Contact

Contact person:

Giuseppe Carlino
Simularia s.r.l.
g.carlino@simularia.it

Contributors

Matteo Paolo Costa