function [Dx,Dy,D_matrix] = image_limitvalue( image1,image2,image3)
%%input:
%image1,image2,image3要比較的三幅圖
%
%output:
% Dx,Dy 三幅圖每26個像素比較后極值的坐標
% D_matrix極值點生成的矩陣

[m,n]=size(image1);
D_matrix=[];
Dy=[];
Dx=[];
 x=0;
 y=0;
for i=2:m-1
   for j=2:n-1
      k1=i-1:i+1;
    k2=j-1:j+1;
    
  max2= max(max(image2(k1,k2)));
   max1= max(max(image1(k1,k2)));
    max3= max(max(image3(k1,k2)));
   B=[max1,max2,max3];
    MAX=max(B);
    
    min2= min(min(image2(k1,k2)));
   min1= min(min(image1(k1,k2)));
    min3= min(min(image3(k1,k2)));
    S=[min1,min2,min3];
   MIN=min(S);
   
       if (image2(i,j)==MAX) | (image2(i,j)==MIN)
           x=x+1
          y=y+1
          Dx(x)=i
           Dy(y)=j;
          
  else
          x=x
           y=y;
   end 
end
end



for x=1:length(Dx)
    D_matrix(Dx(x),Dy(x))=image2(Dx(x),Dy(x));
end