#!/usr/bin/env python

import gdal
import os

a = gdal.Open('/tmp/a.tif')

b = gdal.GetDriverByName('GTiff').Create(
	'/tmp/b.tif',
	a.RasterXSize, a.RasterYSize,
	3, gdal.GDT_Byte,
	[ 'TILED=YES' ]
)

#b.SetDescription(a.GetDescription())
b.SetProjection(a.GetProjection())
b.SetGeoTransform(a.GetGeoTransform())
#b.SetMetadata(a.GetMetadata())

#print a.GetProjection()
#print b.GetProjection()

b = None

os.system('listgeo /tmp/a.tif >/tmp/a.geo')
os.system('listgeo /tmp/b.tif >/tmp/b.geo')
os.system('diff /tmp/a.geo /tmp/b.geo')