view general/numerical/normalise.m @ 61:eff6bddf82e3 tip

Finally implemented perceptual brightness thing.

function x=normalise(x,dim)
% normalise - make x zero mean and unit variance in given dimension
%
% normalise :: [Size], D:natural -> [Size].
%
% The return value has zero mean and unit variance in the
% Dth dimension.

x=divnorm(@(z)std(z,1,dim),addnorm(@(y)mean(y,dim),x));