%%中值

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"); %儲存影象