求大神解救啊！！！
求每句语言是意思~clear all;
close all;
Img = imread('twocells.bmp');% The same cell image in the paper is used here
Img=double(Img(:,:,1));
sigma=1.5; % scale parameter in Gaussian kernel for smoothing.
G=fspecial('gaussian',15,sigma);
Img_smooth=conv2(Img,G,'same');% smooth image by Gaussiin convolution
=gradient(Img_smooth);
f=Ix.^2+Iy.^2;
g=1./(1+f);% edge indicator function.
epsilon=1.5; % the papramater in the definition of smoothed Dirac function
timestep=5;% time step
mu=0.04;% coefficient of the internal (penalizing) energy term P(\phi)
% Note: The product timestep*mu must be less than 0.25 for stable evolution
lambda=5; % coefficient of the weighted length term Lg(\phi)
alf=1.5; % coefficient of the weighted area term Ag(\phi);
% Note: Choose a positive(negative) alf if the initial contour is outside(inside) the object.
% define initial level set function (LSF) as -c0, c0 at points outside and inside of a region R, respectively.
=size(Img);
c0=2; % The constant value used to define binary level set function as initial LSF;
% Using larger value of c0 usually slow down the evolution.
initialLSF=c0*ones(nrow,ncol);
w=8;
initialLSF(w+1:end-w, w+1:end-w)=-c0;
u=initialLSF;
figure;imagesc(Img, );colormap(gray);hold on;
= contour(u,,'r');
title('Initial contour');
% start level set evolution
for n=1:200
u=EVOLUTION_LSD(u, g ,lambda, mu, alf, epsilon, timestep, 1);
if mod(n,20)==0
pause(0.001);
imagesc(Img, );colormap(gray);hold on;
= contour(u,,'r');
iterNum=;
title(iterNum);
hold off;
end
end
imagesc(Img, );colormap(gray);hold on;
= contour(u,,'r');
totalIterNum=;
title(['Final contour, ', totalIterNum]);
figure;
mesh(u);
title('Final level set function');
页:
[1]