acoular 23.6__py3-none-any.whl → 24.3__py3-none-any.whl

acoular/tprocess.py

@@ -40,7 +40,7 @@
 from numpy import array, empty, empty_like, pi, sin, sqrt, zeros, newaxis, unique, \
 int16, nan, concatenate, sum, float64, identity, argsort, interp, arange, append, \
 linspace, flatnonzero, argmin, argmax, delete, mean, inf, asarray, stack, sinc, exp, \
-polymul, arange, cumsum, ceil, split
+polymul, arange, cumsum, ceil, split, array_equal
 
 from numpy.linalg import norm
 from numpy.matlib import repmat
@@ -59,7 +59,7 @@ import numba as nb
 from datetime import datetime
 from os import path
 import wave
-from scipy.signal import butter, lfilter, filtfilt, bilinear
+from scipy.signal import butter, filtfilt, bilinear, tf2sos, sosfilt, sosfiltfilt
 from warnings import warn
 from collections import deque
 from inspect import currentframe
@@ -957,7 +957,7 @@ class SpatialInterpolator(TimeInOut):
         # Interpolation for 1D Arrays 
         if self.array_dimension =='1D' or self.array_dimension =='ring':
             #for rotation add phidelay
-            if not phiDelay == []:
+            if not array_equal(phiDelay,[]):
                 xInterpHelp = repmat(virtNewCoord[0, :], nTime, 1) + repmat(phiDelay, virtNewCoord.shape[1], 1).T
                 xInterp = ((xInterpHelp + pi ) % (2 * pi)) - pi #  shifting phi cootrdinate into feasible area [-pi, pi]
             #if no rotation given
@@ -1001,7 +1001,7 @@ class SpatialInterpolator(TimeInOut):
         # Interpolation for arbitrary 2D Arrays
         elif self.array_dimension =='2D':
             #check rotation
-            if not phiDelay == []:
+            if not array_equal(phiDelay,[]):
                 xInterpHelp = repmat(virtNewCoord[0, :], nTime, 1) + repmat(phiDelay, virtNewCoord.shape[1], 1).T
                 xInterp = ((xInterpHelp + pi ) % (2 * pi)) - pi #shifting phi cootrdinate into feasible area [-pi, pi]
             else:
@@ -1081,7 +1081,7 @@ class SpatialInterpolator(TimeInOut):
         # Interpolation for arbitrary 3D Arrays             
         elif self.array_dimension =='3D':
             #check rotation
-            if not phiDelay == []:
+            if not array_equal(phiDelay,[]):
                 xInterpHelp = repmat(virtNewCoord[0, :], nTime, 1) + repmat(phiDelay, virtNewCoord.shape[1], 1).T
                 xInterp = ((xInterpHelp + pi  ) % (2 * pi)) - pi  #shifting phi cootrdinate into feasible area [-pi, pi]
             else:
@@ -1294,7 +1294,10 @@ class Mixer( TimeInOut ):
         for temp in self.source.result(num):
             sh = temp.shape[0]
             for g in gens:
-                temp1 = next(g)
+                try:
+                    temp1 = next(g)
+                except StopIteration:
+                    return
                 if temp.shape[0] > temp1.shape[0]:
                     temp = temp[:temp1.shape[0]]
                 temp += temp1[:temp.shape[0]]
@@ -1455,10 +1458,10 @@ class Filter(TimeInOut):
     Should not be instanciated by itself
     """
     #: Filter coefficients
-    ba = Property()
+    sos = Property()
 
-    def _get_ba( self ):
-        return [1],[1]
+    def _get_sos( self ):
+        return tf2sos([1],[1])
 
     def result(self, num):
         """ 
@@ -1477,20 +1480,17 @@ class Filter(TimeInOut):
             Delivers the bandpass filtered output of source.
             The last block may be shorter than num.
         """
-        b, a = self.ba
-        zi = zeros((max(len(a), len(b))-1, self.source.numchannels))
+        sos = self.sos
+        zi = zeros((sos.shape[0], 2, self.source.numchannels))
         for block in self.source.result(num):
-            b, a = self.ba # this line is useful in case of changes 
-                            # to self.ba during generator lifetime
-            block, zi = lfilter(b, a, block, axis=0, zi=zi)
+            sos = self.sos # this line is useful in case of changes 
+                            # to self.sos during generator lifetime
+            block, zi = sosfilt(sos, block, axis=0, zi=zi)
             yield block
 
-class FiltFiltOctave( TimeInOut ):
+class FiltOctave( Filter ):
     """
-    Octave or third-octave filter with zero phase delay.
-    
-    This filter can be applied on time signals.
-    It requires large amounts of memory!   
+    Octave or third-octave filter (causal, non-zero phase delay).    
     """
     #: Band center frequency; defaults to 1000.
     band = Float(1000.0, 
@@ -1499,33 +1499,28 @@ class FiltFiltOctave( TimeInOut ):
     #: Octave fraction: 'Octave' or 'Third octave'; defaults to 'Octave'.
     fraction = Trait('Octave', {'Octave':1, 'Third octave':3}, 
         desc = "fraction of octave")
+
+    #: Filter order
+    order = Int(3, desc = "IIR filter order")
         
+    sos = Property( depends_on = ['band', 'fraction', 'source.digest', 'order'])
+
     # internal identifier
     digest = Property( depends_on = ['source.digest', '__class__', \
-        'band', 'fraction'])
+        'band', 'fraction','order'])
 
     @cached_property
     def _get_digest( self ):
         return digest(self)
         
-    def ba(self, order):
-        """ 
-        Internal Butterworth filter design routine.
-        
-        Parameters
-        ----------
-        order : integer
-            The order of the filter.
-        
-        Returns
-        -------
-            b, a : ndarray, ndarray
-                Filter coefficients.
-        """
+    @cached_property
+    def _get_sos( self ):
         # filter design
         fs = self.sample_freq
-        # adjust filter edge frequencies
-        beta = pi/(4*order)
+        # adjust filter edge frequencies for correct power bandwidth (see ANSI 1.11 1987
+        # and Kalb,J.T.: "A thirty channel real time audio analyzer and its applications",
+        # PhD Thesis: Georgia Inst. of Techn., 1975
+        beta = pi/(2*self.order)
         alpha = pow(2.0, 1.0/(2.0*self.fraction_))
         beta = 2 * beta / sin(beta) / (alpha-1/alpha)
         alpha = (1+sqrt(1+beta*beta))/beta
@@ -1534,9 +1529,44 @@ class FiltFiltOctave( TimeInOut ):
             raise ValueError("band frequency too high:%f,%f" % (self.band, fs))
         om1 = fr/alpha 
         om2 = fr*alpha
-#        print om1, om2
-        return butter(order, [om1, om2], 'bandpass') 
-        
+        return butter(self.order, [om1, om2], 'bandpass', output = 'sos') 
+
+class FiltFiltOctave( FiltOctave ):
+    """
+    Octave or third-octave filter with zero phase delay.
+    
+    This filter can be applied on time signals.
+    It requires large amounts of memory!   
+    """
+    #: Filter order (applied for forward filter and backward filter)
+    order = Int(2, desc = "IIR filter half order")
+
+    # internal identifier
+    digest = Property( depends_on = ['source.digest', '__class__', \
+        'band', 'fraction','order'])
+
+    @cached_property
+    def _get_digest( self ):
+        return digest(self)
+ 
+    @cached_property
+    def _get_sos( self ):
+        # filter design
+        fs = self.sample_freq
+        # adjust filter edge frequencies for correct power bandwidth (see FiltOctave)
+        beta = pi/(2*self.order)
+        alpha = pow(2.0, 1.0/(2.0*self.fraction_))
+        beta = 2 * beta / sin(beta) / (alpha-1/alpha)
+        alpha = (1+sqrt(1+beta*beta))/beta
+        # additional bandwidth correction for double-pass
+        alpha = alpha * {6:1.01,5:1.012,4:1.016,3:1.022,2:1.036,1:1.083}.get(self.order,1.0)**(3/self.fraction_)
+        fr = 2*self.band/fs
+        if fr > 1/sqrt(2):
+            raise ValueError("band frequency too high:%f,%f" % (self.band, fs))
+        om1 = fr/alpha 
+        om2 = fr*alpha
+        return butter(self.order, [om1, om2], 'bandpass', output = 'sos')   
+           
     def result(self, num):
         """
         Python generator that yields the output block-wise.
@@ -1554,63 +1584,22 @@ class FiltFiltOctave( TimeInOut ):
             Delivers the zero-phase bandpass filtered output of source.
             The last block may be shorter than num.
         """
-        b, a = self.ba(3) # filter order = 3
+        sos = self.sos 
         data = empty((self.source.numsamples, self.source.numchannels))
         j = 0
         for block in self.source.result(num):
             ns, nc = block.shape
             data[j:j+ns] = block
             j += ns
+        # filter one channel at a time to save memory
         for j in range(self.source.numchannels):
-            data[:, j] = filtfilt(b, a, data[:, j])
+            data[:, j] = sosfiltfilt(sos, data[:, j])
         j = 0
         ns = data.shape[0]
         while j < ns:
             yield data[j:j+num]
             j += num
 
-
-class FiltOctave( Filter ):
-    """
-    Octave or third-octave filter (causal, non-zero phase delay).    
-    """
-    #: Band center frequency; defaults to 1000.
-    band = Float(1000.0, 
-        desc = "band center frequency")
-        
-    #: Octave fraction: 'Octave' or 'Third octave'; defaults to 'Octave'.
-    fraction = Trait('Octave', {'Octave':1, 'Third octave':3}, 
-        desc = "fraction of octave")
-
-    #: Filter order
-    order = Int(3, desc = "IIR filter order")
-        
-    ba = Property( depends_on = ['band', 'fraction', 'source.digest', 'order'])
-
-    # internal identifier
-    digest = Property( depends_on = ['source.digest', '__class__', \
-        'band', 'fraction','order'])
-
-    @cached_property
-    def _get_digest( self ):
-        return digest(self)
-        
-    @cached_property
-    def _get_ba( self ):
-        # filter design
-        fs = self.sample_freq
-        # adjust filter edge frequencies
-        beta = pi/(4*self.order)
-        alpha = pow(2.0, 1.0/(2.0*self.fraction_))
-        beta = 2 * beta / sin(beta) / (alpha-1/alpha)
-        alpha = (1+sqrt(1+beta*beta))/beta
-        fr = 2*self.band/fs
-        if fr > 1/sqrt(2):
-            raise ValueError("band frequency too high:%f,%f" % (self.band, fs))
-        om1 = fr/alpha 
-        om2 = fr*alpha
-        return butter(self.order, [om1, om2], 'bandpass') 
-
 class TimeExpAverage(Filter):
     """
     Computes exponential averaging according to IEC 61672-1
@@ -1622,7 +1611,7 @@ class TimeExpAverage(Filter):
     weight = Trait('F', {'F':0.125, 'S':1.0, 'I':0.035}, 
         desc = "time weighting")    
 
-    ba = Property( depends_on = ['weight', 'source.digest'])
+    sos = Property( depends_on = ['weight', 'source.digest'])
        
     # internal identifier
     digest = Property( depends_on = ['source.digest', '__class__', \
@@ -1633,11 +1622,11 @@ class TimeExpAverage(Filter):
         return digest(self)
         
     @cached_property
-    def _get_ba( self ):
+    def _get_sos( self ):
         alpha = 1-exp(-1/self.weight_/self.sample_freq)
         a = [1, alpha-1]
         b = [alpha]
-        return b,a 
+        return tf2sos(b,a)
 
 class FiltFreqWeight( Filter ):
     """
@@ -1646,17 +1635,18 @@ class FiltFreqWeight( Filter ):
     #: weighting characteristics
     weight = Trait('A',('A','C','Z'), desc="frequency weighting")
 
-    ba = Property( depends_on = ['weight', 'source.digest'])
+    sos = Property( depends_on = ['weight', 'source.digest'])
 
     # internal identifier
-    digest = Property( depends_on = ['source.digest', '__class__'])
+    digest = Property( depends_on = ['source.digest', '__class__', \
+        'weight'])
 
     @cached_property
     def _get_digest( self ):
         return digest(self)
 
     @cached_property
-    def _get_ba( self ):
+    def _get_sos( self ):
         # s domain coefficients
         f1 = 20.598997
         f2 = 107.65265
@@ -1677,7 +1667,7 @@ class FiltFreqWeight( Filter ):
             b = zeros(7)
             b[0] = 1.0
             a = b # 6th order flat response
-        return b,a
+        return tf2sos(b,a)
 
 class FilterBank(TimeInOut):
     """
@@ -1688,7 +1678,7 @@ class FilterBank(TimeInOut):
     """
 
     #: List of filter coefficients for all filters
-    ba = Property()
+    sos = Property()
 
     #: List of labels for bands
     bands = Property()
@@ -1699,8 +1689,8 @@ class FilterBank(TimeInOut):
     #: Number of bands
     numchannels = Property()
 
-    def _get_ba( self ):
-        return [[1]],[[1]]
+    def _get_sos( self ):
+        return [tf2sos([1],[1])]
 
     def _get_bands( self ):
         return ['']
@@ -1729,13 +1719,13 @@ class FilterBank(TimeInOut):
         """
         numbands = self.numbands
         snumch = self.source.numchannels
-        b, a = self.ba
-        zi = [zeros( (max(len(a[0]), len(b[0]))-1, snumch)) for _ in range(numbands)]
+        sos = self.sos
+        zi = [zeros( (sos[0].shape[0],2, snumch)) for _ in range(numbands)]
         res = zeros((num,self.numchannels),dtype='float')
         for block in self.source.result(num):
             bl = block.shape[0]
             for i in range(numbands):
-                res[:,i*snumch:(i+1)*snumch], zi[i] = lfilter(b[i], a[i], block, axis=0, zi=zi[i])
+                res[:,i*snumch:(i+1)*snumch], zi[i] = sosfilt(sos[i], block, axis=0, zi=zi[i])
             yield res
 
 class OctaveFilterBank(FilterBank):
@@ -1780,15 +1770,14 @@ class OctaveFilterBank(FilterBank):
         return len(self.bands)
 
     @cached_property
-    def _get_ba( self ):
+    def _get_sos( self ):
         of = FiltOctave(source=self.source, fraction=self.fraction)
-        b, a = [], []
+        sos = []
         for i in range(self.lband,self.hband,4-self.fraction_):
             of.band = 10**(i/10)
-            b_,a_ = of.ba
-            b.append(b_)
-            a.append(a_)
-        return b, a
+            sos_ = of.sos
+            sos.append(sos_)
+        return sos
 
 class TimeCache( TimeInOut ):
     """

acoular/version.py

@@ -5,5 +5,5 @@
 
 # separate file to find out about version without importing the acoular lib
 __author__ = "Acoular Development Team"
-__date__ = "26 June 2023"
-__version__ = "23.6"
+__date__ = "11 March 2024"
+__version__ = "24.03"

acoular-24.3.dist-info/METADATA

@@ -0,0 +1,181 @@
+Metadata-Version: 2.1
+Name: acoular
+Version: 24.3
+Summary: Python library for acoustic beamforming
+Project-URL: homepage, https://acoular.org
+Project-URL: documentation, https://acoular.org
+Project-URL: repository, https://github.com/acoular/acoular
+Author-email: Acoular Development Team <info@acoular.org>
+Maintainer-email: Adam Kujawski <adam.kujawski@tu-berlin.de>, Art Pelling <a.pelling@tu-berlin.de>, Ennes Sarradj <ennes.sarradj@tu-berlin.de>, Gert Herold <gert.herold@tu-berlin.de>, Mikolaj Czuchaj <mikolaj.czuchaj@tu-berlin.de>, Simon Jekosch <s.jekosch@tu-berlin.de>
+License: Copyright (c) Acoular Development Team.
+        All rights reserved.
+        
+        
+        Redistribution and use in source and binary forms, with or without
+        modification, are permitted provided that the following conditions are met:
+        
+          a. Redistributions of source code must retain the above copyright notice,
+             this list of conditions and the following disclaimer.
+          b. Redistributions in binary form must reproduce the above copyright
+             notice, this list of conditions and the following disclaimer in the
+             documentation and/or other materials provided with the distribution.
+          c. Neither the name of the acoular developers nor the names of
+             its contributors may be used to endorse or promote products
+             derived from this software without specific prior written
+             permission. 
+        
+        
+        THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+        AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+        IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+        ARE DISCLAIMED. IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE FOR
+        ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+        DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
+        SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
+        CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
+        LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
+        OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH
+        DAMAGE.
+License-File: AUTHORS.rst
+License-File: LICENSE
+Keywords: acoustics,beamforming,microphone array
+Classifier: Development Status :: 5 - Production/Stable
+Classifier: Intended Audience :: Education
+Classifier: Intended Audience :: Science/Research
+Classifier: License :: OSI Approved :: BSD License
+Classifier: Programming Language :: Python :: 3.7
+Classifier: Programming Language :: Python :: 3.8
+Classifier: Programming Language :: Python :: 3.9
+Classifier: Programming Language :: Python :: 3.10
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Topic :: Scientific/Engineering :: Physics
+Requires-Python: <=11,>=3.7
+Requires-Dist: numba
+Requires-Dist: numpy
+Requires-Dist: scikit-learn
+Requires-Dist: scipy>=1.1.0
+Requires-Dist: tables>=3.4.4
+Requires-Dist: traits>=6.0
+Provides-Extra: dev
+Requires-Dist: graphviz; extra == 'dev'
+Requires-Dist: ipython; extra == 'dev'
+Requires-Dist: matplotlib; extra == 'dev'
+Requires-Dist: nox; extra == 'dev'
+Requires-Dist: numpydoc; extra == 'dev'
+Requires-Dist: pickleshare; extra == 'dev'
+Requires-Dist: sounddevice; extra == 'dev'
+Requires-Dist: sphinx; extra == 'dev'
+Provides-Extra: full
+Requires-Dist: matplotlib; extra == 'full'
+Requires-Dist: pylops; extra == 'full'
+Requires-Dist: sounddevice; extra == 'full'
+Description-Content-Type: text/markdown
+
+![Acoular Logo](./docs/source/_static/Acoular_logo.png)
+
+[![PyPI](https://img.shields.io/pypi/pyversions/acoular.svg)](https://pypi.org/project/acoular)
+[![PyPI](https://img.shields.io/pypi/v/acoular.svg)](https://pypi.org/project/acoular)
+[![Github](https://github.com/acoular/acoular/actions/workflows/python-package.yml/badge.svg)](https://github.com/acoular/acoular/actions/workflows/python-package.yml)
+
+# Acoular
+Acoular is a Python module for acoustic beamforming that is distributed under the new BSD license. 
+
+It is aimed at applications in acoustic testing. Multichannel data recorded by a microphone array can be processed and analyzed in order to generate mappings of sound source distributions. The maps (acoustic photographs) can then be used to locate sources of interest and to characterize them using their spectra. 
+
+# Features
+- frequency domain beamforming algorithms: delay & sum, Capon (adaptive), MUSIC, functional beamforming, eigenvalue beamforming
+- frequency domain deconvolution algorithms: DAMAS, DAMAS+, Clean, CleanSC, orthogonal deconvolution
+- frequency domain inverse methods: CMF (covariance matrix fitting), general inverse beamforming, SODIX
+- time domain methods: delay & sum beamforming, CleanT deconvolution
+- time domain methods applicable for moving source with arbitrary trajectory (linear, circular, arbitrarily 3D curved), 
+- frequency domain methods for rotating sources via virtual array rotation for arbitrary arrays and with different interpolation techniques
+- 1D, 2D and 3D mapping grids for all methods
+- gridless option for orthogonal deconvolution
+- four different built-in steering vector formulations
+- arbitrary stationary background flow can be considered for all methods
+- efficient cross spectral matrix computation
+- flexible modular time domain processing: n-th octave band filters, fast, slow, and impulse weighting, A-, C-, and Z-weighting, filter bank, zero delay filters
+- time domain simulation of array microphone signals from fixed and arbitrarily moving sources in arbitrary flow
+- fully object-oriented interface
+- lazy evaluation: while processing blocks are set up at any time, (expensive) computations are only performed when needed
+- intelligent and transparent caching: computed results are automatically saved and loaded on the next run to avoid unnecessary re-computation
+- parallel (multithreaded) implementation with Numba for most algorithms
+- easily extendable with new algorithms
+
+# License
+Acoular is licensed under the BSD 3-clause. See [LICENSE](LICENSE)
+
+# Citing
+
+If you use Acoular for academic work, please consider citing our
+[publication](https://doi.org/10.1016/j.apacoust.2016.09.015):
+
+    Ennes Sarradj, Gert Herold,
+    A Python framework for microphone array data processing,
+    Applied Acoustics, Volume 116, 2017, Pages 50-58
+
+# Dependencies
+Acoular runs under Linux, Windows and MacOS and needs Numpy, Scipy, Traits, scikit-learn, pytables, Numba packages available. 
+Matplotlib is needed for some of the examples.
+
+If you want to use input from a soundcard hardware, you will also need to install the [sounddevice](https://python-sounddevice.readthedocs.io/en/0.3.12/installation.html) package. Some solvers for the CMF method need [Pylops](https://pylops.readthedocs.io/en/stable/installation.html).
+
+# Installation
+
+Acoular can be installed via [conda](https://docs.conda.io/en/latest/), which is also part of the [Anaconda Python distribution](https://www.anaconda.com/). It is recommended to install into a dedicated [conda environment](https://docs.conda.io/projects/conda/en/latest/user-guide/tasks/manage-environments.html). After activating this environment, run
+
+    conda install -c acoular acoular
+
+This will install Acoular in your Anaconda Python enviroment and make the Acoular library available from Python. In addition, this will install all dependencies (those other packages mentioned above) if they are not already present on your system. 
+
+A second option is to install Acoular via [pip](https://pip.pypa.io/en/stable/). It is recommended to use a dedicated [virtual environment](https://virtualenv.pypa.io/en/latest/) and then run
+
+    pip install acoular
+
+For more detailed install instructions see the [documentation](http://acoular.org/install/index.html).
+
+# Documentation and help
+Documentation is available [here](http://acoular.org) with a
+[getting started](http://acoular.org/get_started/index.html) section and
+[examples](http://acoular.org/examples/index.html).
+
+The Acoular [blog](https://acoular.github.io/blog/) contains some tutorials.
+
+Problems, suggestions and success using Acoular may be reported via the [acoular-users](https://groups.google.com/forum/#!forum/acoular-users) discussion forum.
+
+# Example
+This reads data from 64 microphone channels and computes a beamforming map for the 8kHz third octave band:
+
+```python
+from os import path
+import acoular
+from matplotlib.pylab import figure, plot, axis, imshow, colorbar, show
+
+# this file contains the microphone coordinates
+micgeofile = path.join(path.split(acoular.__file__)[0],'xml','array_64.xml')
+# set up object managing the microphone coordinates
+mg = acoular.MicGeom( from_file=micgeofile )
+# set up object managing the microphone array data (usually from measurement)
+ts = acoular.TimeSamples( name='three_sources.h5' )
+# set up object managing the cross spectral matrix computation
+ps = acoular.PowerSpectra( time_data=ts, block_size=128, window='Hanning' )
+# set up object managing the mapping grid
+rg = acoular.RectGrid( x_min=-0.2, x_max=0.2, y_min=-0.2, y_max=0.2, z=0.3, \
+increment=0.01 )
+# set up steering vector, implicitely contains also the standard quiescent 
+# environment with standard speed of sound
+st = acoular.SteeringVector( grid = rg, mics=mg )
+# set up the object managing the delay & sum beamformer
+bb = acoular.BeamformerBase( freq_data=ps, steer=st )
+# request the result in the 8kHz third octave band from approriate FFT-Lines
+# this starts the actual computation (data intake, FFT, Welch CSM, beamforming)
+pm = bb.synthetic( 8000, 3 )
+# compute the sound pressure level
+Lm = acoular.L_p( pm )
+# plot the map
+imshow( Lm.T, origin='lower', vmin=Lm.max()-10, extent=rg.extend(), \
+interpolation='bicubic')
+colorbar()
+```
+
+![result](./docs/source/get_started/three_source_py3_colormap.png)