function [Edge,Numberofaera]=pcnn(X)
%X:輸入的灰度圖像，Edge:檢測到的一些邊界點，Numberofaera則表明了在各次迭代時激活的塊區域
X=double(imread('lena128.jpg'));
%X=double(X)
Weight=[0.5 1 0.5;1 0 1;0.5 1 0.5];% quan
Beta=0.4   %
Yuzhi=240
Decay=0.3   %衰減
[a,b]=size(X)
Threshold=zeros(a,b);  %門限值
S=zeros(a+2,b+2);
B=zeros(a,b);%標記樣板，表示該Pixel是否被激活
Y=zeros(a,b);
Edge=zeros(a,b);
Numberofaera=zeros(a,b);
Numberofaera_1=zeros(a,b);
Num_1=0;
Num=0;
n=1

while(sum(sum(B))~=128*128) %若采用128×128圖像
    for i0=2:a+1
         for i1=2:b+1
              V=[S(i0-1,i1-1)  S(i0-1,i1) S(i0-1,i1+1);
                S(i0,i1-1)   S(i0,i1)   S(i0,i1+1);
                S(i0+1,i1-1) S(i0+1,i1) S(i0+1,i1+1)];
              L=sum(sum(V.*Weight));
              F=X(i0-1,i1-1);
              U=double(F)*(1+Beta*double(L));  %double可以不加
                 if U>=Threshold(i0-1,i1-1)|Threshold(i0-1,i1-1)<60
                    T(i0-1,i1-1)=1;
                    Threshold(i0-1,i1-1)=Yuzhi;
                    Y(i0-1,i1-1)=1;
                     if n==1
                        B(i0-1,i1-1)=0;%避免第一次全部激發造成的影響
                        else
                        B(i0-1,i1-1)=1;%已發射過的標記
                        Threshold(i0-1,i1-1)=1000000;%相當于不會第二次激活
                     end
                     else
                       T(i0-1,i1-1)=0;%no use?
                       Y(i0-1,i1-1)=0;
             
                 end
            end
         end
        Threshold(find(B~=1))=exp(-Decay)*Threshold(find(B~=1));
        %被激活過的象素不再參與迭代過程
        if n~=1
           Edge=Edge+judge_edge(Y,n);
           Y(find(Edge<0))=0;%邊界點被置零，Y本來是激發的象素，現在邊界被置0，也不能說白激發了，B矩陣有記錄！當然下次就休想再被激發了
           [Numberofaera_1,Num_1]=bwlabel(Y,4);
          for i=1:a
              for j=1:b
                 if Numberofaera_1(i,j)~=0 
                    Numberofaera_1(i,j)=Numberofaera_1(i,j)+Num;
                    end
               end
               end
       Numberofaera=Numberofaera+Numberofaera_1;
       Num=Num_1;
     end
     if n==1
        S=zeros(a+2,b+2);
        else
        S=Bianhuan(T);
     end
     n=n+1;
     Numberofaera_1=zeros(a,b);
     end