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

piday123

clear all;
: x# b# Y( f& T
! s9 s5 a. R; m, v6 Pclose all;

P=input('inputimage data:');

9 N( ]! A& }( U! Q6 J! b  Vd=size(P);

; D' H2 O. H6 @  k8 E' O- D4 v$ w, Iif(d(3)>1)
! t9 r- T0 |# I6 `
2 |. A+ A. F0 ]) `; X; U, y    P=rgb2gray(P);

+ w' ?3 @4 q4 R- S2 P' w2 ~end

2 I1 M& ~" |$ {- z7 {subplot(131);imshow(P);

8 h2 ?+ c0 W1 b1 f0 O; \8 Ptitle('原始图像');
: ]" H2 L8 f' [8 Y8 ]4 a
" }( ?% [( k$ H  _0 v/ V! O7 l3 n# [P1=imnoise(P,'salt& pepper',0.02);
/ T/ l! h' w9 z
/ V5 t+ ]" P% F4 L1 h# c2 T, @subplot(132)
3 e: s' D3 v: }7 T
imshow (P1);

4 V& M9 T' W$ r# I% o# _5 R9 Dtitle('加入椒盐噪声后的图像');
  I+ K- m4 Z! m( {5 X
; o8 H( V. ]( ^$ Cf=double(P1);

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

g=fftshift(g);%转换数据矩阵fftshift是针对频域的，将FFT的DC分量移到频谱中心
" E: X- a" @* a# y+ n% j2 i3 V/ Y即对频域的图像，
+ N9 a3 l7 F6 \. p; A7 }
[M,N]=size(g);

D0=input('inputnonnegative dhreshold D0=');
- R0 v- U1 `1 D6 b$ d! w7 ?+ H8 G
n=input('input theorder of filtering n=');
& _2 K3 q% N- B  |
$ f9 X- @( j% }* on1=fix(M/2);

! f  j4 {+ T4 q  rn2= fix(N/2);
6 e2 Y, }& n- ~" O( B
# C% ^4 ]$ t2 N5 b2 r# P) A1 hfor i=1:M

7 k7 j! A7 V+ b* p& H    for j=1:N
& E' C% i/ y* D2 Y
        d=sqrt((i-n1)^2+(j-n2)^2);

  H+ s1 e" f3 Rh=1/(1+(d/D0)^(2*n));

        result(i,j)=h*g(i,j);
- z, C  W% R; q2 M! S0 @8 f7 B
0 i- Z8 h9 U2 |4 s7 F! x9 z    end

% T" a& T9 J4 ]end

result=ifftshift(result);
4 H1 F/ q* `% d0 g2 W
X2=ifft2(result);

X3=uint8(real(X2));
  e( y4 ]) J) i5 Q6 Y# @( x+ H( h
' Y9 i) y! E/ R% B1 v7 p! esubplot(133)

imshow (X3);

title('Butterworth低通滤波所得图像');

9 I# V* Z( E3 Z+ {: V  K( [

% J2 i/ A$ T4 r" {

/ g. M( y# `( c

NNNei256

基于matlab的Butterworth低通滤波程序