befproc.m

这个程序能清楚显示灰度图形的灰度直方图

%图像预处理
I=imread('E:\工件图像\裂纹\D10.4Xd5.5X1(1)\1.bmp','bmp');
%I=imread('E:\工件图像\掉角\D10.4Xd5.5X1\1.bmp','bmp');
%I=imread('E:\工件图像\掉角\D10.4Xd5.5X1\11.bmp','bmp');
%I=imread('E:\工件图像\划痕\D10.4Xd5.5X1(1)\1.bmp','bmp');
%I=imread('E:\工件图像\麻坑\D10.4Xd5.5X1\10.bmp','bmp');

figure,imshow(I);
I=medfilt2(I);
figure,imshow(I);
figure,imhist(I);
aw=imhist(I);
%消除直方图噪声
bw=aw;
bw=bw';
x=1:256;
for j=1:5
  for i=1:255
     if abs(bw(i)-bw(i+1))>1
         bw(i+1)=fix(bw(i)+bw(i+1))/2;
     end
  end
end


x=1:1:256;
figure;
plot(x,bw,'-');
ylim([-200,3000]);
xlim([0,250]);
grid on

J=imadjust(I,[0.1 0.9], [0 1],1);
%J=imadjust(I,[0 1], [1 0],1);
%J=imadjust(J,[0 1], [1 0],1);
figure,imshow(J);
%figure,imhist(J);
%%%%%%%%%%%%%%%%%%%%%%%%%555
aw=imhist(J);
%消除直方图噪声
bw=aw;
bw=bw';
x=1:256;
for j=1:5
  for i=1:255
     if abs(bw(i)-bw(i+1))>1
         bw(i+1)=fix(bw(i)+bw(i+1))/2;
     end
  end
end
x=1:1:256;
figure;
plot(x,bw,'-');
ylim([-200,3000]);
xlim([0,250]);
grid on

I1=im2bw(I,70/255);
figure;
imshow(I1);