def main(argv):
try:
file = argv[0]
alks = float(argv[1])
rlks = float(argv[2])
except:
Usage()
sys.exit(1)
ext = os.path.splitext(file)[1]
outName = file.split('.')[0] + '_a' + str(int(alks)) + 'lks_r' + str(
int(rlks)) + 'lks' + ext
if ext == '.int' or ext == '.slc':
a, p, r = readfile.read_complex64(file)
plks = multilook(p, alks, rlks)
alks = multilook(a, alks, rlks)
r['FILE_LENGTH'] = str(dlks.shape[0])
r['WIDTH'] = str(dlks.shape[1])
r['XMAX'] = str(int(r['WIDTH']) - 1)
r['YMAX'] = str(int(r['FILE_LENGTH']) - 1)
try:
r['Y_STEP'] = str(float(r['Y_STEP']) * alks)
r['X_STEP'] = str(float(r['X_STEP']) * rlks)
except:
Geo = 0
f = open(outName + '.rsc', 'w')
for k in r.keys():
f.write(k + ' ' + r[k] + '\n')
f.close()
elif ext == '.unw' or ext == '.cor' or ext == '.hgt':
a, p, r = readfile.read_float32(file)
plks = multilook(p, alks, rlks)
alks = multilook(a, alks, rlks)
writefile.write_float32(plks, outName)
r['FILE_LENGTH'] = str(dlks.shape[0])
r['WIDTH'] = str(dlks.shape[1])
r['XMAX'] = str(int(r['WIDTH']) - 1)
r['YMAX'] = str(int(r['FILE_LENGTH']) - 1)
try:
r['Y_STEP'] = str(float(r['Y_STEP']) * alks)
r['X_STEP'] = str(float(r['X_STEP']) * rlks)
except:
Geo = 0
f = open(outName + '.rsc', 'w')
for k in r.keys():
f.write(k + ' ' + r[k] + '\n')
f.close()
elif ext == ('.dem'):
d, r = readfile.read_dem(file)
dlks = multilook(d, alks, rlks)
print 'writing ' + outName
writefile.write_dem(dlks, outName)
r['FILE_LENGTH'] = str(dlks.shape[0])
r['WIDTH'] = str(dlks.shape[1])
r['XMAX'] = str(int(r['WIDTH']) - 1)
r['YMAX'] = str(int(r['FILE_LENGTH']) - 1)
try:
r['Y_STEP'] = str(float(r['Y_STEP']) * alks)
r['X_STEP'] = str(float(r['X_STEP']) * rlks)
except:
Geo = 0
f = open(outName + '.rsc', 'w')
for k in r.keys():
f.write(k + ' ' + r[k] + '\n')
f.close()
elif ext in ['.jpeg', 'jpg', 'png']:
import Image
im = Image.open(file)
width = im.size[0] / int(rlks)
height = im.size[1] / int(alks)
imlks = im.resize((width, height), Image.NEAREST)
print 'writing ' + outName
imlks.save(outName)
try:
r = readfile.read_rsc_file(file + '.rsc')
except:
sys.exit(1)
r['FILE_LENGTH'] = str(height)
r['WIDTH'] = str(width)
r['XMAX'] = str(int(r['WIDTH']) - 1)
r['YMAX'] = str(int(r['FILE_LENGTH']) - 1)
try:
r['Y_STEP'] = str(float(r['Y_STEP']) * alks)
r['X_STEP'] = str(float(r['X_STEP']) * rlks)
except:
Geo = 0
f = open(outName + '.rsc', 'w')
for k in r.keys():
f.write(k + ' ' + r[k] + '\n')
f.close()
elif ext == ('.h5'):
h5file = h5py.File(file, 'r')
# outName=file.split('.')[0]+'_a'+str(int(alks))+'lks_r'+str(int(rlks))+'lks.h5'
h5file_lks = h5py.File(outName, 'w')
if 'interferograms' in h5file.keys():
print 'Multilooking the interferograms'
gg = h5file_lks.create_group('interferograms')
igramList = h5file['interferograms'].keys()
for igram in igramList:
print igram
unw = h5file['interferograms'][igram].get(igram)
unwlks = multilook(unw, alks, rlks)
group = gg.create_group(igram)
dset = group.create_dataset(igram,
data=unwlks,
compression='gzip')
for key, value in h5file['interferograms'][
igram].attrs.iteritems():
group.attrs[key] = value
group.attrs['WIDTH'] = unwlks.shape[1]
group.attrs['FILE_LENGTH'] = unwlks.shape[0]
try:
group.attrs['Y_STEP'] = alks * float(group.attrs['Y_STEP'])
group.attrs['X_STEP'] = rlks * float(group.attrs['X_STEP'])
except:
group.attrs['AZIMUTH_PIXEL_SIZE'] = alks * float(
group.attrs['AZIMUTH_PIXEL_SIZE'])
group.attrs['RANGE_PIXEL_SIZE'] = rlks * float(
group.attrs['RANGE_PIXEL_SIZE'])
dset1 = h5file['mask'].get('mask')
mask = dset1[0:dset1.shape[0], 0:dset1.shape[1]]
masklks = multilook(mask, alks, rlks)
group = h5file_lks.create_group('mask')
dset = group.create_dataset('mask',
data=masklks,
compression='gzip')
elif 'timeseries' in h5file.keys():
print 'Multilooking the time-series'
group = h5file_lks.create_group('timeseries')
dateList = h5file['timeseries'].keys()
for d in dateList:
print d
unw = h5file['timeseries'].get(d)
unwlks = multilook(unw, alks, rlks)
dset = group.create_dataset(d, data=unwlks, compression='gzip')
for key, value in h5file['timeseries'].attrs.iteritems():
group.attrs[key] = value
group.attrs['WIDTH'] = unwlks.shape[1]
group.attrs['FILE_LENGTH'] = unwlks.shape[0]
try:
group.attrs['Y_STEP'] = alks * float(group.attrs['Y_STEP'])
group.attrs['X_STEP'] = rlks * float(group.attrs['X_STEP'])
except:
group.attrs['AZIMUTH_PIXEL_SIZE'] = alks * float(
group.attrs['AZIMUTH_PIXEL_SIZE'])
group.attrs['RANGE_PIXEL_SIZE'] = rlks * float(
group.attrs['RANGE_PIXEL_SIZE'])
try:
dset1 = h5file['mask'].get('mask')
Mask = dset1[0:dset1.shape[0], 0:dset1.shape[1]]
Masklks = multilook(Mask, alks, rlks)
group = h5file_lks.create_group('mask')
dset = group.create_dataset('mask',
data=Masklks,
compression='gzip')
except:
print 'Multilooked file does not include the maske'
elif 'temporal_coherence' in h5file.keys(
) or 'velocity' in h5file.keys() or 'mask' in h5file.keys():
k = h5file.keys()
print 'multi looking the ' + k[0]
group = h5file_lks.create_group(k[0])
dset1 = h5file[k[0]].get(k[0])
Mask = dset1[0:dset1.shape[0], 0:dset1.shape[1]]
Masklks = multilook(Mask, alks, rlks)
dset = group.create_dataset(k[0], data=Masklks, compression='gzip')
for key, value in h5file[k[0]].attrs.iteritems():
group.attrs[key] = value
try:
group.attrs['Y_STEP'] = alks * float(group.attrs['Y_STEP'])
group.attrs['X_STEP'] = rlks * float(group.attrs['X_STEP'])
except:
group.attrs['AZIMUTH_PIXEL_SIZE'] = alks * float(
group.attrs['AZIMUTH_PIXEL_SIZE'])
group.attrs['RANGE_PIXEL_SIZE'] = rlks * float(
group.attrs['RANGE_PIXEL_SIZE'])
group.attrs['WIDTH'] = Masklks.shape[1]
group.attrs['FILE_LENGTH'] = Masklks.shape[0]
h5file.close()
h5file_lks.close()
def main(argv):
try:
opts, args = getopt.getopt(argv, "h:f:t:p:")
except getopt.GetoptError:
Usage()
sys.exit(1)
if opts == []:
Usage()
sys.exit(1)
for opt, arg in opts:
if opt in ("-h", "--help"):
Usage()
sys.exit()
elif opt == "-f":
file = arg
elif opt == "-t":
filtType = arg
elif opt == "-p":
par = arg
# try:
# file=argv[0]
# alks=float(argv[1])
# rlks=float(argv[2])
# except:
# Usage();sys.exit(1)
ext = os.path.splitext(file)[1]
outName = file.split(".")[0] + "_" + filtType + ext
try:
par
except:
par = []
print "+++++++++++++++++++++++++++"
print "Filter type : " + filtType
print "parameters : " + str(par)
print "+++++++++++++++++++++++++++"
###############################################
if ext == ".int" or ext == ".slc":
a, p, r = readfile.read_complex64(file)
plks = multilook(p, alks, rlks)
alks = multilook(a, alks, rlks)
r["FILE_LENGTH"] = str(dlks.shape[0])
r["WIDTH"] = str(dlks.shape[1])
r["XMAX"] = str(int(r["WIDTH"]) - 1)
r["YMAX"] = str(int(r["FILE_LENGTH"]) - 1)
try:
r["Y_STEP"] = str(float(r["Y_STEP"]) * alks)
r["X_STEP"] = str(float(r["X_STEP"]) * rlks)
except:
Geo = 0
f = open(outName + ".rsc", "w")
for k in r.keys():
f.write(k + " " + r[k] + "\n")
f.close()
elif ext == ".unw" or ext == ".cor" or ext == ".hgt":
a, p, r = readfile.read_float32(file)
plks = multilook(p, alks, rlks)
alks = multilook(a, alks, rlks)
writefile.write_float32(plks, outName)
r["FILE_LENGTH"] = str(dlks.shape[0])
r["WIDTH"] = str(dlks.shape[1])
r["XMAX"] = str(int(r["WIDTH"]) - 1)
r["YMAX"] = str(int(r["FILE_LENGTH"]) - 1)
try:
r["Y_STEP"] = str(float(r["Y_STEP"]) * alks)
r["X_STEP"] = str(float(r["X_STEP"]) * rlks)
except:
Geo = 0
f = open(outName + ".rsc", "w")
for k in r.keys():
f.write(k + " " + r[k] + "\n")
f.close()
elif ext == (".dem"):
d, r = readfile.read_dem(file)
dlks = multilook(d, alks, rlks)
print "writing " + outName
writefile.write_dem(dlks, outName)
r["FILE_LENGTH"] = str(dlks.shape[0])
r["WIDTH"] = str(dlks.shape[1])
r["XMAX"] = str(int(r["WIDTH"]) - 1)
r["YMAX"] = str(int(r["FILE_LENGTH"]) - 1)
try:
r["Y_STEP"] = str(float(r["Y_STEP"]) * alks)
r["X_STEP"] = str(float(r["X_STEP"]) * rlks)
except:
Geo = 0
f = open(outName + ".rsc", "w")
for k in r.keys():
f.write(k + " " + r[k] + "\n")
f.close()
elif ext in [".jpeg", "jpg", "png"]:
import Image
im = Image.open(file)
width = im.size[0] / int(rlks)
height = im.size[1] / int(alks)
imlks = im.resize((width, height), Image.NEAREST)
print "writing " + outName
imlks.save(outName)
try:
r = readfile.read_rsc_file(file + ".rsc")
except:
sys.exit(1)
r["FILE_LENGTH"] = str(height)
r["WIDTH"] = str(width)
r["XMAX"] = str(int(r["WIDTH"]) - 1)
r["YMAX"] = str(int(r["FILE_LENGTH"]) - 1)
try:
r["Y_STEP"] = str(float(r["Y_STEP"]) * alks)
r["X_STEP"] = str(float(r["X_STEP"]) * rlks)
except:
Geo = 0
f = open(outName + ".rsc", "w")
for k in r.keys():
f.write(k + " " + r[k] + "\n")
f.close()
elif ext == (".h5"):
h5file = h5py.File(file, "r")
# outName=file.split('.')[0]+'_a'+str(int(alks))+'lks_r'+str(int(rlks))+'lks.h5'
h5file_lks = h5py.File(outName, "w")
if "interferograms" in h5file.keys():
print "Filtering the interferograms in space"
gg = h5file_lks.create_group("interferograms")
igramList = h5file["interferograms"].keys()
for igram in igramList:
print igram
unwSet = h5file["interferograms"][igram].get(igram)
unw = unwSet[0 : unwSet.shape[0], 0 : unwSet.shape[1]]
unw = filter(unw, filtType, par)
group = gg.create_group(igram)
dset = group.create_dataset(igram, data=unw, compression="gzip")
for key, value in h5file["interferograms"][igram].attrs.iteritems():
group.attrs[key] = value
dset1 = h5file["mask"].get("mask")
mask = dset1[0 : dset1.shape[0], 0 : dset1.shape[1]]
group = h5file_lks.create_group("mask")
dset = group.create_dataset("mask", data=mask, compression="gzip")
elif "timeseries" in h5file.keys():
print "Filtering the time-series"
group = h5file_lks.create_group("timeseries")
dateList = h5file["timeseries"].keys()
for d in dateList:
print d
dset1 = h5file["timeseries"].get(d)
data = dset1[0 : dset1.shape[0], 0 : dset1.shape[1]]
data = filter(data, filtType, par)
dset = group.create_dataset(d, data=data, compression="gzip")
for key, value in h5file["timeseries"].attrs.iteritems():
group.attrs[key] = value
try:
dset1 = h5file["mask"].get("mask")
Mask = dset1[0 : dset1.shape[0], 0 : dset1.shape[1]]
# Masklks=multilook(Mask,alks,rlks)
group = h5file_lks.create_group("mask")
dset = group.create_dataset("mask", data=Mask, compression="gzip")
except:
print "Filterd file does not include the maske"
elif "temporal_coherence" in h5file.keys() or "velocity" in h5file.keys() or "mask" in h5file.keys():
k = h5file.keys()
print "filtering the " + k[0]
group = h5file_lks.create_group(k[0])
dset1 = h5file[k[0]].get(k[0])
data = dset1[0 : dset1.shape[0], 0 : dset1.shape[1]]
data = filter(data, filtType, par)
dset = group.create_dataset(k[0], data=data, compression="gzip")
for key, value in h5file[k[0]].attrs.iteritems():
group.attrs[key] = value
h5file.close()
h5file_lks.close()
def main(argv):
#outName='subsetIgrams.h5'
try:
opts, args = getopt.getopt(argv,"h:f:x:y:o:l:L:")
except getopt.GetoptError:
print 'Error while getting args'
Usage() ; sys.exit(1)
for opt,arg in opts:
if opt in ("-h","--help"):
Usage()
sys.exit()
elif opt == '-f':
File = arg
elif opt=='-y':
ysub=[int(i) for i in arg.split(':')]
ysub.sort()
elif opt=='-x':
xsub = [int(i) for i in arg.split(':')]
xsub.sort()
elif opt=='-o':
outName=arg
elif opt=='-l':
Latsub=[float(i) for i in arg.split(':')]
Latsub.sort()
elif opt=='-L':
Lonsub = [float(i) for i in arg.split(':')]
Lonsub.sort()
#####################################################
try:
File
xsub
ysub
except:
try:
File
Latsub
Lonsub
except:
Usage();sys.exit(1)
try:
outName
except:
outName='subset_'+File
ext = os.path.splitext(File)[1]
if ext == '.h5':
try:
h5file=h5py.File(File,'r')
except:
Usage() ; sys.exit(1)
k=h5file.keys()
# convert LatLon to xy for geocoded file
try:
Latsub
Lonsub
if 'X_FIRST' in h5file[k[0]].attrs.keys():
xsub=[0]*2
ysub=[0]*2
xsub[0]=int((Lonsub[0]-float(h5file[k[0]].attrs['X_FIRST']))/float(h5file[k[0]].attrs['X_STEP']))
xsub[1]=int((Lonsub[1]-float(h5file[k[0]].attrs['X_FIRST']))/float(h5file[k[0]].attrs['X_STEP']))
ysub[0]=int((Latsub[1]-float(h5file[k[0]].attrs['Y_FIRST']))/float(h5file[k[0]].attrs['Y_STEP']))
ysub[1]=int((Latsub[0]-float(h5file[k[0]].attrs['Y_FIRST']))/float(h5file[k[0]].attrs['Y_STEP']))
print 'Subseting geocoded',ext,' file with Latitude and Longitude...'
elif 'X_FIRST' in h5file[k[0]][h5file[k[0]].keys()[0]].attrs.keys(): # for geocoded interferograms/coherence
igramList=h5file[k[0]].keys()
xsub=[0]*2
ysub=[0]*2
xsub[0]=int((Lonsub[0]-float(h5file[k[0]][igramList[0]].attrs['X_FIRST']))/float(h5file[k[0]][igramList[0]].attrs['X_STEP']))
xsub[1]=int((Lonsub[1]-float(h5file[k[0]][igramList[0]].attrs['X_FIRST']))/float(h5file[k[0]][igramList[0]].attrs['X_STEP']))
ysub[0]=int((Latsub[1]-float(h5file[k[0]][igramList[0]].attrs['Y_FIRST']))/float(h5file[k[0]][igramList[0]].attrs['Y_STEP']))
ysub[1]=int((Latsub[0]-float(h5file[k[0]][igramList[0]].attrs['Y_FIRST']))/float(h5file[k[0]][igramList[0]].attrs['Y_STEP']))
print 'Subseting geocoded',ext,' file with Latitude and Longitude...'
else:
print 'Not geocoded file, cannot be subseted with LatLon.'
Usage() ; sys.exit(1)
except:
Geo=0
# k=h5file.keys()
if 'interferograms' in k:
igramList=h5file['interferograms'].keys()
h5out=h5py.File(outName,'w')
gg=h5out.create_group('interferograms')
for igram in igramList:
print igram
dset1=h5file['interferograms'][igram].get(igram)
group=gg.create_group(igram)
dset=group.create_dataset(igram, data=dset1[ysub[0]:ysub[1],xsub[0]:xsub[1]], compression='gzip')
for key, value in h5file['interferograms'][igram].attrs.iteritems():
group.attrs[key] = value
group.attrs['FILE_LENGTH']=shape(dset1[ysub[0]:ysub[1],xsub[0]:xsub[1]])[0]
group.attrs['WIDTH']=shape(dset1[ysub[0]:ysub[1],xsub[0]:xsub[1]])[1]
group.attrs['subset_x0']=xsub[0]
group.attrs['subset_x1']=xsub[1]
group.attrs['subset_y0']=ysub[0]
group.attrs['subset_y1']=ysub[1]
if 'X_FIRST' in h5file['interferograms'][igram].attrs.keys():
group.attrs['X_FIRST']=float(h5file['interferograms'][igram].attrs['X_FIRST']) + xsub[0]*float(h5file['interferograms'][igram].attrs['X_STEP'])
group.attrs['Y_FIRST']=float(h5file['interferograms'][igram].attrs['Y_FIRST']) + ysub[0]*float(h5file['interferograms'][igram].attrs['Y_STEP'])
gm=h5out.create_group('mask')
try:
Mset=h5file['mask'].get('mask')
dset=gm.create_dataset('mask', data=Mset[ysub[0]:ysub[1],xsub[0]:xsub[1]], compression='gzip')
except:
print 'No group for mask found! It may cause problem in other processing steps.'
try:
Cset=h5file['meanCoherence'].get('meanCoherence')
gm=h5out.create_group('meanCoherence')
dset=gm.create_dataset('meanCoherence', data=Cset[ysub[0]:ysub[1],xsub[0]:xsub[1]], compression='gzip')
except:
print 'No average coherence found in the File'
elif k[0] in ('coherence','wrapped'):
corList=h5file[k[0]].keys()
h5out=h5py.File(outName,'w')
gg=h5out.create_group(k[0])
for cor in corList:
print cor
dset1=h5file[k[0]][cor].get(cor)
group=gg.create_group(cor)
dset=group.create_dataset(cor, data=dset1[ysub[0]:ysub[1],xsub[0]:xsub[1]], compression='gzip')
for key, value in h5file[k[0]][cor].attrs.iteritems():
group.attrs[key] = value
group.attrs['FILE_LENGTH']=shape(dset1[ysub[0]:ysub[1],xsub[0]:xsub[1]])[0]
group.attrs['WIDTH']=shape(dset1[ysub[0]:ysub[1],xsub[0]:xsub[1]])[1]
group.attrs['subset_x0']=xsub[0]
group.attrs['subset_x1']=xsub[1]
group.attrs['subset_y0']=ysub[0]
group.attrs['subset_y1']=ysub[1]
if 'X_FIRST' in h5file[k[0]][cor].attrs.keys():
group.attrs['X_FIRST']=float(h5file[k[0]][cor].attrs['X_FIRST']) + xsub[0]*float(h5file[k[0]][cor].attrs['X_STEP'])
group.attrs['Y_FIRST']=float(h5file[k[0]][cor].attrs['Y_FIRST']) + ysub[0]*float(h5file[k[0]][cor].attrs['Y_STEP'])
elif 'timeseries' in h5file.keys():
dateList=h5file['timeseries'].keys()
h5out=h5py.File(outName,'w')
group=h5out.create_group('timeseries')
for d in dateList:
print d
dset1=h5file['timeseries'].get(d)
dset=group.create_dataset(d, data=dset1[ysub[0]:ysub[1],xsub[0]:xsub[1]], compression='gzip')
for key, value in h5file['timeseries'].attrs.iteritems():
group.attrs[key] = value
group.attrs['FILE_LENGTH']=shape(dset1[ysub[0]:ysub[1],xsub[0]:xsub[1]])[0]
group.attrs['WIDTH']=shape(dset1[ysub[0]:ysub[1],xsub[0]:xsub[1]])[1]
group.attrs['subset_x0']=xsub[0]
group.attrs['subset_x1']=xsub[1]
group.attrs['subset_y0']=ysub[0]
group.attrs['subset_y1']=ysub[1]
if 'X_FIRST' in h5file['timeseries'].attrs.keys():
group.attrs['X_FIRST']=float(h5file['timeseries'].attrs['X_FIRST']) + xsub[0]*float(h5file['timeseries'].attrs['X_STEP'])
group.attrs['Y_FIRST']=float(h5file['timeseries'].attrs['Y_FIRST']) + ysub[0]*float(h5file['timeseries'].attrs['Y_STEP'])
h5file.close()
h5out.close()
elif 'temporal_coherence' in h5file.keys() or 'velocity' in h5file.keys() or 'mask' in h5file.keys() or 'rmse' in h5file.keys():
print 'writing >>> ' +outName
dset=h5file[k[0]].get(k[0])
data=dset[ysub[0]:ysub[1],xsub[0]:xsub[1]]
hfout=h5py.File(outName,'w')
group= hfout.create_group(k[0])
group.create_dataset(k[0],data=data,compression='gzip')
for key,value in h5file[k[0]].attrs.iteritems():
group.attrs[key]=value
group.attrs['FILE_LENGTH']=data.shape[0]
group.attrs['WIDTH']=data.shape[1]
group.attrs['XMIN']=0
group.attrs['XMAX']=data.shape[1]-1
group.attrs['YMIN']=0
group.attrs['YMAX']=data.shape[0]-1
group.attrs['subset_x0']=xsub[0]
group.attrs['subset_x1']=xsub[1]
group.attrs['subset_y0']=ysub[0]
group.attrs['subset_y1']=ysub[1]
if 'X_FIRST' in h5file[k[0]].attrs.keys():
group.attrs['X_FIRST']=float(h5file[k[0]].attrs['X_FIRST']) + xsub[0]*float(h5file[k[0]].attrs['X_STEP'])
group.attrs['Y_FIRST']=float(h5file[k[0]].attrs['Y_FIRST']) + ysub[0]*float(h5file[k[0]].attrs['Y_STEP'])
h5file.close()
hfout.close()
elif ext in ['.unw','.cor','.hgt']:
a,p,r = readfile.read_float32(File)
try:
Latsub
Lonsub
try:
r['X_FIRST']
xsub=[0]*2
ysub=[0]*2
xsub[0]=int((Lonsub[0]-float(r['X_FIRST']))/float(r['X_STEP']))
xsub[1]=int((Lonsub[1]-float(r['X_FIRST']))/float(r['X_STEP']))
ysub[0]=int((Latsub[1]-float(r['Y_FIRST']))/float(r['Y_STEP']))
ysub[1]=int((Latsub[0]-float(r['Y_FIRST']))/float(r['Y_STEP']))
print 'Subseting geocoded',ext,' file with Latitude and Longitude...'
except:
print 'Not geocoded file, cannot be subseted with LatLon.'
Usage() ; sys.exit(1)
except:
Geo=0
a=a[ysub[0]:ysub[1],xsub[0]:xsub[1]]
p=p[ysub[0]:ysub[1],xsub[0]:xsub[1]]
print 'writing >>> '+outName
writefile.write_float32(p,outName)
r['FILE_LENGTH']=str(p.shape[0])
r['WIDTH']=str(p.shape[1])
r['XMAX']=str(int(r['WIDTH']) - 1)
r['YMAX']=str(int(r['FILE_LENGTH']) - 1)
r['subset_x0']=str(xsub[0])
r['subset_x1']=str(xsub[1])
r['subset_y0']=str(ysub[0])
r['subset_y1']=str(ysub[1])
try:
r['Y_FIRST']=str(float(r['Y_FIRST'])+ysub[0]*float(r['Y_STEP']))
r['X_FIRST']=str(float(r['X_FIRST'])+xsub[0]*float(r['X_STEP']))
except:
Geo=0
f = open(outName+'.rsc','w')
for k in r.keys():
f.write(k+' '+r[k]+'\n')
f.close()
elif ext== '.dem':
d,r = readfile.read_dem(File)
try:
Latsub
Lonsub
# print Latsub
try:
r['X_FIRST']
xsub=[0]*2
ysub=[0]*2
xsub[0]=int((Lonsub[0]-float(r['X_FIRST']))/float(r['X_STEP']))
xsub[1]=int((Lonsub[1]-float(r['X_FIRST']))/float(r['X_STEP']))
ysub[0]=int((Latsub[1]-float(r['Y_FIRST']))/float(r['Y_STEP']))
ysub[1]=int((Latsub[0]-float(r['Y_FIRST']))/float(r['Y_STEP']))
print 'Subseting',ext,' file with Latitude and Longitude...'
except:
print 'Not geocoded file, cannot be subseted with LatLon.'
Usage() ; sys.exit(1)
except:
Geo=0
d=d[ysub[0]:ysub[1],xsub[0]:xsub[1]]
print 'writing >>> '+outName
writefile.write_dem(d,outName)
r['FILE_LENGTH']=str(d.shape[0])
r['WIDTH']=str(d.shape[1])
r['XMAX']=str(int(r['WIDTH']) - 1)
r['YMAX']=str(int(r['FILE_LENGTH']) - 1)
r['subset_x0']=str(xsub[0])
r['subset_x1']=str(xsub[1])
r['subset_y0']=str(ysub[0])
r['subset_y1']=str(ysub[1])
try:
r['Y_FIRST']=str(float(r['Y_FIRST'])+ysub[0]*float(r['Y_STEP']))
r['X_FIRST']=str(float(r['X_FIRST'])+xsub[0]*float(r['X_STEP']))
except:
Geo=0
f = open(outName+'.rsc','w')
for k in r.keys():
f.write(k+' '+r[k]+'\n')
f.close()
elif ext in ['.jpeg','jpg','png']:
import Image
im = Image.open(File)
try:
r=readfile.read_rsc_file(File+'.rsc')
except:
sys.exit(1)
try:
Latsub
Lonsub
try:
r['X_FIRST']
xsub=[0]*2
ysub=[0]*2
xsub[0]=int((Lonsub[0]-float(r['X_FIRST']))/float(r['X_STEP']))
xsub[1]=int((Lonsub[1]-float(r['X_FIRST']))/float(r['X_STEP']))
ysub[0]=int((Latsub[1]-float(r['Y_FIRST']))/float(r['Y_STEP']))
ysub[1]=int((Latsub[0]-float(r['Y_FIRST']))/float(r['Y_STEP']))
print 'Subseting geocoded',ext,' file with Latitude and Longitude...'
except:
print 'Not geocoded file, cannot be subseted with LatLon.'
Usage() ; sys.exit(1)
except:
Geo=0
box = (xsub[0],ysub[0],xsub[1],ysub[1])
output_img = im.crop(box)
print 'writing >>> '+outName
output_img.save(outName)
# try:
# r=readfile.read_rsc_file(File+'.rsc')
# except:
# sys.exit(1)
r['FILE_LENGTH']=str(ysub[1]-ysub[0])
r['WIDTH']=str(xsub[1]-xsub[0])
r['XMAX']=str(int(r['WIDTH']) - 1)
r['YMAX']=str(int(r['FILE_LENGTH']) - 1)
r['subset_x0']=str(xsub[0])
r['subset_x1']=str(xsub[1])
r['subset_y0']=str(ysub[0])
r['subset_y1']=str(ysub[1])
try:
r['Y_FIRST']=str(float(r['Y_FIRST'])+ysub[0]*float(r['Y_STEP']))
r['X_FIRST']=str(float(r['X_FIRST'])+xsub[0]*float(r['X_STEP']))
except:
Geo=0
f = open(outName+'.rsc','w')
for k in r.keys():
f.write(k+' '+r[k]+'\n')
f.close()
def main(argv):
try:
opts, args = getopt.getopt(argv, "h:f:t:p:")
except getopt.GetoptError:
Usage()
sys.exit(1)
if opts == []:
Usage()
sys.exit(1)
for opt, arg in opts:
if opt in ("-h", "--help"):
Usage()
sys.exit()
elif opt == '-f':
file = arg
elif opt == '-t':
filtType = arg
elif opt == '-p':
par = arg
# try:
# file=argv[0]
# alks=float(argv[1])
# rlks=float(argv[2])
# except:
# Usage();sys.exit(1)
ext = os.path.splitext(file)[1]
outName = file.split('.')[0] + '_' + filtType + ext
try:
par
except:
par = []
print '+++++++++++++++++++++++++++'
print 'Filter type : ' + filtType
print 'parameters : ' + str(par)
print '+++++++++++++++++++++++++++'
###############################################
if ext == '.int' or ext == '.slc':
a, p, r = readfile.read_complex64(file)
plks = multilook(p, alks, rlks)
alks = multilook(a, alks, rlks)
r['FILE_LENGTH'] = str(dlks.shape[0])
r['WIDTH'] = str(dlks.shape[1])
r['XMAX'] = str(int(r['WIDTH']) - 1)
r['YMAX'] = str(int(r['FILE_LENGTH']) - 1)
try:
r['Y_STEP'] = str(float(r['Y_STEP']) * alks)
r['X_STEP'] = str(float(r['X_STEP']) * rlks)
except:
Geo = 0
f = open(outName + '.rsc', 'w')
for k in r.keys():
f.write(k + ' ' + r[k] + '\n')
f.close()
elif ext == '.unw' or ext == '.cor' or ext == '.hgt':
a, p, r = readfile.read_float32(file)
plks = multilook(p, alks, rlks)
alks = multilook(a, alks, rlks)
writefile.write_float32(plks, outName)
r['FILE_LENGTH'] = str(dlks.shape[0])
r['WIDTH'] = str(dlks.shape[1])
r['XMAX'] = str(int(r['WIDTH']) - 1)
r['YMAX'] = str(int(r['FILE_LENGTH']) - 1)
try:
r['Y_STEP'] = str(float(r['Y_STEP']) * alks)
r['X_STEP'] = str(float(r['X_STEP']) * rlks)
except:
Geo = 0
f = open(outName + '.rsc', 'w')
for k in r.keys():
f.write(k + ' ' + r[k] + '\n')
f.close()
elif ext == ('.dem'):
d, r = readfile.read_dem(file)
dlks = multilook(d, alks, rlks)
print 'writing ' + outName
writefile.write_dem(dlks, outName)
r['FILE_LENGTH'] = str(dlks.shape[0])
r['WIDTH'] = str(dlks.shape[1])
r['XMAX'] = str(int(r['WIDTH']) - 1)
r['YMAX'] = str(int(r['FILE_LENGTH']) - 1)
try:
r['Y_STEP'] = str(float(r['Y_STEP']) * alks)
r['X_STEP'] = str(float(r['X_STEP']) * rlks)
except:
Geo = 0
f = open(outName + '.rsc', 'w')
for k in r.keys():
f.write(k + ' ' + r[k] + '\n')
f.close()
elif ext in ['.jpeg', 'jpg', 'png']:
import Image
im = Image.open(file)
width = im.size[0] / int(rlks)
height = im.size[1] / int(alks)
imlks = im.resize((width, height), Image.NEAREST)
print 'writing ' + outName
imlks.save(outName)
try:
r = readfile.read_rsc_file(file + '.rsc')
except:
sys.exit(1)
r['FILE_LENGTH'] = str(height)
r['WIDTH'] = str(width)
r['XMAX'] = str(int(r['WIDTH']) - 1)
r['YMAX'] = str(int(r['FILE_LENGTH']) - 1)
try:
r['Y_STEP'] = str(float(r['Y_STEP']) * alks)
r['X_STEP'] = str(float(r['X_STEP']) * rlks)
except:
Geo = 0
f = open(outName + '.rsc', 'w')
for k in r.keys():
f.write(k + ' ' + r[k] + '\n')
f.close()
elif ext == ('.h5'):
h5file = h5py.File(file, 'r')
# outName=file.split('.')[0]+'_a'+str(int(alks))+'lks_r'+str(int(rlks))+'lks.h5'
h5file_lks = h5py.File(outName, 'w')
if 'interferograms' in h5file.keys():
print 'Filtering the interferograms in space'
gg = h5file_lks.create_group('interferograms')
igramList = h5file['interferograms'].keys()
for igram in igramList:
print igram
unwSet = h5file['interferograms'][igram].get(igram)
unw = unwSet[0:unwSet.shape[0], 0:unwSet.shape[1]]
unw = filter(unw, filtType, par)
group = gg.create_group(igram)
dset = group.create_dataset(igram,
data=unw,
compression='gzip')
for key, value in h5file['interferograms'][
igram].attrs.iteritems():
group.attrs[key] = value
dset1 = h5file['mask'].get('mask')
mask = dset1[0:dset1.shape[0], 0:dset1.shape[1]]
group = h5file_lks.create_group('mask')
dset = group.create_dataset('mask', data=mask, compression='gzip')
elif 'timeseries' in h5file.keys():
print 'Filtering the time-series'
group = h5file_lks.create_group('timeseries')
dateList = h5file['timeseries'].keys()
for d in dateList:
print d
dset1 = h5file['timeseries'].get(d)
data = dset1[0:dset1.shape[0], 0:dset1.shape[1]]
data = filter(data, filtType, par)
dset = group.create_dataset(d, data=data, compression='gzip')
for key, value in h5file['timeseries'].attrs.iteritems():
group.attrs[key] = value
try:
dset1 = h5file['mask'].get('mask')
Mask = dset1[0:dset1.shape[0], 0:dset1.shape[1]]
# Masklks=multilook(Mask,alks,rlks)
group = h5file_lks.create_group('mask')
dset = group.create_dataset('mask',
data=Mask,
compression='gzip')
except:
print 'Filterd file does not include the maske'
elif 'temporal_coherence' in h5file.keys(
) or 'velocity' in h5file.keys() or 'mask' in h5file.keys():
k = h5file.keys()
print 'filtering the ' + k[0]
group = h5file_lks.create_group(k[0])
dset1 = h5file[k[0]].get(k[0])
data = dset1[0:dset1.shape[0], 0:dset1.shape[1]]
data = filter(data, filtType, par)
dset = group.create_dataset(k[0], data=data, compression='gzip')
for key, value in h5file[k[0]].attrs.iteritems():
group.attrs[key] = value
h5file.close()
h5file_lks.close()
def main(argv):
#outName='subsetIgrams.h5'
try:
opts, args = getopt.getopt(argv, "h:f:x:y:o:l:L:")
except getopt.GetoptError:
print 'Error while getting args'
Usage()
sys.exit(1)
for opt, arg in opts:
if opt in ("-h", "--help"):
Usage()
sys.exit()
elif opt == '-f':
File = arg
elif opt == '-y':
ysub = [int(i) for i in arg.split(':')]
ysub.sort()
elif opt == '-x':
xsub = [int(i) for i in arg.split(':')]
xsub.sort()
elif opt == '-o':
outName = arg
elif opt == '-l':
Latsub = [float(i) for i in arg.split(':')]
Latsub.sort()
elif opt == '-L':
Lonsub = [float(i) for i in arg.split(':')]
Lonsub.sort()
#####################################################
try:
File
xsub
ysub
except:
try:
File
Latsub
Lonsub
except:
Usage()
sys.exit(1)
try:
outName
except:
outName = 'subset_' + File
ext = os.path.splitext(File)[1]
if ext == '.h5':
try:
h5file = h5py.File(File, 'r')
except:
Usage()
sys.exit(1)
k = h5file.keys()
# convert LatLon to xy for geocoded file
try:
Latsub
Lonsub
if 'X_FIRST' in h5file[k[0]].attrs.keys():
xsub = [0] * 2
ysub = [0] * 2
xsub[0] = int(
(Lonsub[0] - float(h5file[k[0]].attrs['X_FIRST'])) /
float(h5file[k[0]].attrs['X_STEP']))
xsub[1] = int(
(Lonsub[1] - float(h5file[k[0]].attrs['X_FIRST'])) /
float(h5file[k[0]].attrs['X_STEP']))
ysub[0] = int(
(Latsub[1] - float(h5file[k[0]].attrs['Y_FIRST'])) /
float(h5file[k[0]].attrs['Y_STEP']))
ysub[1] = int(
(Latsub[0] - float(h5file[k[0]].attrs['Y_FIRST'])) /
float(h5file[k[0]].attrs['Y_STEP']))
print 'Subseting geocoded', ext, ' file with Latitude and Longitude...'
elif 'X_FIRST' in h5file[k[0]][h5file[k[0]].keys()[0]].attrs.keys(
): # for geocoded interferograms/coherence
igramList = h5file[k[0]].keys()
xsub = [0] * 2
ysub = [0] * 2
xsub[0] = int(
(Lonsub[0] -
float(h5file[k[0]][igramList[0]].attrs['X_FIRST'])) /
float(h5file[k[0]][igramList[0]].attrs['X_STEP']))
xsub[1] = int(
(Lonsub[1] -
float(h5file[k[0]][igramList[0]].attrs['X_FIRST'])) /
float(h5file[k[0]][igramList[0]].attrs['X_STEP']))
ysub[0] = int(
(Latsub[1] -
float(h5file[k[0]][igramList[0]].attrs['Y_FIRST'])) /
float(h5file[k[0]][igramList[0]].attrs['Y_STEP']))
ysub[1] = int(
(Latsub[0] -
float(h5file[k[0]][igramList[0]].attrs['Y_FIRST'])) /
float(h5file[k[0]][igramList[0]].attrs['Y_STEP']))
print 'Subseting geocoded', ext, ' file with Latitude and Longitude...'
else:
print 'Not geocoded file, cannot be subseted with LatLon.'
Usage()
sys.exit(1)
except:
Geo = 0
# k=h5file.keys()
if 'interferograms' in k:
igramList = h5file['interferograms'].keys()
h5out = h5py.File(outName, 'w')
gg = h5out.create_group('interferograms')
for igram in igramList:
print igram
dset1 = h5file['interferograms'][igram].get(igram)
group = gg.create_group(igram)
dset = group.create_dataset(igram,
data=dset1[ysub[0]:ysub[1],
xsub[0]:xsub[1]],
compression='gzip')
for key, value in h5file['interferograms'][
igram].attrs.iteritems():
group.attrs[key] = value
group.attrs['FILE_LENGTH'] = shape(dset1[ysub[0]:ysub[1],
xsub[0]:xsub[1]])[0]
group.attrs['WIDTH'] = shape(dset1[ysub[0]:ysub[1],
xsub[0]:xsub[1]])[1]
group.attrs['subset_x0'] = xsub[0]
group.attrs['subset_x1'] = xsub[1]
group.attrs['subset_y0'] = ysub[0]
group.attrs['subset_y1'] = ysub[1]
if 'X_FIRST' in h5file['interferograms'][igram].attrs.keys():
group.attrs['X_FIRST'] = float(
h5file['interferograms']
[igram].attrs['X_FIRST']) + xsub[0] * float(
h5file['interferograms'][igram].attrs['X_STEP'])
group.attrs['Y_FIRST'] = float(
h5file['interferograms']
[igram].attrs['Y_FIRST']) + ysub[0] * float(
h5file['interferograms'][igram].attrs['Y_STEP'])
gm = h5out.create_group('mask')
try:
Mset = h5file['mask'].get('mask')
dset = gm.create_dataset('mask',
data=Mset[ysub[0]:ysub[1],
xsub[0]:xsub[1]],
compression='gzip')
except:
print 'No group for mask found! It may cause problem in other processing steps.'
try:
Cset = h5file['meanCoherence'].get('meanCoherence')
gm = h5out.create_group('meanCoherence')
dset = gm.create_dataset('meanCoherence',
data=Cset[ysub[0]:ysub[1],
xsub[0]:xsub[1]],
compression='gzip')
except:
print 'No average coherence found in the File'
elif k[0] in ('coherence', 'wrapped'):
corList = h5file[k[0]].keys()
h5out = h5py.File(outName, 'w')
gg = h5out.create_group(k[0])
for cor in corList:
print cor
dset1 = h5file[k[0]][cor].get(cor)
group = gg.create_group(cor)
dset = group.create_dataset(cor,
data=dset1[ysub[0]:ysub[1],
xsub[0]:xsub[1]],
compression='gzip')
for key, value in h5file[k[0]][cor].attrs.iteritems():
group.attrs[key] = value
group.attrs['FILE_LENGTH'] = shape(dset1[ysub[0]:ysub[1],
xsub[0]:xsub[1]])[0]
group.attrs['WIDTH'] = shape(dset1[ysub[0]:ysub[1],
xsub[0]:xsub[1]])[1]
group.attrs['subset_x0'] = xsub[0]
group.attrs['subset_x1'] = xsub[1]
group.attrs['subset_y0'] = ysub[0]
group.attrs['subset_y1'] = ysub[1]
if 'X_FIRST' in h5file[k[0]][cor].attrs.keys():
group.attrs['X_FIRST'] = float(
h5file[k[0]][cor].attrs['X_FIRST']) + xsub[0] * float(
h5file[k[0]][cor].attrs['X_STEP'])
group.attrs['Y_FIRST'] = float(
h5file[k[0]][cor].attrs['Y_FIRST']) + ysub[0] * float(
h5file[k[0]][cor].attrs['Y_STEP'])
elif 'timeseries' in h5file.keys():
dateList = h5file['timeseries'].keys()
h5out = h5py.File(outName, 'w')
group = h5out.create_group('timeseries')
for d in dateList:
print d
dset1 = h5file['timeseries'].get(d)
dset = group.create_dataset(d,
data=dset1[ysub[0]:ysub[1],
xsub[0]:xsub[1]],
compression='gzip')
for key, value in h5file['timeseries'].attrs.iteritems():
group.attrs[key] = value
group.attrs['FILE_LENGTH'] = shape(dset1[ysub[0]:ysub[1],
xsub[0]:xsub[1]])[0]
group.attrs['WIDTH'] = shape(dset1[ysub[0]:ysub[1],
xsub[0]:xsub[1]])[1]
group.attrs['subset_x0'] = xsub[0]
group.attrs['subset_x1'] = xsub[1]
group.attrs['subset_y0'] = ysub[0]
group.attrs['subset_y1'] = ysub[1]
if 'X_FIRST' in h5file['timeseries'].attrs.keys():
group.attrs['X_FIRST'] = float(
h5file['timeseries'].attrs['X_FIRST']) + xsub[0] * float(
h5file['timeseries'].attrs['X_STEP'])
group.attrs['Y_FIRST'] = float(
h5file['timeseries'].attrs['Y_FIRST']) + ysub[0] * float(
h5file['timeseries'].attrs['Y_STEP'])
h5file.close()
h5out.close()
elif 'temporal_coherence' in h5file.keys(
) or 'velocity' in h5file.keys() or 'mask' in h5file.keys(
) or 'rmse' in h5file.keys():
print 'writing >>> ' + outName
dset = h5file[k[0]].get(k[0])
data = dset[ysub[0]:ysub[1], xsub[0]:xsub[1]]
hfout = h5py.File(outName, 'w')
group = hfout.create_group(k[0])
group.create_dataset(k[0], data=data, compression='gzip')
for key, value in h5file[k[0]].attrs.iteritems():
group.attrs[key] = value
group.attrs['FILE_LENGTH'] = data.shape[0]
group.attrs['WIDTH'] = data.shape[1]
group.attrs['XMIN'] = 0
group.attrs['XMAX'] = data.shape[1] - 1
group.attrs['YMIN'] = 0
group.attrs['YMAX'] = data.shape[0] - 1
group.attrs['subset_x0'] = xsub[0]
group.attrs['subset_x1'] = xsub[1]
group.attrs['subset_y0'] = ysub[0]
group.attrs['subset_y1'] = ysub[1]
if 'X_FIRST' in h5file[k[0]].attrs.keys():
group.attrs['X_FIRST'] = float(
h5file[k[0]].attrs['X_FIRST']) + xsub[0] * float(
h5file[k[0]].attrs['X_STEP'])
group.attrs['Y_FIRST'] = float(
h5file[k[0]].attrs['Y_FIRST']) + ysub[0] * float(
h5file[k[0]].attrs['Y_STEP'])
h5file.close()
hfout.close()
elif ext in ['.unw', '.cor', '.hgt']:
a, p, r = readfile.read_float32(File)
try:
Latsub
Lonsub
try:
r['X_FIRST']
xsub = [0] * 2
ysub = [0] * 2
xsub[0] = int(
(Lonsub[0] - float(r['X_FIRST'])) / float(r['X_STEP']))
xsub[1] = int(
(Lonsub[1] - float(r['X_FIRST'])) / float(r['X_STEP']))
ysub[0] = int(
(Latsub[1] - float(r['Y_FIRST'])) / float(r['Y_STEP']))
ysub[1] = int(
(Latsub[0] - float(r['Y_FIRST'])) / float(r['Y_STEP']))
print 'Subseting geocoded', ext, ' file with Latitude and Longitude...'
except:
print 'Not geocoded file, cannot be subseted with LatLon.'
Usage()
sys.exit(1)
except:
Geo = 0
a = a[ysub[0]:ysub[1], xsub[0]:xsub[1]]
p = p[ysub[0]:ysub[1], xsub[0]:xsub[1]]
print 'writing >>> ' + outName
writefile.write_float32(p, outName)
r['FILE_LENGTH'] = str(p.shape[0])
r['WIDTH'] = str(p.shape[1])
r['XMAX'] = str(int(r['WIDTH']) - 1)
r['YMAX'] = str(int(r['FILE_LENGTH']) - 1)
r['subset_x0'] = str(xsub[0])
r['subset_x1'] = str(xsub[1])
r['subset_y0'] = str(ysub[0])
r['subset_y1'] = str(ysub[1])
try:
r['Y_FIRST'] = str(
float(r['Y_FIRST']) + ysub[0] * float(r['Y_STEP']))
r['X_FIRST'] = str(
float(r['X_FIRST']) + xsub[0] * float(r['X_STEP']))
except:
Geo = 0
f = open(outName + '.rsc', 'w')
for k in r.keys():
f.write(k + ' ' + r[k] + '\n')
f.close()
elif ext == '.dem':
d, r = readfile.read_dem(File)
try:
Latsub
Lonsub
# print Latsub
try:
r['X_FIRST']
xsub = [0] * 2
ysub = [0] * 2
xsub[0] = int(
(Lonsub[0] - float(r['X_FIRST'])) / float(r['X_STEP']))
xsub[1] = int(
(Lonsub[1] - float(r['X_FIRST'])) / float(r['X_STEP']))
ysub[0] = int(
(Latsub[1] - float(r['Y_FIRST'])) / float(r['Y_STEP']))
ysub[1] = int(
(Latsub[0] - float(r['Y_FIRST'])) / float(r['Y_STEP']))
print 'Subseting', ext, ' file with Latitude and Longitude...'
except:
print 'Not geocoded file, cannot be subseted with LatLon.'
Usage()
sys.exit(1)
except:
Geo = 0
d = d[ysub[0]:ysub[1], xsub[0]:xsub[1]]
print 'writing >>> ' + outName
writefile.write_dem(d, outName)
r['FILE_LENGTH'] = str(d.shape[0])
r['WIDTH'] = str(d.shape[1])
r['XMAX'] = str(int(r['WIDTH']) - 1)
r['YMAX'] = str(int(r['FILE_LENGTH']) - 1)
r['subset_x0'] = str(xsub[0])
r['subset_x1'] = str(xsub[1])
r['subset_y0'] = str(ysub[0])
r['subset_y1'] = str(ysub[1])
try:
r['Y_FIRST'] = str(
float(r['Y_FIRST']) + ysub[0] * float(r['Y_STEP']))
r['X_FIRST'] = str(
float(r['X_FIRST']) + xsub[0] * float(r['X_STEP']))
except:
Geo = 0
f = open(outName + '.rsc', 'w')
for k in r.keys():
f.write(k + ' ' + r[k] + '\n')
f.close()
elif ext in ['.jpeg', 'jpg', 'png']:
import Image
im = Image.open(File)
try:
r = readfile.read_rsc_file(File + '.rsc')
except:
sys.exit(1)
try:
Latsub
Lonsub
try:
r['X_FIRST']
xsub = [0] * 2
ysub = [0] * 2
xsub[0] = int(
(Lonsub[0] - float(r['X_FIRST'])) / float(r['X_STEP']))
xsub[1] = int(
(Lonsub[1] - float(r['X_FIRST'])) / float(r['X_STEP']))
ysub[0] = int(
(Latsub[1] - float(r['Y_FIRST'])) / float(r['Y_STEP']))
ysub[1] = int(
(Latsub[0] - float(r['Y_FIRST'])) / float(r['Y_STEP']))
print 'Subseting geocoded', ext, ' file with Latitude and Longitude...'
except:
print 'Not geocoded file, cannot be subseted with LatLon.'
Usage()
sys.exit(1)
except:
Geo = 0
box = (xsub[0], ysub[0], xsub[1], ysub[1])
output_img = im.crop(box)
print 'writing >>> ' + outName
output_img.save(outName)
# try:
# r=readfile.read_rsc_file(File+'.rsc')
# except:
# sys.exit(1)
r['FILE_LENGTH'] = str(ysub[1] - ysub[0])
r['WIDTH'] = str(xsub[1] - xsub[0])
r['XMAX'] = str(int(r['WIDTH']) - 1)
r['YMAX'] = str(int(r['FILE_LENGTH']) - 1)
r['subset_x0'] = str(xsub[0])
r['subset_x1'] = str(xsub[1])
r['subset_y0'] = str(ysub[0])
r['subset_y1'] = str(ysub[1])
try:
r['Y_FIRST'] = str(
float(r['Y_FIRST']) + ysub[0] * float(r['Y_STEP']))
r['X_FIRST'] = str(
float(r['X_FIRST']) + xsub[0] * float(r['X_STEP']))
except:
Geo = 0
f = open(outName + '.rsc', 'w')
for k in r.keys():
f.write(k + ' ' + r[k] + '\n')
f.close()
def main(argv):
try:
file=argv[0]
alks=float(argv[1])
rlks=float(argv[2])
except:
Usage();sys.exit(1)
ext = os.path.splitext(file)[1]
outName=file.split('.')[0]+'_a'+str(int(alks))+'lks_r'+str(int(rlks))+'lks'+ext
if ext == '.int' or ext == '.slc':
a,p,r = readfile.read_complex64(file)
plks=multilook(p,alks,rlks)
alks=multilook(a,alks,rlks)
r['FILE_LENGTH']=str(dlks.shape[0])
r['WIDTH']=str(dlks.shape[1])
r['XMAX']=str(int(r['WIDTH']) - 1)
r['YMAX']=str(int(r['FILE_LENGTH']) - 1)
try:
r['Y_STEP']=str(float(r['Y_STEP'])*alks)
r['X_STEP']=str(float(r['X_STEP'])*rlks)
except:
Geo=0
f = open(outName+'.rsc','w')
for k in r.keys():
f.write(k+' '+r[k]+'\n')
f.close()
elif ext == '.unw' or ext == '.cor' or ext == '.hgt':
a,p,r = readfile.read_float32(file)
plks=multilook(p,alks,rlks)
alks=multilook(a,alks,rlks)
writefile.write_float32(plks,outName)
r['FILE_LENGTH']=str(dlks.shape[0])
r['WIDTH']=str(dlks.shape[1])
r['XMAX']=str(int(r['WIDTH']) - 1)
r['YMAX']=str(int(r['FILE_LENGTH']) - 1)
try:
r['Y_STEP']=str(float(r['Y_STEP'])*alks)
r['X_STEP']=str(float(r['X_STEP'])*rlks)
except:
Geo=0
f = open(outName+'.rsc','w')
for k in r.keys():
f.write(k+' '+r[k]+'\n')
f.close()
elif ext == ('.dem'):
d,r = readfile.read_dem(file)
dlks=multilook(d,alks,rlks)
print 'writing '+outName
writefile.write_dem(dlks,outName)
r['FILE_LENGTH']=str(dlks.shape[0])
r['WIDTH']=str(dlks.shape[1])
r['XMAX']=str(int(r['WIDTH']) - 1)
r['YMAX']=str(int(r['FILE_LENGTH']) - 1)
try:
r['Y_STEP']=str(float(r['Y_STEP'])*alks)
r['X_STEP']=str(float(r['X_STEP'])*rlks)
except:
Geo=0
f = open(outName+'.rsc','w')
for k in r.keys():
f.write(k+' '+r[k]+'\n')
f.close()
elif ext in ['.jpeg','jpg','png']:
import Image
im = Image.open(file)
width = im.size[0] / int(rlks)
height = im.size[1] / int(alks)
imlks = im.resize((width, height), Image.NEAREST)
print 'writing ' + outName
imlks.save(outName)
try:
r=readfile.read_rsc_file(file+'.rsc')
except:
sys.exit(1)
r['FILE_LENGTH']=str(height)
r['WIDTH']=str(width)
r['XMAX']=str(int(r['WIDTH']) - 1)
r['YMAX']=str(int(r['FILE_LENGTH']) - 1)
try:
r['Y_STEP']=str(float(r['Y_STEP'])*alks)
r['X_STEP']=str(float(r['X_STEP'])*rlks)
except:
Geo=0
f = open(outName+'.rsc','w')
for k in r.keys():
f.write(k+' '+r[k]+'\n')
f.close()
elif ext == ('.h5'):
h5file=h5py.File(file,'r')
# outName=file.split('.')[0]+'_a'+str(int(alks))+'lks_r'+str(int(rlks))+'lks.h5'
h5file_lks=h5py.File(outName,'w')
if 'interferograms' in h5file.keys():
print 'Multilooking the interferograms'
gg = h5file_lks.create_group('interferograms')
igramList=h5file['interferograms'].keys()
for igram in igramList:
print igram
unw = h5file['interferograms'][igram].get(igram)
unwlks=multilook(unw,alks,rlks)
group = gg.create_group(igram)
dset = group.create_dataset(igram, data=unwlks, compression='gzip')
for key, value in h5file['interferograms'][igram].attrs.iteritems():
group.attrs[key] = value
group.attrs['WIDTH']=unwlks.shape[1]
group.attrs['FILE_LENGTH']=unwlks.shape[0]
try:
group.attrs['Y_STEP']=alks*float(group.attrs['Y_STEP'])
group.attrs['X_STEP']=rlks*float(group.attrs['X_STEP'])
except:
group.attrs['AZIMUTH_PIXEL_SIZE']=alks*float(group.attrs['AZIMUTH_PIXEL_SIZE'])
group.attrs['RANGE_PIXEL_SIZE']=rlks*float(group.attrs['RANGE_PIXEL_SIZE'])
dset1=h5file['mask'].get('mask')
mask=dset1[0:dset1.shape[0],0:dset1.shape[1]]
masklks=multilook(mask,alks,rlks)
group=h5file_lks.create_group('mask')
dset = group.create_dataset('mask', data=masklks, compression='gzip')
elif 'timeseries' in h5file.keys():
print 'Multilooking the time-series'
group = h5file_lks.create_group('timeseries')
dateList=h5file['timeseries'].keys()
for d in dateList:
print d
unw = h5file['timeseries'].get(d)
unwlks=multilook(unw,alks,rlks)
dset = group.create_dataset(d, data=unwlks, compression='gzip')
for key,value in h5file['timeseries'].attrs.iteritems():
group.attrs[key] = value
group.attrs['WIDTH']=unwlks.shape[1]
group.attrs['FILE_LENGTH']=unwlks.shape[0]
try:
group.attrs['Y_STEP']=alks*float(group.attrs['Y_STEP'])
group.attrs['X_STEP']=rlks*float(group.attrs['X_STEP'])
except:
group.attrs['AZIMUTH_PIXEL_SIZE']=alks*float(group.attrs['AZIMUTH_PIXEL_SIZE'])
group.attrs['RANGE_PIXEL_SIZE']=rlks*float(group.attrs['RANGE_PIXEL_SIZE'])
try:
dset1 = h5file['mask'].get('mask')
Mask = dset1[0:dset1.shape[0],0:dset1.shape[1]]
Masklks=multilook(Mask,alks,rlks)
group=h5file_lks.create_group('mask')
dset = group.create_dataset('mask', data=Masklks, compression='gzip')
except:
print 'Multilooked file does not include the maske'
elif 'temporal_coherence' in h5file.keys() or 'velocity' in h5file.keys() or 'mask' in h5file.keys():
k=h5file.keys()
print 'multi looking the '+ k[0]
group=h5file_lks.create_group(k[0])
dset1 = h5file[k[0]].get(k[0])
Mask = dset1[0:dset1.shape[0],0:dset1.shape[1]]
Masklks=multilook(Mask,alks,rlks)
dset = group.create_dataset(k[0], data=Masklks, compression='gzip')
for key , value in h5file[k[0]].attrs.iteritems():
group.attrs[key]=value
try:
group.attrs['Y_STEP']=alks*float(group.attrs['Y_STEP'])
group.attrs['X_STEP']=rlks*float(group.attrs['X_STEP'])
except:
group.attrs['AZIMUTH_PIXEL_SIZE']=alks*float(group.attrs['AZIMUTH_PIXEL_SIZE'])
group.attrs['RANGE_PIXEL_SIZE']=rlks*float(group.attrs['RANGE_PIXEL_SIZE'])
group.attrs['WIDTH']=Masklks.shape[1]
group.attrs['FILE_LENGTH']=Masklks.shape[0]
h5file.close()
h5file_lks.close()
