function image = rescalingFunction(image,i_lower,i_upper);

image=double(image);

maxValue=double(i_upper);

i_upper=0.9*double(i_upper);i_lower=double(i_lower);

image(image>i_upper) = maxValue;

image(image<i_lower) = 0;

indeces = find(image>i_lower & image<i_upper);

image(indeces) = (double(image(indeces))-i_lower)/(i_upper-i_lower)*maxValue;

%image = medfilt2 ( image, [5,5] );

image=uint16(image);