from osgeo import gdal
import numpy as np

data = np.zeros((10, 10), dtype=np.uint8)
driver = gdal.GetDriverByName('GTiff')
ds = driver.Create('tiny.tif', 10, 10, 1, gdal.GDT_Byte)
ds.GetRasterBand(1).WriteArray(data)
ds.SetGeoTransform([0, 10, 0, 100, 0, -10])
ds.SetProjection('GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]]')
ds = None