数字图像处理 matlab实现图像分割

代码

Image1=im2double(imread('1.jpg'));
gray=rgb2gray(Image1);
T=graythresh(gray);%使用 Otsu 方法获取阈值，T 被归一化到[0,1]区间。
BW=im2bw(gray,T);%以 T 为阈值把灰度图像 I 转变为二值图像。
figure,imshow(BW),title('二值化图像');
SE=strel('square',3);%创建一个由square指定的结构元素，参数3控制 SE 的大小。
Morph=imopen(BW,SE); %对图像BW进行开运算，返回图像为Morph，SE是结构元素
Morph=imclose(Morph,SE);
figure,imshow(Morph),title('形态学滤波');
[B L]=bwboundaries(1-Morph);%搜索二值图像BW的外边界和内边界
figure,imshow(L),title('划分的区域');
hold on;
for i=1:length(B)
 boundary=B{i};
 plot(boundary(:,2),boundary(:,1),'r','LineWidth',2);
end
M=zeros(length(B)); 
for k=1:length(B) 
 N=length(B{k}); 
 if N/2~=round(N/2) 
 B{k}(end+1,:)=B{k}(end,:); 
 N=N+1;
 end
 M(k)=[N*3/4]; 
end
S=zeros(size(Morph)); 
figure,imshow(S);
hold on;
for k=1:length(B) 
 z=B{k}(:,2)+1i*B{k}(:,1); 
 Z=fft(z); 
 [Y I]=sort(abs(Z)); 
 for count=1:M(k) 
 Z(I(count))=0; 
 end
 zz=ifft(Z);%对Z进行IDFT 运算。 
 plot(real(zz),imag(zz),'w'); 
end

原图：

运行结果：