图像处理代码共享

Colbie

图像处理代码共享
; N+ Y& ^) s! O
---------------------------------我是分割线--------------------------------------
3 f/ d5 u4 R1 A* p0 ~3 nh=fspecial('average',[3 3])%3*3的均值滤波器
h=fspecial('disk',5)%半径为5的圆形平均模板
h=fspecial('gaussian',[5 5],0.5)%5*5,标准差为0.5的高斯低通滤波器
5 ?# {- Y, p* u7 `: V$ `- Q( Ch=fspecial('sobel')%加强水平边缘的竖直梯度算子
3 q( g& D4 q( H; \, A6 ~
---------------------------------我是分割线--------------------------------------
I=imread('baby_noise.bmp');
5 x1 v) S6 `) K( J  bfigure,imshow(I);
$ X; x0 f9 ^& w( ]h3_5=fspecial('gaussian',3,0.5);
, v4 G' q/ K$ n+ u6 yI3_5=imfilter(I,h3_5);
figure,imshow(I3_5);
/ m9 A) V! c. hh3_18=fspecial('gaussian',3,1.8);
. p7 v$ K1 V" gI3_18=imfilter(I,h3_18);
0 l" O) V- S9 J5 B3 }, wfigure,imshow(I3_18);
h5_8=fspecial('gaussian',5,0.8);
I5_8=imfilter(I,h5_8);
+ [" Q+ w2 K1 Y& f. @! c& wfigure,imshow(I5_8);
imwrite(I5_8,'baby5_8.bmp');
$ A% }6 }. d1 jh7_12=fspecial('gaussian',7,1.2);
" [! G; ~3 u; @2 fI7_12=imfilter(I,h7_12);
figure,imshow(I7_12);
imwrite(I7_12,'baby7_12.bmp');
: C# K, G" T8 a% ^$ j/ d
---------------------------------我是分割线--------------------------------------
& X4 J/ m. t- v- z; `; t" h" Bf=imread('cameraman.tif');
" r3 `4 o; Y! s' i; @imshow(f);
w=[0.1 0.1 0.1;0.1 0.1 0.1;0.1 0.1 0.1]
+ `- J. a5 [0 r0 x( N. Zw
4 ^! Z8 s% M( ?8 e. m; W0 @) Ig=imfilter(f,w,'corr','replicate');
figure,imshow(g);

( ?* O; j& o, L% l; m7 x2 q7 E---------------------------------我是分割线--------------------------------------

I=imread('baby_noise.bmp');
figure,imshow(I);
- V, k  ?( N3 @h=fspecial('average',3);
! `1 H! d; q7 J2 m' M3 jI3=imfilter(I,h,'corr','replicate');
figure,imshow(I3);
h=fspecial('average',5);
I5=imfilter(I,h,'corr','replicate');
figure,imshow(I5);
3 V& M) F9 }* E% K7 P# w
3 {. Y% n$ I2 o. d7 ]* a0 e. `  G* {6 g- r---------------------------------我是分割线--------------------------------------
%Robert 交叉梯度
I=imread('bacteria.bmp');
: d& I) e+ {/ d3 timshow(I);
w1=[-1 0;0 1];
: O3 x0 R, U" V2 ^- zw2=[0 -1;1 0];
G1=imfilter(I,w1,'corr','replicate');
: e( n3 q  U% CG2=imfilter(I,w2,'corr','replicate');
+ b0 \# g7 D% M4 aG=abs(G1)+abs(G2);
figure,imshow(G,[]);
$ A' N8 ^" L- Q& Q3 i! _figure,imshow(abs(G1),[]);
" w, F0 ^$ B) [- M8 }figure,imshow(abs(G2),[]);
4 W2 t, X% g0 _8 [6 F" o, f
---------------------------------我是分割线--------------------------------------
0 S$ p% X) V5 Q5 L9 L* K! B%三种拉普拉斯模板
: b: D6 S! k  [1 V/ L: t$ U8 sI=imread('bacteria.bmp');
* R/ K; P! Q0 U8 ?% Ximshow(I);
4 I  o( W% f6 d1 i; ?* n) a7 kI=double(I);
/ y3 t" R$ N# d7 I" |2 Sw1=[0 -1 0;-1 4 -1;0 -1 0];
4 |+ L* T& m; P) H' XL1=imfilter(I,w1,'corr','replicate');
7 G. n8 I$ y9 H! d% G$ Dw2=[-1 -1 -1;-1 8 -1;-1 -1 -1];
L2=imfilter(I,w2,'corr','replicate');
4 g( [4 F1 X: `% f( m  @: gfigure,imshow(abs(L1),[]);
figure,imshow(abs(L2),[]);
; m( K! ~  y4 U, p. ?  k/ xw3=[1 4 1;4 -20 4;1 4 1];
L3=imfilter(I,w3,'corr','replicate');
7 C3 i9 c! t; f  nfigure,imshow(abs(L3),[]);

( G' x; j9 S* k4 r$ t$ E' D! u---------------------------------我是分割线--------------------------------------
6 Q) _' K/ W+ |. e) E! M- |
4 _# G. x1 A5 _3 _%sobel梯度
I=imread('bacteria.bmp');
w1=fspecial('sobel');
0 o5 W5 h. T) N+ D$ ~w2=w1';
G1=imfilter(I,w1);
. k+ L# H) o9 n2 wG2=imfilter(I,w2);
G=abs(G1)+abs(G2);
. T: V$ M- X$ Q1 [" l/ yfigure,imshow(G1,[]);
- k1 Q+ D  `' w) P9 O/ a% B! xfigure,imshow(G2,[]);
figure,imshow(G,[]);

---------------------------------我是分割线--------------------------------------
I=imread('baby.bmp');
figure,imshow(I,[]);
Id=double(I);
h_lap=[-1 -1 -1;-1 8 -1;-1 -1 -1];
I_lap=imfilter(Id,h_lap,'corr','replicate');
$ Y4 k4 u! z) S1 T. jfigure,imshow(uint8(abs(I_log)),[]);
! D2 E" d& q$ Y1 c5 q! _; Ch_log=fspecial('log',5,0.5);
I_log=imfilter(Id,h_log,'corr','replicate');
figure,imshow(uint8(abs(I_log)),[]);
h_log=fspecial('log',5,2);
6 M/ @% w6 O# ^- h% Y. _$ gI_log=imfilter(Id,h_log,'corr','replicate');
' L' ^, {! ^* t" A& A6 ~figure,imshow(uint8(abs(I_log)),[]);

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

aplerll

看看都有哪些代码。