def decompress_image(im1_blocks_q):
im1_blocks_dq = dequantize(im1_blocks_q)
im1_blocks_idct = inverse_transform(im1_blocks_dq)
im1_recovered = merge_blocks(im1_blocks_idct)
return im1_recovered
def decompress_image_with_dual_reference(im1_offsets, im1_residuals_q, reference_image1, reference_image2):
im1_residuals_dq = dequantize(im1_residuals_q)
im1_residuals_idct = inverse_transform(im1_residuals_dq)
im1_blocks_recovered = apply_dual_reference(im1_offsets, im1_residuals_idct, reference_image1, reference_image2)
im1_recovered = merge_blocks(im1_blocks_recovered)
return im1_recovered
from our_mpeg import transform, inverse_transform
from our_mpeg import quantize, dequantize
from our_mpeg import zigzag, inverse_zigzag
print('loading...', flush=True)
image_data = load(sys.argv[1])
print('finished loading')
sys.stdout.flush()
im1 = image_data[0,0,0]
im2 = image_data[1,0,0]
print('encoding...')
im1_blocks = break_blocks(im1)
im1_offsets, im1_residuals = compute_reference(im1_blocks, im2)
im1_residuals_dct = transform(im1_residuals)
im1_residuals_q = quantize(im1_residuals_dct)
im1_residuals_q_zig = zigzag(im1_residuals_q)
print('finished encoding')
print('decoding...')
im1_residuals_q_dzig = inverse_zigzag(im1_residuals_q_zig)
im1_residuals_dq = dequantize(im1_residuals_q_dzig)
im1_residuals_idct = inverse_transform(im1_residuals_dq)
im1_blocks_recovered = apply_reference(im1_offsets, im1_residuals_idct, im2)
im1_recovered = merge_blocks(im1_blocks_recovered)
print('finished decoding')
print(np.linalg.norm(im1_recovered - im1))
