直方图规范化.docVIP

加入噪声 I=imread(dzw.jpg);%读取图像 K=rgb2gray(I); J=imnoise(K,gaussian,0,0.01);%给图像加高斯噪声 subplot(1,2,1);imshow(J); title(加入高斯噪声的图像); H = IMNOISE(K,salt pepper,0.02);%给图像加椒盐噪声 subplot(1,2,2);imshow(H); title(加入椒盐噪声的图像); clear all; f=imread(dzw.jpg); subplot(2,2,1) imshow(f);%显示原图像% title(处理前图像); f1=rgb2gray(f); %转换成灰度图 %绘制直方图% subplot(2,2,2) imhist(f1);%显示直方图% title(原图直方图) g=histeq(f1);%均衡化后的图像% subplot(2,2,3); imshow(g);%显示均衡后的图像% title(均衡后的图像); subplot(2,2,4); imhist(g);%均衡后的直方图% title(均衡后的直方图); I=imread(dzw.jpg);%读取图像 K=rgb2gray(I); J=imnoise(K,gaussian,0,0.01);%给图像加高斯噪声 subplot(1,2,1);imshow(J); title(加入高斯噪声的图像); H = IMNOISE(K,salt pepper,0.02);%给图像加椒盐噪声 subplot(1,2,2);imshow(H); title(加入椒盐噪声的图像); M1=(1/4)*[0,1,0;1,0,1;0,1,0]; M2=(1/8)*[1,1,1;1,0,1;1,1,1]; L=imfilter(J,M1); G=imfilter(J,M2); subplot(1,2,1);imshow(L); title(4邻域平均) subplot(1,2,2);imshow(G); title(8邻域平均) I=imread(dzw.jpg);%读取图像 K=rgb2gray(I); J=imnoise(K,gaussian,0,0.01);%给图像加高斯噪声 subplot(1,2,1);imshow(J); title(加入高斯噪声的图像); H = IMNOISE(K,salt pepper,0.02);%给图像加椒盐噪声 subplot(1,2,2);imshow(H); title(加入椒盐噪声的图像); M3=(1/5)*[0,1,0;1,1,1;0,1,0]; M4=(1/16)*[1,2,1;2,4,2;1,2,1]; L=imfilter(J,M3); G=imfilter(J,M4); subplot(1,2,1);imshow(L); title(1/5加权平均) subplot(1,2,2);imshow(G); title(1/16加权平均) clear all; f=imread(dzw.jpg); subplot(1,3,1) imshow(f);%显示原图像% title(原图); %w=fspecial(laplacian,0); h=fspecial(laplacian); h1=[0 -1 0;-1 5 -1;0 -1 0]; h2=[0 -2 0;-2 9 -2;0 -2 0]; g1=imfilter(f,h1,replicate); g2=imfilter(f,h2,replicate); subplot(1,3,2) imshow(g1,[]); title(锐化后的图像a=1); subplot(1,3,3); imshow(g2,[]); title(锐化后的图像a=2); figure(2); BW1=edge(f,sobel); %用SOBEL算子进行边缘检测 BW2=edge(f,roberts);%用Roberts算子进行边缘检测 BW3=edge(f,prewitt); %用prewitt算子进行边缘检测 subplot(1,3,1), imshow(BW1); title(sobel edge check); subplot(1,3,2), imshow(BW2); title(roberts edge check); subplot(1,3,3), imshow(BW3); title(prewitt edge check); clear all; I=imread(dzw.jpg); %读取图像， subplot(2,2,1);imshow(I);title(原图); %显示