% Test the routine MyBitDepth.m.
%
% Samuli Siltanen November 2014

% Load image 
im = imread('Tokyo2048.tif','tif');
im = double(im);

% Normalize image: pixel values between zero and one
im = im-min(im(:));
im = im/max(im(:));

% Downsample image
im = MyDS2(im); % 1024x1024
im = MyDS2(im); % 512x512
im = MyDS2(im); % 256x256
%im = MyDS2(im); % 128x128
%im = MyDS2(im); % 64x64
%im = MyDS2(im); % 32x32

% Discretize (quantize) the pixel values
B   = 8;
imB = MyBitDepth(im,B);

% View the result
figure(1)
clf
imagesc([im,imB])
axis equal
colormap gray
axis off