tests for spectral_coverage and filter_width_factor

Author: DrMarc

.github/workflows/python-package.yml

@@ -36,7 +36,7 @@ jobs:
         pip install -e .[testing]
     - name: Lint with flake8
       run: |
-        flake8 . --count --select=E9,F63,F7,F82 --show-source --statistics
+        flake8 . --count --select=E9,F63,F7 --show-source --statistics
         flake8 . --count --exit-zero --max-complexity=10 --max-line-length=127 --statistics
     - name: Test with pytest
       run: |

slab/__init__.py

@@ -2,7 +2,7 @@
 import os
 import pathlib
 
-__version__ = '1.7.0'
+__version__ = '1.8.0'
 
 # The variable _in_notebook is used to enable audio playing in Jupiter notebooks
 # and on Google colab (see slab.sound.play())

tests/test_filter.py

@@ -67,9 +67,13 @@ def test_cos_filterbank():
         length = numpy.random.randint(1000, 5000)
         low_cutoff = numpy.random.randint(0, 500)
         high_cutoff = numpy.random.choice([numpy.random.randint(5000, 15000), None])
+        filter_width_factor = numpy.random.uniform(0.1, 2.5)
         n_filters = []
         for bandwidth in numpy.linspace(0.1, 0.9, 9):
-            fbank = slab.Filter.cos_filterbank(length=length, bandwidth=bandwidth, low_cutoff=low_cutoff, high_cutoff=high_cutoff, pass_bands=False, samplerate=sound.samplerate)
+            fbank = slab.Filter.cos_filterbank(length=length, bandwidth=bandwidth,
+                                               filter_width_factor=filter_width_factor,
+                                               low_cutoff=low_cutoff, high_cutoff=high_cutoff,
+                                               pass_bands=False, samplerate=sound.samplerate)
             n_filters.append(fbank.n_filters)
             filtsound = fbank.apply(sound)
             assert filtsound.n_channels == fbank.n_filters

tests/test_sounds.py

@@ -30,6 +30,7 @@ def test_sound_generation():
     # test if .name attribute of loaded sound object is the path string
     assert loaded1.name == f"{tmpdir / 'sound.wav'}"
 
+
 def test_read_write():
     for _ in range(20):
         for normalize in [True, False]:
@@ -99,3 +100,11 @@ def test_frames():
             center1 = window[frame_dur][0]
             center2 = sound[numpy.where(sound.times == center)[0][0]][0]
             numpy.testing.assert_almost_equal(center1, center2, decimal=1)
+
+
+def test_spectral_coverage():
+    sig1 = slab.Sound.harmoniccomplex(f0=50, amplitude=[0,0,0])
+    sig2 = slab.Sound.harmoniccomplex(f0=50, amplitude=[0,0,0,0,0])
+    sig3 = slab.Sound.pinknoise()
+    assert sig1.spectral_coverage(threshold='otsu') < sig2.spectral_coverage(threshold='otsu')
+    assert sig2.spectral_coverage() < sig3.spectral_coverage()