<< Chapter < Page Chapter >> Page >
This module describes how a matched filter algorithm is used to identify birdcall audio files.

Matched filter implementation

Matched filter

Matched filters do an excellent job of identifying sound samples, so we decided to apply the method here to identify birdcall audiofiles. A matched filter searches for a sample clip, the filter, within a longer audio recording. Convolution compares the filter to the longer signal at eachpossible offset. The greater the maximum amplitude of the convolution result, the stronger the match. By having a different filter for each birdcall, we cansearch an audio file to identify which birdcall it contains.

The matched filter algorithm is as follows:

  • Reverse the filters in the time-domain.
  • Normalize the energy of each of the filters.
  • Convolve each filter with the input signal and take the maximum amplitude of the resulting convolution signals.
  • The filter that gives us the greatest maximum value indicates which birdcall the signal contains.

Filter library creation

Our first step in implementing the matched filter algorithm was to create a library of birdcall filters. To do this, we looked atthe spectrograms of a few sample audio files of the same birdcall and selected a portion that looked representative of the call.

For each birdcall, one of the representative audio segments was saved as a filter. Because the first two steps of the above matchedfilter algorithm affect only the filter library and are independent of the input signal, we reversed the filters and normalized their energy before saving themto wave files.

Matlab implementation

The following MATLAB script performed our matched filter algorithm. When given a wave file as input, it would tell us how well theaudio sample matched against each of the 6 birdcall filters.

We used circular convolution instead of linear convolution because it computed much faster. MATLAB's built-in cconv functionzero-pads the two signals before multiplying their FFTs, generating the convolution result we are looking for.
function result = birdcheck(file) [sig, fs, nbits]= wavread(file); signal=sig(:,1);signal=signal/max(abs(signal)); filters{1}=wavread('filters/bob.wav');filters{2}=wavread('filters/lt.wav'); filters{3}=wavread('filters/lw.wav');filters{4}=wavread('filters/pygmy.wav'); filters{5}=wavread('filters/red.wav');filters{6}=wavread('filters/redcry.wav'); for i=1:6filter=filters{i}; result(i) = max(abs(cconv(signal,filter(end:-1:1))));end end

The script was able to correctly identify several birdcalls. It did fail to correctly identify four cases in twocategories:

  • Two of our loon tremolo files registered as pygmy owl common songs.
  • Two of our red-tailed hawk shriek files registered as red-tailed hawk cries.

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Birdcall identification project. OpenStax CNX. Dec 17, 2008 Download for free at http://cnx.org/content/col10616/1.2
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Birdcall identification project' conversation and receive update notifications?

Ask