%圖片的預處理
I=imread（‘D:\matlab\radar.png‘）;
I=rgb2gray（I）;
figure（1）imshow（I）title（‘原圖‘）;

%對于高斯噪聲均值濾波效果好，反之，對椒鹽噪聲中值濾波效果好;
%對本例更多是椒鹽噪聲

J1=medfilt2（I）;?%缺省的【3?3】窗口對矩陣I進行二維中值濾波
figure（2）imshow（J1）title（‘中值濾波‘）;
%中值濾波后，回波圖像較為集中了，但丟失了部分面積較小的回波區域。

h1=fspecial（‘average‘3）;%fspecial（‘‘average‘3）用來創建一個3*3的均值濾波器
J2=filter2（h1I）/255;
figure（3）imshow（J2）title（‘均值濾波‘）;
%均值濾波邊緣模糊，但去除了孤立點，回波較為集中

J3=wiener2（I）;
figure（4）imshow（J3）title（‘維納濾波‘）;
%維納濾波與均值濾波類似，邊緣模糊厲害，但去除了孤立點，回波較為集中

h2=fspecial（‘gaussian‘）;
J4=filter2（h2I）/255;
figure（5）imshow（J4）title（‘高斯濾波‘）;
%高斯濾波還有孤立點，不好

%數學形態學濾波
%B=[0?1?0;1?1?1;0?1?0];
B1=ones（22）;%考慮到氣象是所有附近的點都需考慮，采用8領域算子。
J5=imopen（IB1）;
J6=imclose（J5B1）;
figure（6）imshow（J6）titl

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

?????文件???????1918??2013-04-10?09:25??quzao.m

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

?????????????????1918????????????????????1