try:
	from osgeo import gdal
except:
	import gdal
from osgeo import osr
import numpy as np
import os,sys
import cv2
import csv
import firebase_admin
from firebase_admin import credentials
from firebase_admin import db
from firebase_admin import messaging
import time
import threading
import os
from numba import jit, cuda
import traceback
cred = credentials.Certificate('servicekey.json')
try:
    firebase_admin.initialize_app(cred, {'databaseURL': 'https://farmbase-b2f7e-31c0c.firebaseio.com/','storageBucket': 'farmbase-b2f7e.appspot.com'})
except:
    print('fire running')
from firebase_admin import storage

bucket = storage.bucket()

uid = 'vPM6EKWHfHRxmRfc0YJ7PZ8J9tR2'
uid = 'zkBZGiByagbnwFvjOmtLnSBVLeu1'
uid = 'KzkCnSmHLpOkhnseAhhRWrI35BD2'
uid = 'JgbkYJkMVsRXpSndLsDlabDomOp2'
uid = 'VWoPJ2nfWDRlrYQo9dmctdLHeJN2'
#uid = 'wEJQVaxiZZaaId2aJe6lRJ8Zd0S2'
uid = 'xk6kws4YMkdeBAstSWDNsd42w9m2'

def create_geotiff(uid, fieldid, sensedday,imagetype):
	isdir = os.path.isdir(uid)
	if isdir != True:
		os.mkdir(uid)

	fieldobj = db.reference('PaidMonitoredFields').child('PMF').child(uid).child(fieldid).get()

	try:
		source_blob_name = 'PaidMonitoredFields/' + uid + '/' + fieldid + '/' + str(sensedday) + '/' + imagetype
		blob = bucket.blob(source_blob_name)
		print(source_blob_name)
		file_name = uid + '/'+ fieldid + '_' + imagetype + '.png'
		blob.download_to_filename(file_name)
		f_name= file_name
		f_tiff_name = uid + '/'+ fieldid + '_' + imagetype + '.tif'
		new_string_json = fieldobj
		print(f_name)
		print(f_tiff_name)
		n_file = cv2.imread(f_name)		
		orig_h,orig_w,orig_rgb = n_file.shape
		w_new = orig_w
		h_new = orig_h
		# factor = 14000
		# w = round(factor*w_new/h_new)
		# w_new = w
		# h_new = factor
		print('new_wh')
		print(w_new)
		print(h_new)
		n_file = cv2.resize(n_file,(w_new,h_new))
		image_size = (w_new,h_new)
		newfieldminlat =new_string_json["FieldMinLat"]
		newfieldmaxlat =new_string_json["FieldMaxLat"]
		newfieldminlong =new_string_json["FieldMinLong"]
		newfieldmaxlong =new_string_json["FieldMaxLong"]
		lat = [newfieldminlat, newfieldmaxlat]
		lon = [newfieldminlong, newfieldmaxlong]
		b_pixels,g_pixels,r_pixels = cv2.split(n_file)
		nx = image_size[1]
		ny = image_size[0]
		xmin,ymin,xmax,ymax = [newfieldminlong, newfieldminlat, newfieldmaxlong, newfieldmaxlat]
		xres = (float(xmax)-float(xmin))/float(ny)
		yres = (float(ymax)-float(ymin))/float(nx)
		xmin = float(xmin)
		ymin = float(ymin)
		xmax = float(xmax)
		ymax = float(ymax)
		xres = float(xres)
		yres = float(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(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
	except:
		aw = 1
		print(traceback.format_exc())