分形几何

开发平台：

Matlab

tztqex.txt：源码内容

							  %读取待处理图像
 I=imread('4-1.jpg'); 
 imshow(I);
 %对输入的图像I进行形态学开操作，去掉不完全包含在圆盘中的对象
   bw=imopen(I,strel('line',10,45));
   %反白处理
 bw=imcomplement(bw);
 %将图像转化为二值图像
 bw=im2bw(I,0.75);
figure,imshow(bw);
 %细化（骨架化）
bw=bwperim(bw);   
 figure,imshow(bw);
%去掉象素小于10的细小物体
bw = bwareaopen(bw, 10);
 figure,imshow(bw);
 %清除与边界连通的物体
 bw=imclearborder(bw,4);
 figure,imshow(bw);
 %确定图像中纤维的个数，标示二值图像bw中所有相关部分，并且返回在图像中找到的对象的个数
 [labeled,numObjects]=bwlabel(bw,8);
%测量对象的属性
 graindata=regionprops(labeled,'Area');
 allgrains=[graindata.Area]
 allgrains =266  
            183  
            200   
            148
  M=max(allgrains)
 N=min(allgrains)
 mean(allgrains)
num=size(allgrains)
%绘制一个包含500个柱的直方图来说明纤维大小的分布情况
 hist(allgrains,500)
xlabel('Pixel');
ylabel('Distribution');
 title (['Total Fibre Numbers are:',num2str(num(2))],'color','r')
 %打印图表
 %print-depsc-tiff-r200myplot 
%保存命令窗口数据
 diary filename.m  

						