%close?all
clear?all;clc
h=helpdlg（‘按左鍵選擇分割區(qū)域雙擊‘）;
uiwait（h）

%程序中對(duì)圖像特征最敏感的兩個(gè)參數(shù)
%T=5;%鄰域特征的域值，決定聚類中心的一個(gè)關(guān)鍵參數(shù)
r=300;%歸為該聚類中心的最小距離

prompt={‘最大循環(huán)次數(shù):‘‘鄰域灰度閾值:‘};
name=‘蟻群聚類分割‘;
numlines=1;
defaultanswer={‘10‘‘10‘};
options.Resize=‘on‘;
options.Windowstyle=‘normal‘;?
options.Interpreter=‘tex‘;
itermax1=inputdlg（promptnamenumlinesdefaultansweroptions）;
itermax=str2num（itermax1{1}）;
T=str2num（itermax1{2}）;
%itermax=10;?%



%load?aa;%aa.mat里面儲(chǔ)存了一幅灰度圖象的矩陣數(shù)據(jù)
%bb=imcrop（bb）;%從該圖像中選中一塊矩形區(qū)域，程序?qū)x中區(qū)域進(jìn)行處理，這樣可以減少運(yùn)行時(shí)間?

%[xmap]=imread（‘Dorsal?10.4.gif‘）;%這幾句讀入索引類型的圖象
%grayimg=ind2gray（xmap）;%把索引類型圖像轉(zhuǎn)換為灰度圖象
dicom1=dicomread（‘33091247‘）;
grayimg=dicom2gray（dicom1）;
%grayimg=imadjust（dicom1）;
figure;imshow（grayimg）

%grayimg=imresize（grayimg0.5）;%把圖片縮小為原來大小的0.5倍
%figure;imshow（grayimg）;%把讀入的圖片顯示出來
set（gcf‘outerposition‘get（0‘screensize‘））;%設(shè)置figure全屏顯示
[ffimgrect]=imcrop（grayimg）;

set（gcf‘outerposition‘get（0‘screensize‘））;%設(shè)置figure全屏顯示
title（‘原始圖像‘）;
%h1?=?zoom;
%set（h1‘Enable‘‘on‘）
%figure;
%imshow（ffimg）;

%w8=[111;1-81;111];
%ffimg=ffimg-imfilter（ffimgw8‘replicate‘）;

%ffimg=grayimg;
%ff=imread（‘LENA256.BMP‘）;%讀入另一幅灰度圖象對(duì)算法進(jìn)行效果測(cè)試

%ffimg=imread（‘cameraman.tif‘）;%處理256*256的圖片
%ffimg=imresize（ffimg0.5）;

%cannyedgeimg=edge（ff‘canny‘）;%對(duì)讀入的待處理灰度圖象進(jìn)行canny邊緣檢測(cè)
%figure;imshow（cannyedgeimg）%顯示邊緣檢測(cè)輸出的二值邊緣圖像?

%width=70;%設(shè)定一個(gè)正方形的區(qū)域的邊長(zhǎng)
%另一種選取圖像中一個(gè)正方形區(qū)域的方法
%ff=ff（floor（a/2）-width:floor（a/2）+width-1floor（b/2）-width:floor（b/2）+width-1）;

%[ab]=size（ff）;%輸出圖像的大小
%figure;imshow（ff）;%顯示待分割的圖像
ff=double（ffimg）;%把圖像數(shù)據(jù)類型由unit8轉(zhuǎn)化為double類型，便于后面進(jìn)行數(shù)值處理

ff?=?medfilt2（ff）;%進(jìn)行中值濾波處理，效果比原來在連續(xù)性這一點(diǎn)上要好
%figure;imshow（uint8（ff））;


%w1=fspecial（‘laplacian‘）;w2=w1‘;%采用拉普拉斯梯度算子?
%w1=fspecial（‘log‘3）;w2=w1‘;%采用log梯度算子
%w1=fspecial（‘prewitt‘）;w2=w1‘;%采用prewitt梯度算子
w1=[-1-2-1;000;121];w2=[-101;-202;-101];%采用sobel梯度算子
g=sqrt（imfilter（ffw1）.^2+imfilter（ffw2）.^2）;%計(jì)算梯度圖像
gmax=max（sum（g（:3:size（ff2）-2）））/（size（ff1）-4）;%取出最大梯度列的均值??注意:?由于邊界采用了人為填充所以剃度很大應(yīng)該忽略?

fp=padarray（ff[11]‘replicate‘‘both‘）;%對(duì)邊界進(jìn)行填充?目的是為了進(jìn)行3*3的8鄰域特征提取
?for?i=2:size（ff1）+1
??????for?j=2:size（ff2）+1
??????????a=[?fp（i-1j-1）?fp（i-1j）?fp（i-1j+1）;fp（ij-1）??fp（ij）?fp（ij+1）;?fp（i+1j-1）??fp（i+1j）?fp（i+1j+1）];
??????????ne（i-1j-1）=length（find（（a-fp（ij））??????end
?end

%下面是測(cè)試ne中小于等于6對(duì)應(yīng)的圖象中的點(diǎn)所組成的圖像邊緣
%??as1=size（ff1）;as2=size（ff2）;
%??m=zeros（as1as2）;
%??m=reshape（m1as1*as2）;
%??m（find（ne<=5））=1;
%??m=reshape（mas1as2）;
%??figure;imshow（m）

?
?%??初始化聚類中心
?f=reshape（ff1size（ff1）*size（ff2））;
?h=hist（f256）;%灰度直方圖圖象點(diǎn)統(tǒng)計(jì)
?
?counter=0;%計(jì)算出灰度直方圖上的極大值點(diǎn)，以它們作為聚類中心
?for?i=2:length（h）-1
??if?h（i-1）<=h（i）&h（i）>=h（i+1）
???counter=counter+1;
???V（counter）=i;?
??end
end
h1=sum（h（V））/length（h（V））;%取整體均值作為區(qū)分背景目標(biāo)與邊緣的區(qū)域值，即像素?cái)?shù)遠(yuǎn)遠(yuǎn)大于該整體均值的灰度作為背景或者目標(biāo)的灰度特征
h2=sum（h（find（h（V）??for?i=1:counter
??????if?h（V（i））>5*h1
?

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件??????10878??2009-06-04?15:23??LCL_ANT_COLONY\Colony_partition\Automatic\ACO_edge_detection.m

?????文件?????118946??2009-05-20?17:38??LCL_ANT_COLONY\Colony_partition\Automatic\Dorsal?10.4.gif

?????文件????????215??2009-05-20?17:38??LCL_ANT_COLONY\Colony_partition\Automatic\fname.mat

?????文件???????8356??2009-05-21?19:30??LCL_ANT_COLONY\Colony_partition\Automatic\waitbar.m

?????文件???????3431??2009-05-22?12:31??LCL_ANT_COLONY\Colony_partition\Automatic\drawedgelist.m

?????文件???????7646??2009-05-20?17:38??LCL_ANT_COLONY\Colony_partition\Automatic\edgelink.m

?????文件?????189886??2009-05-20?10:53??LCL_ANT_COLONY\Colony_partition\Automatic\33091247

?????文件????????267??2009-05-24?22:45??LCL_ANT_COLONY\Colony_partition\Automatic\a_dicom.m

?????文件????????501??2009-05-25?17:02??LCL_ANT_COLONY\Colony_partition\Automatic\dicom2gray.m

?????文件?????118946??2009-05-20?17:38??LCL_ANT_COLONY\Colony_partition\FCM\Dorsal?10.4.gif

?????文件?????????41??2009-05-20?17:38??LCL_ANT_COLONY\Colony_partition\FCM\fdd.m

?????文件???????8217??2009-05-20?17:43??LCL_ANT_COLONY\Colony_partition\FCM\waitbar.m

?????文件???????8572??2009-06-04?15:28??LCL_ANT_COLONY\Colony_partition\FCM\MYFCM.m

?????文件???????4040??2009-05-20?17:43??LCL_ANT_COLONY\Colony_partition\FCM\fcm.m

?????文件?????189886??2009-05-20?10:53??LCL_ANT_COLONY\Colony_partition\FCM\33091247

?????文件???????3388??2009-05-25?08:29??LCL_ANT_COLONY\Colony_partition\FCM\nlfilter.m

?????文件????????376??2003-01-26?01:00??LCL_ANT_COLONY\Colony_partition\FCM\mkconstarray.m

?????文件???????3329??2009-05-25?08:55??LCL_ANT_COLONY\Colony_partition\FCM\rand.asv

?????文件???????8395??2009-05-30?09:28??LCL_ANT_COLONY\Colony_partition\FCM\MYFCM.asv

?????文件??????41938??2009-05-30?09:31??LCL_ANT_COLONY\Colony_partition\FCM\lcl.mat

?????文件?????118946??2009-05-20?17:38??LCL_ANT_COLONY\Colony_partition\Semi-automatic\Dorsal?10.4.gif

?????文件???????3472??2009-05-22?12:55??LCL_ANT_COLONY\Colony_partition\Semi-automatic\drawedgelist.m

?????文件???????7646??2009-05-20?17:38??LCL_ANT_COLONY\Colony_partition\Semi-automatic\edgelink.m

?????文件???????2601??2009-05-22?09:09??LCL_ANT_COLONY\Colony_partition\Semi-automatic\edge_linking.m

?????文件???????9515??2009-05-20?17:38??LCL_ANT_COLONY\Colony_partition\Semi-automatic\source_MY_ACO_edge_detection.gif

?????文件???????8217??2009-05-21?08:33??LCL_ANT_COLONY\Colony_partition\Semi-automatic\waitbar.m

?????文件????????500??2009-05-22?08:41??LCL_ANT_COLONY\Colony_partition\Semi-automatic\READ?ME.txt

?????目錄??????????0??2009-11-12?10:17??LCL_ANT_COLONY\Colony_partition\Automatic

?????目錄??????????0??2009-11-12?10:17??LCL_ANT_COLONY\Colony_partition\FCM

?????目錄??????????0??2009-11-12?10:17??LCL_ANT_COLONY\Colony_partition\Semi-automatic

............此處省略5個(gè)文件信息