def plot_file(file):
name = basename(file)
fname = 'plots/%s.pdf' % name
#fname = 'plots/%s.png'%name
if exists(fname):
print('Skip %s!' % fname)
return
print(fname)
with h5py.File(file, 'r') as fid:
slip = nres = fid['Bm'][...]
mplt = MapPlotSlab()
mplt.plot_top()
mplt = MapPlotFault(fault_file)
mplt.plot_slip(slip)
mplt = vj.MapPlotDisplacement()
mplt.plot_sites_seafloor(sites_seafloor=sites_seafloor)
mo, mw = vj.MomentCalculator(fault_file, earth_file).moment(slip)
plt.title('Mo=%g, Mw=%.2f' % (mo, mw))
#plt.show()
plt.savefig(fname)
plt.close()
def plot_file(file):
name = basename(file)
fname = 'plots/%s.pdf'%name
#fname = 'plots/%s.png'%name
if exists(fname):
print('Skip %s!'%fname)
return
print(fname)
with h5py.File(file,'r') as fid:
slip = nres = fid['Bm'][...]
mplt = MapPlotSlab()
mplt.plot_top()
mplt = MapPlotFault(fault_file)
mplt.plot_slip(slip)
mplt = vj.MapPlotDisplacement()
mplt.plot_sites_seafloor(sites_seafloor = sites_seafloor)
mo, mw = vj.MomentCalculator(fault_file, earth_file).moment(slip)
plt.title('Mo=%g, Mw=%.2f'%(mo, mw))
#plt.show()
plt.savefig(fname)
plt.close()
def plot_slip(self, slip, fn):
mplt = MapPlotFault(
'/home/zy/workspace/viscojapan/inversions/static_inversion/coseismic_inversion_wider/fault_model/fault_He50km_east.h5'
)
mplt.plot_slip(slip)
plt.savefig(join(self.outs_dir, fn))
plt.close()
def test_dip(self):
dips = self.read_pars('dip')
MapPlotFault(fault_file=self.original_fault_file).plot_slip(
dips.reshape([-1, 1]))
plt.title('dip')
plt.savefig(join(self.outs_dir, 'dip.png'))
plt.close()
def test_lonlat_grids(self):
lons = self.read_pars('lons')
lats = self.read_pars('lats')
self.plot_original_fault()
mplt = MapPlot()
mplt.basemap.plot(lons,lats,
latlon=True,marker='.', ls='', color='red')
MapPlotSlab().plot_top()
plt.savefig(join(self.outs_dir, 'lonlat.png'))
plt.close()
def test_lonlat_grids(self):
lons = self.read_pars('lons')
lats = self.read_pars('lats')
self.plot_original_fault()
mplt = MapPlot()
mplt.basemap.plot(lons,
lats,
latlon=True,
marker='.',
ls='',
color='red')
MapPlotSlab().plot_top()
plt.savefig(join(self.outs_dir, 'lonlat.png'))
plt.close()
def plot_file(file):
name = basename(file)
fname = 'plots/%s.png'%name
if exists(fname):
print('Skip %s!'%fname)
return
print(fname)
with h5py.File(file,'r') as fid:
slip = nres = fid['Bm'][...]
mplt = MapPlotSlab()
mplt.plot_top()
mplt = MapPlotFault('../fault_bott40km.h5')
mplt.plot_slip(slip)
mplt = vj.MapPlotDisplacement()
mplt.plot_sites_seafloor(sites_seafloor = sites_seafloor)
#plt.show()
plt.savefig(fname)
plt.close()
def plot_file(file):
name = basename(file)
fname = 'plots/%s.png' % name
if exists(fname):
print('Skip %s!' % fname)
return
print(fname)
with h5py.File(file, 'r') as fid:
slip = nres = fid['Bm'][...]
mplt = MapPlotSlab()
mplt.plot_top()
mplt = MapPlotFault('../fault_bott40km.h5')
mplt.plot_slip(slip)
mplt = vj.MapPlotDisplacement()
mplt.plot_sites_seafloor(sites_seafloor=sites_seafloor)
#plt.show()
plt.savefig(fname)
plt.close()
for file in files:
print(file)
bn = os.path.basename(file).split('.')[0]
plot_dir = 'plots/%s'%bn
if not exists(plot_dir):
makedirs(plot_dir)
with h5py.File(file) as fid:
Bm = fid['Bm'][...]
d_pred = fid['d_pred'][...]
for nth, epoch in enumerate(epochs[1:]):
print(epoch)
mplt = MapPlotFault(fault_file)
mplt.plot_slip(Bm[nth*num_subflts:
(nth+1)*num_subflts])
## mplt = vj.MapPlotDisplacement()
## mplt.plot_disp(d_pred[nth*num_obs:
## (nth+1)*num_obs],sites, scale=scale)
## sites = [ii.decode() for ii in d_ep.sites]
## mplt.plot_disp(d_ep[epoch], sites, color='red', scale=scale)
mplt = MapPlotSlab()
mplt.plot_top()
savefig(join(plot_dir, 'incr_slip_%04d.png'%epoch))
# plt.show()
plt.close()
for file in files:
print(file)
bn = os.path.basename(file).split('.')[0]
plot_dir = 'plots/%s'%bn
if not exists(plot_dir):
makedirs(plot_dir)
with h5py.File(file) as fid:
Bm = fid['Bm'][...]
d_pred = fid['d_pred'][...]
for nth, epoch in enumerate(epochs[1:]):
print(epoch)
mplt = MapPlotFault(fault_file)
mplt.plot_slip(Bm[nth*num_subflts:
(nth+1)*num_subflts])
## mplt = vj.MapPlotDisplacement()
## mplt.plot_disp(d_pred[nth*num_obs:
## (nth+1)*num_obs],sites, scale=scale)
## sites = [ii.decode() for ii in d_ep.sites]
## mplt.plot_disp(d_ep[epoch], sites, color='red', scale=scale)
mplt = MapPlotSlab()
mplt.plot_top()
savefig(join(plot_dir, 'incr_slip_%04d.png'%epoch))
# plt.show()
plt.close()
def plot_slip(self, slip, fn):
fault_file = '/home/zy/workspace/viscojapan/tests/share/fault_bott80km.h5'
mplt = MapPlotFault(fault_file)
mplt.plot_slip(slip)
plt.savefig(join(self.outs_dir, fn))
plt.close()
def test_dep_top(self):
dips = self.read_pars('top')
MapPlotFault(fault_file=self.original_fault_file).pcolor_on_fault(dips)
plt.title('top')
plt.savefig(join(self.outs_dir, 'dep_top.png'))
plt.close()
def plot_slip(self, slip, fn):
fault_file = '/home/zy/workspace/viscojapan/tests/share/fault_bott80km.h5'
mplt = MapPlotFault(fault_file)
mplt.plot_slip(slip)
plt.savefig(join(self.outs_dir,fn))
plt.close()
def plot_slip(self, slip, fn):
mplt = MapPlotFault('/home/zy/workspace/viscojapan/inversions/static_inversion/coseismic_inversion_wider/fault_model/fault_He50km_east.h5')
mplt.plot_slip(slip)
plt.savefig(join(self.outs_dir,fn))
plt.close()
def test_dip(self):
dips = self.read_pars('dip')
MapPlotFault(fault_file = self.original_fault_file).plot_slip(dips.reshape([-1,1]))
plt.title('dip')
plt.savefig(join(self.outs_dir, 'dip.png'))
plt.close()
def test_dep_top(self):
dips = self.read_pars('top')
MapPlotFault(fault_file = self.original_fault_file).pcolor_on_fault(dips)
plt.title('top')
plt.savefig(join(self.outs_dir, 'dep_top.png'))
plt.close()
