Skip to content
This repository was archived by the owner on Feb 7, 2024. It is now read-only.

[SIGNAL] Added peak calibration #222

Closed
adrian-stepien wants to merge 4 commits into from
Closed

[SIGNAL] Added peak calibration #222

Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
ba46cda
[SIGNAL] Added peak calibration
adrian-stepien Dec 7, 2018
e9fabe9
[SIGNAL] Renamed calibrate_peak to calibrate_with_peak
adrian-stepien Feb 17, 2020
67c5b75
Merge branch 'python-acoustics:master' into dev/calibrate_peak
adrian-stepien Apr 21, 2022
cc2be50
Merge branch 'python-acoustics:master' into dev/calibrate_peak
adrian-stepien Sep 15, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
18 changes: 17 additions & 1 deletion acoustics/_signal.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -5,7 +5,7 @@
from scipy.signal import detrend, lfilter, bilinear, spectrogram, filtfilt, resample, fftconvolve
import acoustics

from acoustics.standards.iso_tr_25417_2007 import REFERENCE_PRESSURE
from acoustics.standards.iso_tr_25417_2007 import REFERENCE_PRESSURE, sound_pressure_level
from acoustics.standards.iec_61672_1_2013 import WEIGHTING_SYSTEMS
from acoustics.standards.iec_61672_1_2013 import (NOMINAL_OCTAVE_CENTER_FREQUENCIES,
NOMINAL_THIRD_OCTAVE_CENTER_FREQUENCIES)
Expand Down Expand Up @@ -119,6 +119,22 @@ def calibrate_with(self, other, decibel, inplace=False):
gain = decibel - other.leq()
return self.gain(gain, inplace=inplace)

def calibrate_with_peak(self, decibel, inplace=False):
"""Calibrate signal knowing its SPL for peak FS.

:param decibel: Peak FS value to calibrate to.
:param inplace: Whether to perform inplace or not.
:returns: Calibrated signal.
:rtype: :class:`Signal`

Values of `decibel` are broadcasted. To set a value per channel, use `decibel[...,None]`.
"""
# Calculate the gain so that samples represent the pressure values
# We know that peak value equals 1. 1 Pa is roughly 94 dB SPL.
# In that case, we have to boost the signal by decibel - ~94.
gain = decibel * np.ones(self.shape) - sound_pressure_level(1)
return self.gain(gain, inplace=inplace)

def decimate(self, factor, zero_phase=False, ftype='iir', order=None):
"""Decimate signal by integer `factor`. Before downsampling a low-pass filter is applied.

Expand Down