def main():
fig, axes = setup_figure()
functions = [in_phase_traces, out_of_phase_traces]
labels = ['Different Amplitude', 'Shifted Traces']
for i, (func, label) in enumerate(zip(functions, labels)):
row = axes[i]
wiggles(row[0], *func())
cross_corr(row[1], func)
for ax in row[2:]:
setup_trace_xy(ax, *func())
semblance(row[2], func)
eigenstructure(row[3], func)
row[0].set_ylabel(label, size=14)
# Locatable axes with shared axis's don't observe the adjustable param??
# Setting limits manually... Shouldn't be necessary...
for ax in axes[0,2:]:
ax.axis([-2.1, 2.1, -2.1, 2.1])
for ax in axes[1,2:]:
ax.axis([-1.1, 1.1, -1.1, 1.1])
axes[0, 0].set_title('Input Traces')
bottom_label(axes[1, 0], 'A')
axes[0, 1].set_title('Bahorich & Farmer (1995)\nCross Correlation')
bottom_label(axes[1, 1], 'B')
axes[0, 2].set_title('Marfurt, et al (1998)\nSemblance')
bottom_label(axes[1, 2], 'C')
axes[0, 3].set_title('Gersztenkorn & Marfurt (1999)\nEigenstructure')
bottom_label(axes[1, 3], 'D')
fig.savefig(image_dir('figure_2.png'), dpi=200, bbox_inches='tight')
def run(coherence_func, title, template):
trace_funcs = [figure_2.identical_traces, figure_2.in_phase_traces, figure_2.out_of_phase_traces]
ylabels = ["Identical Traces", "Different Amplitude", "Shifted Traces"]
for trace_func, ylabel in zip(trace_funcs, ylabels):
filename = template.format(ylabel.replace(" ", "_"))
fig = compare(trace_func, coherence_func, title, ylabel)
fig.savefig(image_dir(filename), bbox_inches="tight", dpi=80)
def run(coherence_func, title, template):
trace_funcs = [figure_2.identical_traces, figure_2.in_phase_traces,
figure_2.out_of_phase_traces]
ylabels = ['Identical Traces', 'Different Amplitude', 'Shifted Traces']
for trace_func, ylabel in zip(trace_funcs, ylabels):
filename = template.format(ylabel.replace(' ', '_'))
fig = compare(trace_func, figure_2.semblance, title, ylabel)
fig.savefig(image_dir(filename), bbox_inches='tight', dpi=80)
def run(coherence_func, title, template):
trace_funcs = [
figure_2.identical_traces, figure_2.in_phase_traces,
figure_2.out_of_phase_traces
]
ylabels = ['Identical Traces', 'Different Amplitude', 'Shifted Traces']
for trace_func, ylabel in zip(trace_funcs, ylabels):
filename = template.format(ylabel.replace(' ', '_'))
fig = compare(trace_func, coherence_func, title, ylabel)
fig.savefig(image_dir(filename), bbox_inches='tight', dpi=80)
def main():
seismic = data.load_seismic()
fig, axes = setup_figure()
plot(axes[:,0], 'Seismic Data', '', seismic, 'A')
plot(axes[:,1], 'Bahorich & Farmer (1995)', 'Cross-correlation',
basic_methods.bahorich_coherence(seismic, 21), 'B')
plot(axes[:,2], 'Marfurt, et al, (1998)', 'Semblance-based',
moving_window(seismic, (3,3,9), basic_methods.marfurt_semblance), 'C')
plot(axes[:,3], 'Gersztenkorn and Marfurt (1999)', 'Eigenstructure-based',
moving_window(seismic, (3,3,9), basic_methods.eig), 'D')
plot(axes[:,4], 'Randen, et al, (2000)', 'Gradient Structure Tensor Based',
gst_coherence(seismic, (3, 3, 9), sigma=2), 'E')
fig.savefig(image_dir('figure_1.png'), dpi=200, bbox_inches='tight')
plt.show()
def main():
seismic = data.load_seismic()
fig, axes = setup_figure()
plot(axes[:, 0], 'Seismic Data', '', seismic, 'A')
plot(axes[:, 1], 'Bahorich & Farmer (1995)', 'Cross-correlation',
basic_methods.bahorich_coherence(seismic, 21), 'B')
plot(axes[:, 2], 'Marfurt, et al, (1998)', 'Semblance-based',
moving_window(seismic, (3, 3, 9), basic_methods.marfurt_semblance),
'C')
plot(axes[:, 3], 'Gersztenkorn and Marfurt (1999)', 'Eigenstructure-based',
moving_window(seismic, (3, 3, 9), basic_methods.eig), 'D')
plot(axes[:, 4], 'Randen, et al, (2000)',
'Gradient Structure Tensor Based',
gst_coherence(seismic, (3, 3, 9), sigma=2), 'E')
fig.savefig(image_dir('figure_1.png'), dpi=200, bbox_inches='tight')
plt.show()