% Usage: RIM = im_remap(IM,IN_FROM,IN_TO,OUT_FROM,OUT_TO,GAMMA)
% IM is clipped and remapped so that IN_FROM:IN_TO remaps to OUT_FROM:OUT_TO.
% GAMMA is optional (default 1).

function rim = im_remap(im,in_from,in_to,out_from,out_to,gamma)

if nargin < 6
    gamma = 1;
end

rim = im;
rim = max(rim,in_from);
rim = min(rim,in_to);
rim = (rim-in_from)/(in_to-in_from);
rim = rim.^gamma;
rim = out_from + rim*(out_to-out_from);