annotate demo/demo_aim2006_dcgc.m @ 4:537f939baef0 tip

various bug fixes and changed copyright message

% Script to compare the AIM2006 output for four files (listed below)
% This scrips creates six figure windows, one for each stage of the AIM
% processing model and plots four panes in each window, one for each file

% (c) 2006-2008, University of Cambridge, Medical Research Council 
% Tom Walters ([e-mail])
% http://www.pdn.cam.ac.uk/cnbh/aim2006
% $Date: 2008-06-10 18:00:16 +0100 (Tue, 10 Jun 2008) $
% $Revision: 585 $

filenames{1}='aa110p122s100t.wav';
filenames{2}='aa256p122s100t.wav';
filenames{3}='aa110p089s100t.wav';
filenames{4}='aa256p089s100t.wav';


analysis_start_time=0;
analysis_end_time=0.2;

display_start_time=0.13;
display_end_time=0.154;

bmm_scaling=0.25;
nap_scaling=1;
sai_scaling=0.25;
mellin_scaling=16;

for ii=1:length(filenames)
    name=filenames{ii};
    
    % Find the pitch and scale of the input sound from the filename
    inp=strfind(name,'p');
    pitch=str2num(name(inp-3:inp-1));
    ins=strfind(name,'s');
    scale=str2num(name(ins-3:ins-1));
    
    % Set the parameters for analysis of the sound
    dcgc_options=setparams_dcgc(filenames{ii}, analysis_start_time, analysis_end_time);
    
    %dcgc_output{ii}.filename=name;
    %dcgc_output{ii}.scale=scale;
    %dcgc_output{ii}.pitch=pitch;
    
    % Analyse with aim2006
    dcgc_output{ii}=aim(dcgc_options);
    
    % dcgc_output{ii}.filename=name;
    % dcgc_output{ii}.scale=scale;
    % dcgc_output{ii}.pitch=pitch;
    
    options.minimum_time=display_start_time;
    options.maximum_time=display_end_time;
    
    title_string=['Pitch:' num2str(pitch) ' Scale:' num2str(scale)];
    
    figure(1)
    set(gcf, 'Name', ['PCP: ' dcgc_output{ii}.info.current_pcp_module]);
    subplot(2,2,ii);
    plot(dcgc_output{ii}.data.pcp, [display_start_time display_end_time]);
    
    figure(2)
    set(gcf, 'Name', ['BMM: ' dcgc_output{ii}.info.current_bmm_module]);
    subplot(2,2,ii);
    plot_bmm(dcgc_output{ii}, options, bmm_scaling,title_string);
    
    figure(3)
    set(gcf, 'Name', ['NAP: ' dcgc_output{ii}.info.current_nap_module]);
    subplot(2,2,ii);
    plot_nap(dcgc_output{ii},options, nap_scaling,title_string);
    
    figure(4)
    set(gcf, 'Name', ['SP: ' dcgc_output{ii}.info.current_strobes_module]);
    subplot(2,2,ii);
    %plot(dcgc_output{ii}.data.nap);
    %hold on;
    plot_strobes(dcgc_output{ii}, options, nap_scaling,title_string);
    
    max_frame=length(dcgc_output{ii}.data.sai);
    
    figure(5)
    set(gcf, 'Name', ['SAI: ' dcgc_output{ii}.info.current_sai_module]);
    subplot(2,2,ii);
    plot_sai(dcgc_output{ii},max_frame, sai_scaling,title_string);
    
    figure(6)
    set(gcf, 'Name', ['Usermodule: ' dcgc_output{ii}.info.current_usermodule_module]);
    subplot(2,2,ii);
    plot_mellin(dcgc_output{ii},max_frame, mellin_scaling, title_string);
    
end

for ii=1:length(filenames)
    mellin_images(:,:,ii)=getvalues(dcgc_output{ii}.data.usermodule{1,max_frame});
    compressed_image=mellin_images.^(0.5);
    mi_mean=mean(compressed_image, 3);
end


mi_mean_fr=dcgc_output{ii}.data.usermodule{max_frame};
mi_mean_fr=setvalues(mi_mean_fr, mi_mean);
input.data.usermodule{1}=mi_mean_fr;
input.all_options.usermodule=dcgc_output{ii}.all_options.usermodule;
title_string='Mean mellin image';

figure;
plot_mellin(input,1, mellin_scaling, title_string);