view Code/Collation/get_Pitch.m @ 4:92ca03a8fa99 tip

Update to ICASSP 2013 benchmark

function [] = get_Pitch( dirName, statsFileID )

DEBUG=0;
% this function collates the results of all HNR values

% identify the speaker in the stats file
fprintf( statsFileID, '%s ', dirName );

metricName = '_YIN_pitch';
metricFileName = [ dirName metricName '.txt'];
[ metrics, nanValues ] = readValuesFromFile( statsFileID, metricFileName, metricName );

%---------------- GET THE VOICED FRAME VALUES -------------------

% only wish to consider values from voiced frames.

% going to disregard all unvoiced frames as these produce outliers
[vuv] = detect_VoicedUnvoiced( [ dirName '.wav' ], 0 );
% if vuv is longer than the pitch array, truncate.
if( length(vuv) ~= length(nanValues) )
    if( (length(vuv) - length(nanValues)) > 1 )
        disp('EEEK!');
    end
    if( length(vuv) > length(nanValues) )
        %increase metrics
        nanValues = [ nanValues; 0 ];
    else
        vuv = [vuv; 0];
    end
end


% remove any NaN pitch values from the voicing information
i=1;
while( i <= length(vuv) )

    if( nanValues(i) == 1 )
        if( i == 1 )
            vuv = vuv(2:length(vuv));
            nanValues = nanValues(2:length(nanValues));
        elseif(i == length(vuv) )
            vuv = vuv(1:length(vuv)-1);
            nanValues = nanValues(1:length(nanValues)-1);
        else
            vuv = [ vuv(1: i-1); vuv(i+1:end) ];
            nanValues = [ nanValues(1: i-1); nanValues(i+1:end) ];
        end
    else
        i=i+1;
    end
    
end
    
% remove the non-voiced frames (includes unvoiced and silent frames) from the metrics

voicedFrames = find( vuv == 1 );
metrics = metrics( voicedFrames );


if(DEBUG)
    subplot(211);hold on;
    if( strfind( dirName, 'neg') ) 
        plot(metrics,'b');
    elseif( strfind( dirName, 'pos') ) 
        plot(metrics,'r');
    end
end

% get basic metric set
[ metricSet ] = basicMetricSet( metrics, statsFileID );

if(DEBUG)
    subplot(212); 
    if( strfind( dirName, 'neg') )
        plot(metrics, 'xb'); hold on;
    elseif( strfind( dirName, 'pos') )
        plot(metrics, 'xr'); hold on;
    end
end


% get basic metric set for the first derivative
firstDerivative = gradient( metrics );
[ firstDerivativeMetricSet ] = basicMetricSet( firstDerivative, statsFileID );

% get basic metric set for the second derivative
secondDerivative = gradient( firstDerivative );
[ secondDerivativeMetricSet ] = basicMetricSet( secondDerivative, statsFileID );

fprintf( statsFileID, '\n');