comparison Code/Collation/get_AudioPower.m @ 4:92ca03a8fa99 tip

Update to ICASSP 2013 benchmark

+function [] = get_AudioPower( 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 = '_AP';
+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 );
+unvoicedFrames = find( vuv == 2 );
+
+% find the average energies for voiced and unvoiced speech
+voicedEnergy = sum( metrics( voicedFrames )) / length( voicedFrames );
+unvoicedEnergy = sum( metrics( unvoicedFrames )) / length( unvoicedFrames );
+
+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 );
+
+% print the average voiced and unvoiced energy to file
+fprintf( statsFileID, '\t %f \t %f ', voicedEnergy, unvoicedEnergy);
+
+fprintf( statsFileID, '\n');