Skip to content

Reproduce graphical soundscape analysis for a sampling site.

License

Notifications You must be signed in to change notification settings

PEM-Humboldt/graph_soundscapes

Repository files navigation

Graphical soundscape analysis

Visualization of the patterns of acoustic activity displayed throughout the course of a day, or "graphical soundscape", is a popular way to characterize environmental audio recordings. This repository provides a simple step by step procedure to apply this analysis using the approach proposed by Campos-Cerqueira et al., 2017, Forumo et al., 2019, and Campos-Cerqueira et al., 2017.

Getting started

A simple step by step that walks you through the process is provided in the script simple_step_by_step.R.

Dependencies

To run the scripts you will need R version 3.6 or higher, and the following R packages:

  • seewave
  • tuneR
  • vegan: Community Ecology Package
  • ggplot
  • viridis
  • reshape2

Install the packages and load the libraries

install.packages(c("seewave", "tuneR", "vegan", "ggplot2", "viridis", "reshape2"))

1. Load packages and functions

library(viridis)
library(tuneR)
library(seewave)
library(reshape2)
library(ggplot2)
source('graph_soundscape_fcns.R')

2. Get audio metadata

Set the location of the audio recordings. For this simple example, the audio directory should have all the files of a single recording plot.

path_files = './audio/V6DA/'  # location of audio dataset

Get a file list of all .wav samples

flist = list.files(path_files, recursive = T, pattern = '.WAV', ignore.case = T)

Get metadata from raw recordings. Songmeter and Audiomoth are the currently supported audio recorders.

df = metadata_audio(flist, path_files, verbose = T, rec_model = 'SM')

Plot sampling scheme to verify the recording scheme.

ggplot(df, aes(y=hour)) + geom_bar(width=0.3, alpha=0.5) + theme_minimal()
ggplot(df, aes(y=day)) + geom_bar(width=0.3, alpha=0.5)

3. Compute and plot graphical soundscapes

gs = graphical_soundscape(df, spec_wl=256, fpeaks_th=0.003, fpeaks_f=0, verbose=T)
plot_graphical_soundscape(gs)

Example of a 24 hour graphical soundscape computed from 1108 audio recordings collected in a tropical dry forest patch in Dagua, Valle del Cauca, Colombia.

Batch process

Further scripts are provided to batch analyse data from multiple recording plots:

  • read_audio_metadata.R
  • batch_compute_graph_soundscape.R

Author

This script has been developed by Juan Sebastián Ulloa at the Instituto Humboldt, Colombia, following a previously proposed approach from other researchers (see References).

License

This project is licensed under the MIT License - see the LICENSE.txt file for details

References

  • Campos‐Cerqueira, M., et al., 2020. How does FSC forest certification affect the acoustically active fauna in Madre de Dios, Peru? Remote Sensing in Ecology and Conservation 6, 274–285. https://doi.org/10.1002/rse2.120
  • Furumo, P.R., Aide, T.M., 2019. Using soundscapes to assess biodiversity in Neotropical oil palm landscapes. Landscape Ecology 34, 911–923.
  • Campos-Cerqueira, M., Aide, T.M., 2017. Changes in the acoustic structure and composition along a tropical elevational gradient. JEA 1, 1–1. https://doi.org/10.22261/JEA.PNCO7I

About

Reproduce graphical soundscape analysis for a sampling site.

Topics

Resources

License

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages