view aim-mat/modules/nap/hl/gen_hl.m @ 4:537f939baef0 tip

various bug fixes and changed copyright message

% generating function for 'aim-mat'
% 
%   INPUT VALUES:
%  
%   RETURN VALUE:
%
% 
% (c) 2011, University of Southampton
% Maintained by Stefan Bleeck ([e-mail])
% download of current version is on the soundsoftware site: 
% http://code.soundsoftware.ac.uk/projects/aimmat
% documentation and everything is on http://www.acousticscale.org



function nap=gen_hl(bmm,options)

waithand=waitbar(0,'generating NAP');

nap=halfwayrectify(bmm);

vals=getvalues(nap);

%gtfb output values are: 0< gt_vals <1
%therefore we will scale for 16bit values
%operationally only 15 bits are used as we
%half wave rectify

% TCW AIM 2006 - 07/03/06 autoscaling now fixed
%vals=vals.*2.^15;
% Note: tcw Nov '05 - only re-added this part (from the log code) since otherwise it really
% screws up the display code later on. See do_aim_autoscale, cases {4,5}
% for why - the scaling seems a bit arbitrary to me!

nap=setvalues(nap,vals);

% lowpassfiltering in the end
if isfield(options,'do_lowpassfiltering') % this is multiplied to the threshold_time_const
    do_lowpassfiltering=options.do_lowpassfiltering;
else
    do_lowpassfiltering=1;
end

% cut off frequency lowpassfiltering in the end
if isfield(options,'lowpass_cutoff_frequency') % this is multiplied to the threshold_time_const
    lowpass_cutoff_frequency=options.lowpass_cutoff_frequency;
else
    lowpass_cutoff_frequency=1200;
end



% order of the lowpassfiltering in the end
if isfield(options,'lowpass_order') % this is multiplied to the threshold_time_const
    lowpass_order=options.lowpass_order;
else
    lowpass_order=2;
end


nr_chan=getnrchannels(nap);
if options.do_lowpassfiltering==1
    for ii=1:nr_chan  % through all channels: prepare working variable
		if mod(ii,10)==0
			waitbar(ii/nr_chan);
		end
        sig=getsinglechannel(nap,ii);
        newsig=leakyintegrator(sig,lowpass_cutoff_frequency,lowpass_order);
        nap=setsinglechannel(nap,ii,newsig);
    end
end



timesum=getsum(nap);
maxt=max(timesum);
nap=setscalesumme(nap,maxt);
%save nap1.mat nap;
close(waithand);