function imo = CS4640_scale(im,a,b)
%

[M,N] = size(im);

imo = zeros(M,N);
im_min = min(im(:));
im_max = max(im(:));
im_range = im_max - im_min;

for r = 1:M
    for c = 1:N
        imo(r,c) = a + (b-a)*(im(r,c) - im_min)/im_range;
    end
end
