基于matlab的Butterworth低通滤波程序分享

piday123

clear all;
6 ?/ h6 z3 ?9 C. j5 Z1 h
' Q2 Z8 W& C- p& c0 D: E/ f# U$ n% _close all;
0 m- k; K: L6 [0 F
P=input('inputimage data:');
$ D0 ~; F! }+ d/ x; k6 z
d=size(P);

4 `; _* [3 c7 |5 j/ H! c2 Iif(d(3)>1)
; ^$ V+ z- \$ u7 v4 n, U6 O
    P=rgb2gray(P);

end
" G2 O; G( k9 c: ^3 F
- i5 b6 z! {: Y5 H4 @subplot(131);imshow(P);

title('原始图像');

P1=imnoise(P,'salt& pepper',0.02);
$ V/ }& s" v% X, S4 U5 z, F4 a4 _
subplot(132)
" h5 E; o; i4 H  n* }
imshow (P1);
# u3 |2 Z$ h/ Y' j) }
title('加入椒盐噪声后的图像');
5 D# k. O  x7 T
f=double(P1);

g=fft2(f);%傅里叶变换

g=fftshift(g);%转换数据矩阵fftshift是针对频域的，将FFT的DC分量移到频谱中心
4 l" W  C7 o3 p" S1 t& A即对频域的图像，
3 n. X9 c! u2 C" F" ]1 O" n! `
: ]0 k, m2 t2 M( l9 n0 j6 U- q [M,N]=size(g);
/ j5 Y  L0 x4 ^! r
D0=input('inputnonnegative dhreshold D0=');
& N; Q9 A7 D8 ^* n& q
n=input('input theorder of filtering n=');
4 o, ^5 [9 j, Q$ ~5 C
* U! _, g& W# P# F6 Q/ w- Jn1=fix(M/2);

n2= fix(N/2);
, f$ k% s* I2 R. ?3 L0 b1 k9 b
% S  j7 k: w( F& O4 l1 n0 ^, Efor i=1:M

    for j=1:N

0 X/ R1 ]. u( B  W5 I        d=sqrt((i-n1)^2+(j-n2)^2);
3 u& q/ c( [% f: T/ E
! }% V  i; n' j" W$ _2 Ah=1/(1+(d/D0)^(2*n));
! m& @0 h9 ?' M4 K+ p* a
8 Q& O! l, V& ?* q. Y! n        result(i,j)=h*g(i,j);

    end
6 N! c& o- j$ M6 V" M
end
8 v' E9 b. e4 a! h  B- F2 `0 m
result=ifftshift(result);
; c& t: e) z2 Y, F! ]8 F. G
X2=ifft2(result);

X3=uint8(real(X2));
6 a# X  O; _* x; \3 I' [& t
' `- ]2 K6 ^* k: `- E& c( p& s6 csubplot(133)

. P* V  W9 R, f0 x" {imshow (X3);
" W  p; g. u% I
6 L9 a; h9 f/ ?title('Butterworth低通滤波所得图像');
9 O0 y" K/ g6 u8 S; i
, O) W/ n6 t3 y$ C) A9 e$ G
) P" B- d9 ]) d5 @8 e; i
: H! `  R5 Q/ \' ]

: {. z" R) I1 H/ J" A

NNNei256

基于matlab的Butterworth低通滤波程序