clear all, close all;

% choisir une des images
Ima=imread('lettres.tif');
Ima=imread('boat.tif');
Ima=imread('D93.tif');


figure, imshow(Ima);

[Ima2,dctIma2,QdctIma2]=compress(double(Ima),8);
figure, imshow(Ima2);
figure, imagesc(dctIma2), colorbar; title('tableau des coefs de DCT'), colormap(gray)
figure, imagesc(QdctIma2), colorbar; title('tableau des coefs de DCT quantifiés'), colormap(gray)
					     
% proportion de coefficient de dct à 0:
numel(find(dctIma2==0))/numel(Ima) % avant compression
numel(find(QdctIma2==0))/numel(Ima) % après compression

% proportion de coefficents mis à 0 par compression:
(numel(find(QdctIma2==0))-numel(find(dctIma2==0)))/numel(Ima) % après compression

figure, imagesc(QdctIma2+1000*(QdctIma2==0)), colorbar; title('tableau des coefs de DCT quantifiés (valeurs à 0 en évidence)') 
