The hardware and bandwidth for this mirror is donated by dogado GmbH, the Webhosting and Full Service-Cloud Provider. Check out our Wordpress Tutorial.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]dogado.de.

simulariatools

R-CMD-check CRAN status CRAN downloads DOI

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

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 = "./test/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()
#> Warning: Removed 1225 rows containing non-finite values
#> (`stat_contour_filled()`).

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()
#> Warning: Removed 280 rows containing missing values (`geom_raster()`).

List of functions

Available functions are listed below:

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

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.
Health stats visible at Monitor.