FCMLSM 基于FCM和水平集的图像分割

基于FCM和水平集的图像分割，以脑颅图像为例，分割结果良好

function [imglsoptssls]=FCLSM（imgimgfcmbeta）

img=double（img）;

se=5;       %template radius for spatial filtering
sigma=2;    %spatial filter weight
d0=.5;      %fuzzy thresholding
epsilon=1.5;    %Dirac regulator

%adaptive definition of penalizing item mu
u=d0<=imgfcm;
bwa=bwarea（u）;  %area of initial contour
bw2=bwperim（u）;
bwp=sum（sum（bw2））;  %peripherium of initial contour
mu=bwp/bwa;     %Coefficient of the internal （penalizing） energy term P（\phi）;
opts.mu=mu;
timestep=0.2/mu; %The product timestep*mu must be less than 0.25 for stability
opts.timestep=timestep;
%end

% fs=fspecial（‘disk‘5）;
% imgs=imfilter（imgfs‘replicate‘）;
fs=fspecial（‘gaussian‘sesigma）;
img_smooth=conv2（double（img）double（fs）‘same‘）;
[IxIy]=gradient（img_smooth）;
f=Ix.^2+Iy.^2;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2011-04-16 23:04  FCMLSM\
     文件     1150012  2011-04-16 01:35  FCMLSM\CBM2011.pdf
     文件       51080  2009-10-30 05:41  FCMLSM\ctlivertumor.bmp
     文件        3102  2009-12-04 15:16  FCMLSM\FCLSM.m
     文件        3564  2011-04-16 23:03  FCMLSM\fuzzyLSM.m
     文件       51078  2009-06-22 01:39  FCMLSM\mrihead.bmp
     文件         980  2011-04-16 23:04  FCMLSM\runtest.m
     文件        2455  2009-10-30 20:23  FCMLSM\SFCM2D.m
     文件        1332  2011-07-16 01:02  license.txt