echo?on
?clc

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%???讀取圖像
?f=imread（‘F:\myimage\noise?experiment\test.BMP‘）;
?f=f（::1）;
?figure（1）
?imshow（f）
?pause?
?close?all
?clc
?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%????給圖像添加噪聲
?g=imnoise（f‘salt?&?pepper‘0.09）;??????????
?figure（2）
?imshow（g）
?title（‘椒鹽噪聲圖像‘）
?pause
?close?all
?clc
?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%???中值濾波
?J=medfilt2（g）;
?figure（3）
?imshow（J）
?title（‘中值濾波‘）????????????????????????????%%%%%%%%%%%%%%?????????中值濾波
?
?
?pause
?close?all
?clc
?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%????形態(tài)學(xué)濾波
?
w3=strel（‘square‘3）;
?L=imopen（Jw3）;
?figure（6）
?imshow（L）
?title（‘開運(yùn)算后的圖像‘）???????????????????????%%%%%%%%%%%%%%%%%%%??開運(yùn)算后的圖像
?
?
?m=imclose（Lw3）;
?figure（7）
?imshow（m）
title（‘開運(yùn)算后再閉運(yùn)算后的圖像‘）??????????????%%%%%%%%%%%%%%??開運(yùn)算后再閉運(yùn)算后的圖像
pause
close?all
clc
?
?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%??????????形態(tài)學(xué)灰度不均勻校正
?w5=strel（‘square‘105）;
?n11=imclose（mw5）;
?figure（28）
?imshow（n11）
?title（‘經(jīng)閉運(yùn)算后的圖像150‘）
?
?NNN=imsubtract（n11m）;
?figure（29）
?
?imshow（NNN）
?title（‘頂帽變換150‘）???
?
?pause?
?close?all
?clc
?

level?=?graythresh（m）;
BW?=?im2bw（mlevel）;
figure?imshow（BW）
pause
close?all
clc

?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1448??2008-09-30?10:28??noise5.m

-----------?---------??----------?-----??----

?????????????????1666????????????????????2