area_between_zeros.m

来自「The aip file contains few Matlab routine」· M 代码 · 共 42 行

function [str_index,end_index ]  = Area_Between_Zeros(input_arr);

%%%%%%%% Finding the area between each zero crossings for COL proj hist %%%%%%%%%%%%%%%%
% vect_value = find(input_arr ==0);
% 
% for j = 1 :size(vect_value,2)-1
%     
%     zero_str =  vect_value(1,j); 
%     zero_end =  vect_value(1,j+1);
%     
%     area_bw_zeros = input_arr(zero_str:zero_end);
%     
%     mean_value(j) = mean(area_bw_zeros);
%     sigma_value(j) = std(area_bw_zeros);
%     
%     weighted_value(j) = 0.75 * ( 2 * sigma_value(j)) + 0.25 * ( mean_value(j));
% end
% 
% [max_value,max_index] = max(weighted_value) ; 
% 
% str_index = vect_value(max_index);
% end_index = vect_value(max_index+1);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

k = 1;
for i = 1 : length(input_arr)

    if( input_arr(i) > 0 )
      non_zero_index_array(k) = i;
      k = k+1;
    end;
    
end;

str_index = min(non_zero_index_array)
end_index = max(non_zero_index_array)

keyboard;