close all;
clear all;

[I1, lut1] = imread('lena.tif');
I2=imnoise(I1, 'gaussian', 0, 0.005);
figure
image(I1)
colormap(requant(256));

figure
image(I1)
colormap(requant(8));

figure
image(I2)
colormap(requant(8));

figure
image(I1)
colormap(requant(16));