Skip to content
Snippets Groups Projects

Resolve "Run snow annual map refactoring"

Merged Resolve "Run snow annual map refactoring"
54-run-snow-annual-map-refactoring into develop
Merged Aurore Dupuis requested to merge 54-run-snow-annual-map-refactoring into develop
Viewing commit eab10194
Prev Next
Show latest version
1 file
+ 2
− 2
Compare changes
  • Side-by-side
  • Inline
  • eab10194
    Aurore Dupuis
    #54 - Fix arguments order in convert_mask_line · eab10194
    Aurore Dupuis authored
python/s2snow/snow_annual_map.py
+ 262
− 377
@@ -19,86 +19,284 @@
# See the License for the specific language governing permissions and
# limitations under the License.
#
import logging
import os
import os.path as op
from os.path import basename
import shutil
import logging
import multiprocessing
from xml.dom import minidom
from datetime import timedelta
from lxml import etree
import gdal
from gdalconst import GA_ReadOnly
# OTB Applications
import otbApplication as otb
# Import python decorators for the different needed OTB applications
from s2snow.app_wrappers import band_math, get_app_output, super_impose, band_mathX
from s2snow.utils import str_to_datetime, datetime_to_str
from s2snow.utils import write_list_to_file, read_list_from_file
from s2snow.app_wrappers import band_math, get_app_output, super_impose, band_mathX, gap_filling
from s2snow.snow_product_parser import load_snow_product
from s2snow.utils import datetime_to_str
from s2snow.utils import write_list_to_file, read_list_from_file
# Build gdal option to generate maks of 1 byte using otb extended filename
# syntaxx
GDAL_OPT = "?&gdal:co:NBITS=1&gdal:co:COMPRESS=DEFLATE"
def parse_xml(filepath):
""" Parse an xml file to return the zs value of a snow product
"""
logging.debug("Parsing " + filepath)
xmldoc = minidom.parse(filepath)
group = xmldoc.getElementsByTagName('Global_Index_List')[0]
zs = group.getElementsByTagName("QUALITY_INDEX")[0].firstChild.data
LABEL_NO_SNOW = "0"
LABEL_SNOW = "100"
LABEL_CLOUD = "205"
LABEL_NO_DATA = "255"
LABEL_NO_DATA_OLD = "254"
def compute_SCD(binary_snow_mask_list, multitemp_cloud_vrt, input_dates_filename, output_dates_filename, output_dates,
snow_coverage_duration, multitemp_snow_vrt, gapfilled_timeserie, multitemp_snow100,
multitemp_snow100_gapfilled, ram):
logging.info("Start compute_SCD")
# build snow mask vrt
logging.info("Building multitemp snow mask vrt")
logging.info("snow vrt: {}", multitemp_snow_vrt)
gdal.BuildVRT(multitemp_snow_vrt,
binary_snow_mask_list,
separate=True,
srcNodata='None')
# multiply by 100 for the temporal interpolation
logging.info("Scale by 100 multitemp snow mask vrt")
bandMathXApp = band_mathX([multitemp_snow_vrt],
multitemp_snow100,
"im1 mlt 100",
ram,
otb.ImagePixelType_uint8)
bandMathXApp.ExecuteAndWriteOutput()
bandMathXApp = None
# gap filling the snow timeserie
logging.info("Gap filling the snow timeserie")
app_gap_filling = gap_filling(multitemp_snow100,
multitemp_cloud_vrt,
multitemp_snow100_gapfilled + "?&gdal:co:COMPRESS=DEFLATE",
input_dates_filename,
output_dates_filename,
ram,
otb.ImagePixelType_uint8)
# @TODO the mode is for now forced to DEBUG in order to generate img on disk
# img_in = get_app_output(app_gap_filling, "out", mode)
# if mode == "DEBUG":
# shutil.copy2(gapfilled_timeserie, path_out)
# app_gap_filling = None
img_in = get_app_output(app_gap_filling, "out", "DEBUG")
app_gap_filling = None
# threshold to 0 or 1
logging.info("Round to binary series of snow occurrence")
bandMathXApp = band_mathX([img_in],
gapfilled_timeserie + GDAL_OPT,
"(im1 mlt 2) dv 100",
ram,
otb.ImagePixelType_uint8)
bandMathXApp.ExecuteAndWriteOutput()
bandMathXApp = None
# generate the annual map
logging.info("Generate Snow coverate duration: {}", snow_coverage_duration)
band_index = list(range(1, len(output_dates) + 1))
logging.debug("Bande index: {}", band_index)
expression = "+".join(["im1b" + str(i) for i in band_index])
logging.debug("expression: {}", expression)
bandMathApp = band_math([gapfilled_timeserie],
snow_coverage_duration,
expression,
ram,
otb.ImagePixelType_uint16)
bandMathApp.ExecuteAndWriteOutput()
bandMathApp = None
logging.info("End compute_SCD")
def compute_CCD(binary_cloud_mask_list, ccd_file_path, multitemp_cloud_vrt, ram):
logging.info("Start compute_CCD")
# build cloud mask vrt
logging.info("Building multitemp cloud mask vrt")
logging.info("cloud vrt: {}".format(multitemp_cloud_vrt))
gdal.BuildVRT(multitemp_cloud_vrt,
binary_cloud_mask_list,
separate=True,
srcNodata='None')
# generate the summary map
logging.info("Generate Cloud_Occurence: {}", ccd_file_path)
band_index = list(range(1, len(binary_cloud_mask_list) + 1))
logging.debug("bande index: {}", band_index)
expression = "+".join(["im1b" + str(i) for i in band_index])
logging.debug("expression: {}", expression)
bandMathApp = band_math([multitemp_cloud_vrt],
ccd_file_path,
expression,
ram,
otb.ImagePixelType_uint16)
bandMathApp.ExecuteAndWriteOutput()
bandMathApp = None
logging.info("End compute_CCD")
def convert_snow_masks_into_binary_cloud_masks(path_out, ram, product_dict):
logging.info("Start convert_snow_masks_into_binary_cloud_masks")
# convert the snow masks into binary cloud masks
expression = "im1b1==" + LABEL_CLOUD + "?1:(im1b1==" + LABEL_NO_DATA + "?1:(im1b1==" + LABEL_NO_DATA_OLD + "?1:0))"
binary_cloud_mask_list = convert_mask_list(path_out, product_dict, expression, "cloud", ram,
mask_format=GDAL_OPT)
logging.debug("Binary cloud mask list:")
logging.debug(binary_cloud_mask_list)
logging.info("Start convert_snow_masks_into_binary_cloud_masks")
return binary_cloud_mask_list
def convert_snow_masks_into_binary_snow_masks(path_out, ram, product_dict):
logging.info("Start convert_snow_masks_into_binary_snow_masks")
# convert the snow masks into binary snow masks
expression = "(im1b1==" + LABEL_SNOW + ")?1:0"
binary_snow_mask_list = convert_mask_list(path_out, product_dict, expression, "snow", ram,
mask_format=GDAL_OPT)
logging.debug("Binary snow mask list:")
logging.debug(binary_snow_mask_list)
logging.info("End convert_snow_masks_into_binary_snow_masks")
return binary_snow_mask_list
def merge_product_at_same_date(path_out, product_dict, ram):
logging.info("Start merge_product_at_same_date")
merge_product_dict = {}
for key in list(product_dict.keys()):
if len(product_dict[key]) > 1:
merged_mask = op.join(path_out, key + "_merged_snow_product.tif")
merge_masks_at_same_date(product_dict[key],
merged_mask,
LABEL_SNOW,
ram)
merge_product_dict[key] = merged_mask
else:
merge_product_dict[key] = product_dict[key][0].get_snow_mask()
logging.info("End merge_product_at_same_date")
return merge_product_dict
#TODO move this function in app_wrappers.py along other otb applications
def gap_filling(img_in, mask_in, img_out, input_dates_file=None,
output_dates_file=None, ram=None, out_type=None):
""" Create and configure the ImageTimeSeriesGapFilling application
using otb.Registry.CreateApplication("ImageTimeSeriesGapFilling")
def compute_output_dates(date_start, date_stop, output_dates_file_path):
logging.info("Start compute_output_dates")
output_dates = []
if op.exists(output_dates_file_path):
logging.debug("Read output_date_file : {}", output_dates_file_path)
output_dates = read_list_from_file(output_dates_file_path)
else:
logging.debug("Compute output_dates from {} to {}", date_start, date_stop)
tmp_date = date_start
while tmp_date <= date_stop:
output_dates.append(datetime_to_str(tmp_date))
tmp_date += timedelta(days=1)
write_list_to_file(output_dates_file_path, output_dates)
logging.info("End compute_output_dates")
return output_dates
def load_densification_products(date_margin, date_start, date_stop, densification_path_list, path_tmp, product_dict,
ram):
logging.info("Start load_densification_products")
# load densification snow products
densification_product_dict = load_products(date_start, date_stop, date_margin, densification_path_list, None,
None)
logging.debug("Densification product dict:")
logging.debug(densification_product_dict)
# Get the footprint of the first snow product
s2_footprint_ref = product_dict[list(product_dict.keys())[0]][0].get_snow_mask()
if densification_product_dict:
# Reproject the densification products on S2 tile before going further
for densifier_product_key in list(densification_product_dict.keys()):
for densifier_product in densification_product_dict[densifier_product_key]:
original_mask = densifier_product.get_snow_mask()
reprojected_mask = op.join(path_tmp,
densifier_product.product_name + "_reprojected.tif")
if not os.path.exists(reprojected_mask):
super_impose_app = super_impose(s2_footprint_ref,
original_mask,
reprojected_mask,
"nn",
int(LABEL_NO_DATA),
ram,
otb.ImagePixelType_uint8)
super_impose_app.ExecuteAndWriteOutput()
super_impose_app = None
densifier_product.snow_mask = reprojected_mask
logging.debug(densifier_product.snow_mask)
# Add the products to extend the product_dict
if densifier_product_key in list(product_dict.keys()):
product_dict[densifier_product_key].extend(densification_product_dict[densifier_product_key])
else:
product_dict[densifier_product_key] = densification_product_dict[densifier_product_key]
Keyword arguments:
img_in -- the input timeserie image
mask_in -- the input masks
img_out -- the output image
ram -- the ram limitation (not mandatory)
out_type -- the output image pixel type (not mandatory)
"""
if img_in and mask_in and img_out:
logging.info("Processing ImageTimeSeriesGapFilling with args:")
logging.info("img_in = " + img_in)
logging.info("mask_in = " + mask_in)
logging.info("img_out = " + img_out)
gap_filling_app = otb.Registry.CreateApplication("ImageTimeSeriesGapFilling")
gap_filling_app.SetParameterString("in", img_in)
gap_filling_app.SetParameterString("mask", mask_in)
gap_filling_app.SetParameterString("out", img_out)
gap_filling_app.SetParameterInt("comp", 1)
gap_filling_app.SetParameterString("it", "linear")
if input_dates_file is not None:
logging.info("input_dates_file = " + input_dates_file)
gap_filling_app.SetParameterString("id", input_dates_file)
if output_dates_file is not None:
logging.info("output_dates_file = " + output_dates_file)
gap_filling_app.SetParameterString("od", output_dates_file)
if ram is not None:
logging.info("ram = " + str(ram))
gap_filling_app.SetParameterString("ram", str(ram))
if out_type is not None:
logging.info("out_type = " + str(out_type))
gap_filling_app.SetParameterOutputImagePixelType("out", out_type)
return gap_filling_app
logging.debug(product_dict[densifier_product_key])
else:
logging.error("Parameters img_in, img_out and mask_in are required")
logging.warning("No Densifying candidate product found!")
logging.info("End load_densification_products")
def load_products(date_start, date_stop, date_margin, snow_products_list, tile_id=None, product_type=None):
logging.info("Start load_products")
logging.info("Parsing provided snow products list")
product_dict = {}
search_start_date = date_start - date_margin
search_stop_date = date_stop + date_margin
for product_path in snow_products_list:
try:
product = load_snow_product(str(product_path))
logging.info(str(product))
current_day = datetime_to_str(product.acquisition_date)
test_result = True
if search_start_date > product.acquisition_date or \
search_stop_date < product.acquisition_date:
test_result = False
if (tile_id is not None) and (tile_id not in product.tile_id):
test_result = False
if (product_type is not None) and (product_type not in product.platform):
test_result = False
if test_result:
if current_day not in list(product_dict.keys()):
product_dict[current_day] = [product]
else:
product_dict[current_day].append(product)
logging.info("Keeping: {}", str(product))
else:
logging.warning("Discarding: {}", str(product))
except Exception:
logging.error("Unable to load product: {}", product_path)
logging.debug("Product dictionnary:")
logging.debug(product_dict)
logging.info("End load_products")
return product_dict
def convert_mask_list(path_tmp, resulting_snow_mask_dict, expression, type_name, ram, mask_format=""):
binary_mask_list = []
for mask_date in sorted(resulting_snow_mask_dict):
binary_mask = op.join(path_tmp,
mask_date + "_" + type_name + "_binary.tif")
binary_mask = extract_binary_mask(resulting_snow_mask_dict[mask_date],
binary_mask,
expression,
ram,
mask_format)
binary_mask_list.append(binary_mask)
return binary_mask_list
def extract_binary_mask(mask_in, mask_out, expression, ram, mask_format=""):
bandMathApp = band_math([mask_in],
mask_out + mask_format,
expression,
ram,
otb.ImagePixelType_uint8)
bandMathApp.ExecuteAndWriteOutput()
return mask_out
def merge_masks_at_same_date(snow_product_list, merged_snow_product, threshold=100, ram=None):
""" This function implement the fusion of multiple snow mask
@@ -109,7 +307,7 @@ def merge_masks_at_same_date(snow_product_list, merged_snow_product, threshold=1
threshold -- the threshold between valid <= invalid data
ram -- the ram limitation (not mandatory)
"""
logging.info("Merging products into " + merged_snow_product)
logging.info("Merging products into {}", merged_snow_product)
# the merging is performed according the following selection:
# if img1 < threshold use img1 data
@@ -118,12 +316,12 @@ def merge_masks_at_same_date(snow_product_list, merged_snow_product, threshold=1
# the order of the images in the input list is important:
# we expect to have first the main input products
# and then the densification products
img_index = list(range(1, len(snow_product_list)+1))
img_index = list(range(1, len(snow_product_list) + 1))
expression_merging = "".join(["(im" + str(i) + "b1<=" + str(threshold) + "?im" + str(i) + "b1:" for i in img_index])
expression_merging += "im"+str(img_index[-1])+"b1"
expression_merging += "im" + str(img_index[-1]) + "b1"
expression_merging += "".join([")" for i in img_index])
img_list= [i.get_snow_mask() for i in snow_product_list]
img_list = [i.get_snow_mask() for i in snow_product_list]
bandMathApp = band_math(img_list,
merged_snow_product,
expression_merging,
@@ -131,316 +329,3 @@ def merge_masks_at_same_date(snow_product_list, merged_snow_product, threshold=1
otb.ImagePixelType_uint8)
bandMathApp.ExecuteAndWriteOutput()
bandMathApp = None
""" This module provide the implementation of the snow annual map """
class snow_annual_map():
def __init__(self, params):
logging.info("Init snow_multitemp")
self.tile_id = params.get("tile_id")
self.date_start = str_to_datetime(params.get("date_start"), "%d/%m/%Y")
self.date_stop = str_to_datetime(params.get("date_stop"), "%d/%m/%Y")
self.date_margin = timedelta(days=params.get("date_margin", 0))
self.output_dates_filename = params.get("output_dates_filename", None)
self.mode = params.get("mode", "RUNTIME")
# Compute an id like T31TCH_20170831_20180901 to label the map
self.processing_id = str(self.tile_id + "_" + \
datetime_to_str(self.date_start) + "_" + \
datetime_to_str(self.date_stop))
# Retrive the input_products_list
self.input_path_list = params.get("input_products_list", [])
# @TODO an available path_tmp must be provide or the TMPDIR variable must be avaible
self.path_tmp = str(params.get("path_tmp", os.environ.get('TMPDIR')))
if not os.path.exists(self.path_tmp):
logging.error(self.path_tmp + ", the target does not exist and can't be used for processing")
self.path_out = op.join(str(params.get("path_out")), self.processing_id)
if not os.path.exists(self.path_out):
os.mkdir(self.path_out)
self.ram = params.get("ram", 512)
self.nbThreads = params.get("nbThreads", None)
self.use_densification = params.get("use_densification", False)
if self.use_densification:
self.densification_path_list = params.get("densification_products_list", [])
# Define label for output snow product (cf snow product format)
self.label_no_snow = "0"
self.label_snow = "100"
self.label_cloud = "205"
self.label_no_data = "255"
self.label_no_data_old = "254"
# Build useful paths
self.input_dates_filename = op.join(self.path_tmp, "input_dates.txt")
if not self.output_dates_filename:
self.output_dates_filename = op.join(self.path_tmp, "output_dates.txt")
self.multitemp_snow_vrt = op.join(self.path_tmp, "multitemp_snow_mask.vrt")
self.multitemp_cloud_vrt = op.join(self.path_tmp, "multitemp_cloud_mask.vrt")
self.gapfilled_timeserie = op.join(self.path_tmp, "DAILY_SNOW_MASKS_" + self.processing_id + ".tif")
self.annual_snow_map = op.join(self.path_tmp, "SCD_" + self.processing_id + ".tif")
self.cloud_occurence_img = op.join(self.path_tmp, "CLOUD_OCCURENCE_" + self.processing_id +".tif")
self.metadata_path = op.join(self.path_out, "LIS_METADATA.XML")
def run(self):
logging.info("Run snow_annual_map")
# Set maximum ITK threads
if self.nbThreads:
os.environ["ITK_GLOBAL_DEFAULT_NUMBER_OF_THREADS"] = str(self.nbThreads)
# search matching snow product
self.product_dict = self.load_products(self.input_path_list, self.tile_id, None)
logging.debug("Product dictionnary:")
logging.debug(self.product_dict)
# Exiting with error if none of the input products were loaded
if not self.product_dict:
logging.error("Empty product list!")
return
# Do the loading of the products to densify the timeserie
if self.use_densification:
# load densification snow products
densification_product_dict = self.load_products(self.densification_path_list, None, None)
logging.info("Densification product dict:")
logging.info(densification_product_dict)
# Get the footprint of the first snow product
s2_footprint_ref = self.product_dict[list(self.product_dict.keys())[0]][0].get_snow_mask()
if densification_product_dict:
# Reproject the densification products on S2 tile before going further
for densifier_product_key in list(densification_product_dict.keys()):
for densifier_product in densification_product_dict[densifier_product_key]:
original_mask = densifier_product.get_snow_mask()
reprojected_mask = op.join(self.path_tmp,
densifier_product.product_name + "_reprojected.tif")
if not os.path.exists(reprojected_mask):
super_impose_app = super_impose(s2_footprint_ref,
original_mask,
reprojected_mask,
"nn",
int(self.label_no_data),
self.ram,
otb.ImagePixelType_uint8)
super_impose_app.ExecuteAndWriteOutput()
super_impose_app = None
densifier_product.snow_mask = reprojected_mask
logging.debug(densifier_product.snow_mask)
# Add the products to extend the self.product_dict
if densifier_product_key in list(self.product_dict.keys()):
self.product_dict[densifier_product_key].extend(densification_product_dict[densifier_product_key])
else:
self.product_dict[densifier_product_key] = densification_product_dict[densifier_product_key]
else:
logging.warning("No Densifying candidate product found!")
# re-order products according acquisition date
input_dates = sorted(self.product_dict.keys())
write_list_to_file(self.input_dates_filename, input_dates)
# compute or retrieve the output dates
output_dates = []
if op.exists(self.output_dates_filename):
output_dates = read_list_from_file(self.output_dates_filename)
else:
tmp_date = self.date_start
while tmp_date <= self.date_stop:
output_dates.append(datetime_to_str(tmp_date))
tmp_date += timedelta(days=1)
write_list_to_file(self.output_dates_filename, output_dates)
shutil.copy2(self.input_dates_filename, self.path_out)
shutil.copy2(self.output_dates_filename, self.path_out)
# merge products at the same date
self.resulting_snow_mask_dict={}
for key in list(self.product_dict.keys()):
if len(self.product_dict[key]) > 1:
merged_mask = op.join(self.path_tmp, key + "_merged_snow_product.tif")
merge_masks_at_same_date(self.product_dict[key],
merged_mask,
self.label_snow,
self.ram)
self.resulting_snow_mask_dict[key] = merged_mask
else:
self.resulting_snow_mask_dict[key] = self.product_dict[key][0].get_snow_mask()
# convert the snow masks into binary snow masks
expression = "(im1b1==" + self.label_snow + ")?1:0"
self.binary_snowmask_list = self.convert_mask_list(expression, "snow", GDAL_OPT)
logging.debug("Binary snow mask list:")
logging.debug(self.binary_snowmask_list)
# convert the snow masks into binary cloud masks
expression = "im1b1=="+self.label_cloud+"?1:(im1b1=="+self.label_no_data+"?1:(im1b1=="+self.label_no_data_old+"?1:0))"
self.binary_cloudmask_list = self.convert_mask_list(expression, "cloud", GDAL_OPT)
logging.debug("Binary cloud mask list:")
logging.debug(self.binary_cloudmask_list)
# build cloud mask vrt
logging.info("Building multitemp cloud mask vrt")
logging.info("cloud vrt: " + self.multitemp_cloud_vrt)
gdal.BuildVRT(self.multitemp_cloud_vrt,
self.binary_cloudmask_list,
separate=True,
srcNodata='None')
shutil.copy2(self.multitemp_cloud_vrt, self.path_out)
# generate the summary map
band_index = list(range(1, len(self.binary_cloudmask_list)+1))
expression = "+".join(["im1b" + str(i) for i in band_index])
bandMathApp = band_math([self.multitemp_cloud_vrt],
self.cloud_occurence_img,
expression,
self.ram,
otb.ImagePixelType_uint16)
bandMathApp.ExecuteAndWriteOutput()
bandMathApp = None
logging.info("Copying outputs from tmp to output folder")
shutil.copy2(self.cloud_occurence_img, self.path_out)
# build snow mask vrt
logging.info("Building multitemp snow mask vrt")
logging.info("snow vrt: " + self.multitemp_snow_vrt)
gdal.BuildVRT(self.multitemp_snow_vrt,
self.binary_snowmask_list,
separate=True,
srcNodata='None')
# multiply by 100 for the temporal interpolation
logging.info("Scale by 100 multitemp snow mask vrt")
multitemp_snow100 = op.join(self.path_tmp, "multitemp_snow100.tif")
bandMathXApp = band_mathX([self.multitemp_snow_vrt],
multitemp_snow100,
"im1 mlt 100",
self.ram,
otb.ImagePixelType_uint8)
bandMathXApp.ExecuteAndWriteOutput()
bandMathXApp = None
# gap filling the snow timeserie
multitemp_snow100_gapfilled = op.join(self.path_tmp, "multitemp_snow100_gapfilled.tif")
app_gap_filling = gap_filling(multitemp_snow100,
self.multitemp_cloud_vrt,
multitemp_snow100_gapfilled+"?&gdal:co:COMPRESS=DEFLATE",
self.input_dates_filename,
self.output_dates_filename,
self.ram,
otb.ImagePixelType_uint8)
# @TODO the mode is for now forced to DEBUG in order to generate img on disk
#img_in = get_app_output(app_gap_filling, "out", self.mode)
#if self.mode == "DEBUG":
#shutil.copy2(self.gapfilled_timeserie, self.path_out)
#app_gap_filling = None
img_in = get_app_output(app_gap_filling, "out", "DEBUG")
app_gap_filling = None
# threshold to 0 or 1
logging.info("Round to binary series of snow occurrence")
bandMathXApp = band_mathX([img_in],
self.gapfilled_timeserie+GDAL_OPT,
"(im1 mlt 2) dv 100",
self.ram,
otb.ImagePixelType_uint8)
bandMathXApp.ExecuteAndWriteOutput()
bandMathXApp = None
# generate the annual map
band_index = list(range(1, len(output_dates)+1))
expression = "+".join(["im1b" + str(i) for i in band_index])
bandMathApp = band_math([self.gapfilled_timeserie],
self.annual_snow_map,
expression,
self.ram,
otb.ImagePixelType_uint16)
bandMathApp.ExecuteAndWriteOutput()
bandMathApp = None
logging.info("Moving outputs from tmp to output folder")
shutil.copy2(self.annual_snow_map, self.path_out)
shutil.copy2(self.gapfilled_timeserie, self.path_out)
os.remove(self.annual_snow_map)
os.remove(self.gapfilled_timeserie)
logging.info("End of snow_annual_map")
if self.mode == "DEBUG":
dest_debug_dir = op.join(self.path_out, "tmpdir")
if op.exists(dest_debug_dir):
shutil.rmtree(dest_debug_dir)
shutil.copytree(self.path_tmp, dest_debug_dir)
def create_metadata(self):
# Compute and create the content for the product metadata file.
root = etree.Element("INPUTS_LIST")
for product_path in self.input_path_list:
product_name = basename(str(product_path))
etree.SubElement(root, "PRODUCT_NAME").text = product_name
et = etree.ElementTree(root)
et.write(self.metadata_path, pretty_print=True)
def load_products(self, snow_products_list, tile_id=None, product_type=None):
logging.info("Parsing provided snow products list")
product_dict = {}
search_start_date = self.date_start - self.date_margin
search_stop_date = self.date_stop + self.date_margin
for product_path in snow_products_list:
try:
product = load_snow_product(str(product_path))
logging.info(str(product))
current_day = datetime_to_str(product.acquisition_date)
test_result = True
if search_start_date > product.acquisition_date or \
search_stop_date < product.acquisition_date:
test_result = False
if (tile_id is not None) and (tile_id not in product.tile_id):
test_result = False
if (product_type is not None) and (product_type not in product.platform):
test_result = False
if test_result:
if current_day not in list(product_dict.keys()):
product_dict[current_day] = [product]
else:
product_dict[current_day].append(product)
logging.info("Keeping: " + str(product))
else:
logging.warning("Discarding: " + str(product))
except Exception:
logging.error("Unable to load product :" + product_path)
return product_dict
def convert_mask_list(self, expression, type_name, mask_format=""):
binary_mask_list = []
for mask_date in sorted(self.resulting_snow_mask_dict):
binary_mask = op.join(self.path_tmp,
mask_date + "_" + type_name + "_binary.tif")
binary_mask = self.extract_binary_mask(self.resulting_snow_mask_dict[mask_date],
binary_mask,
expression,
mask_format)
binary_mask_list.append(binary_mask)
return binary_mask_list
def extract_binary_mask(self, mask_in, mask_out, expression, mask_format=""):
bandMathApp = band_math([mask_in],
mask_out + mask_format,
expression,
self.ram,
otb.ImagePixelType_uint8)
bandMathApp.ExecuteAndWriteOutput()
return mask_out