close all; %关闭当前所有图形窗口，清空工作空间变量，清除工作空间所有变量
clc;
I=imread('D:\迅雷下载\weiqi.jpg');
imshow(I);
J=rgb2gray(I);
J=medfilt2(J,[3 3]);
figure,imshow(J);title('灰度图');
K=imhist(J);%计算直方图
bar(K);%画直方图
[Y,X]=imhist(J);
yy1=smooth(Y,5);
%plot(X,yy1);%得到直方图曲线
[width,height]=size(J);
T1=100;
for i=1:width
for j=1:height
if(J(i,j)<T1)
BW1(i,j)=255;
else
BW1(i,j)=0;
end
end
end
BW1=imfill(BW1,'holes');
figure;imshow(BW1),title('人工阈值进行分割');
se=strel('disk',2);
BW1=imerode(BW1,se);
Q=~BW1;
Q=imfill(Q,'holes');
Q=imclearborder(Q,4);%抑制和图像边界相连的亮对象
figure,imshow(Q);title('最终图像');
T2=graythresh(J);
%BW2=im2bw(J,T2);%Otus阈值进行分割
%figure;imshow(BW2),title('Otus阈值进行分割');
image = imread('D:\迅雷下载\weiqi.jpg');
imager = image(:,:,1);
Z=imager;
imshow(Z);
M=imhist(Z);%计算直方图
bar(M);%画直方图
[width,height]=size(Z);
T1=225;
for i=1:width
for j=1:height
if(J(i,j)>T1)
BW2(i,j)=255;
else
BW2(i,j)=0;
end
end
end
figure;imshow(BW2),title('人工阈值进行分割');
BW1 = bwmorph(BW1,'close'); %闭运算
imshow(BW1)
BW2 = bwmorph(BW2,'open'); %开运算
figure, imshow(BW2)
[m,n]=size(BW2);
% 求垂直投影
for y=1:n
S(y)=sum(I(1:m,y));
end
y=1:n;
figure
subplot(211),plot(y,S(y));
title('垂直投影');
% 求水平投影
for x=1:m
S(x)=sum(I(x,:));
end
x=1:m;
subplot(212),plot(x,S(x));
title('水平投影');
N=imadd(BW1,BW2);%将2个图像相加
figure;imshow(N);
N= bwmorph(N,'open'); %开运算
figure;imshow(N);
%N = bwmorph(N,'close'); %闭运算
%figure;imshow(N);
[m,n]=size(N);
% 求垂直投影
for y=1:n
S(y)=sum(I(1:m,y));
end
y=1:n;
S(y)=smooth(S(y),10);
figure;plot(y,S(y));
title('垂直投影');
% 求水平投影
for x=1:m
S(x)=sum(I(x,:));
end
x=1:m;
S(x)=smooth(S(x),10);
figure;plot(x,S(x));
title('水平投影');
se = strel('disk',2);
N= imerode(N,se);
figure;imshow(N);title('腐蚀图像');
评论0