#!/usr/bin/env python

import gdal
import gdalconst
import glob
import posixpath

scale = 50

#for current_filename in glob.glob('chunked/scale_%03d/*.tif' % (scale)):
for current_filename in glob.glob('chunked/scale_%03d/chunk_0_4.tif' % (scale)):
	basename = posixpath.basename(current_filename)
	old_filename = 'chunked-old/scale_%03d/%s' % (scale, basename)
	diff_filename = '../tmp/diff-%s' % (basename)
	print old_filename, current_filename
	old_image = gdal.Open(old_filename, gdalconst.GA_ReadOnly)
	current_image = gdal.Open(current_filename, gdalconst.GA_ReadOnly)
	diff_image = gdal.GetDriverByName('GTiff').Create(
	        diff_filename,
		old_image.RasterXSize, old_image.RasterYSize,
		3, gdal.GDT_Byte, 
		[ 'TILED=YES' ]
	)


	old_array = old_image.:wq