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()