path(path, '../');

im = jpeg_read('lenalo.jpg');

jpeg_write(im, 'lena1.jpg')

qt = im.quant_tables{1};
c = dequantize(im.coef_arrays{1}, qt);

om = im;

oqt = ones(8,8) * sum(qt(:)) / 64 /2;
oqt(oqt>255)=255;

om.coef_arrays{1} = quantize( c, oqt);
om.quant_tables{1} = oqt;

jpeg_write(om, 'lena2.jpg')

l1=imread('lena1.jpg');
l2=imread('lena2.jpg');

fx=200;
fy=180;
fw=190;
saveEPS(l1(fx:fx+fw,fy:fy+fw), 'lenaLo.eps');
saveEPS(l2(fx:fx+fw,fy:fy+fw), 'lenaLoUQ.eps');
saveEPS(abs(c(fx:fx+fw,fy:fy+fw)), 'lenaDCTface.eps');