I=imread('lena.bmp'); %讀原圖
J1=imnoise(I,'salt & pepper',0.02); %加均值為0,方差為0.02的椒鹽噪聲
J2=imnoise(I,'gaussian',0.02); %加均值為0,方差為0.02的高斯噪聲。
subplot(2,2,1),imshow(J1); %顯示有椒鹽噪聲圖像
subplot(2,2,2),imshow(J2); %顯示有高斯噪聲圖像
I1= medfilt2(J1,[5,5]); %對有椒鹽噪聲圖像進行5×5方形窗口中值濾波
I2= medfilt2(J2,[5,5]); %對有高斯噪聲圖像進行5×5方形窗口中值濾波
subplot(2,2,3),imshow(I1); %顯示有椒鹽噪聲圖像的濾波結果
subplot(2,2,4),imshow(I2); %顯示有高斯噪聲圖像的濾波結果
%%均衡
I=imread('lena.bmp'); %將圖讀入到I
imshow(I); %顯示圖像
figure,imhist(I); %直方圖
I2=histeq(I); %均衡化
figure;imshow(I2); %處理後圖像顯示
figure;imhist(I2); %均衡化後直方圖
imwrite(I2,'lena2.bmp'); %保存圖像