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)