close?all;%最大值濾波
inputimage=imread（‘0.bmp‘）;
figure（1）imshow（inputimage）;title（‘原始圖‘）;
inputimage1=inputimage（::）;
Y1=ordfilt2（inputimage11ones（33））;
%?Z（::1）=Y1;
?Z=Y1;
figure（2）imshow（Z）;title（‘均值濾波‘）;

[M_CUTN_CUT]=size（Z）;
kth=2.1;
q=Z（::）;
q=q（~isnan（q））;
mu=mean（q）;
sigma=std（d