www.gusucode.com > emg_functions > emg_functions/image_processing/im_jointEntropy.m

    % Betthauser 2016 -- Calculate joint entropy of 2 rvs in a 2d joint
% histogram
function jointEntropy = im_jointEntropy(jointHistogram)
    jointEntropy = 0; hsizey = size(jointHistogram,1); hsizex = size(jointHistogram,2);
    for i = 1:hsizey
        for j = 1:hsizex
            p = jointHistogram(i,j);
            if p > 0
                jointEntropy = jointEntropy - p*log2(p);        
            end
        end
    end
end