def main(argv):
try:
file = argv[0]
except:
usage()
sys.exit(1)
print '\n*************** Save to GRD file for GMT ****************'
print 'Input file: ' + file
ext = os.path.splitext(file)[1]
########## PySAR HDF5 Files ################
if ext == '.h5':
h5 = h5py.File(file, 'r')
k = h5.keys()
if 'interferograms' in k: k[0] = 'interferograms'
elif 'coherence' in k: k[0] = 'coherence'
elif 'timeseries' in k: k[0] = 'timeseries'
if k[0] in ('interferograms', 'coherence', 'wrapped'):
atr = h5[k[0]][h5[k[0]].keys()[0]].attrs
elif k[0] in ('dem', 'velocity', 'mask', 'temporal_coherence', 'rmse',
'timeseries'):
atr = h5[k[0]].attrs
if 'timeseries' in k:
try:
d = sys.argv[2]
except:
print 'No input date... continue to convert the last date of timeseries'
dateList = h5['timeseries'].keys()
d = dateList[-1]
print 'reading ' + d + ' ... '
dset = h5['timeseries'].get(d)
elif k[0] in ['interferograms', 'coherence', 'wrapped']:
print 'interferograms is not supported currently.'
sys.exit(1)
else:
dset = h5[k[0]].get(k[0])
z = dset[0:dset.shape[0], 0:dset.shape[1]]
########## ROI_PAC Files ##################
elif ext in ['.unw', '.cor', '.hgt', '.dem']:
if ext == '.dem': z, atr = readfile.read_real_int16(file)
if ext in ['.unw', '.cor', '.hgt']:
a, z, atr = readfile.read_float32(file)
else:
print 'Unsupported file extension: ' + ext
sys.exit(1)
lats, lons = get_geo_lat_lon(atr)
#fname='velocity.grd'
outName = os.path.basename(file).replace(ext, '') + '.grd'
print 'Output (GMT grd file): ' + outName
_gmt.write_gmt_simple(lons, flipud(lats), flipud(z), outName, title='default', name='velocity',\
scale=1.0, offset=0, units='meters')
def main(argv):
try: file = argv[0]
except: Usage(); sys.exit(1)
print '\n*************** Save to GRD file for GMT ****************'
print 'Input file: '+file
ext = os.path.splitext(file)[1]
########## PySAR HDF5 Files ################
if ext == '.h5':
h5=h5py.File(file,'r')
k=h5.keys()
if 'interferograms' in k: k[0] = 'interferograms'
elif 'coherence' in k: k[0] = 'coherence'
elif 'timeseries' in k: k[0] = 'timeseries'
if k[0] in ('interferograms','coherence','wrapped'):
atr = h5[k[0]][h5[k[0]].keys()[0]].attrs
elif k[0] in ('dem','velocity','mask','temporal_coherence','rmse','timeseries'):
atr = h5[k[0]].attrs
if 'timeseries' in k:
try:
d=sys.argv[2]
except:
print 'No input date... continue to convert the last date of timeseries'
dateList=h5['timeseries'].keys()
d=dateList[-1]
print 'reading '+d + ' ... '
dset=h5['timeseries'].get(d)
elif k[0] in ['interferograms','coherence','wrapped']:
print 'interferograms is not supported currently.'; sys.exit(1)
else: dset=h5[k[0]].get(k[0])
z=dset[0:dset.shape[0],0:dset.shape[1]]
########## ROI_PAC Files ##################
elif ext in ['.unw','.cor','.hgt','.dem']:
import pysar._readfile as readfile
if ext == '.dem' : z,atr = readfile.read_real_int16(file)
if ext in ['.unw','.cor','.hgt']: a,z,atr = readfile.read_float32(file)
else: print 'Unsupported file extension: '+ext; sys.exit(1)
lats,lons=get_geo_lat_lon(atr)
#fname='velocity.grd'
outName=os.path.basename(file).replace(ext,'')+'.grd'; print 'Output (GMT grd file): '+outName
_gmt.write_gmt_simple(lons, flipud(lats), flipud(z), outName, title='default', name='velocity',\
scale=1.0, offset=0, units='meters')
def write_grd_file(data, atr, fname_out=None):
'''Write GMT .grd file for input data matrix, using giant._gmt module.
Inputs:
data - 2D np.array in int/float, data matrix to write
atr - dict, attributes of input data matrix
fname_out - string, output file name
Output:
fname_out - string, output file name
'''
if not fname_out:
fname_out = os.path.splitext(atr['FILE_PATH'])[0] + '.grd'
# Get 1D array of lats and lons
lats, lons = get_geo_lat_lon(atr)
# writing
print 'writing >>> ' + fname_out
gmt.write_gmt_simple(lons, np.flipud(lats), np.flipud(data), fname_out,\
title='default', name=atr['FILE_TYPE'], scale=1.0, offset=0, units=atr['UNIT'])
return fname_out
lats,lons=get_geo_lat_lon(h5)
if 'timeseries' in k:
try:
d=sys.argv[2]
except:
print 'No input date... continue to convert the last date of timeseries'
dateList=h5['timeseries'].keys()
d=dateList[-1]
print 'reading '+d + ' ... '
dset=h5['timeseries'].get(d)
elif 'interferograms' in k:
print 'interferograms'
else:
dset=h5[k[0]].get(k[0])
z=dset[0:dset.shape[0],0:dset.shape[1]]
#fname='velocity.grd'
outName=os.path.basename(file).replace('.h5','')
fname = outName + '.grd'
print 'Output (GMT grd file): '+ fname
_gmt.write_gmt_simple(lons, flipud(lats), flipud(z), fname, title='default', name='velocity', scale=1.0, offset=0, units='meters')
#_gmt.write_gmt_simple(lons, lats, res, fname, title=title, name=name, units=units)
print 'No input date... continue to convert the last date of timeseries'
dateList = h5['timeseries'].keys()
d = dateList[-1]
print 'reading ' + d + ' ... '
dset = h5['timeseries'].get(d)
elif 'interferograms' in k:
print 'interferograms'
else:
dset = h5[k[0]].get(k[0])
z = dset[0:dset.shape[0], 0:dset.shape[1]]
#fname='velocity.grd'
outName = os.path.basename(file).replace('.h5', '')
fname = outName + '.grd'
print 'Output (GMT grd file): ' + fname
_gmt.write_gmt_simple(lons,
flipud(lats),
flipud(z),
fname,
title='default',
name='velocity',
scale=1.0,
offset=0,
units='meters')
#_gmt.write_gmt_simple(lons, lats, res, fname, title=title, name=name, units=units)
