#Get current directory
currdir = os.getcwd()
#Check if GIAnT dir exists. Create it if not.
if os.path.isdir(inps.prepDir):
print "{0} directory already exists".format(inps.prepDir)
else:
os.mkdir(inps.prepDir)
pairs = []
#####Get list of IFGS
if inps.ilist in ['', None]:
pairs = SS.getPairDirs(dirname=inps.srcDir)
else:
pairs = SS.pairDirs_from_file(inps.ilist, base=inps.srcDir)
#Create ifg.list
ifglist = os.path.join(inps.prepDir, 'ifg.list')
xObj = None
if (not os.path.exists(ifglist)) or inps.force:
fid = open(ifglist, 'w')
for pair in pairs:
dates = os.path.basename(pair).split('_')
print pair
xmlFile = os.path.join(pair, 'insarProc.xml')
xObj = OB.fromstring(open(xmlFile, 'r').read())
try:
bTop = xObj.baseline.perp_baseline_top
bBot = xObj.baseline.perp_baseline_bottom
inps = parser.parse_args()
return inps
if __name__ == '__main__':
'''
The main driver.
'''
inps = parse()
currDir = os.getcwd()
if inps.subset in ['', None]:
pairDirs = SS.getPairDirs(dirname=inps.dirname)
else:
pairDirs = SS.pairDirs_from_file(inps.subset, base=inps.dirname)
print 'Number of IFGs : ', len(pairDirs)
print 'Start step: ', inps.start
print 'End step: ', inps.end
for pair in pairDirs:
unwfile = os.path.join(pair, 'filt_topophase.unw')
if os.path.exists(unwfile) and (not inps.force) and (inps.start !=
'geocode'):
print 'Interferogram {0} already processed upto unwrapping'.format(
pair)
else:
print 'Processing: {0}'.format(pair)
os.chdir(pair)
#Get current directory
currdir = os.getcwd()
#Check if GIAnT dir exists. Create it if not.
if os.path.isdir(inps['prepDir']):
print("{0} directory already exists".format(inps['prepDir']))
else:
os.mkdir(inps['prepDir'])
pairs = []
#####Get list of IFGS
if inps['list'] in ['', None]:
pairs = SS.getPairDirs(dirname=inps['insarDir'])
else:
pairs = SS.pairDirs_from_file(inps['list'], base=inps['insarDir'])
#Create ifg.list
ifglist = os.path.join(inps['prepDir'], 'ifg.list')
xObj = None
if (not os.path.exists(ifglist)) or inps['force']:
fid = open(ifglist, 'w')
for pair in pairs:
dates = os.path.basename(pair).split('_')
xmlFile = os.path.join(pair, 'insarProc.xml')
xObj = OB.fromstring(open(xmlFile, 'r').read())
try:
bTop = xObj.baseline.perp_baseline_top
bBot = xObj.baseline.perp_baseline_bottom
bPerp = 0.5 * (bTop + bBot)