comparison Code/Descriptors/Matlab/MPEG7/FromWeb/VoiceSauce/func_SnackPitch.m @ 4:92ca03a8fa99 tip

Update to ICASSP 2013 benchmark

+function [F0, V, err] = func_SnackPitch(wavfile, windowsize, frameshift, maxF0, minF0)
+% [F0, V, err] = func_SnackPitch(wavfile, windowsize, frameshift, maxF0, minF0)
+% Input:  wavfile - input wav file
+%         windowlength - windowlength in seconds
+%         frameshift - in seconds
+%         maxF0/minF0 - max and min thresholds
+% Output: F0, V - F0 and voicing vectors
+%         err - error flag
+% Notes:  If on Windows/PC, this function will call up the Snack executable
+% in the Windows folder. Otherwise, it will try to call up Snack - this
+% require Snack to be installed beforehand.
+%
+% Author: Yen-Liang Shue, Speech Processing and Auditory Perception Laboratory, UCLA
+% Copyright UCLA SPAPL 2009
+
+
+%settings
+iwantfilecleanup = 1;  %delete tcl,f0,frm files when done
+tclfilename = 'tclforsnackpitch.tcl'; %produce tcl file to call snack
+    
+% make the wavfile acceptable to tcl in Snack
+wavfile = strrep(wavfile, '[', '\[');
+wavfile = strrep(wavfile, ']', '\]');
+wavfile = ['"' wavfile '"'];
+%wavfile = strrep(wavfile, '\', '\\');   % for PC
+
+if (ispc)  % pc can run snack.exe
+    wavfile = strrep(wavfile, '\', '\\');    
+    err = system(['Windows\snack.exe pitch ' wavfile ' -method esps -framelength ' num2str(frameshift) ' -windowlength ' num2str(windowsize) ' -maxpitch ' num2str(maxF0) ' -minpitch ' num2str(minF0)]);
+
+else % for macs and possibly others
+    %cmd = 'tclsh'; % for older Macs and 'nixes
+    cmd = 'wish8.4'; % seems to work for OSX Snow Leopard
+    %cmd = 'wish84';  % for pc
+    
+    fid = fopen(tclfilename, 'wt');
+
+    fprintf(fid, '#!/bin/sh\n');
+    fprintf(fid, '# the next line restarts with wish \\\n');
+    fprintf(fid, 'exec wish8.4 "$0" "$@"\n\n');
+    fprintf(fid, 'package require snack\n\n');
+    fprintf(fid, 'snack::sound s\n\n');
+    fprintf(fid, 's read %s\n\n', wavfile);
+    fprintf(fid, 'set fd [open [file rootname %s].f0 w]\n', wavfile);
+    fprintf(fid, 'puts $fd [join [s pitch -method esps -framelength %f -windowlength %f -maxpitch %d -minpitch %d] \\n]\n', frameshift, windowsize, maxF0, minF0);
+    fprintf(fid, 'close $fd\n\n');
+    fprintf(fid, 'exit');
+    fclose(fid);
+    
+    err = system([cmd ' ' tclfilename]);
+    
+end
+
+% change back into original file
+wavfile = strrep(wavfile, '\[', '[');
+wavfile = strrep(wavfile, '\]', ']');
+pitchfile = [wavfile(2:end-4) 'f0'];  % this should be where the F0's are stored
+
+% oops, exit now
+if (err ~= 0)
+    F0 = NaN;
+    V = NaN;
+    if (iwantfilecleanup)
+        pitchfile = strrep(pitchfile, '\\', '\');
+        if (exist(pitchfile, 'file') ~= 0)
+            delete(pitchfile);
+        end
+        
+        if (exist(tclfilename, 'file') ~= 0)
+            delete(tclfilename);
+        end
+    end
+    return;
+end
+    
+% snack executed successfully, read out the F0 values
+[F0, V, dummy1, dummy2] = textread(pitchfile, '%f %f %f %f');
+
+if (iwantfilecleanup)
+    pitchfile = strrep(pitchfile, '\\', '\');
+    if (exist(pitchfile, 'file') ~= 0)
+        delete(pitchfile);
+    end
+
+    if (exist(tclfilename, 'file') ~= 0)
+        delete(tclfilename);
+    end
+end
+
+