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