function v = variance(image)
% image intensity variance

m = mean(mean(image));
msqr = mean(mean(image.*image));

v = msqr - m*m;
