def main():
fig, ax = plt.subplots()
ax.set_xlabel('Eastward Slip (m)')
ax.set_ylabel('Northward Slip (m)')
offsets = [[0, 0]]
heaves = [[0, 0, 0]]
for hor, name in zip(data.horizons[::-1], data.gulick_names[::-1]):
results = get_result(hor.name)
discarded = 200 - results.shape[0]
print 'Discarded {} points from {}'.format(discarded, hor.name)
# Plot covariance ellipse...
dx, dy, slip = plot(ax, results, '# ' + name)
offsets.append([dx, dy])
heaves.append(utilities.calculate_heave([dx, dy], hor))
# Plot plate motion over 200kyr
utilities.plot_plate_motion(xy=offsets[3], time=2e5)
# Plot lines connecting the horizons...
ax.plot(*zip(*offsets), marker='o', color='darkred')
# Plot heaves...
heaves = np.array(heaves)
ax.plot(*heaves[:, :2].T, marker='o', color='green')
# Set aspect ratio of plot to 1 so that azimuths are properly represented
ax.axis('equal')
plt.show()
def main():
fig, ax = plt.subplots()
ax.set_xlabel('Eastward Slip (m)')
ax.set_ylabel('Northward Slip (m)')
offsets = [[0, 0]]
heaves = [[0,0,0]]
for hor, name in zip(data.horizons[::-1], data.gulick_names[::-1]):
results = get_result(hor.name)
discarded = 200 - results.shape[0]
print 'Discarded {} points from {}'.format(discarded, hor.name)
# Plot covariance ellipse...
dx, dy, slip = plot(ax, results, '# ' + name)
offsets.append([dx, dy])
heaves.append(utilities.calculate_heave([dx, dy], hor))
# Plot plate motion over 200kyr
utilities.plot_plate_motion(xy=offsets[3], time=2e5)
# Plot lines connecting the horizons...
ax.plot(*zip(*offsets), marker='o', color='darkred')
# Plot heaves...
heaves = np.array(heaves)
ax.plot(*heaves[:,:2].T, marker='o', color='green')
# Set aspect ratio of plot to 1 so that azimuths are properly represented
ax.axis('equal')
plt.show()
import geoprobe
import data
import utilities
from fault_kinematics import homogeneous_simple_shear
from process_bootstrap_results import get_result, load
fault = data.to_world(data.to_xyz(data.fault))
alpha = data.alpha
f, group = load()
for hor in data.horizons:
print hor.name
xyz = data.to_world(data.to_xyz(hor))
slip = get_result(hor.name, group)
restored = homogeneous_simple_shear.inclined_shear(fault, xyz, slip, alpha)
print 'Resampling...'
restored = data.to_model(restored)
restored = utilities.grid_xyz(restored)
new_hor = geoprobe.horizon(*restored.T)
new_hor.write('restored_horizons/' + hor.name + '.hzn')
f.close()
