Perceptual Analysis, Visualization and Organization of Spectral Color Data in R

A cohesive framework for parsing, analyzing and organizing color from spectral data.

cran version cran downloads Build Status Coverage status

An R package for the spectral and spatial analysis of color patterns

Currently maintained by Rafael Maia, Thomas White, and Hugo Gruson.


pavo is an R package developed with the goal of establishing a flexible and integrated workflow for working with spectral and spatial colour data. It includes functions that take advantage of new data classes to work seamlessly from importing raw spectra and images, to visualisation and analysis. It provides flexible ways to input spectral data from a variety of equipment manufacturers, process these data, extract variables, and produce publication-quality figures.

pavo was written with the following workflow in mind:

  • Organise data by importing and processing spectra and images (e.g., to remove noise, negative values, smooth curves, etc.).
  • Analyse the resulting files, using spectral analyses of shape (hue, saturation, brightness), visual models based on perceptual data, and/or spatial adjacency and boundary strength analyses.
  • Visualise the output, with multiple options provided for exploration, presentation, and analysis.

Need more information, or help with the package?

  • Read the Package Vignettes (or via browseVignettes('pavo')) for detailed examples and discussion.
  • Check out the Latest News for changes and updates.
  • Can't find what you're looking for? Send an email to the mailing list: [email protected]

Citing pavo

The manuscripts describing the package has been published and are free to access:

> v. 2.0

Maia R., Gruson H., Endler J.A. and White T.E. 2018 pavo 2: new tools for the spectral and spatial analysis of colour in R. bioRxiv. doi: 10.1101/427658

< v. 2.0

Maia R., Eliason C.M., Bitton P.-P., Doucet S.M. and Shawkey M.D. 2013. pavo: an R Package for the analysis, visualization and organization of spectral data. Methods in Ecology and Evolution 4(10):609-613. doi: 10.1111/2041-210X.12069


This is the development page for pavo. The stable release is available from CRAN. Simply use install.packages('pavo') to install.

If you want to install the bleeding edge version of pavo, you can:

# install.packages('remotes')
  • download files from GitHub and install using $R CMD INSTALL or, from within R:
install.packages(path, type='source', repos=NULL)


pavo 2.1.0


  • added the argument reclass to procimg(), which allows users to interactively correct areas within images that have been misclassified
  • added the rod sensitivity of Canis familiaris
  • peakshape() uses a completely different algorithm to find the FWHM. It now works as expected for spectra with multiple peaks. See PR #137 for a detailed overview of the changes.
  • data used internally by pavo (bgandilum, transmissiondata, ttvertex, vissyst) is no longer exposed to users


  • new functions is.vismodel() and is.colspace() are exported to test whether an object is of class vismodel or colspace, respectively
  • fixed a bug where images would sometimes be wrongly detected as user-classified in as.rimg()
  • the UV-sensitive cone is now only always named "u", even for VS species (such as pfowl and avg.v in vismodel() and sensdata()). This removes an unnecessary but harmless warning when colspace() was used to place quantum catches of such species in the tetrahedral colour space.
  • the achro argument in coldist() has been changed for achromatic to better match the arguments from vismodel(). Older scripts that use achro should not be affected and still work as before.
  • the package imager is no longer a dependency, and is only loaded if using some features of procimg().
  • the package mapproj is no longer a dependency, and is only loaded if using projplot().
  • added the argument labels.stack to plot.rspec, which allows the use of custom spectra labels in stacked plots.
  • users now receive a warning when interpolating beyond the limits of the data using as.rspec, and can control the behaviour with the new argument exceed.limits.
  • all deprecated functions and arguments have now been fully removed.
  • as.rspec() now accepts both numeric and character vectors to identify the wavelength column using whichwl (eg. whichwl = "wl").
  • Reference images in classify() can now be specified using either a numeric vector (to identify by image position in a list) or character vector (to identify by image name).
  • fixed a bug in aggspec() when wavelength column was previously removed by the user.
  • fixed a bug where cocplot() would failed whenever type graphical parameter was specified.
  • spec2rgb() has been simplified to rely more on vismodel(). As a result, output values may be slightly different but upon testing, we found that differences between the old and the new version were barely noticeable.
  • the vignette have been split into three smaller parts, which should help new users to get started with pavo
  • numerous under-the-hood changes for stability and speed, with thanks to three reviewers and an associate editor at MEE.

pavo 2.0.0


  • image-based workflow for the combined analysis of colour and pattern geometry
  • added the visual systems of the (trichromatic) jumping spider Habronattus pyrrithrix and the (trichromatic) triggerfish Rhinecanthus aculeatus
  • getspec() can now read Avantes binary files (.TRM files)


  • Carotenoid chroma (S9) in summary.rspec() has been fixed to (R700 - R450)/R700. This gives the same result as before but with a flipped sign, and better reflects the original formula in the literature.
  • cieLAB values have been rescaled, and are expressed in the more standard range: L [0,100], ab [-128,127]
  • getspec() has an additional argument set to TRUE by default to ignore case in file extension matching
  • fix a bug where getspec() would sometimes fail with files including numbers in scientific format
  • add a new option in tetraplot() to add cone names (u,s,m,l)

pavo 1.4.0


  • getspec() can now read OceanOptics .ProcSpec files
  • added the visual system of Ctenophorus ornatus, the (trichromatic) ornate dragon lizard


  • getspecf() (and the argument fast = TRUE in getspec()) have been deprecated
  • summary.rspec() returned incorrect values for S7. If you use S7, please re-run your analyses


  • summary.rspec() now properly outputs NA for monotonically decreasing spectra
  • fixed warning when subset.rspec() was provided with a logical vector
  • fixed harmless warning when summary.colspace() was used on a tcs object
  • by argument in merge.rspec() is no longer ignored
  • fixed bug in voloverlap() when plot = TRUE
  • fixed bug in vismodel() when transmission has more than one column
  • fixed bug in vismodel() that applied von Kries correction to achromatic channel
  • added argument fill=FALSE in voloverlap()
  • fixed bug in jndplot() when suppressing the plotting of arrows
  • better handling of subset data when using summary.colspace() and summary.vismodel()
  • fixed bug in coldist() when noise = "quantum" and achro = TRUE were used
  • fixed bug in jndplot() when arrow = "none" and achro = TRUE
  • spec2rgb() now takes into account the 390-400 nm wavelength range into account when possible
  • as.rspec() no longer fails with tibbles
  • bin option of procspec() now works for all values of bins
  • non-relative quantum catches from dataframe object were not correctly scaled in "di", "tri", "categorical" and "coc" colorspaces
  • fixed a bug in colspace() where it would incorrectly infer a preference for a general trichromatic space, when a cie model is more appropriate
  • fixed a bug so that cie color matching functions can be more easily be used in a general trichromatic space (i.e. maxwell triangle)

pavo 1.3.1


  • bootcoldist() uses bootstrap to calculate the confidence intervals on the mean color distance between samples


  • fixed bug in coldist() when fewer samples than the number of photoreceptors
  • fixed bug in getspec() that would cause it to crash with errors
  • fixed bug in tetraplot() when saving images

pavo 1.3.0


  • jnd2xyz() converts distances (in JND, resulting from a coldist() call) into cartesian coordinates
  • plot() methods for objects resulting from jnd2xyz()
  • jndrot() produces rotations of Cartesian coordinates resulting from jnd2xyz()
  • coldist2mat() converts coldist() result from a pairwise data.frame to a distance matrix
  • sensdata() function for retrieving and/or visualising pavo's in-build spectral data


  • tetraplot() and cieplot() have been completely rewritten to allow finer viewing control
  • tetraplot() allows forced perspective using size to denote distance
  • voloverlap() and vol() have also been changed to work with the new tetraplot() options
  • getspec() has been rewritten to be faster, more general, and allow parallel processing
  • subset functions now allow more than one argument to be used, and allow further attributes to be passed onto grep (e.g. invert = TRUE)


  • fixed bug in coldist() on log-transformation when object was neither of class vismodel nor colspace
  • fixed bug in dL calculation when input is a colspace object
  • fixed bug in vismodel() when a data frame, matrix or rspec object was passed as the background
  • fixed bug in colspace() models when using non-standard receptor names or ordering
  • fixed bug in hexagon() model when calculating location & metrics for achromatic stimuli
  • fixed location of red vertex in tetraplot()
  • fixed bug in the argument names for expanding text labels in colspace plots
  • removed na.rm argument from aggspec() that was causing a bug when the error function did not have that argument. User should pass it as an argument to the function if necessary.
  • changed default to achro=FALSE in coldist()
  • replaced the modelled receptor sensitivities of the honeybee Apis mellifera with the empirical sensitivities from Peitsch et al (1992)
  • the built-in 'green' background spectrum is no longer normalized
  • removed wavelength limitations in the calculation of H3 from summary.rspec
  • all visual systems (except CIE) have been normalized to have an integral of 1

pavo 1.2.0


  • added the CIELch model accessed via colspace(space = 'cielch')
  • added the sensdata() function for retrieving and/or visualising pavo's in-build spectral data


  • vignettes have been amalgamated & the single, main vignette is now up-to-date
  • added more informative labels for the segment analysis plot

pavo 1.1.0


  • segspace() replaces the deprecated segclass(), and is accessed via the colspace() argument space = 'segment'. The results of segspace() are also now compatible with coldist() for the estimation of Euclidean colour-distances.
  • segplot() is a plot for Endler's (1990) segment analysis, and is accessed — along with all other 2d plots — via plot.colspace()


  • the use of relative quantum catches is now optional in the categorical colorspace (though still produces a warning), for greater flexibility
  • updated several functions to work when rspec object has only one spectrum
  • fixed bug in voloverlap where interactive plots would result in error
  • fixed incorrect labels in the maxwell triangle plot
  • fixed a bug in as.rspec() in which lim was not applied when interpolate = FALSE
  • fixed bug in aggplot() which resulted in error when using lty, lwd arguments
  • warning if ocular media is being used in both vismodel() and sensmodel()
  • added an 'all' option to the achromatic argument in vismodel()
  • added the ability to calculate dL for cielab models in coldist()
  • added some more informative messages and warnings

pavo 1.0

  • See vignette for detailed description of changes.


  • coldist() arguments have been changed. Now the empirically estimated value for the Weber fraction must be entered, instead of the noise-to-signal ratio. The noise-to-signal ratio is then calculated based on the empirically estimated Weber fraction for the reference cone type, and applied to the remaining cone types. This should avoid confusion between empirically estimated values for the Weber fraction and the noise-to-signal ratio, which are currently prevalent in the literature.
  • coldist() now has an additional argument, weber.achro, so that the value for the Weber fraction to be used to calculate achromatic contrast can be input independently of the cone ratios.
  • tcs() is deprecated, replaced by colspace().


  • colspace() replaces tcs() and introduces several new colorspaces
  • plot() methods for several colspace() outputs, including a static tetrahedral colorspace
  • projpoints() allows the plotting of points in a projplot() figure
  • vol() draws volume polygons in static tetrahedral plots
  • axistetra() draws reference x, y and z axis arrows in static tetrahedral plots
  • legendtetra() adds legends to static tetrahedral plots


  • summary.colspace() for tcs spaces now returns relative color volume as well as absolute
  • tcsvol() and voloverlap() now allow control for line width
  • procspec() fixed error when attempting to smooth rspec object without column names
  • procspec() handles smoothing before fixing negatives to avoid re-adding negatives when smoothing
  • procspec(), aggplot() accept additional arguments to summary functions (e.g., na.rm=TRUE)
  • peakshape() default wavelength limits (lim) now taken from rspec object rather than 300-700
  • peakshape() returns warning if a spectrum contains duplicate reflectance values
  • summary.rspec() works with single spectra
  • aggplot() and aggspec() fixed bug on ordering of levels when they don't match the sequence in the rspec object
  • aggplot() added logical argument "legend" for automatically adding a legend to the plot
  • vismodel() returns error if bkg=NULL
  • getspec() patched to stop returning warnings in Yosemite
  • getspec() has a faster (~5-10X), but less flexible, algorithm used when all input files are from the same source.

pavo 0.5-6


  • fixed bug in calculation of dichromat contrast in coldist()

pavo 0.5-5


  • fixed bug in calculation of H3 in summary.rspec()

pavo 0.5-4


  • changed default values for coldist()

pavo 0.5-2


  • fixes to the blue tit visual system, changed vismodel() argument to "bluetit"

pavo 0.5-1


  • vismodel() accepts matrix, data.frame or rspec objects for the illuminant, updated warning messages associated with this use
  • vismodel() accepts user-defined achromatic receptors
  • tcsplot(), tcsvol() & tcspoints(): transparency control passed to user
  • getspec() works with OceanView files

pavo 0.5

  • updated citation()


  • vismodel() vonkries = TRUE does not return a NULL result
  • vismodel() works with a single spectrum object
  • tcsplot() allows greater control of tetrahedron appearance
  • summary.rspec() allows for user-defined minimum wavelength (for calculation of UV variables)

pavo 0.3-1


  • irrad2flux() and flux2irrad() to convert illuminant measurements


  • vismodel() less cryptic error messages
  • as_rspec() fix message pertaining wavelength column
  • getspec() now removes empty columns (generated by bad tabulation)
  • plot.rspec() fix color labelling issue (previously, when user specified fewer colors than number of spectra, 'stack' and 'overlay' colored spectra differently)

pavo 0.3


  • subset() class methods for rspec, vismodel and tcs
  • summary() method for vismodel, returns attributes used in visual model


  • vismodel() output includes only Qi or fi (as selected by qcatch argument)
  • vismodel() von Kries correction is now an optional argument
  • coldist() subset argument for partial filtering
  • Updated sensitivity curves


  • getspec() works with Avasoft 8 output
  • aggspec() if no "by" argument is supplied, applies function to all spectra


  • tcs() accepts receptor names other than usml, issues warning
  • aggplot() allows control for different line types per spectra
  • voloverlap() now assigns darkgrey color for overlap if color vector is of length 1 or 2

pavo 0.1-2


  • merge.rspec() combines rspec objects in a single object


  • Changed vismodel() output: von Kries correction is now an optional argument, output includes only Qi and fi
  • vismodel() accepts a scale parameter (multiplies illuminant in order to make it in flux units)
  • coldist() includes option for receptor noise calculation
  • Updated sensitivity curves


  • sensmodel() accepts user-defined ocular media transmission
  • voloverlap() allows greater user control of plotting options
  • voloverlap() includes Monte Carlo option for overlap calculation
  • as.rspec() accepts "whichwl" argument for user-defined wavelength column selection
  • as.rspec() includes "interp" argument; updated automatic search for wavelength column
  • aggspec() aggregates spectra by multiple vectors (e.g. average spectra by species and sex using by=list(species,sex) )
  • peakshape() gives plot titles, allows ask=TRUE and mfrow


  • explorespec() "free" argument fixed
  • summary.rspec() S5 variable fixed; segments now divided by B1 (brightness-independent measure of saturation)
  • summary.rspec() checks for positive values when bmaxneg is calculated, returns NA
  • tcs() accepts "v" cone class
  • aggspec() fixed matching of names
  • procspec() works with rspec objects that include a single reflectance spectrum
  • peakshape() minor fixes

Reference manual

It appears you don't have a PDF plugin for this browser. You can click here to download the reference manual.


2.1.0 by Rafael Maia, 15 days ago,

Report a bug at

Browse source code at

Authors: Rafael Maia [aut, cre] , Thomas White [aut] , Hugo Gruson [aut] , John Endler [aut] , Chad Eliason [aut] , Pierre-Paul Bitton [aut]

Documentation:   PDF Manual  

GPL (>= 2) license

Imports geometry, pbmcapply, plot3D, readbitmap, rcdd, sp, viridisLite

Suggests digest, imager, knitr, mapproj, rgl, rmarkdown, spelling, testthat, NbClust, xml2

Suggested by photobiologyInOut.

See at CRAN