|
Daniel@0
|
1 /* Part of DML (Digital Music Laboratory)
|
|
Daniel@0
|
2 Copyright 2014-2015 Samer Abdallah, University of London
|
|
Daniel@0
|
3
|
|
Daniel@0
|
4 This program is free software; you can redistribute it and/or
|
|
Daniel@0
|
5 modify it under the terms of the GNU General Public License
|
|
Daniel@0
|
6 as published by the Free Software Foundation; either version 2
|
|
Daniel@0
|
7 of the License, or (at your option) any later version.
|
|
Daniel@0
|
8
|
|
Daniel@0
|
9 This program is distributed in the hope that it will be useful,
|
|
Daniel@0
|
10 but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
Daniel@0
|
11 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
Daniel@0
|
12 GNU General Public License for more details.
|
|
Daniel@0
|
13
|
|
Daniel@0
|
14 You should have received a copy of the GNU General Public
|
|
Daniel@0
|
15 License along with this library; if not, write to the Free Software
|
|
Daniel@0
|
16 Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
|
|
Daniel@0
|
17 */
|
|
Daniel@0
|
18
|
|
Daniel@0
|
19 :- module(dsp,
|
|
Daniel@0
|
20 [ spectrogram/2
|
|
Daniel@0
|
21 , spectrogram/3
|
|
Daniel@0
|
22 , spectrum/2
|
|
Daniel@0
|
23 , sum/2
|
|
Daniel@0
|
24 ]).
|
|
Daniel@0
|
25
|
|
Daniel@0
|
26 :- use_module(library(semweb/rdf_db)).
|
|
Daniel@0
|
27 :- use_module(library(sandbox)).
|
|
Daniel@0
|
28 :- use_module(library(rdfutils)).
|
|
Daniel@0
|
29 :- use_module(library(memo)).
|
|
Daniel@0
|
30 :- use_module(library(mlserver)).
|
|
Daniel@0
|
31
|
|
Daniel@0
|
32 :- set_prolog_flag(double_quotes,string).
|
|
Daniel@0
|
33
|
|
Daniel@0
|
34 :- setting(specgram_dynamic_range,number,90,"Dynamic range in dB for spectrograms").
|
|
Daniel@0
|
35
|
|
Daniel@0
|
36 spectrogram(URI,X) :-
|
|
Daniel@0
|
37 rdf_text(URI,beets:path,P),
|
|
Daniel@0
|
38 X===specgrm( monofile(P), hanning(2048),512).
|
|
Daniel@0
|
39
|
|
Daniel@0
|
40 :- volatile_memo spectrum(+atom,-ground).
|
|
Daniel@0
|
41 spectrum(URI,X) :-
|
|
Daniel@0
|
42 spectrogram(URI,Y),
|
|
Daniel@0
|
43 X1===sum(Y,2),
|
|
Daniel@0
|
44 persist_item(X1,X).
|
|
Daniel@0
|
45
|
|
Daniel@0
|
46 spectrogram(URI,Offset,Length) :-
|
|
Daniel@0
|
47 setting(specgram_dynamic_range,DBs),
|
|
Daniel@0
|
48 ( cp_audio:audio_file(URI,P,Fmt),
|
|
Daniel@0
|
49 ( Fmt=just(aac)
|
|
Daniel@0
|
50 -> MSignal=sndfile(q(P),q(enc),q(aac))
|
|
Daniel@0
|
51 ; MSignal=sndfile(q(P))
|
|
Daniel@0
|
52 )
|
|
Daniel@0
|
53 ; cp_audio:audio_link(URI,P,just(mp3)),
|
|
Daniel@0
|
54 MSignal=mp3file(q(P))
|
|
Daniel@0
|
55 ), !,
|
|
Daniel@0
|
56 ?? (specgrm( taket(Length,dropt(Offset,mixdown(MSignal))),
|
|
Daniel@0
|
57 hanning(2048),512,"range",DBs,"offset",Offset);
|
|
Daniel@0
|
58 ylim([0,15])).
|
|
Daniel@0
|
59
|
|
Daniel@0
|
60 sum(Items,Sum) :-
|
|
Daniel@0
|
61 once(member(X,Items)),
|
|
Daniel@0
|
62 D===numdims(X)+1,
|
|
Daniel@0
|
63 Sum===sum(cellcat(D,cell(Items)),D).
|
|
Daniel@0
|
64
|
|
Daniel@0
|
65 sandbox:safe_primitive(dsp:spectrogram(_,_,_)).
|
