def geocode_file(inps):
"""
Geocodes the input file
:param inps: input name space
"""
if 'stripmap' in inps.prefix:
sys.path.append(os.path.join(os.getenv('ISCE_STACK'), 'stripmapStack'))
else:
sys.path.append(os.path.join(os.getenv('ISCE_STACK'), 'topsStack'))
import geocodeGdal as gg
inps.cropbox = [val for val in inps.cropbox.split()]
if len(inps.cropbox) != 4:
raise Exception('Bbox should contain 4 floating point values')
scp_arg = {'latFile': os.path.abspath(os.path.join(inps.geo_reference_dir, inps.lat_file)),
'lonFile': os.path.abspath(os.path.join(inps.geo_reference_dir, inps.lon_file)),
'xOff': 0,
'yOff': 0,
'prodlist': [inps.input_file + '.slc.ml'],
'resamplingMethod': 'near'}
scp_arg = argparse.Namespace(**scp_arg)
#scp_arg.prodlist = [inps.prod_list + '.slc.ml']
WSEN = str(inps.cropbox[2]) + ' ' + str(inps.cropbox[0]) + ' ' + str(inps.cropbox[3]) + ' ' + str(inps.cropbox[1])
latFile, lonFile = gg.prepare_lat_lon(scp_arg)
gg.getBound(latFile, float(inps.cropbox[0]), float(inps.cropbox[1]), 'lat')
gg.getBound(lonFile, float(inps.cropbox[2]), float(inps.cropbox[3]), 'lon')
infile = os.path.abspath(inps.input_file + '.slc.ml')
print('geocoding ' + infile)
outFile = os.path.join(os.path.dirname(infile), "geo_" + os.path.basename(infile))
gg.writeVRT(infile, latFile, lonFile)
cmd = 'gdalwarp -of ENVI -geoloc -te ' + WSEN + ' -tr ' + str(inps.lat_step) + ' ' + \
str(inps.lon_step) + ' -srcnodata 0 -dstnodata 0 ' + ' -r ' + inps.resampling_method + \
' -co INTERLEAVE=BIL ' + infile + '.vrt ' + outFile
print(cmd)
os.system(cmd)
# write_xml(outFile)
cmd = "gdal2isce_xml.py -i " + outFile
os.system(cmd)
return
def geocode_file(inps):
"""
Geocodes the input file
:param inps: input name space
"""
import geocodeGdal as gg
inps.cropbox = [val for val in inps.cropbox.split()]
if len(inps.cropbox) != 4:
raise Exception('Bbox should contain 4 floating point values')
inps.lat_file = os.path.abspath(
os.path.join(inps.geo_master_dir, inps.lat_file))
inps.lon_file = os.path.abspath(
os.path.join(inps.geo_master_dir, inps.lon_file))
inps.prod_list = [inps.prod_list + '.slc.ml']
WSEN = str(inps.cropbox[2]) + ' ' + str(inps.cropbox[0]) + ' ' + str(
inps.cropbox[3]) + ' ' + str(inps.cropbox[1])
latFile, lonFile = gg.prepare_lat_lon(inps)
gg.getBound(latFile, float(inps.cropbox[0]), float(inps.cropbox[1]), 'lat')
gg.getBound(lonFile, float(inps.cropbox[2]), float(inps.cropbox[3]), 'lon')
for infile in inps.prod_list:
infile = os.path.abspath(infile)
print('geocoding ' + infile)
outFile = os.path.join(os.path.dirname(infile),
"geo_" + os.path.basename(infile))
gg.writeVRT(infile, latFile, lonFile)
cmd = 'gdalwarp -of ENVI -geoloc -te ' + WSEN + ' -tr ' + str(inps.lat_step) + ' ' + \
str(inps.lon_step) + ' -srcnodata 0 -dstnodata 0 ' + ' -r ' + inps.resampling_method + \
' -co INTERLEAVE=BIL ' + infile + '.vrt ' + outFile
print(cmd)
os.system(cmd)
# write_xml(outFile)
cmd = "gdal2isce_xml.py -i " + outFile
os.system(cmd)
return