www.gusucode.com > mean-shift用于图像分割 Matlab程序 > Ms_segmenter/M_shift1.m
function [Mode, Number_values] = M_shift1 (Values, Window_radius, Initial) % assumes that Values contains cumulative information, i.e. % the number of elements of the particular value in each Values cell Current_val = round(Initial+2); Next_val = round(Initial); L = size(Values, 2); % global plot_number plot_xmin plot_xmax plot_ymin plot_ymax % if plot_number == 1 % figure; % end % II = find(Values > 0); % if plot_number == 1 % plot_xmin = min(II)-Window_radius/2; % plot_xmax = max(II)+Window_radius/2; % plot_ymin = min(Values)-Window_radius/2; % plot_ymax = max(Values)+Window_radius/2; % end % subplot(2, 2, plot_number); % plot(Values, '.'); % xlabel('feature space (intensity)'); % ylabel('density'); % axis([plot_xmin plot_xmax plot_ymin plot_ymax]) % text(Next_val+Window_radius, Values(Next_val)-Window_radius/4, 'Initial'); while abs(Next_val-Current_val) > 1 % Draw_rect(1.8*Window_radius, Window_radius, [Next_val, Values(Next_val)], 'magenta'); Current_val = Next_val; Range = [max(1, Current_val-Window_radius): ... min(L, Current_val+Window_radius)]; Window = Values(Range); Sum_values = full(Window * Range'); Number_values = full(sum(Window)); if Number_values > 0 Next_val = round(Sum_values/Number_values); end end Mode = Next_val; % text(Next_val+Window_radius, Values(Next_val)+Window_radius/4, 'Final'); % plot_number = plot_number + 1;