def calc_trmm_da_monthly(year, month, file_prefix="trmm", calc=False): try: file_input_month = output_folder + file_prefix + "_" + month + "_" + year + ".tif" file_input_avg = output_folder + "avg/" + file_prefix + "_" + month + "_avg.tif" file_output = output_folder + "da/" + file_prefix + "_" + month + "_" + year + "_da.tif" file_output_processed = output_folder + "output/" + file_prefix + "_" + month + "_" + year + "_da.tif" files = [file_input_month, file_input_avg] if calc: if os.path.isfile(file_input_month) and os.path.isfile( file_input_avg): try: calc_layers(files, file_output, "ratio") except Exception, e: log.error(e) pass process_layers(file_output, file_output_processed) # TODO: reprojection except Exception, e: log.error(e)
def calc_trmm_monthly(year, month, file_prefix="trmm", calc=False): try: files_path = input_folder + year + "/" + month + "/*.tif" output_filename = file_prefix + "_" + month + "_" + year + ".tif" output_file = output_folder + output_filename default_workspace = manager.geoserver.get_default_workspace() # calculate layers if calc: calc_layers(files_path, output_file, "sum") if os.path.isfile(output_file): # covert to geotiff 3857 the file if os.path.isdir(output_folder + "/output"): pass else: os.mkdir(output_folder + "/output") output_processed_file = output_folder + "/output/" + output_filename process_layers(output_file, output_processed_file) else: log.error("No file named %s " % output_file) except Exception, e: print e pass
def calc_trmm_da_monthly(year, month, file_prefix="trmm", calc=False): try: file_input_month = output_folder + file_prefix + "_" + month +"_" + year + ".tif" file_input_avg = output_folder + "avg/" + file_prefix + "_" + month +"_avg.tif" file_output = output_folder + "da/" + file_prefix + "_" + month + "_" + year + "_da.tif" file_output_processed = output_folder + "output/" + file_prefix + "_" + month + "_" + year + "_da.tif" files = [file_input_month, file_input_avg] if calc: if os.path.isfile(file_input_month) and os.path.isfile(file_input_avg): try: calc_layers(files, file_output, "ratio") except Exception, e: log.error(e) pass process_layers(file_output, file_output_processed) # TODO: reprojection except Exception, e: log.error(e)
def calc_trmm(): # take folders folders = glob.glob(input_folder) for folder in folders: print input_folder # create output folder /output # covert to geotiff 3857 the file if os.path.isdir(folder + "/output"): remove(folder + "/output/*.tif") else: os.mkdir(folder + "/output") folder_name = folder[folder.rindex("/")+1:] print folder_name file_output = folder + "/3b42_" + folder_name + "_4326.tif" print file_output # sum *.tif to output folder input_files = glob.glob(folder +"/*.tif") calc_layers(input_files, file_output, "sum") file_output_processed = folder + "/output/3b42_" + folder_name + ".tif" process_layers(file_output, file_output_processed)
def calc_trmm_monthly(year, month, file_prefix="trmm", calc=False): try: files_path = input_folder + year+ "/" + month + "/*.tif" output_filename = file_prefix+ "_"+ month +"_"+ year +".tif" output_file = output_folder + output_filename default_workspace = manager.geoserver.get_default_workspace() # calculate layers if calc: calc_layers(files_path, output_file, "sum") if os.path.isfile(output_file): # covert to geotiff 3857 the file if os.path.isdir(output_folder + "/output"): pass else: os.mkdir(output_folder + "/output") output_processed_file = output_folder + "/output/" + output_filename process_layers(output_file, output_processed_file) else: log.error("No file named %s " % output_file) except Exception, e: print e pass
def process(): if os.path.isdir(output_folder): log.info("already exists") else: os.mkdir(output_folder) output_file = output_folder + get_filename(input_file) + ".tif" process_layers(input_file, output_file) print output_file
def process_earthstat(): if os.path.isdir(output_folder): log.info("already exists") else: os.mkdir(output_folder) dir = glob.glob(input_folder + "*") for d in dir: if os.path.isdir(d): input_files = glob.glob(d + "/*.tif") for input_file in input_files: output_file = output_folder + get_filename(input_file) + ".tif" process_layers(input_file, output_file) print output_file
def calc_trmm_avg_monthly(month, file_prefix="trmm", calc=False): try: files_path = output_folder + file_prefix + "_" + month + "*" file_output = output_folder + "avg/" + file_prefix + "_" + month + "_avg.tif" file_output_processed = output_folder + "output/" + file_prefix + "_" + month + "_avg.tif" log.info(files_path) if calc: calc_layers(files_path, file_output, "avg") process_layers(file_output, file_output_processed) # TODO: reproject except Exception, e: log.error(e)
def calc_trmm_avg_monthly(month, file_prefix="trmm", calc=False): try: files_path = output_folder + file_prefix + "_" + month +"*" file_output = output_folder + "avg/" + file_prefix + "_" + month + "_avg.tif" file_output_processed = output_folder + "output/" + file_prefix + "_" + month + "_avg.tif" log.info(files_path) if calc: calc_layers(files_path, file_output, "avg") process_layers(file_output, file_output_processed) # TODO: reproject except Exception, e: log.error(e)
def calc(): # take folders folders = glob.glob(input_folder) for folder in folders: print input_folder # create output folder /output # covert to geotiff 3857 the file output_folder = folder + "/output" if os.path.isdir(output_folder): remove(output_folder + "/*.tif") else: os.mkdir(output_folder) input_files = glob.glob(folder +"/*.tif") for input_file in input_files: output_filename = output_folder + "/" + get_filename(input_file) + ".tif" print(output_filename) # create a geotiff + overviews process_layers(input_file, output_filename)
def calc_trmm(): # take folders folders = glob.glob(input_folder) for folder in folders: print input_folder # create output folder /output # covert to geotiff 3857 the file output_folder = folder + "/output" if os.path.isdir(output_folder): remove(output_folder + "/*.tif") else: os.mkdir(output_folder) input_files = glob.glob(folder + "/*.tif") for input_file in input_files: output_filename = output_folder + "/" + get_filename( input_file) + ".tif" print(output_filename) # create a geotiff + overviews process_layers(input_file, output_filename)