clear;
clc;
load?T72_mag.mat
I=mobandata{150};
subplot（231）
imshow（I）;
title（‘原圖像‘）
%求圖像背景的均值和方差
[mn]=size（I）;
N=m*n;
s=0;
v=0;
for?i=1:40
????for?j=1:128
????????s=s+I（ij）;
????????v=v+I（ij）*I（ij）;
????end????
end

for?i=80:128
????for?j=1:128
????????s=s+I（ij）;
????????v=v+I（ij）*I（ij）;
????end
end

for?i=1:128
????for?j=1:44
????????s=s+I（ij）;
????????v=v+I（ij）*I（ij）;
????end
end

for?i=1:128
????for?j=86:128
????????s=s+I（ij）;
????????v=v+I（ij）*I（ij）;
????end
end
mean=s/（2*40*128+2*48*128）;
var=v/（2*40*128+48*128）-（s/（2*40*128+2*48*128））^2;
sigma=sqrt（var）;
T=raylinv（0.999sigma）;
for?i=1:m
????for?j=1:n
????????if?I（ij）>T
????????????I（ij）=1;
????????else?I（ij）=0;
????????end
????end
end
subplot（232）
imshow（I）;
title（‘閾值判別后的圖象‘）
for?i=3:m-2
????for?j=3:n-2
????????count=0;
????????for?x=i-2:i+2
????????????for?y=j-2:j+2
??????????????count=count+I（xy）;
????????????end
????????end
????????if?count<9