matlab PCNN 显微图像处理

%------Pulse Coupled Neural Network-----------

%------名称：脉冲耦合神经网络（PCNN）

function H = PCNN2(I)

I = imread(‘BanLanGen.png‘);

figure,imshow(I);

title(‘原始图像‘);

th = graythresh(I);

I = im2bw(I,th);

%H = Pcnn(I,30);

%************************************************************

% Initialize Function

%************************************************************

[m,n] = size (I);

figure,imshow(I);

title(‘二值图像‘);

X = im2double(I);

%************************************************************

% Initialize PCNN Parameters

%************************************************************

% P ：L E F VF VL VE beta

% V: 1.0 1.0 0.1 0.5 0.2 20 0.1

al = 1.0; ae = 1.0; af = 0.1; vf = 0.5; vl = 0.2; ve = 20; B = 0.1;

W =[0.5 1 0.5;...

1 0 1;...

0.5 1 0.5];

M = W; Y = zeros(m,n); F = Y; L = Y; U = Y; E = Y;

%************************************************************

% PCNN 点火过程

%************************************************************

for i = 1:30

wk = conv2(Y,M,‘same‘);

F = exp(-af).* F + vf.* wk + X;

L = exp(-al).* L + vl.* wk;

U = F.*(1 + B.* L);

Y = double(U>E);

E = exp(-ae).* E + ve.* Y;

end

figure,imshow(Y)

title(‘PCNN处理后图像‘);

fs = fspecial(‘sobel‘);

filter_Y = imfilter(Y, fs);

figure,imshow(filter_Y);

title(‘滤波图像‘);