def setUp(self):
'''setup'''
self.input_xml = os.path.join(self.validation_path,
'LT50420342011119PAC01.xml')
# Specify the XML metadata file defining the data to process
self.processor = AuxNARRGribProcessor(self.input_xml)
# Process the associated AUX data
self.processor.extract_aux_data()
def generate_lst(xml_filename,
only_extract_aux_data=False,
keep_lst_temp_data=False,
keep_intermediate_data=False,
debug=False):
'''
Description:
Provides the glue code for generating LST products.
'''
# Get the logger
logger = logging.getLogger(__name__)
# Retrieval and initial processing of the required AUX data
try:
logger.info('Extracting LST AUX data')
current_processor = AuxNARRGribProcessor(xml_filename)
current_processor.extract_aux_data()
except Exception:
logger.error('Failed processing auxillary NARR data')
raise
if only_extract_aux_data:
logger.info('Stopping - User requested to stop after extracting'
' LST AUX data')
return
# Generate the thermal, upwelled, and downwelled radiance bands as well as
# the atmospheric transmittance band
cmd = ['lst_intermediate_data',
'--xml', xml_filename,
'--verbose']
if debug:
cmd.append('--debug')
cmd = ' '.join(cmd)
output = ''
try:
logger.info('Calling [{0}]'.format(cmd))
output = util.System.execute_cmd(cmd)
except Exception:
logger.error('Failed creating intermediate data')
raise
finally:
if len(output) > 0:
logger.info(output)
# Generate Estimated Landsat Emissivity band
try:
current_processor = (
estimate_landsat_emissivity.EstimateLandsatEmissivity(
xml_filename, keep_intermediate_data))
current_processor.generate_product()
except Exception:
logger.error('Failed creating Estimated Landsat Emissivity data')
raise
# Generate Land Surface Temperature band
try:
current_processor = build_lst_data.BuildLSTData(xml_filename)
current_processor.generate_data()
except Exception:
logger.error('Failed processing Land Surface Temperature')
raise
# Cleanup
if not keep_intermediate_data:
# Remove the grib extraction directories
shutil.rmtree('HGT_1', ignore_errors=True)
shutil.rmtree('HGT_2', ignore_errors=True)
shutil.rmtree('SPFH_1', ignore_errors=True)
shutil.rmtree('SPFH_2', ignore_errors=True)
shutil.rmtree('TMP_1', ignore_errors=True)
shutil.rmtree('TMP_2', ignore_errors=True)
# Remove the point directories generated during the core processing
remove_dirs = set()
point_filename = 'point_list.txt'
with open(point_filename, 'r') as point_list_fd:
remove_dirs = set(list([line.strip()
for line in point_list_fd.readlines()]))
for dirname in remove_dirs:
shutil.rmtree(dirname, ignore_errors=False)
# Finally remove the file
os.unlink(point_filename)
if not keep_lst_temp_data:
util.Metadata.remove_products(xml_filename, ['lst_temp'])
