function y=entropy(x)
%-----------------------------------------------
[count,xi]=imhist(x);
p=count./65536+eps;
y=-sum(sum(p.*log2(p)));%entropy is 7.5696.
%------------------------------------------------