图像处理代码共享

Colbie

图像处理代码共享
0 O, b: z2 {) v% J& e+ g
---------------------------------我是分割线--------------------------------------
1 \2 g6 L- m" _h=fspecial('average',[3 3])%3*3的均值滤波器
/ g9 X$ V, T; nh=fspecial('disk',5)%半径为5的圆形平均模板
h=fspecial('gaussian',[5 5],0.5)%5*5,标准差为0.5的高斯低通滤波器
h=fspecial('sobel')%加强水平边缘的竖直梯度算子

---------------------------------我是分割线--------------------------------------
2 p: O1 \; F* C* y8 hI=imread('baby_noise.bmp');
+ `+ g1 @! K0 l  w! Xfigure,imshow(I);
h3_5=fspecial('gaussian',3,0.5);
: c& m+ w4 E/ g3 \0 mI3_5=imfilter(I,h3_5);
) z% P* H2 a9 Y# [  K, w  \2 pfigure,imshow(I3_5);
% _( p; q2 O1 _h3_18=fspecial('gaussian',3,1.8);
; w1 @9 f. C  N4 \5 E- y* M0 n$ Z1 SI3_18=imfilter(I,h3_18);
3 R: D9 @+ _3 o6 e4 y. ^0 `! K7 H0 [figure,imshow(I3_18);
& y/ p8 G  N& Z) S; N8 Q7 j  Ah5_8=fspecial('gaussian',5,0.8);
6 q$ ~/ ~7 ^- i1 qI5_8=imfilter(I,h5_8);
figure,imshow(I5_8);
imwrite(I5_8,'baby5_8.bmp');
h7_12=fspecial('gaussian',7,1.2);
8 h! K7 M: l% Y7 K0 v8 e2 g; ^I7_12=imfilter(I,h7_12);
* x# @; N, I/ I; Y( z8 \  }. Gfigure,imshow(I7_12);
imwrite(I7_12,'baby7_12.bmp');

---------------------------------我是分割线--------------------------------------
8 T; X0 P) L( D: C) v( C/ f" Nf=imread('cameraman.tif');
imshow(f);
+ z- c; n# `! s  p; r, \4 Dw=[0.1 0.1 0.1;0.1 0.1 0.1;0.1 0.1 0.1]
w
g=imfilter(f,w,'corr','replicate');
/ T3 E7 h: J' m& G4 w/ g" Ofigure,imshow(g);

0 c- ~, P* x/ ~4 T---------------------------------我是分割线--------------------------------------
8 g+ f  a- G# [6 x
2 C* |; _) f. W- DI=imread('baby_noise.bmp');
3 s9 o* q" y, \5 ?; G; ofigure,imshow(I);
h=fspecial('average',3);
, {8 J" e- D0 _5 h; i, q! ]I3=imfilter(I,h,'corr','replicate');
figure,imshow(I3);
; W9 q0 H8 Z5 th=fspecial('average',5);
% B; U6 Z. `* |8 QI5=imfilter(I,h,'corr','replicate');
+ q6 ]3 O1 t5 U! Pfigure,imshow(I5);

6 L1 n. i) u1 f/ w---------------------------------我是分割线--------------------------------------
%Robert 交叉梯度
8 e' g( g  |. AI=imread('bacteria.bmp');
imshow(I);
w1=[-1 0;0 1];
9 X$ _0 h6 D3 y. B  }w2=[0 -1;1 0];
% n' [6 d$ R7 m% R) wG1=imfilter(I,w1,'corr','replicate');
G2=imfilter(I,w2,'corr','replicate');
! ~) R: C0 m; |1 e, RG=abs(G1)+abs(G2);
6 Z/ ?; D; Z7 ]# p; G& Gfigure,imshow(G,[]);
+ }; ]' J4 f  e6 ?3 e, D2 T% ifigure,imshow(abs(G1),[]);
figure,imshow(abs(G2),[]);

---------------------------------我是分割线--------------------------------------
, l: n: y' X* O5 `: A9 |%三种拉普拉斯模板
6 c+ f# _/ ~4 S7 E5 I" I; bI=imread('bacteria.bmp');
imshow(I);
  r# _+ N' l6 x1 rI=double(I);
w1=[0 -1 0;-1 4 -1;0 -1 0];
- k/ ?$ n; l$ r/ Q! ?8 `: q8 v( IL1=imfilter(I,w1,'corr','replicate');
/ {5 V. l! @6 Nw2=[-1 -1 -1;-1 8 -1;-1 -1 -1];
L2=imfilter(I,w2,'corr','replicate');
figure,imshow(abs(L1),[]);
- ?6 E$ E2 O! t. Ofigure,imshow(abs(L2),[]);
w3=[1 4 1;4 -20 4;1 4 1];
L3=imfilter(I,w3,'corr','replicate');
figure,imshow(abs(L3),[]);

3 c' s2 b; N7 M( d- |1 F8 o  o$ v---------------------------------我是分割线--------------------------------------

%sobel梯度
6 \1 m7 W( V3 l5 a4 H0 L# iI=imread('bacteria.bmp');
) j' \& |7 s' s& Q5 ?6 Dw1=fspecial('sobel');
( ?( j! H) t  e! ]4 o' Tw2=w1';
# Y/ H" z  B* `) z* u. L( ]1 EG1=imfilter(I,w1);
G2=imfilter(I,w2);
- W+ T& H- ^/ c. V( a2 eG=abs(G1)+abs(G2);
figure,imshow(G1,[]);
figure,imshow(G2,[]);
2 Z4 \9 f/ N& t8 @, S4 kfigure,imshow(G,[]);
( B( Z" y$ C9 p- Z  z
---------------------------------我是分割线--------------------------------------
I=imread('baby.bmp');
) }: P6 s' h! Y! z9 `figure,imshow(I,[]);
Id=double(I);
h_lap=[-1 -1 -1;-1 8 -1;-1 -1 -1];
$ g6 S. I$ I& h$ Q& P6 XI_lap=imfilter(Id,h_lap,'corr','replicate');
figure,imshow(uint8(abs(I_log)),[]);
& }! H+ H( x3 Y) U1 n# @3 D* Xh_log=fspecial('log',5,0.5);
I_log=imfilter(Id,h_log,'corr','replicate');
' E. G+ o9 ~- }$ J7 P  ~+ [figure,imshow(uint8(abs(I_log)),[]);
7 z( e5 Z. j  A4 R8 a; K$ Qh_log=fspecial('log',5,2);
I_log=imfilter(Id,h_log,'corr','replicate');
figure,imshow(uint8(abs(I_log)),[]);

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

aplerll

看看都有哪些代码。