def extract_metadata4dem_geo(fname):
"""Read/extract attribute for .dem file from Gamma to ROI_PAC
For example, it read input file, sim_150911-150922.utm.dem,
find its associated par file, sim_150911-150922.utm.dem.par, read it, and
convert to ROI_PAC style and write it to an rsc file, sim_150911-150922.utm.dem.rsc
"""
atr = {}
atr['PROCESSOR'] = 'gamma'
atr['FILE_TYPE'] = os.path.splitext(fname)[1]
atr['Y_UNIT'] = 'degrees'
atr['X_UNIT'] = 'degrees'
par_file = fname+'.par'
par_dict = readfile.read_gamma_par(par_file)
atr.update(par_dict)
# Write to .rsc file
rsc_file = fname+'.rsc'
try:
atr_orig = readfile.read_roipac_rsc(rsc_file)
except:
atr_orig = dict()
if not set(atr.items()).issubset(set(atr_orig.items())):
atr_out = {**atr_orig, **atr}
print('writing >>> '+os.path.basename(rsc_file))
writefile.write_roipac_rsc(atr_out, out_file=rsc_file)
return rsc_file
def extract_metadata4dem_radar(fname):
"""Read/extract attribute for .hgt_sim file from Gamma to ROI_PAC
Input:
sim_150911-150922.hgt_sim
sim_150911-150922.rdc.dem
Search for:
sim_150911-150922.diff_par
Output:
sim_150911-150922.hgt_sim.rsc
sim_150911-150922.rdc.dem.rsc
"""
atr = {}
atr['PROCESSOR'] = 'gamma'
atr['FILE_TYPE'] = os.path.splitext(fname)[1]
# Get basename of file
fname_base = os.path.splitext(fname)[0]
for i in range(5):
fname_base = os.path.splitext(fname_base)[0]
par_file = fname_base+'.diff_par'
par_dict = readfile.read_gamma_par(par_file)
atr.update(par_dict)
# Write to .rsc file
rsc_file = fname+'.rsc'
try:
atr_orig = readfile.read_roipac_rsc(rsc_file)
except:
atr_orig = dict()
if not set(atr.items()).issubset(set(atr_orig.items())):
atr_out = {**atr_orig, **atr}
print('writing >>> '+os.path.basename(rsc_file))
writefile.write_roipac_rsc(atr_out, out_file=rsc_file)
return rsc_file
def extract_metadata4dem_geo(fname):
"""Read/extract attribute for .dem file from Gamma to ROI_PAC
For example, it read input file, sim_150911-150922.utm.dem,
find its associated par file, sim_150911-150922.utm.dem.par, read it, and
convert to ROI_PAC style and write it to an rsc file, sim_150911-150922.utm.dem.rsc
"""
atr = {}
atr['PROCESSOR'] = 'gamma'
atr['FILE_TYPE'] = os.path.splitext(fname)[1]
atr['Y_UNIT'] = 'degrees'
atr['X_UNIT'] = 'degrees'
par_file = fname + '.par'
#print('read '+os.path.basename(par_file))
par_dict = readfile.read_gamma_par(par_file)
atr.update(par_dict)
# Write to .rsc file
rsc_file = fname + '.rsc'
try:
atr_orig = readfile.read_roipac_rsc(rsc_file)
except:
atr_orig = None
#keyList = [i for i in atr_orig.keys() if i in atr.keys()]
if any((i not in atr_orig.keys() or atr_orig[i] != atr[i])
for i in atr.keys()):
print('writing >>> ' + os.path.basename(rsc_file))
writefile.write_roipac_rsc(atr, out_file=rsc_file)
return rsc_file
def extract_metadata4geometry_geo(fname):
"""Read/extract attribute for *.dem / *.UTM_TO_RDC file from Gamma to ROI_PAC
Inputs:
sim_20070813_20080310.utm.dem
sim_20070813_20080310.UTM_TO_RDC
Search for:
sim_20070813_20080310.utm.dem.par
Outputs:
sim_20070813_20080310.utm.dem.rsc
sim_20070813_20080310.UTM_TO_RDC.rsc
"""
# Get/read GAMMA par file
ext = os.path.splitext(fname)[1]
if ext in ['.UTM_TO_RDC']:
par_file = os.path.splitext(fname)[0] + '.utm.dem.par'
elif fnames[0].endswith('.utm.dem'):
par_file = fname + '.par'
par_dict = readfile.read_gamma_par(par_file)
# initiate ROIPAC dict
atr = {}
atr['PROCESSOR'] = 'gamma'
atr['FILE_TYPE'] = ext
atr['Y_UNIT'] = 'degrees'
atr['X_UNIT'] = 'degrees'
atr.update(par_dict)
# Write to .rsc file
rsc_file = fname + '.rsc'
try:
atr_orig = readfile.read_roipac_rsc(rsc_file)
except:
atr_orig = dict()
if not set(atr.items()).issubset(set(atr_orig.items())):
atr_out = {**atr_orig, **atr}
print('writing >>> ' + os.path.basename(rsc_file))
writefile.write_roipac_rsc(atr_out, out_file=rsc_file)
return rsc_file
def extract_metadata4dem_radar(fname):
"""Read/extract attribute for .hgt_sim file from Gamma to ROI_PAC
Input:
sim_150911-150922.hgt_sim
sim_150911-150922.rdc.dem
Search for:
sim_150911-150922.diff_par
Output:
sim_150911-150922.hgt_sim.rsc
sim_150911-150922.rdc.dem.rsc
"""
atr = {}
atr['PROCESSOR'] = 'gamma'
atr['FILE_TYPE'] = os.path.splitext(fname)[1]
# Get basename of file
fname_base = os.path.splitext(fname)[0]
for i in range(5):
fname_base = os.path.splitext(fname_base)[0]
par_file = fname_base + '.diff_par'
#print('read '+os.path.basename(par_file))
par_dict = readfile.read_gamma_par(par_file)
atr.update(par_dict)
# Write to .rsc file
rsc_file = fname + '.rsc'
try:
atr_orig = readfile.read_roipac_rsc(rsc_file)
except:
atr_orig = None
#keyList = [i for i in atr_orig.keys() if i in atr.keys()]
if any((i not in atr_orig.keys() or atr_orig[i] != atr[i])
for i in atr.keys()):
print('writing >>> ' + os.path.basename(rsc_file))
writefile.write_roipac_rsc(atr, out_file=rsc_file)
return rsc_file
def extract_metadata4interferogram(fname):
"""Read/extract attributes for PySAR from Gamma .unw, .cor and .int file
Parameters: fname : str, Gamma interferogram filename or path,
i.e. /PopoSLT143TsxD/diff_filt_HDR_130118-130129_4rlks.unw
Returns: atr : dict, Attributes dictionary
"""
file_dir = os.path.dirname(fname)
file_basename = os.path.basename(fname)
rsc_file = fname+'.rsc'
# if os.path.isfile(rsc_file):
# return rsc_file
atr = {}
atr['PROCESSOR'] = 'gamma'
atr['FILE_TYPE'] = os.path.splitext(fname)[1]
# Get info: date12, num of loooks
try:
date12 = str(re.findall('\d{8}[-_]\d{8}', file_basename)[0])
except:
date12 = str(re.findall('\d{6}[-_]\d{6}', file_basename)[0])
m_date, s_date = date12.replace('-', '_').split('_')
atr['DATE12'] = ptime.yymmdd(m_date)+'-'+ptime.yymmdd(s_date)
lks = os.path.splitext(file_basename.split(date12)[1])[0]
# Read .off and .par file
off_file = file_dir+'/*'+date12+lks+'.off'
m_par_file = [file_dir+'/*'+m_date+lks+i for i in ['.amp.par', '.ramp.par']]
s_par_file = [file_dir+'/*'+s_date+lks+i for i in ['.amp.par', '.ramp.par']]
try:
off_file = ut.get_file_list(off_file)[0]
except:
print('\nERROR: Can not find .off file, it supposed to be like: '+off_file)
try:
m_par_file = ut.get_file_list(m_par_file)[0]
except:
print('\nERROR: Can not find master date .par file, it supposed to be like: '+m_par_file)
try:
s_par_file = ut.get_file_list(s_par_file)[0]
except:
print('\nERROR: Can not find slave date .par file, it supposed to be like: '+s_par_file)
par_dict = readfile.read_gamma_par(m_par_file)
off_dict = readfile.read_gamma_par(off_file)
atr.update(par_dict)
atr.update(off_dict)
# Perp Baseline Info
atr = get_perp_baseline(m_par_file, s_par_file, off_file, atr)
# LAT/LON_REF1/2/3/4
atr = get_lalo_ref(m_par_file, atr)
# Write to .rsc file
try:
atr_orig = readfile.read_roipac_rsc(rsc_file)
except:
atr_orig = dict()
if not set(atr.items()).issubset(set(atr_orig.items())):
atr_out = {**atr_orig, **atr}
print('merge %s, %s and %s into %s' % (os.path.basename(m_par_file),
os.path.basename(s_par_file),
os.path.basename(off_file),
os.path.basename(rsc_file)))
writefile.write_roipac_rsc(atr_out, out_file=rsc_file)
return rsc_file
def extract_metadata4geometry_radar(fname):
"""Read/extract attribute for .hgt_sim file from Gamma to ROI_PAC
Input:
sim_20070813_20080310.hgt_sim
sim_20070813_20080310.rdc.dem
Search for:
sim_20070813_20080310.diff_par
Output:
sim_20070813_20080310.hgt_sim.rsc
sim_20070813_20080310.rdc.dem.rsc
"""
# Get/read GAMMA par file
# for loop to get rid of multiple dot in filename
fname_base = os.path.splitext(fname)[0]
for i in range(5):
fname_base = os.path.splitext(fname_base)[0]
par_file = fname_base + '.diff_par'
par_dict = readfile.read_gamma_par(par_file)
# Get/read LAT/LON_REF1/2/3/4
msg = 'grab LAT/LON_REF1/2/3/4 from par file: '
# get date of one acquisition
try:
m_date = str(re.findall('\d{8}', fname_base)[0])
except:
m_date = str(re.findall('\d{6}', fname_base)[0])
# search existing par file
geom_dir = os.path.dirname(fname) #PROJECT_DIR/geom_master
ifg_dir = os.path.join(geom_dir, '../*/{}_*'.format(
m_date)) #PROJECT_DIR/interferograms/{m_date}_20141225
m_par_files = [
os.path.join(geom_dir, '*{}*{}'.format(m_date, ext))
for ext in PAR_EXT_LIST
]
m_par_files += [
os.path.join(ifg_dir, '*{}*{}'.format(m_date, ext))
for ext in PAR_EXT_LIST
]
m_par_files = ut.get_file_list(m_par_files)
# read par file
if len(m_par_files) > 0:
m_par_file = m_par_files[0]
msg += m_par_file
par_dict = get_lalo_ref(m_par_file, par_dict)
else:
msg += ' no par file found with date: {}'.format(m_date)
print(msg)
# initiate ROIPAC dict
atr = {}
atr['PROCESSOR'] = 'gamma'
atr['FILE_TYPE'] = os.path.splitext(fname)[1]
atr.update(par_dict)
# Write to .rsc file
rsc_file = fname + '.rsc'
try:
atr_orig = readfile.read_roipac_rsc(rsc_file)
except:
atr_orig = dict()
if not set(atr.items()).issubset(set(atr_orig.items())):
atr_out = {**atr_orig, **atr}
print('writing >>> ' + os.path.basename(rsc_file))
writefile.write_roipac_rsc(atr_out, out_file=rsc_file)
return rsc_file