winw = 1500;
winl = -700;
imin = winl - winw/2;
imax = winl + winw/2;
disp([' window level: ', num2str(winl)]);
disp([' window width: ', num2str(winw)]);
disp('processing the window level and width ...');
imgdata(imgdata < imin) = imin;
imgdata(imgdata > imax) = imax;
for z = 1 : zmax
for y = 1 : ymax
for x = 1 : xmax
value = imgdata(x, y, z);
if value <= imin
value = 0;
elseif value < imax
value = (value - imin) / winw * 255;
elseif value >= imax
value = 255;
end
imgdata(x, y, z) = value;
end
end
disp([' processing the slice: ', num2str(z), '/', num2str(zmax)]);
end