a b/Thoracic Organs Segmentation code/rescalingFunction.m 

  1
  





  
  

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





  
  

  3
  





  
  

  4
  
image=double(image);





  
  

  5
  
maxValue=double(i_upper);





  
  

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





  
  

  7
  
image(image>i_upper) = maxValue;





  
  

  8
  
image(image<i_lower) = 0;





  
  

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





  
  

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





  
  

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





  
  

  12
  
image=uint16(image);