def plot_z_slice(data,
qx=None,
qy=None,
qz=None,
iz=0,
fignum=None,
if_faults=False,
cscale=None):
#Plot a map-view slice through the 3D matrix data at index iz
# qx,qy,qz are the coordinates along the dimensions of the matrix data
if qx is None:
qx = arange(0, data.shape[1])
if qy is None:
qy = arange(0, data.shape[0])
if qz is None:
qz = arange(0, data.shape[2])
figid = plt.figure(fignum)
figid.clf()
if if_faults:
from misc_tools import load_faults
faultlon, faultlat = load_faults()
plt.plot(faultlon, faultlat, 'k')
plt.axis([min(qx), max(qx), min(qy), max(qy)])
pc_id = plt.pcolor(qx, qy, data[:, :, iz])
if cscale is not None:
plt.clim(cscale)
figid.colorbar(pc_id)
figid.show()
def plot_z_slice(data,qx=None,qy=None,qz=None,iz=0,fignum=None,if_faults=False,cscale=None):
#Plot a map-view slice through the 3D matrix data at index iz
# qx,qy,qz are the coordinates along the dimensions of the matrix data
if qx is None:
qx=arange(0,data.shape[1])
if qy is None:
qy=arange(0,data.shape[0])
if qz is None:
qz=arange(0,data.shape[2])
figid=plt.figure(fignum)
figid.clf();
if if_faults:
from misc_tools import load_faults
faultlon,faultlat=load_faults()
plt.plot(faultlon,faultlat,'k')
plt.axis([min(qx),max(qx),min(qy),max(qy)])
pc_id=plt.pcolor(qx,qy,data[:,:,iz] )
if cscale is not None:
plt.clim(cscale)
figid.colorbar(pc_id)
figid.show()
arrsta = [] #a list of station names
for arrival in ev.arrivals:
if arrival.phase is 'P':
arrvec.append(arrival.ttime)
absvec.append(arrival.epoch)
arrsta.append(arrival.staname)
if not os.path.isfile(arrival.staname + 'traveltime'):
continue
absvec = asarray(absvec)
arrvec = asarray(arrvec)
#Return the entire misfit matrix
dstep = int(loc_params['dstep2'])
dx, dy, dz = nlon / dstep, nlat / dstep, nz / dstep
dx, dy, dz = 1, 1, 1
qx, qy, qz = range(0, nlon, dx), range(0, nlat, dy), range(0, nz, dz)
#minx, miny, minz, orgmin ,orimean,oristd= misc_tools.exp_grid_search(arrsta, qx, qy,qz, absvec, li,ev.epoch)
#orimean,oristd= misc_tools._grid_search_traveltimes_origin(arrsta, qx, qy,qz, absvec, li)
minx, miny, minz, orimean, oristd = misc_tools.exp_grid_search(
arrsta, qx, qy, qz, absvec, li)
#Plot
idepth = misc_tools.find_nearest_index(earth_rad - qdep, ev.depth)
faultlon, faultlat = misc_tools.load_faults()
figid = plt.figure(1)
figid.clf()
plt.plot(faultlon, faultlat, 'k')
plt.axis([min(qlon), max(qlon), min(qlat), max(qlat)])
pc_id = plt.pcolor(qlon, qlat, oristd[:, :, idepth])
figid.colorbar(pc_id)
figid.show()
arrvec=[]
arrsta=[] #a list of station names
for arrival in ev.arrivals:
if arrival.phase is 'P':
arrvec.append(arrival.ttime)
absvec.append(arrival.epoch)
arrsta.append(arrival.staname)
if not os.path.isfile(arrival.staname+'traveltime'):
continue
absvec=asarray(absvec)
arrvec=asarray(arrvec)
#Return the entire misfit matrix
dstep = int(loc_params['dstep2'])
dx, dy, dz = nlon / dstep, nlat / dstep, nz / dstep
dx,dy,dz=1,1,1
qx, qy, qz = range(0, nlon, dx), range(0, nlat, dy), range(0, nz, dz);
#minx, miny, minz, orgmin ,orimean,oristd= misc_tools.exp_grid_search(arrsta, qx, qy,qz, absvec, li,ev.epoch)
#orimean,oristd= misc_tools._grid_search_traveltimes_origin(arrsta, qx, qy,qz, absvec, li)
minx,miny,minz,orimean,oristd=misc_tools.exp_grid_search(arrsta,qx,qy,qz,absvec,li)
#Plot
idepth=misc_tools.find_nearest_index(earth_rad-qdep,ev.depth)
faultlon,faultlat=misc_tools.load_faults()
figid=plt.figure(1)
figid.clf();
plt.plot(faultlon,faultlat,'k')
plt.axis([min(qlon),max(qlon),min(qlat),max(qlat)])
pc_id=plt.pcolor(qlon,qlat,oristd[:,:,idepth] )
figid.colorbar(pc_id)
figid.show()
