图像处理代码共享

Colbie

图像处理代码共享

---------------------------------我是分割线--------------------------------------
h=fspecial('average',[3 3])%3*3的均值滤波器
% n8 X' [7 a  Z$ _h=fspecial('disk',5)%半径为5的圆形平均模板
* w' K8 t( u& |/ ~: N/ D8 X. A  Nh=fspecial('gaussian',[5 5],0.5)%5*5,标准差为0.5的高斯低通滤波器
h=fspecial('sobel')%加强水平边缘的竖直梯度算子

---------------------------------我是分割线--------------------------------------
1 K. f$ v" o) s( N! o* t- lI=imread('baby_noise.bmp');
3 j4 x$ }/ w- H! d$ F& z. zfigure,imshow(I);
4 ^5 m$ K: |9 _. [* Jh3_5=fspecial('gaussian',3,0.5);
I3_5=imfilter(I,h3_5);
figure,imshow(I3_5);
1 X1 v! ]* L$ k4 w5 Q5 B2 F9 Qh3_18=fspecial('gaussian',3,1.8);
$ e/ U- _5 c2 w; II3_18=imfilter(I,h3_18);
% W- X$ Y2 N( F  W% P8 u* d& ifigure,imshow(I3_18);
  c1 j& A0 K, T3 ih5_8=fspecial('gaussian',5,0.8);
I5_8=imfilter(I,h5_8);
4 g: O$ }  }7 A- w  sfigure,imshow(I5_8);
imwrite(I5_8,'baby5_8.bmp');
h7_12=fspecial('gaussian',7,1.2);
/ z: C' M7 t0 Y9 Z4 RI7_12=imfilter(I,h7_12);
figure,imshow(I7_12);
imwrite(I7_12,'baby7_12.bmp');
! y6 I( p1 A/ Q: F/ ^
---------------------------------我是分割线--------------------------------------
f=imread('cameraman.tif');
) ?$ Z  [8 v8 r  I8 fimshow(f);
. L6 `5 r3 H$ Rw=[0.1 0.1 0.1;0.1 0.1 0.1;0.1 0.1 0.1]
: h/ e8 P% H6 q4 tw
  b3 H% @" z" V2 \$ dg=imfilter(f,w,'corr','replicate');
( c$ x7 C. ?: Z! Rfigure,imshow(g);

---------------------------------我是分割线--------------------------------------

/ [; c) \: m1 i) w6 p  MI=imread('baby_noise.bmp');
figure,imshow(I);
9 ~+ H% L: z/ X7 C) v5 S1 T- k6 X$ Th=fspecial('average',3);
9 \; z( U$ g6 m) XI3=imfilter(I,h,'corr','replicate');
' p  Z6 D, m3 k. ^- v( Rfigure,imshow(I3);
: i: C9 {. @; Yh=fspecial('average',5);
I5=imfilter(I,h,'corr','replicate');
0 \) M- y+ o' [; J, t8 Ofigure,imshow(I5);
4 v2 q+ b2 Q& u7 }6 m( Z. P
: F# k9 C" K$ b% \% v---------------------------------我是分割线--------------------------------------
%Robert 交叉梯度
I=imread('bacteria.bmp');
imshow(I);
w1=[-1 0;0 1];
w2=[0 -1;1 0];
G1=imfilter(I,w1,'corr','replicate');
2 v% a% p4 b3 S$ E' KG2=imfilter(I,w2,'corr','replicate');
6 ]% [, Y1 J! K6 c- ?+ P  `# _G=abs(G1)+abs(G2);
# T/ Q' K  C$ U) t( O; n2 s# ^# I& Efigure,imshow(G,[]);
3 p+ _  E; B: Kfigure,imshow(abs(G1),[]);
6 h+ W9 U' ?! d$ x* A% Yfigure,imshow(abs(G2),[]);

, A: b  ^8 x8 ]& ^+ s5 m, W---------------------------------我是分割线--------------------------------------
%三种拉普拉斯模板
! [. o- U. `; w! t( m& B7 VI=imread('bacteria.bmp');
imshow(I);
  C" E# X& F: a8 WI=double(I);
w1=[0 -1 0;-1 4 -1;0 -1 0];
L1=imfilter(I,w1,'corr','replicate');
$ _6 j0 [! F. {( m% P3 Ow2=[-1 -1 -1;-1 8 -1;-1 -1 -1];
/ I# @5 S+ q3 p& t, lL2=imfilter(I,w2,'corr','replicate');
/ X2 I2 @& {9 q. m1 tfigure,imshow(abs(L1),[]);
6 `  B. Y$ |8 u' ~( m$ j# Cfigure,imshow(abs(L2),[]);
w3=[1 4 1;4 -20 4;1 4 1];
L3=imfilter(I,w3,'corr','replicate');
! j9 b2 d0 R3 R& l! Gfigure,imshow(abs(L3),[]);
: g8 n6 f$ A; A/ V
/ [2 \) ]5 `7 I6 K8 V0 z% l: R9 w) \---------------------------------我是分割线--------------------------------------
9 ]) _  \+ ]5 v- K
%sobel梯度
I=imread('bacteria.bmp');
w1=fspecial('sobel');
) C2 {- n  F6 m' zw2=w1';
G1=imfilter(I,w1);
; Q6 m3 Q0 N$ D- \& Y( KG2=imfilter(I,w2);
G=abs(G1)+abs(G2);
( ^5 M# J& k& tfigure,imshow(G1,[]);
figure,imshow(G2,[]);
3 x) \8 q; V8 O; Q2 \" qfigure,imshow(G,[]);
) }& n3 b0 q# ?/ W
---------------------------------我是分割线--------------------------------------
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');
figure,imshow(uint8(abs(I_log)),[]);
h_log=fspecial('log',5,0.5);
I_log=imfilter(Id,h_log,'corr','replicate');
* S* r6 O( T; `9 n: K% ifigure,imshow(uint8(abs(I_log)),[]);
h_log=fspecial('log',5,2);
I_log=imfilter(Id,h_log,'corr','replicate');
/ C4 ?$ E6 J9 F$ b! qfigure,imshow(uint8(abs(I_log)),[]);

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

aplerll

看看都有哪些代码。