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()
def dip_corrected(seismic, window, func):
surface = data.load_horizon()
flat = basic_methods.flatten(seismic, surface, seismic.shape[-1])
sembl = moving_window(flat, window, func)
return basic_methods.unflatten(sembl, surface, seismic.shape)
def dip_corrected(seismic, window, func):
surface = data.load_horizon()
flat = basic_methods.flatten(seismic, surface, seismic.shape[-1])
sembl = moving_window(flat, window, func)
return basic_methods.unflatten(sembl, surface, seismic.shape)