errorgraphs.m

jpg compression and decompression~

% first do S-CIELAB errors
dcTuneError(1)=scieError('parrot.tif', 'parrot-2k.jpg');
pspError(1)=scieError('parrot.tif', 'parrot-q99.jpg');

dcTuneError(2)=scieError('parrot.tif', 'parrot-6k.jpg');
pspError(2)=scieError('parrot.tif', 'parrot-q75.jpg');
edgeError(1)=scieError('parrot.tif', 'parrot-bj-6k.jpg');

dcTuneError(3)=scieError('parrot.tif', 'parrot-9k.jpg');
pspError(3)=scieError('parrot.tif', 'parrot-q50.jpg');
edgeError(2)=scieError('parrot.tif', 'parrot-bj-9k.jpg');

dcTuneError(4)=scieError('parrot.tif', 'parrot-13k.jpg');
pspError(4)=scieError('parrot.tif', 'parrot-q25.jpg');
edgeError(3)=scieError('parrot.tif', 'parrot-bj-13k.jpg');

% now do SNR errors
dcTuneSNR(1)=imageSNR('parrot.tif', 'parrot-2k.jpg');
pspSNR(1)=imageSNR('parrot.tif', 'parrot-q99.jpg');

dcTuneSNR(2)=imageSNR('parrot.tif', 'parrot-6k.jpg');
pspSNR(2)=imageSNR('parrot.tif', 'parrot-q75.jpg');
edgeSNR(1)=imageSNR('parrot.tif', 'parrot-bj-6k.jpg');

dcTuneSNR(3)=imageSNR('parrot.tif', 'parrot-9k.jpg');
pspSNR(3)=imageSNR('parrot.tif', 'parrot-q50.jpg');
edgeSNR(2)=imageSNR('parrot.tif', 'parrot-bj-9k.jpg');

dcTuneSNR(4)=imageSNR('parrot.tif', 'parrot-13k.jpg');
pspSNR(4)=imageSNR('parrot.tif', 'parrot-q25.jpg');
edgeSNR(3)=imageSNR('parrot.tif', 'parrot-bj-13k.jpg');

% bit rate vectors
bitRates=[2000*8/(384*256) 6000*8/(384*256) 9000*8/(384*256) ...
      13000*8/(384*256)];
bitRateEdge=[6000*8/(384*256) 9000*8/(384*256) 13000*8/(384*256)];

% Plot the stuff
figure;
subplot(2, 1, 1);
plot(bitRates, dcTuneError, 'o-',bitRates, pspError,'x:', ...
    bitRateEdge, edgeError, 's--');
xlabel('Bit Rate');
ylabel('S-CIELAB Error');
title('S-CIELAB Error vs. Bit Rate');

subplot(2, 1, 2);
plot(bitRates, dcTuneSNR, 'o-',bitRates, pspSNR, 'x:', ...
    bitRateEdge, edgeSNR, 's--');
xlabel('Bit Rate');
ylabel('PSNR (dB)');
title('PSNR vs. Bit Rate');

% Now calculate the error values for the other three pictures

% Buildings
dcBuildingError=scieError('buildings.tif', 'buildings-13k.jpg');
dcBuildingsSNR=imageSNR('buildings.tif', 'buildings-13k.jpg');

pspBuildingError=scieError('buildings.tif', 'buildings-q75.jpg');
pspBuildingsSNR=imageSNR('buildings.tif', 'buildings-q75.jpg');

edgeBuildingError=scieError('buildings.tif', 'buildings-bj-13k.jpg');
edgeBuildingsSNR=imageSNR('buildings.tif', 'buildings-bj-13k.jpg');

% Hats
dcHatsError=scieError('hats.tif', 'hats-6k.jpg');
dcHatssSNR=imageSNR('hats.tif', 'hats-6k.jpg');

pspHatsError=scieError('hats.tif', 'hats-q75.jpg');
pspHatsSNR=imageSNR('hats.tif', 'hats-q75.jpg');

edgeHatsError=scieError('hats.tif', 'hats-bj-6k.jpg');
edgeHatsSNR=imageSNR('hats.tif', 'hats-bj-6k.jpg');

% Island
dcIslandError=scieError('island.tif', 'island-6k.jpg');
dcIslandsSNR=imageSNR('island.tif', 'island-6k.jpg');

pspIslandError=scieError('island.tif', 'island-q75.jpg');
pspIslandSNR=imageSNR('island.tif', 'island-q75.jpg');

edgeIslandError=scieError('island.tif', 'island-bj-7k.jpg');
edgeIslandSNR=imageSNR('island.tif', 'island-bj-7k.jpg');