function[Edge,Numberofaera] = pcnn(X)
% X:輸入的灰度圖像，Edge:為檢測到的一些邊界點
% Numberofaera:表明了在各次迭代激活的塊區域
X = double(imread('C:\shi.BMP'));
ADD = imnoise(X,'salt & pepper', 0.1);
%X=ADD./max(max(ADD))
figure,imshow(ADD);
Weight = [0.5 1 0.5;1 0 1;0.5 1 0.5];
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))~=256*256)
    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-1,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));
            if U>=Threshold(i0-1,i1-1)|Threshold(i0-1,i1-1)<60
                T(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;
                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本身是激發的象素，現在邊界被置零，也不能說是白激發了
        %B矩陣有記錄！當然下次不能再被激發了
        Y(find(Edge<0)) = 0;
        [ 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
figure,imshow(Numberofaera_1);