Example #1
0
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)
Example #5
0
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()
Example #6
0
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()