def bperp_btemp_master_slave(drsFiles=None):
from dateutil import parser
mresfiles = []
sresfiles = []
iresfiles = []
for f in drsFiles:
d = adore.drs2dict(f)
mresfiles.append(d['general']['m_resfile'].strip())
sresfiles.append(d['general']['s_resfile'].strip())
iresfiles.append(d['general']['i_resfile'].strip())
bperp = []
btemp = []
master = []
slave = []
for f in N.r_[0:len(mresfiles)]:
mobj = adore.dict2obj(adore.res2dict(mresfiles[f]))
sobj = adore.dict2obj(adore.res2dict(sresfiles[f]))
ires = adore.res2dict(iresfiles[f])
iobj = adore.dict2obj(ires)
master.append(parser.parse(mobj.readfiles.First_pixel_azimuth_time))
slave.append(parser.parse(sobj.readfiles.First_pixel_azimuth_time))
bperp.append(iobj.coarse_orbits.Bperp)
btemp.append(iobj.coarse_orbits.Btemp)
return bperp, btemp, master, slave
def bperp_btemp_master_slave(drsFiles=None):
from dateutil import parser
mresfiles=[]
sresfiles=[]
iresfiles=[]
for f in drsFiles:
d=adore.drs2dict(f);
mresfiles.append(d['general']['m_resfile'].strip());
sresfiles.append(d['general']['s_resfile'].strip());
iresfiles.append(d['general']['i_resfile'].strip());
bperp=[]
btemp=[]
master=[]
slave=[]
for f in N.r_[0:len(mresfiles)]:
mobj=adore.dict2obj(adore.res2dict(mresfiles[f]))
sobj=adore.dict2obj(adore.res2dict(sresfiles[f]))
ires=adore.res2dict(iresfiles[f])
iobj=adore.dict2obj(ires)
master.append( parser.parse(mobj.readfiles.First_pixel_azimuth_time) )
slave.append( parser.parse(sobj.readfiles.First_pixel_azimuth_time) )
bperp.append( iobj.coarse_orbits.Bperp )
btemp.append( iobj.coarse_orbits.Btemp )
return bperp, btemp, master, slave
def main(argv=None):
if argv is None:
argv = sys.argv
drsfile = argv[1]
try:
ct = float(argv[2])
except:
ct = 0.0
if drsfile is not None:
d = adore.drs2dict(drsfile)
mresfile = d['general']['m_resfile'].strip()
sresfile = d['general']['s_resfile'].strip()
iresfile = d['general']['i_resfile'].strip()
iobj = adore.dict2obj(adore.res2dict(iresfile))
mobj = adore.dict2obj(adore.res2dict(mresfile))
sobj = adore.dict2obj(adore.res2dict(sresfile))
A = iobj.fine_coreg.results
#ires['fine_coreg']['results'];
A = A[A[:, 5] > ct, :]
#figure
pylab.figure()
pylab.title('Range: Fine Correlation - Coregistration Polynomial')
rg = adore.basic.rescale(
A[:, 2], [-2, 2], arrlim=[mobj.crop.First_pixel, mobj.crop.Last_pixel])
az = adore.basic.rescale(
A[:, 1], [-2, 2], arrlim=[mobj.crop.First_line, mobj.crop.Last_line])
pixelOffset = adore.polyval(rg, az,
iobj.comp_coregpm.Estimated_coefficientsP)
lineOffset = adore.polyval(rg, az,
iobj.comp_coregpm.Estimated_coefficientsL)
Q = pylab.quiver(A[:, 2], A[:, 1], A[:, 4] - pixelOffset,
A[:, 3] - lineOffset, A[:, 5])
azStd = round(pylab.np.std(A[:, 3] - lineOffset))
rgStd = round(pylab.np.std(A[:, 4] - pixelOffset))
qk = pylab.quiverkey(Q,
0.33,
0.92,
azStd, ('%d [px az]' % (azStd)),
labelpos='W',
fontproperties={'weight': 'bold'})
qk = pylab.quiverkey(Q,
0.66,
0.92,
rgStd, ('%d [px rg]' % (rgStd)),
labelpos='W',
fontproperties={'weight': 'bold'})
pylab.colorbar()
pylab.show()
def average_coherence_adore(iresfiles=None, drsFiles=None):
""" avg_coherence=average_coherence_adore(iresfiles=None, drsFiles=None):
"""
from dateutil import parser
if drsFiles is not None:
iresfiles = []
for f in drsFiles:
d = adore.drs2dict(f)
iresfiles.append(d['general']['i_resfile'].strip())
coh = []
for f in N.r_[0:len(iresfiles)]:
ires = adore.res2dict(iresfiles[f])
iobj = adore.dict2obj(ires)
coh.append(adore.getProduct(ires, 'coherence'))
coh = N.rollaxis(N.dstack(coh), 2, 0)
return coh.mean(0)
def average_coherence_adore(iresfiles=None, drsFiles=None):
""" avg_coherence=average_coherence_adore(iresfiles=None, drsFiles=None):
"""
from dateutil import parser
if drsFiles is not None:
iresfiles=[]
for f in drsFiles:
d=adore.drs2dict(f);
iresfiles.append(d['general']['i_resfile'].strip());
coh=[]
for f in N.r_[0:len(iresfiles)]:
ires=adore.res2dict(iresfiles[f])
iobj=adore.dict2obj(ires)
coh.append( adore.getProduct(ires,'coherence') );
coh=N.rollaxis(N.dstack(coh),2,0)
return coh.mean(0);
def bperp_btemp(drsFiles):
from dateutil import parser
mresfiles = []
sresfiles = []
iresfiles = []
for f in drsFiles:
d = adore.drs2dict(f)
mresfiles.append(d['general']['m_resfile'].strip())
sresfiles.append(d['general']['s_resfile'].strip())
iresfiles.append(d['general']['i_resfile'].strip())
bperp = []
btemp = []
for f in N.r_[0:len(mresfiles)]:
mobj = adore.dict2obj(adore.res2dict(mresfiles[f]))
sobj = adore.dict2obj(adore.res2dict(sresfiles[f]))
ires = adore.res2dict(iresfiles[f])
iobj = adore.dict2obj(ires)
bperp.append(iobj.coarse_orbits.Bperp)
btemp.append(iobj.coarse_orbits.Btemp)
return bperp, btemp
def bperp_btemp(drsFiles):
from dateutil import parser
mresfiles=[]
sresfiles=[]
iresfiles=[]
for f in drsFiles:
d=adore.drs2dict(f);
mresfiles.append(d['general']['m_resfile'].strip());
sresfiles.append(d['general']['s_resfile'].strip());
iresfiles.append(d['general']['i_resfile'].strip());
bperp=[]
btemp=[]
for f in N.r_[0:len(mresfiles)]:
mobj=adore.dict2obj(adore.res2dict(mresfiles[f]))
sobj=adore.dict2obj(adore.res2dict(sresfiles[f]))
ires=adore.res2dict(iresfiles[f])
iobj=adore.dict2obj(ires)
bperp.append( iobj.coarse_orbits.Bperp )
btemp.append( iobj.coarse_orbits.Btemp )
return bperp, btemp