clear all
I=imread('E:\圖像處理\DIP\dip9\CT.jpg');
[M,N]=size(I);
figure(1);imshow(I);title('原始圖像');
y=mean(x,1);
yz1=(y(1)+y(2))/2 ; %求均值,單閾值分割
I1=I;
for i=1:M
for j=1:N
if I1(i,j)<yz1
I2(i,j)=I1(i,j);
else
I3(i,j)=I1(i,j);
end
end
end
figure(2);imshow(I2,[]);title('小於均值的圖像');
figure(3);imshow(I3,[]);title('大於均值的圖像');