二維-高斯低通濾波器matlab程序設(shè)計

1、精選優(yōu)質(zhì)文檔-傾情為你奉上二維高斯濾波器圖像濾波MTATLAB程序設(shè)計用MATLAB設(shè)計一個3×3模板標(biāo)準(zhǔn)差為0.5，1.5，2.5的二維高斯低通濾波器(Gaussion low pass filter分別對灰度圖像，真彩色圖像，偽彩色圖像進(jìn)行濾波處理%二維高斯低通濾波器%文件為glpf.mfunction glpf(J)%處理索引圖像if isind(J) I,map=imread(J); I=imnoise(I,'gaussian',0.005); subplot(2,2,1) imshow(I,map); title('a.原始圖像'); K1=

2、filter2(fspecial('gaussian',3*3,0.5),I); K2=filter2(fspecial('gaussian',3*3,1.5),I); K3=filter2(fspecial('gaussian',3*3,2.5),I); subplot(2,2,2) imshow(K1,map); title('b.濾波器3*3,sigma=0.5') subplot(2,2,3) imshow(K2,map); title('c.濾波器3*3,sigma=1.5') subplot(2,2,4

3、) imshow(K3,map); title('d.濾波器3*3,sigma=2.5')end%處理灰度,真彩色圖像if isind(J)J=imread(J);%讀入圖像I=imnoise(J,'gaussian');%加入高斯噪聲%輸出原始圖像subplot(2,2,1)imshow(I);title('a.原始圖像')h1 = fspecial('gaussian',3,3, 0.5);%用預(yù)定義的gaussian函數(shù)m n p=size(I);%m n p m,n代表像素點的位置,p代表空間分量%p=1,代表亮度分量,即

4、灰度圖像%p=3，代表亮度分量(Y),色差分量(Cr,Cb)if p=1%處理灰度圖像 I=double(I); I=conv2(I,h1,'same');%I與h的二維離散卷積endif p=3%處理真彩色 I=double(I); I(:,:,1)=conv2(I(:,:,1),h1,'same'); I(:,:,2)=conv2(I(:,:,2),h1,'same'); I(:,:,3)=conv2(I(:,:,3),h1,'same');endI=uint8(I);%8位無符號整型%經(jīng)過3*3，sigma=0.5二維高斯低

5、通濾波器濾波后的圖像subplot(2,2,2)imshow(I);title('b.濾波器3*3,sigma=0.5')h2 = fspecial('gaussian',3,3,1.5);m n p=size(I);if p=1 I=double(I); I=conv2(I,h2,'same');%二維離散卷積endif p=3 I=double(I); I(:,:,1)=conv2(I(:,:,1),h2,'same'); I(:,:,2)=conv2(I(:,:,2),h2,'same'); I(:,:,3)

6、=conv2(I(:,:,3),h2,'same');endI=uint8(I);%經(jīng)過3*3，sigma=1.5二維高斯低通濾波器濾波后的圖像subplot(2,2,3)imshow(I);title('c.濾波器3*3,sigma=1.5')h3 = fspecial('gaussian',3,3, 2.5);m n p=size(I);if p=1 I=double(I); I=conv2(I,h3,'same');endif p=3 I=double(I); I(:,:,1)=conv2(I(:,:,1),h3,'same'); I(:,:,2)=conv2(I(:,:,2),h3,'same'); I(:,:,3)=conv2(I(:,:,3),h3,'same');endI=uint8(I);%經(jīng)過3*3，sigma=2.5二維高斯低通濾波器濾波后的圖像subplot(2,2,4)imshow(I);title('d.濾波器3*3,sigma=2.5')end程序調(diào)用%glpfUse.m%調(diào)用function glpf()處理灰度，真彩色，索引(偽彩色)圖像并繪圖glpf('填寫你需要處理圖像的路徑'