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 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([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()])) process_saveroi(inps) multitrack_utilities.delete_tmpgeo(inps.outdir, 'geo_', '.h5')
def main(iargs=None): inps = cmd_line_parse(iargs) print('single track!') inps.startDate, inps.endDate = multitrack_utilities.find_start_end_date( os.getcwd(), inps) print(inps) if str.find(inps.file, 'ifgramStack') != -1: process_ifgramStack(inps) elif str.find(inps.file, 'velocity') != -1: process_geocode(inps) process_velocity(inps) elif str.find(inps.file, 'timeseries') != -1: process_geocode(inps) process_timeseries(inps) else: process_HDFEOS(inps) # rename *.rsc1 to *.rsc outfile = multitrack_utilities.seperate_str_byspace(['srtm.dem' + '.rsc']) write_rsc_file( inps, outfile, multitrack_utilities.seperate_str_byspace(['srtm.dem' + '.rsc1'])) os.remove(outfile) print('rename *.rsc1 to *.rsc') os.rename( multitrack_utilities.seperate_str_byspace(['srtm.dem' + '.rsc1']), outfile) # generate dem.jpeg dem_jpeg('srtm.dem')
def process_geocode(inps): """process temporalCoherence.h5 and geometryRadar.h5 file""" # process cor and dem dataset if os.path.exists("".join(inps.outdir)) == False: os.makedirs("".join(inps.outdir)) # geocode corname = 'temporalCoherence.h5' cmd_args = [ corname, '-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()) 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())
def process_saveroi(inps): #save_roipac cmd_args = [ 'geo_temporalCoherence.h5', '-o', "".join(['geo_', inps.startDate, '_', inps.endDate, '.cor']) ] 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_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()]))
def process_velocity(inps): """process velocity.h5 file""" atr_asc = inps.file #velocity 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()) os.chdir("".join(inps.outdir)) process_saveroi(inps) print('save unw file') multitrack_utilities.velo_disp(inps) multitrack_utilities.delete_tmpgeo(inps.outdir, 'geo_', '.h5')
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 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')