CT 影像调整窗宽窗位代码

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