def read_geometry(inps):
ts_obj = timeseries(inps.timeseries_file)
ts_obj.open(print_msg=False)
# 2D / 3D geometry
if inps.geom_file:
geom_obj = geometry(inps.geom_file)
geom_obj.open()
print(('read 2D incidenceAngle,slantRangeDistance from {} file:'
' {}').format(geom_obj.name, os.path.basename(geom_obj.file)))
inps.incAngle = geom_obj.read(datasetName='incidenceAngle', print_msg=False).flatten()
inps.rangeDist = geom_obj.read(datasetName='slantRangeDistance', print_msg=False).flatten()
if 'bperp' in geom_obj.datasetNames:
print('read 3D bperp from {} file: {} ...'.format(geom_obj.name, os.path.basename(geom_obj.file)))
inps.pbase = geom_obj.read(datasetName='bperp', print_msg=False).reshape((geom_obj.numDate, -1))
inps.pbase -= inps.pbase[ts_obj.refIndex]
else:
print('read mean bperp from {} file'.format(ts_obj.name))
inps.pbase = ts_obj.pbase.reshape((-1, 1))
# 0D geometry
else:
print('read mean incidenceAngle,slantRangeDistance,bperp value from {} file'.format(ts_obj.name))
inps.incAngle = ut.incidence_angle(ts_obj.metadata, dimension=0)
inps.rangeDist = ut.range_distance(ts_obj.metadata, dimension=0)
inps.pbase = ts_obj.pbase.reshape((-1, 1))
inps.sinIncAngle = np.sin(inps.incAngle * np.pi / 180.)
return inps
def get_slant_range_distance(self, box=None):
if not self.extraMetadata or 'Y_FIRST' in self.extraMetadata.keys():
return None
data = ut.range_distance(self.extraMetadata,
dimension=2,
print_msg=False)
if box is not None:
data = data[box[1]:box[3], box[0]:box[2]]
return data
def main(argv):
try:
File = argv[0]
atr = readfile.read_attribute(File)
except:
usage()
sys.exit(1)
try:
outFile = argv[1]
except:
outFile = 'rangeDistance.h5'
# Calculate look angle
range_dis = ut.range_distance(atr, dimension=2)
# Geo coord
if 'Y_FIRST' in atr.keys():
print(
'Input file is geocoded, only center range distance is calculated: '
)
print(range_dis)
length = int(atr['LENGTH'])
width = int(atr['WIDTH'])
range_dis_mat = np.zeros((length, width), np.float32)
range_dis_mat[:] = range_dis
range_dis = range_dis_mat
print('writing >>> ' + outFile)
atr['FILE_TYPE'] = 'mask'
atr['UNIT'] = 'm'
try:
atr.pop('REF_DATE')
except:
pass
writefile.write(range_dis, out_file=outFile, metadata=atr)
return outFile