图像处理代码共享

Colbie

图像处理代码共享
5 n! C$ _4 P7 X, _( f, w" e0 m& ?' E
) G. w; m; U, x9 @3 {---------------------------------我是分割线--------------------------------------
h=fspecial('average',[3 3])%3*3的均值滤波器
6 m4 [8 K. P  t: |- H* v2 |h=fspecial('disk',5)%半径为5的圆形平均模板
h=fspecial('gaussian',[5 5],0.5)%5*5,标准差为0.5的高斯低通滤波器
( ?' Q* F7 Z3 Q% m2 A* V% Fh=fspecial('sobel')%加强水平边缘的竖直梯度算子
# C% Y2 V2 @: c& Z/ S3 `- F! b0 V
---------------------------------我是分割线--------------------------------------
# t% s: ^4 A% U4 HI=imread('baby_noise.bmp');
figure,imshow(I);
2 }  i5 M% a- ], x% @2 L2 ah3_5=fspecial('gaussian',3,0.5);
I3_5=imfilter(I,h3_5);
6 Y9 ~$ m. G7 M$ N) d' Zfigure,imshow(I3_5);
& X, @, ~" z) h0 V$ sh3_18=fspecial('gaussian',3,1.8);
I3_18=imfilter(I,h3_18);
! {7 p) {8 a  qfigure,imshow(I3_18);
h5_8=fspecial('gaussian',5,0.8);
1 C. q8 n* F0 P# h  oI5_8=imfilter(I,h5_8);
figure,imshow(I5_8);
imwrite(I5_8,'baby5_8.bmp');
h7_12=fspecial('gaussian',7,1.2);
I7_12=imfilter(I,h7_12);
figure,imshow(I7_12);
imwrite(I7_12,'baby7_12.bmp');
* K  {+ [* o- R( e
---------------------------------我是分割线--------------------------------------
f=imread('cameraman.tif');
0 L. K: C. q( D1 a6 d- yimshow(f);
* c$ f5 Q/ L6 J2 N+ Q! o  mw=[0.1 0.1 0.1;0.1 0.1 0.1;0.1 0.1 0.1]
, H, H+ T1 s3 e, yw
) @" K7 V2 ~& H! W3 v; S' |0 ^g=imfilter(f,w,'corr','replicate');
6 ^5 D: c+ X3 S7 t8 O* N8 ~figure,imshow(g);
$ |) T) \9 p+ v
---------------------------------我是分割线--------------------------------------

( \: ]/ M* D; s' bI=imread('baby_noise.bmp');
figure,imshow(I);
9 P" k0 z3 C; A7 f5 I5 h, uh=fspecial('average',3);
) _) J9 H* Y( o8 ^9 H% [; k2 YI3=imfilter(I,h,'corr','replicate');
# p$ C/ q8 o" i: gfigure,imshow(I3);
h=fspecial('average',5);
' ?1 X; C" z& N* F/ Y7 F8 \I5=imfilter(I,h,'corr','replicate');
6 @" F9 u3 t3 n" \figure,imshow(I5);

---------------------------------我是分割线--------------------------------------
%Robert 交叉梯度
. l/ ?- ~$ I8 V5 }I=imread('bacteria.bmp');
imshow(I);
w1=[-1 0;0 1];
w2=[0 -1;1 0];
G1=imfilter(I,w1,'corr','replicate');
8 U( t" ?7 K/ L' n$ A: T$ d: CG2=imfilter(I,w2,'corr','replicate');
G=abs(G1)+abs(G2);
4 Y0 f* {6 M) ~& U% wfigure,imshow(G,[]);
1 Z9 ^/ ]( y9 Kfigure,imshow(abs(G1),[]);
figure,imshow(abs(G2),[]);

5 g7 w/ f: A; G% P* h3 f---------------------------------我是分割线--------------------------------------
# G2 p+ e' _* d! |# j%三种拉普拉斯模板
& ?9 m* K: C( o- L* B7 `; iI=imread('bacteria.bmp');
imshow(I);
: r% I" s. u* X& v! B- tI=double(I);
' H" e' R9 Z4 e) w# g! Ww1=[0 -1 0;-1 4 -1;0 -1 0];
+ H2 e+ ]$ w% |7 a7 V; [6 SL1=imfilter(I,w1,'corr','replicate');
w2=[-1 -1 -1;-1 8 -1;-1 -1 -1];
+ y$ n# W# {' Q3 A* {L2=imfilter(I,w2,'corr','replicate');
+ T8 _7 D' {9 }* n( ]: [7 Q1 rfigure,imshow(abs(L1),[]);
( K' n2 q! v' yfigure,imshow(abs(L2),[]);
- [) z4 X; }4 Z# sw3=[1 4 1;4 -20 4;1 4 1];
7 Z& {5 c) a1 ]L3=imfilter(I,w3,'corr','replicate');
figure,imshow(abs(L3),[]);

3 i; P3 K! z9 V, o! a3 T---------------------------------我是分割线--------------------------------------

/ o# C& C8 v4 N/ c0 K" n% ~; i# ^%sobel梯度
I=imread('bacteria.bmp');
w1=fspecial('sobel');
$ i: A) b/ y4 Q  r4 t# d+ Z6 Xw2=w1';
G1=imfilter(I,w1);
6 ^6 e8 s. d) \5 ~$ j1 YG2=imfilter(I,w2);
# o4 k: a! S$ I& WG=abs(G1)+abs(G2);
! l# V6 V6 {* e  x6 ~figure,imshow(G1,[]);
figure,imshow(G2,[]);
figure,imshow(G,[]);

$ u* {6 K8 I: K" I( ?2 S---------------------------------我是分割线--------------------------------------
9 i6 J4 C3 A, R8 }$ GI=imread('baby.bmp');
1 }' B: ]8 K+ }figure,imshow(I,[]);
: H+ Q: k: f, HId=double(I);
. c) i" R- _; w) w8 `' j. L1 eh_lap=[-1 -1 -1;-1 8 -1;-1 -1 -1];
1 ~/ a# S0 A% J. A+ o6 U  cI_lap=imfilter(Id,h_lap,'corr','replicate');
figure,imshow(uint8(abs(I_log)),[]);
h_log=fspecial('log',5,0.5);
' d% Z( O; g6 F( Z4 `' h0 zI_log=imfilter(Id,h_log,'corr','replicate');
figure,imshow(uint8(abs(I_log)),[]);
h_log=fspecial('log',5,2);
I_log=imfilter(Id,h_log,'corr','replicate');
figure,imshow(uint8(abs(I_log)),[]);

游客，如果您要查看本帖隐藏内容请回复

aplerll

看看都有哪些代码。