def process_ifgramStack(inps, track, startDate, endDate, outdir, tmpdir):
"""process ifgramStack.h5 file"""
if os.path.exists("".join(tmpdir)) == False:
os.makedirs("".join(tmpdir))
atr_asc = ['./inputs/' + str(inps.DataType) + '.h5']
# mask
if inps.maskfile and inps.maskfile != 'None':
date_msk1 = re.split(r'[_.]', str(inps.maskfile))[0]
date_msk1f = multitrack_utilities.find_nearest_date(
os.getcwd(), str(date_msk1))
date_msk2 = re.split(r'[_.]', str(inps.maskfile))[1]
date_msk2f = multitrack_utilities.find_nearest_date(
os.getcwd(), str(date_msk2))
maskfile = date_msk1f + '_' + date_msk2f + '.cor'
cmd_args = [atr_asc, '-m', maskfile]
print('mask.py', cmd_args)
args_str = multitrack_utilities.seperate_str_byspace(cmd_args)
os.system(
multitrack_utilities.seperate_str_byspace(
['mask.py', args_str.split()]))
atr_asc = './inputs/' + multitrack_utilities.seprate_filename_extension(
"".join(atr_asc))[1] + '_msk.h5'
#geocode ifgramStack file
cmd_args = [
atr_asc, '-b', inps.SNWE, '-y', inps.latStep, '-x', inps.lonStep,
'--outdir', "".join(tmpdir)
]
print("geocode.py", cmd_args)
args_str = multitrack_utilities.seperate_str_byspace(cmd_args)
mintpy.geocode.main(args_str.split())
#save dataset of unw
os.chdir("".join(tmpdir))
filename, extension = multitrack_utilities.seprate_filename_extension(
"".join(atr_asc))[1:3]
cmd_args = [
'geo_' + filename + extension,
"".join(['unwrapPhase-', startDate, '_', endDate])
]
print("save_roipac.py", cmd_args)
args_str = multitrack_utilities.seperate_str_byspace(cmd_args)
os.system(
multitrack_utilities.seperate_str_byspace(
['save_roipac.py', args_str.split()]))
copy_file('geo', 'unw', outdir)
shutil.rmtree(tmpdir)
def process_HDFEOS(inps):
"""process S1*.h5 file"""
atr_asc = inps.file
#save dataset of unw and geometry
filename, extension = multitrack_utilities.seprate_filename_extension(
"".join(atr_asc))[1:3]
cmd_args = [
filename + extension,
"".join(['displacement-', inps.startDate, '_', inps.endDate]), '-o',
"".join(['geo_', inps.startDate, '_', inps.endDate, '.unw'])
]
print("save_roipac.py", cmd_args)
asct_str = multitrack_utilities.seperate_str_byspace(cmd_args)
os.system(
multitrack_utilities.seperate_str_byspace(
['save_roipac.py', asct_str.split()]))
#mv
if os.path.exists(inps.outdir) == False:
os.makedirs(inps.outdir)
else:
shutil.rmtree(inps.outdir)
os.makedirs(inps.outdir)
key1 = 'geo_'
for file in os.listdir(os.getcwd()):
if str.find(file, key1) != -1:
shutil.move(file, inps.outdir)
def process_ifgramStack(inps):
"""process ifgramStack.h5 file"""
if os.path.exists("".join(inps.outdir)) == False:
os.makedirs("".join(inps.outdir))
#ifgramStack file
atr_asc = ['./inputs/' + inps.file]
cmd_args = [
atr_asc, '-b', inps.SNWE, '-y', inps.latStep, '-x', inps.lonStep,
'--outdir', "".join(inps.outdir)
]
print("geocode.py", cmd_args)
args_str = multitrack_utilities.seperate_str_byspace(cmd_args)
mintpy.geocode.main(args_str.split())
#save dataset of unw
os.chdir("".join(inps.outdir))
filename, extension = multitrack_utilities.seprate_filename_extension(
"".join(atr_asc))[1:3]
cmd_args = [
'geo_' + filename + extension,
"".join(['unwrapPhase-', inps.startDate, '_', inps.endDate])
]
print("save_roipac.py", cmd_args)
asct_str = multitrack_utilities.seperate_str_byspace(cmd_args)
os.system(
multitrack_utilities.seperate_str_byspace(
['save_roipac.py', asct_str.split()]))
def process_timeseries(inps):
"""geocode timeseries**.h5 file and get the deformation field of two time periods"""
atr_asc = inps.file
#unw file
cmd_args = [
atr_asc, '-b', inps.SNWE, '-y', inps.latStep, '-x', inps.lonStep,
'--outdir', "".join(inps.outdir)
]
print("geocode.py", cmd_args)
args_str = multitrack_utilities.seperate_str_byspace(cmd_args)
mintpy.geocode.main(args_str.split())
#save dataset of unw
os.chdir("".join(inps.outdir))
filename, extension = multitrack_utilities.seprate_filename_extension(
"".join(atr_asc))[1:3]
cmd_args = [
'geo_' + filename + extension,
"".join([inps.startDate, '_', inps.endDate])
]
print("save_roipac.py", cmd_args)
asct_str = multitrack_utilities.seperate_str_byspace(cmd_args)
os.system(
multitrack_utilities.seperate_str_byspace(
['save_roipac.py', asct_str.split()]))
def process_velocity(inps, track, startDate, endDate, outdir, tmpdir):
"""process velocity.h5 file"""
if os.path.exists("".join(tmpdir)) == False:
os.makedirs("".join(tmpdir))
atr_asc = "".join([str(inps.DataType) + '.h5'])
# mask
if inps.maskfile and inps.maskfile != 'None':
cmd_args = [atr_asc, '-m', inps.maskfile]
print('mask.py', cmd_args)
args_str = multitrack_utilities.seperate_str_byspace(cmd_args)
os.system(
multitrack_utilities.seperate_str_byspace(
['mask.py', args_str.split()]))
atr_asc = multitrack_utilities.seprate_filename_extension(
atr_asc)[1] + '_msk.h5'
print(atr_asc)
#geocode velocity file
cmd_args = [
atr_asc, '-b', inps.SNWE, '-y', inps.latStep, '-x', inps.lonStep,
'--outdir', "".join(tmpdir)
]
print("geocode.py", cmd_args)
args_str = multitrack_utilities.seperate_str_byspace(cmd_args)
mintpy.geocode.main(args_str.split())
# save dataset of unw
os.chdir("".join(tmpdir))
print('save unw file')
atr_asc = 'geo_' + atr_asc
velocity_displacement(atr_asc, startDate, endDate)
copy_file('geo', 'unw', outdir)
return atr_asc
def process_HDFEOS(inps, track, startDate, endDate, outdir, tmpdir):
"""process S1*.h5 file"""
atr_asc = multitrack_utilities.find_HDFEOS_fullname("".join(
[os.getenv('SCRATCHDIR') + '/' + track + '/mintpy/']))
#save dataset of unw
filename, extension = multitrack_utilities.seprate_filename_extension(
"".join(atr_asc))[1:3]
cmd_args = [
filename + extension,
"".join(['displacement-', startDate, '_', endDate]), '-o',
"".join(['geo_', startDate, '_', endDate, '.unw'])
]
print("save_roipac.py", cmd_args)
args_str = multitrack_utilities.seperate_str_byspace(cmd_args)
os.system(
multitrack_utilities.seperate_str_byspace(
['save_roipac.py', args_str.split()]))
# mv
if os.path.exists(tmpdir) == False:
os.makedirs(tmpdir)
else:
shutil.rmtree(tmpdir)
os.makedirs(tmpdir)
key1 = 'geo_'
for file in os.listdir(os.getcwd()):
if str.find(file, key1) != -1:
shutil.move(file, tmpdir)
copy_file('geo', 'unw', outdir)
shutil.rmtree(tmpdir)
def process_timeseries(inps, track, startDate, endDate, outdir, tmpdir):
"""geocode timeseries**.h5 file and get the deformation field of two time periods"""
if os.path.exists("".join(tmpdir)) == False:
os.makedirs("".join(tmpdir))
atr_asc = multitrack_utilities.find_timeseries_horzvert("".join(
[os.getenv('SCRATCHDIR') + '/' + track + '/mintpy/']))
# mask
if inps.maskfile and inps.maskfile != 'None':
cmd_args = [atr_asc, '-m', inps.maskfile]
print('mask.py', cmd_args)
args_str = multitrack_utilities.seperate_str_byspace(cmd_args)
os.system(
multitrack_utilities.seperate_str_byspace(
['mask.py', args_str.split()]))
atr_asc = multitrack_utilities.seprate_filename_extension(
atr_asc)[1] + '_msk.h5'
#geocode timeseries file
cmd_args = [
atr_asc, '-b', inps.SNWE, '-y', inps.latStep, '-x', inps.lonStep,
'--outdir', tmpdir
]
print("geocode.py", cmd_args)
args_str = multitrack_utilities.seperate_str_byspace(cmd_args)
mintpy.geocode.main(args_str.split())
#save dataset of unw
os.chdir("".join(tmpdir))
filename, extension = multitrack_utilities.seprate_filename_extension(
"".join(atr_asc))[1:3]
cmd_args = [
'geo_' + filename + extension, "".join([startDate, '_', endDate])
]
print("save_roipac.py", cmd_args)
args_str = multitrack_utilities.seperate_str_byspace(cmd_args)
os.system(
multitrack_utilities.seperate_str_byspace(
['save_roipac.py', args_str.split()]))
copy_file('geo', 'unw', outdir)
shutil.rmtree(tmpdir)
def process_HDFEOS(inps):
"""process *.he5 file"""
atr_asc = inps.file
#save dataset of unw cor and dem
filename, extension = multitrack_utilities.seprate_filename_extension(
"".join(atr_asc))[1:3]
cmd_args = [
filename + extension,
"".join(['displacement-', inps.startDate, '_', inps.endDate]), '-o',
"".join(['geo_', inps.startDate, '_', inps.endDate, '.unw'])
]
print("save_roipac.py", cmd_args)
asct_str = multitrack_utilities.seperate_str_byspace(cmd_args)
os.system(
multitrack_utilities.seperate_str_byspace(
['save_roipac.py', asct_str.split()]))
cmd_args = [
filename + extension, 'temporalCoherence', '-o',
"".join(['geo_', inps.startDate, '_', inps.endDate, '.cor'])
]
print("save_roipac.py", cmd_args)
asct_str = multitrack_utilities.seperate_str_byspace(cmd_args)
completion_status = os.system(
multitrack_utilities.seperate_str_byspace(
['save_roipac.py', asct_str.split()]))
cmd_args = [filename + extension, 'height', '-o', 'srtm.dem']
print("save_roipac.py", cmd_args)
asct_str = multitrack_utilities.seperate_str_byspace(cmd_args)
os.system(
multitrack_utilities.seperate_str_byspace(
['save_roipac.py', asct_str.split()]))
# mv
if os.path.exists(inps.outdir) == False:
os.makedirs(inps.outdir)
else:
shutil.rmtree(inps.outdir)
os.makedirs(inps.outdir)
key1 = 'geo_'
key2 = 'srtm'
for file in os.listdir(os.getcwd()):
if str.find(file, key1) != -1 or str.find(file, key2) != -1:
shutil.move(file, inps.outdir)
os.chdir("".join(inps.outdir))
multitrack_utilities.delete_tmpgeo(inps.outdir, 'geo_', '.h5')
def run_software(inpsdict, inps):
"""run save_geodmod/gbis.py in proper directory"""
if inpsdict['DataSet'][0] == 'None':
tempfilename = inps.templateFile
folders = multitrack_utilities.find_folder(
multitrack_utilities.seprate_filename_extension(tempfilename)[1])
print(folders)
else:
folders = inpsdict['DataSet']
print(folders)
# if datatype is HDFEOS,first judge whether they have same lat_step and lon_step
if inpsdict['DataType'] == 'HDFEOS':
flag = multitrack_utilities.check_X_Y_step(folders)
if flag:
multitrack_run_software(inpsdict, folders, inps)
else:
multitrack_run_software(inpsdict, folders, inps)
def process_HDFEOS(inps):
"""process *.he5 file"""
atr_asc = inps.file
#save dataset of unw cor and dem
filename, extension = mu.seprate_filename_extension("".join(atr_asc))[1:3]
output_unw = "".join(['geo_', inps.startDate, '_', inps.endDate, '.unw'])
cmd_args = [
'../' + filename + extension,
"".join(['displacement-', inps.startDate, '_',
inps.endDate]), '-o', output_unw
]
print("save_roipac.py", cmd_args)
asct_str = mu.seperate_str_byspace(cmd_args)
os.system(mu.seperate_str_byspace(['save_roipac.py', asct_str.split()]))
output_cor = "".join(['geo_', inps.startDate, '_', inps.endDate, '.cor'])
cmd_args = [
'../' + filename + extension, 'temporalCoherence', '-o', output_cor
]
print("save_roipac.py", cmd_args)
asct_str = mu.seperate_str_byspace(cmd_args)
completion_status = os.system(
mu.seperate_str_byspace(['save_roipac.py',
asct_str.split()]))
output_dem = 'srtm.dem'
cmd_args = ['../' + filename + extension, 'height', '-o', output_dem]
print("save_roipac.py", cmd_args)
asct_str = mu.seperate_str_byspace(cmd_args)
os.system(mu.seperate_str_byspace(['save_roipac.py', asct_str.split()]))
mask_filter(inps, output_unw)
mask_filter(inps, output_cor)
#mask_filter(inps,output_dem)
if inps.SNWE:
print('Subset data based on bbox')
subset_data_based_bbox(inps, output_unw)
subset_data_based_bbox(inps, output_cor)
subset_data_based_bbox(inps, output_dem)
def process_ifgramStack(inps):
"""process ifgramStack.h5 file"""
if os.path.exists("".join(inps.outdir)) == False:
os.makedirs("".join(inps.outdir))
# dem file
demname = 'geometryRadar.h5'
if not os.path.isfile(demname):
demname_f = './inputs/geometryRadar.h5'
else:
demname_f = 'geometryRadar.h5'
cmd_args = [
demname_f, '-b', inps.SNWE, '-y', inps.latStep, '-x', inps.lonStep,
'--outdir', "".join(inps.outdir)
]
print("geocode.py", cmd_args)
args_str = multitrack_utilities.seperate_str_byspace(cmd_args)
mintpy.geocode.main(args_str.split())
#ifgramStack file
atr_asc = ['./inputs/' + inps.file]
cmd_args = [
atr_asc, '-b', inps.SNWE, '-y', inps.latStep, '-x', inps.lonStep,
'--outdir', "".join(inps.outdir)
]
print("geocode.py", cmd_args)
args_str = multitrack_utilities.seperate_str_byspace(cmd_args)
mintpy.geocode.main(args_str.split())
#save dataset of unw cor and dem
os.chdir("".join(inps.outdir))
filename, extension = multitrack_utilities.seprate_filename_extension(
"".join(atr_asc))[1:3]
cmd_args = [
'geo_' + filename + extension,
"".join(['unwrapPhase-', inps.startDate, '_', inps.endDate])
]
print("save_roipac.py", cmd_args)
asct_str = multitrack_utilities.seperate_str_byspace(cmd_args)
os.system(
multitrack_utilities.seperate_str_byspace(
['save_roipac.py', asct_str.split()]))
cmd_args = [
'geo_' + filename + extension,
"".join(['coherence-', inps.startDate, '_', inps.endDate])
]
print("save_roipac.py", cmd_args)
asct_str = multitrack_utilities.seperate_str_byspace(cmd_args)
completion_status = os.system(
multitrack_utilities.seperate_str_byspace(
['save_roipac.py', asct_str.split()]))
cmd_args = ['geo_geometryRadar.h5', 'height', '-o', 'srtm.dem']
print("save_roipac.py", cmd_args)
asct_str = multitrack_utilities.seperate_str_byspace(cmd_args)
os.system(
multitrack_utilities.seperate_str_byspace(
['save_roipac.py', asct_str.split()]))
multitrack_utilities.delete_tmpgeo(inps.outdir, 'geo_', '.h5')