def main(argv):
try:
igramsFile = argv[0]
except:
Usage()
sys.exit(1)
try:
inversion_method = argv[1]
except:
inversion_method = 'L2'
#try:
h5file = h5py.File(igramsFile, 'r')
if not 'interferograms' in h5file.keys():
print '**********************************************************************'
print 'ERROR:'
print ' ' + igramsFile + ' was not found or the file is not readable!'
print '**********************************************************************'
Usage()
sys.exit(1)
numIfgrams = len(h5file['interferograms'].keys())
if numIfgrams == 0.:
print "load interferograms first by running: load_data.py TEMPLATEFILE"
sys.exit(1)
h5timeseries = h5py.File('timeseries.h5', 'w')
if inversion_method in ['L2', 'l2']:
ut.timeseries_inversion(h5file, h5timeseries)
elif inversion_method in ['L1', 'l1']:
ut.timeseries_inversion_L1(h5file, h5timeseries)
else:
print '******************************'
print 'WARNING:'
print ' Inversion method not recognized!'
print ' L2 norm minimization is used.'
print '******************************'
ut.timeseries_inversion(h5file, h5timeseries)
dset1 = h5file['mask'].get('mask')
Mask = dset1[0:dset1.shape[0], 0:dset1.shape[1]]
group = h5timeseries.create_group('mask')
dset = group.create_dataset('mask', data=Mask, compression='gzip')
h5file.close()
h5timeseries.close()
def main(argv):
try:
igramsFile = argv[0]
except:
Usage()
sys.exit(1)
try:
inversion_method = argv[1]
except:
inversion_method = "L2"
# try:
h5file = h5py.File(igramsFile, "r")
if not "interferograms" in h5file.keys():
print "**********************************************************************"
print "ERROR:"
print " " + igramsFile + " was not found or the file is not readable!"
print "**********************************************************************"
Usage()
sys.exit(1)
numIfgrams = len(h5file["interferograms"].keys())
if numIfgrams == 0.0:
print "load interferograms first by running: load_data.py TEMPLATEFILE"
sys.exit(1)
h5timeseries = h5py.File("timeseries.h5", "w")
if inversion_method in ["L2", "l2"]:
ut.timeseries_inversion(h5file, h5timeseries)
elif inversion_method in ["L1", "l1"]:
ut.timeseries_inversion_L1(h5file, h5timeseries)
else:
print "******************************"
print "WARNING:"
print " Inversion method not recognized!"
print " L2 norm minimization is used."
print "******************************"
ut.timeseries_inversion(h5file, h5timeseries)
dset1 = h5file["mask"].get("mask")
Mask = dset1[0 : dset1.shape[0], 0 : dset1.shape[1]]
group = h5timeseries.create_group("mask")
dset = group.create_dataset("mask", data=Mask, compression="gzip")
h5file.close()
h5timeseries.close()
def main(argv):
inps = cmdLineParse()
# Input file info
atr = readfile.read_attribute(inps.ifgram_file)
k = atr['FILE_TYPE']
if not k == 'interferograms':
sys.exit('ERROR: only interferograms file supported, input is ' + k +
' file!')
# Network Inversion
if not inps.inverse_method == 'L1':
print('Inverse time-series using L2 norm minimization')
ut.timeseries_inversion(inps.ifgram_file, inps.timeseries_file)
else:
print('Inverse time-series using L1 norm minimization')
ut.timeseries_inversion_L1(inps.ifgram_file, inps.timeseries_file)
return inps.timeseries_file
def main(argv):
inversion_method = 'L2'
try:
opts, args = getopt.getopt(argv,"h:f:m:l:")
except getopt.GetoptError:
Usage() ; sys.exit(1)
if len(sys.argv)>2:
for opt,arg in opts:
if opt in ("-h","--help"): Usage(); sys.exit()
elif opt == '-f': igramsFile = arg
elif opt == '-m': maskFile = arg
elif opt == '-l': inversion_method = arg
elif len(sys.argv)==2:
if os.path.isfile(argv[0]): igramsFile = argv[0]
else: Usage(); sys.exit(1)
elif len(sys.argv)<2: Usage(); sys.exit(1)
#try: igramsFile = argv[0]
#except: Usage() ; sys.exit(1)
#try: inversion_method = argv[1]
#except: inversion_method = 'L2'
#try:
h5file = h5py.File(igramsFile,'r')
if not 'interferograms' in h5file.keys():
print '**********************************************************************'
print 'ERROR:'
print ' '+igramsFile+ ' was not found or the file is not readable!'
print '**********************************************************************'
Usage();sys.exit(1)
numIfgrams = len(h5file['interferograms'].keys())
if numIfgrams == 0.:
print "load interferograms first by running: load_data.py TEMPLATEFILE"
sys.exit(1)
h5timeseries = h5py.File('timeseries.h5','w')
if inversion_method in ['L2','l2']:
ut.timeseries_inversion(h5file,h5timeseries)
elif inversion_method in ['L1','l1']:
ut.timeseries_inversion_L1(h5file,h5timeseries)
else:
print '******************************'
print 'WARNING:'
print ' Inversion method not recognized!'
print ' L2 norm minimization is used.'
print '******************************'
ut.timeseries_inversion(h5file,h5timeseries)
## generate 'mask' for timeseries.h5
try:
maskFile
h5filem = h5py.File(maskFile,'r')
dset1 = h5filem['mask'].get('mask')
Mask = dset1[0:dset1.shape[0],0:dset1.shape[1]]
group=h5timeseries.create_group('mask')
dset = group.create_dataset('mask', data=Mask, compression='gzip')
h5filem.close()
except:
dset1 = h5file['mask'].get('mask')
Mask = dset1[0:dset1.shape[0],0:dset1.shape[1]]
group=h5timeseries.create_group('mask')
dset = group.create_dataset('mask', data=Mask, compression='gzip')
h5file.close()
h5timeseries.close()
def main(argv):
inversion_method = 'L2'
try:
opts, args = getopt.getopt(argv, "h:f:m:l:")
except getopt.GetoptError:
Usage()
sys.exit(1)
if len(sys.argv) > 2:
for opt, arg in opts:
if opt in ("-h", "--help"):
Usage()
sys.exit()
elif opt == '-f':
igramsFile = arg
elif opt == '-m':
maskFile = arg
elif opt == '-l':
inversion_method = arg
elif len(sys.argv) == 2:
if os.path.isfile(argv[0]): igramsFile = argv[0]
else:
Usage()
sys.exit(1)
elif len(sys.argv) < 2:
Usage()
sys.exit(1)
#try: igramsFile = argv[0]
#except: Usage() ; sys.exit(1)
#try: inversion_method = argv[1]
#except: inversion_method = 'L2'
#try:
h5file = h5py.File(igramsFile, 'r')
if not 'interferograms' in h5file.keys():
print '**********************************************************************'
print 'ERROR:'
print ' ' + igramsFile + ' was not found or the file is not readable!'
print '**********************************************************************'
Usage()
sys.exit(1)
numIfgrams = len(h5file['interferograms'].keys())
if numIfgrams == 0.:
print "load interferograms first by running: load_data.py TEMPLATEFILE"
sys.exit(1)
h5timeseries = h5py.File('timeseries.h5', 'w')
if inversion_method in ['L2', 'l2']:
ut.timeseries_inversion(h5file, h5timeseries)
elif inversion_method in ['L1', 'l1']:
ut.timeseries_inversion_L1(h5file, h5timeseries)
else:
print '******************************'
print 'WARNING:'
print ' Inversion method not recognized!'
print ' L2 norm minimization is used.'
print '******************************'
ut.timeseries_inversion(h5file, h5timeseries)
## generate 'mask' for timeseries.h5
try:
maskFile
h5filem = h5py.File(maskFile, 'r')
dset1 = h5filem['mask'].get('mask')
Mask = dset1[0:dset1.shape[0], 0:dset1.shape[1]]
group = h5timeseries.create_group('mask')
dset = group.create_dataset('mask', data=Mask, compression='gzip')
h5filem.close()
except:
dset1 = h5file['mask'].get('mask')
Mask = dset1[0:dset1.shape[0], 0:dset1.shape[1]]
group = h5timeseries.create_group('mask')
dset = group.create_dataset('mask', data=Mask, compression='gzip')
h5file.close()
h5timeseries.close()
