%%?讀取圖像并灰度轉換
RGB_fig=imread（‘C:\Users\ylxhm\Desktop\圖像分析\圖片處理的圖像\166.jpg‘）;%讀取圖像
gray_fig=rgb2gray（RGB_fig）;%得到灰度圖像
figure;imshow（RGB_fig）;title（‘原圖‘）;
figure;imshow（gray_fig）;title（‘灰度圖像‘）;

%%?改進的中值濾波去噪
%改進的中值濾波
outcome1_figure?=gray_fig;
outcome1_figure（:）=0;
alreadyProcessed=false（size（gray_fig））;%生成邏輯非的矩陣
%?迭代
noise_figure=imnoise（gray_fig‘salt?&?pepper‘0.3）;%添加椒鹽噪聲后的圖像
figure;imshow（noise_figure）;title（‘加噪圖像‘）;

Smax=7;
for?k=3:2:Smax
???zmin=ordfilt2（noise_figure1ones（k?k）‘symmetric‘）;
???zmax=ordfilt2（noise_figurek*kones（k?k）‘symmetric‘）;
???zmed=medfilt2（noise_figure[k?k]‘symmetric‘）;
???processUsingLevelB=（zmed>zmin）&（zmax>zmed）&?...
???????~alreadyProcessed;?
???zB=（noise_figure>zmin）&（zmax>noise_figure）;
???outputZxy=processUsingLevelB&zB;
???outputZmed=processUsingLevelB&~zB;
???outcome1_figure（outputZxy）=noise_figure（outputZxy）;
???outcome1_figure（outputZmed）=zmed（outputZmed）;