def gen_samp_tif(rand_temp): from osgeo import gdal from osgeo import osr import numpy as np import os,sys import cv2 import csv f_name='TCI.png' f_tiff_name = 'TCI.tif' n_file = cv2.imread(f_name) orig_h,orig_w,orig_rgb = n_file.shape image_size = (orig_w,orig_h) b_pixels,g_pixels,r_pixels = cv2.split(n_file) nx = image_size[0] ny = image_size[1] xmin,ymin,xmax,ymax = [newfieldminlong, newfieldminlat, newfieldmaxlong, newfieldmaxlat] xres = (float(xmax)-float(xmin))/float(nx) yres = (float(ymax)-float(ymin))/float(ny) xmin = float(xmin) ymin = float(ymin) xmax = float(xmax) ymax = float(ymax) xres = float(xres) yres = float(yres) print('res') #print(xres) #print(yres) geotransform = (xmin,xres,0,ymax,0,-yres) dst_ds = gdal.GetDriverByName('GTiff').Create(f_tiff_name,ny,nx,3,gdal.GDT_Byte) dst_ds.SetGeoTransform(geotransform) srs = osr.SpatialReference() #srs.ImportFromEPSG(3857) srs.ImportFromEPSG(4326) dst_ds.SetProjection(srs.ExportToWkt()) dst_ds.GetRasterBand(1).WriteArray(r_pixels) dst_ds.GetRasterBand(2).WriteArray(g_pixels) dst_ds.GetRasterBand(3).WriteArray(b_pixels) dst_ds.FlushCache() dst_ds=None