دو تا اولی گلوبال و دوتای آخری هم لوکال
I= double(I);
variance = var (I(:));
miangin = mean (I(:));
[r c] = size(I);
img = zeros(r+2,c+2);
img(2:(r + 1) ,2: (c + 1)) = I;
img = [ img(1:r, 1:c)(:), img(1:r, 2: (c + 1))(:),img(1:r, 3: (c + 2 ))(:), img(2:(r + 1),1: c)(:), img(2:(r + 1),2: (c + 1))(:), img(2:(r + 1),3: (c + 2))(:), img(3:(r + 2),1: c)(:), img(3:(r + 2),2: (c + 1))(:), img(3:(r + 2),3: (c + 2))(:)];
temp = zeros(r+2,c+2);
temp(2:(r + 1) ,2: (c + 1)) = 1;
temp = [ temp(1:r, 1:c)(:), temp(1:r, 2: (c + 1))(:),temp(1:r, 3: (c + 2 ))(:), temp(2:(r + 1),1: c)(:), temp(2:(r + 1),2: (c + 1))(:), temp(2:(r + 1),3: (c + 2))(:), temp(3:(r + 2),1: c)(:), temp(3:(r + 2),2: (c + 1))(:), temp(3:(r + 2),3: (c + 2))(:)];
sumpix = sum(temp,2);
miangin_local = sum(img, 2) ./sumpix;
variance_local = (((I(:) -sum(img, 2) ./sumpix).^2) ./ sumpix);
miangin_local = reshape(miangin_local,r,c);
variance_local = reshape(variance_local,r,c);