clear;
clc;
I=im2double（imread（‘F:\matlab\an\bin\tupian\wuhan3.tif‘））;
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+