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 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)