rar~.zip_firm6od_matlab_matlab图像分割_图像分割_遗传算法图像

%%%利用最佳直方图熵法（KSW熵法）及传统遗传算法实现灰度图像二阈值分割 %%%主程序 %%初始部分，读取图像及计算相关信息 clear; close all; clc; I=imread('cman.bmp'); figure imshow(I); I=rgb2gray(I); % I=imread('Lenna.bmp'); hist=imhist(I);%显示图像数据柱状图 total=0; for i=0:255 total=total+hist(i+1); end hist1=hist/total;%求每点的归一化//求像素为i的概率Pi %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %程序主干部分 %种群随机初始化，种群数取20，染色体二进制编码取16位 t0=clock; population=20; X00=round(rand(1,population)*255); X01=round(rand(1,population)*255); for i=1:population X0(i,:)=[X00(i) X01(i)]; end for i=1:population if X0(i,1)>X0(i,2) temp=X0(i,1); X0(i,2)=temp; X0(i,1)=temp; else end adapt_value0(i)=ksw_2(X0(i,1),X0(i,2),0,255,hist1); end adapt_average0=mean(adapt_value0); X1=X0; adapt_value1=adapt_value0; adapt_average1=adapt_average0; %循环搜索，搜索代数取100 generation=100; for k=1:generation s1=Selection_2d(X1,adapt_value1); s_code10=dec2bin(s1(:,1),8); s_code11=dec2bin(s1(:,2),8); [c10,c11]=Crossover_2d(s_code10,s_code11); [v10,v11]=mutation_2d(c10,c11); X20=(bin2dec(v10))'; X21=(bin2dec(v11))'; for i=1:population X2(i,:)=[X20(i) X21(i)]; end for i=1:population adapt_value2(i)=ksw_2(X2(i,1),X2(i,2),0,255,hist1); end adapt_average2=mean(adapt_value2); if abs(adapt_average2-adapt_average1)<=0.03 break; else X1=X2; adapt_value1=adapt_value2; adapt_average1=adapt_average2; end end max_value=max(adapt_value2); number=find(adapt_value2==max_value); opt=X2(number(1),:); t1=clock; search_time=etime(t1,t0); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%阈值分割及显示部分 I_temp1=I; [height,width]=size(I_temp1) for i=1:height for j=1:width if I_temp1(i,j)<opt(1); I_temp1(i,j)=0; else if I_temp1(i,j)>opt(2); I_temp1(i,j)=255; else I_temp1(i,j)=180; end end end end I1= I_temp1; disp('灰度图像阈值分割的效果如图所示：'); %disp('源图为：Fifure No.1'); %disp('最佳直方图熵法及传统遗传算法阈二值分割后的图像为：Fifure No.2'); figure(1); imshow(I); title('源图'); figure(2); imshow(I1); title('ksw遗传算法二阈值分割后的图像'); disp('最佳直方图熵法及传统遗传算法二阈值为(s,t)：'); disp(opt(1)); disp(opt(2)); disp('最佳直方图熵法及传统遗传算法二阈值搜索所用时间（s）：'); disp(search_time); %%程序结束