window_width_level.m

本函数实现对16位位深医学灰度图像的窗宽窗位调节

%
function [out]=Window_Width_Level(in,W,L)




% clear all
% M=256;
% N=512;
% NewM=256;
% NewN=256;
% W=4096;
% L=16384;
% 
% 
% 
%  p=1;% repeation of same point.
% %%offset=4096+8*N*p;
% offset=0;
% 
% fname='D:\My Codes\Matlab_old_2005.11.13\Sense_Own\full_01.raw';
% fid = fopen('data.txt','w');
%         
% %%coil01=READ_DATA(fname,offset,M,N)
% coil0=FFT_RECONSTRUCT(fname,M,N,NewM,NewN);    
% in=abs(coil0);
% figure(1);
% imshow(in);
% 

 [M,N]=size(in);
 if((W<1)&(L<1))
  MaxGray=1;
 else 
    
     MaxGray=65535;
     
 end    
inmode=in*MaxGray/Maximum(in);

out=zeros(M,N);

low=L-W/2;
if low<0 
    low=0;
    end
high=L+W/2;

if high>MaxGray
    high=MaxGray;
end  

for j=1:M
    for k=1:N
       
     %%  inmode(j,k)=inmode(j,k)*32767/pmax;
        if (inmode(j,k)<low)
            inmode(j,k)=low;
        
        elseif (inmode(j,k)>high)
                inmode(j,k)=high;
        end

        out(j,k)=(inmode(j,k)-low)*MaxGray/(high-low);
        
        
        
    end
       

end