|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
clear all;, }. q) q) j5 n# F: q& f
7 \) @; M0 d7 m2 G& w1 b
close all;
% B9 l! |/ }* E$ v6 H/ e/ t$ y# p7 H" v k
P=input('inputimage data:');
6 h6 A; r3 X' f- {+ P5 ]/ J" \1 k- W5 I- ]( K
d=size(P);; h9 }( O: |& t( C& j( J7 Q6 K6 C' K
' X& w K$ |/ s5 M6 H/ qif(d(3)>1)" B% O' V* i% y% ]
: ^( K( j) r$ i2 O# u
P=rgb2gray(P);
2 K `3 l3 x9 w t* F" }- E& U' a
/ O8 U. x# [% ]: N* @9 Cend8 [4 s5 U8 O: z* \9 N" u
# i# G% T3 G Z( y5 J: A; v/ }subplot(131);imshow(P);
, o" S' u0 S' }- g( C" D* ?% I' k& {: r, P+ v* n3 ]
title('原始图像');
* _3 A1 |9 G' {' M) Z, z. K: w, G9 U% @: Q- f. e
P1=imnoise(P,'salt& pepper',0.02);; V8 k# }& P1 p) k) |7 w
9 V& H7 s7 u4 P" Y( i
subplot(132)( E! V$ [- x A7 O9 B
; v) s& [% `8 g1 U( n, timshow (P1);" u; R$ j+ w8 r" }7 [( i
3 O5 h' g3 z7 G* ~/ o3 U
title('加入椒盐噪声后的图像');
( D) Z9 ^8 ?& u, S3 Q& v1 \3 j; x* a
0 N3 o$ @4 Y0 v m" qf=double(P1);
7 C% F1 Z( m2 w1 p9 b& @; j2 O1 x" H! g4 Q# C) r6 W3 n$ T
g=fft2(f);%傅里叶变换3 f& ~0 p4 L9 n L, n. ^
2 ^8 a, {2 Q4 s0 T) O8 l9 R0 R
g=fftshift(g);%转换数据矩阵fftshift是针对频域的,将FFT的DC分量移到频谱中心
) ^5 P6 a* O) S j! X& H即对频域的图像,
0 A8 W& U3 B, n: W% n# r" h( v' W3 R
[M,N]=size(g);: _. E) R4 m+ V8 G8 R! K
3 f; v5 Q3 Z( u: O1 U7 Z; V% wD0=input('inputnonnegative dhreshold D0=');: m4 _- c7 z7 F0 m+ ]9 b5 R
5 L3 L$ `8 x$ H7 n0 j( u! w# cn=input('input theorder of filtering n=');8 M! j$ z: h% `+ B9 z; K, c) ?! E
/ {' L Y5 g6 X$ u$ d5 g
n1=fix(M/2);
) r6 s8 [: s# A: C1 P3 M" _) y4 Q
n2= fix(N/2);
, C8 k+ m/ J9 `+ I* s0 t, Y* s2 B- J0 H* v0 a$ @( a) P1 G) H4 l
for i=1:M6 a/ Y; z* n* O. s
8 t: ?# c9 w4 G* A
for j=1:N
. ^) V$ c) M6 Q! P( T
# h& V. j- |1 w" A$ a: ~9 @) D d=sqrt((i-n1)^2+(j-n2)^2);( z4 K2 L+ h. B- V% w) Z! q
2 o* b; A# m. o# O
h=1/(1+(d/D0)^(2*n));
3 ], w+ C6 v6 l! h6 k' N" h& T1 G/ x- w) o# _4 ^# g
result(i,j)=h*g(i,j);
, z3 ^0 {/ Y6 I% Z9 D& b+ N$ _6 r. W# v+ \: ~' M7 A% F
end
& E, a" ]+ W1 a7 Z! Y5 i$ @; Q
' ]- L2 A7 Q; \# Y. g, Y4 ]end
: T6 X2 X% F& J6 x; q1 Q! H
K9 p& F8 o* C _2 xresult=ifftshift(result);7 C% Y$ C9 ]4 v9 y" }* p
' k# G/ Q5 G# @" M1 P+ YX2=ifft2(result);. k& a; Z$ N" G+ {
$ J, v" A2 E! T$ n" [X3=uint8(real(X2));& G% y3 l8 \3 m% C. q- P4 g
) k; T8 Z8 ^3 gsubplot(133)+ z O! O% v( W7 q
. ~* |- p7 u9 X3 S; m, D. Gimshow (X3);
/ q6 m- e# k0 i8 }7 A9 X5 g, O: P4 S# m3 W; G, t
title('Butterworth低通滤波所得图像');
; v$ q9 D! `3 _- L+ d w) M' w) U
7 X( X9 f$ d" G4 M/ y
9 }/ K1 @2 v( `, b+ {/ Q
, i/ z# K+ T! n1 V
: R5 P }$ N/ Z8 i4 U |
|