提取峰值

郁闷了好几天，怎样求峰值呢，
昨天回去的路上才子说：可以用三个点。
晕，突然想到啦，这么弱智的问题居然困扰了我好久。
下面是求峰值的代码，可能有点乱。
求出峰值，下面的步骤就比较简单啦。

a = imread( ' c:oat.png');   % 读入图像
figure,imshow(a);  %显示图像
b = a(:);
k = size(b);
%  下述操作可以实现像素的显示
c = linspace( 0 , 0 , 256 );
for  i = 1 :k
     for  j = 1 : 256
         if (b(i) == j)
            c(j) = c(j) + 1 ;
         end
     end
end
%  以下步骤简单实现读取整个图像的峰值点

max = c( 1 );
maxi = 0 ;
for  i = 1 : 256
     if (c(i) > max)
        max = c(i);
        maxi = i;
     end
end
max
maxi




%  下面步骤实现对各个峰值的读取
d = c;  %作为临时操作数组，以不破坏原数组
emax = linspace( 0 , 0 , 256 );
emaxi = linspace( 0 , 0 , 256 );
a = 1 ;
for  j = 2 : 256
         if (d(j) > d(j - 1 )  & d(j) > d(j + 1 ))
            emax(a) = d(j);
            emaxi(a) = j;
            d(j) = d(j - 1 );
            a = a + 1 ;   %开始没有想明白，用额外的循环变量
         end
end


%  下面的程序读取四个高峰值，并进行显示
fmax = linspace( 0 , 0 , 4 );
fmaxi = linspace( 0 , 0 , 4 );
b = 1 ;
for  i = 1 : 4
    fmax(b) = emax( 1 );
     for  j = 2 : 256
         if (emax(j) > fmax(b))
            fmax(b) = emax(j);
            fmaxi(b) = emaxi(j);
            t = j;
         end
     end
    emax(t) = 0 ;
    b = b + 1 ;
end


for  i = 1 : 4
    fmax(i)
    fmaxi(i)
end