%
%
% Mutual Information (MI) Thresholding
%
% [T1, T2, mi, imT1, imT2, imF, quality, miscore, mii1, mii2] = mutualinfoThreshold(im1, im2, mask, sigma, porder)
%
% Input:
%   im1/im2 : input data pair to be thresholded. Must be the same size
%
%   mask : binary mask the same size as im1/im2. Zero values will be removed from consideration in the data.
%      default = 1
%
%   sigma : smoothing factor for MI surface
%      default = no smoothing
%
%   porder : peak to choose (1=highest peak, 2=second highest, 3=third, etc..  -1=select peak with highest quality)
%      default = 1
%
%
%
% Output:
%
% T1,T2 : threshold for im1/im2
% mi : 256x256 MI surface
% imT1 : im1 thresholded with T1
% imT2 : im2 thresholded with T2
% imF : fusion of imT1 and imT2 using a simple neighbourhood support weighting
% quality : quality score at selected peak
% miscore : MI value at selected peak
% mii1,mii2 : (y,x) position of selected peak on MI surface
%
%
% author: Ciar醤