def extract_multilook_number( geom_dir, meta=dict(), fext_list=['.rdr', '.geo', '.rdr.full', '.geo.full']): for fbase in ['hgt', 'lat', 'lon', 'los']: fbase = os.path.join(geom_dir, fbase) for fext in fext_list: fnames = glob.glob(fbase + fext) if len(fnames) > 0: break if len(fnames) > 0: fullXmlFile = '{}.full.xml'.format(fnames[0]) if os.path.isfile(fullXmlFile): fullXmlDict = readfile.read_isce_xml(fullXmlFile) xmlDict = readfile.read_attribute(fnames[0]) meta['ALOOKS'] = int( int(fullXmlDict['LENGTH']) / int(xmlDict['LENGTH'])) meta['RLOOKS'] = int( int(fullXmlDict['WIDTH']) / int(xmlDict['WIDTH'])) break # default value for key in ['ALOOKS', 'RLOOKS']: if key not in meta: meta[key] = 1 # NCORRLOOKS for coherence calibration rgfact = float(meta['rangeResolution']) / float(meta['rangePixelSize']) azfact = float(meta['azimuthResolution']) / float(meta['azimuthPixelSize']) meta['NCORRLOOKS'] = meta['RLOOKS'] * meta['ALOOKS'] / (rgfact * azfact) return meta
def extract_multilook_number(geom_dir, metadata=dict()): for fbase in ['hgt','lat','lon','los']: fbase = os.path.join(geom_dir, fbase) fnames = glob.glob('{}*.rdr'.format(fbase)) + glob.glob('{}*.geo'.format(fbase)) if len(fnames) > 0: fullXmlFile = '{}.full.xml'.format(fnames[0]) if os.path.isfile(fullXmlFile): fullXmlDict = readfile.read_isce_xml(fullXmlFile) xmlDict = readfile.read_attribute(fnames[0]) metadata['ALOOKS'] = int(int(fullXmlDict['LENGTH']) / int(xmlDict['LENGTH'])) metadata['RLOOKS'] = int(int(fullXmlDict['WIDTH']) / int(xmlDict['WIDTH'])) break # default for key in ['ALOOKS', 'RLOOKS']: if key not in metadata: metadata[key] = 1 return metadata
def extract_multilook_number(geom_dir, metadata=dict()): for fbase in ['hgt','lat','lon','los']: fbase = os.path.join(geom_dir, fbase) fnames = glob.glob('{}*.rdr'.format(fbase)) + glob.glob('{}*.geo'.format(fbase)) if len(fnames) > 0: fullXmlFile = '{}.full.xml'.format(fnames[0]) if os.path.isfile(fullXmlFile): fullXmlDict = readfile.read_isce_xml(fullXmlFile) xmlDict = readfile.read_attribute(fnames[0]) metadata['ALOOKS'] = int(int(fullXmlDict['LENGTH']) / int(xmlDict['LENGTH'])) metadata['RLOOKS'] = int(int(fullXmlDict['WIDTH']) / int(xmlDict['WIDTH'])) break # default for key in ['ALOOKS', 'RLOOKS']: if key not in metadata: metadata[key] = 1 return metadata