seg/tutorials-2014

Geophysical Tutorials for 2014

tutorials-2014's Issues

issue with syntheic data generation

hi
i have an issue if you plz help me out here
when i use the code

area = [5000, 25000, 5000, 25000]
bounds = [0, 30000, 0, 30000, 0, 5000]
mpl_fig = plt.figure()

mpl.axis('scaled')
mpl.square(area)
dike_verts = mpl.draw_polygon(bounds[:4], mpl.gca(), xy2ne=True)
dike_verts
dike = PolygonalPrism(dike_verts, 0, 5000, {'magnetization': 10})
myv.figure(size=(600, 400))
myv.polyprisms([dike], linewidth=2)
myv.axes(myv.outline(bounds), ranges=[b*0.001 for b in bounds], nlabels=3, fmt='%.1f')
myv.wall_north(bounds)
myv.wall_bottom(bounds)
myv.savefig('tmp/model_dike.png')
myv.show()
Image(filename='tmp/model_dike.png')

i receive this error

C:\ProgramData\Anaconda2\lib\site-packages\fatiando\vis\mpl.py:70: UserWarning: This module will be removed in v0.6. We recommend the use of matplotlib.pyplot module directly. Some of the fatiando specific functions will remain.
"specific functions will remain.")
C:\ProgramData\Anaconda2\lib\site-packages\fatiando\vis\myv.py:53: UserWarning: This module will be removed in v0.7.
warnings.warn("This module will be removed in v0.7.")

AttributeError Traceback (most recent call last)
in ()
8 from fatiando.vis import myv, mpl
9 import fatiando
---> 10 print("Using Fatiando a Terra version {0}".format((fatiando.version)))
11 area = [5000, 25000, 5000, 25000]
12 bounds = [0, 30000, 0, 30000, 0, 5000]

AttributeError: 'module' object has no attribute 'version'

i am stuck here would you please help me thanks

Amplitude extraction

Hello! I've been working for a few days with the wedge model code (1412_Tuning_and_AVO/tuning_wedge_v2.py). Thanks for sharing it. I'm not a programmer and it is really helping me.

I would like to share a question and hear your opinion about it: in the algorithm, the graph "thickness vs amplitude" is constructed by extracting the amplitude over the upper time of the modeled wedge layer, that is, a time-slice extraction. However, as the algorithm itself shows, sometimes the seismic reflector mismatches the modeled interface.

So, my question is: wouldn't it be more appropriate to make te code extract the maximum amplitude of the synthetic reflector? By doing this we would be simulating how the amplitude would actually be if the reflector was mapped.

To illustrate, I made this change in the code and generated two figures, comparing this graph to the scenarios:

First, as the code is now:

And with the suggested amplitude extraction to the maximum:

Thanks!

New Script for Python version 3+ >> example in the 1406_Make_a_synthetic folder

Changes to "las.py" and "how_to_make_synthetic.ipynb" files, to meet Python 3+ requirements.

https://gist.github.com/marcfreir/fd4bea2921a2c77e3ad68266fbaf2799

LASReader gives invalid syntax error

The code for the LASReader https://scipy-cookbook.readthedocs.io/items/LASReader.html , which is used in this notebook, gives an invalid syntax error on line 102 and 174. Can someone explain why? and how can it be solved.
Thanks.