fteikpy
|License| |Stars| |Pyversions| |Version| |Downloads| |Code style: black| |Codacy Badge| |Codecov| |Build| |Docs| |DOI|
fteikpy is a Python library that computes accurate first arrival traveltimes in 2D and 3D heterogeneous isotropic velocity models. The algorithm handles properly the curvature of wavefronts close to the source which can be placed without any problem between grid points.
The code is based on FTeik <https://github.com/Mark-Noble/FTeik-Eikonal-Solver>
__ implemented in Python and compiled just-in-time <https://en.wikipedia.org/wiki/Just-in-time_compilation>
__ with numba <https://numba.pydata.org/>
__.
.. figure:: https://raw.githubusercontent.com/keurfonluu/fteikpy/master/.github/sample.gif
:alt: sample-marmousi
:width: 100%
:align: center
Computation of traveltimes and ray-tracing on smoothed Marmousi velocity model.
Features
Forward modeling:
- Compute traveltimes in 2D and 3D Cartesian grids with the possibility to use a different grid spacing in Z, X and Y directions,
- Compute traveltime gradients at runtime or a posteriori,
- A posteriori 2D and 3D ray-tracing.
Parallel:
- Traveltime grids are seemlessly computed in parallel for different sources,
- Raypaths from a given source to different locations are also evaluated in parallel.
Installation
The recommended way to install fteikpy and all its dependencies is through the Python Package Index:
.. code:: bash
pip install fteikpy --user
Otherwise, clone and extract the package, then run from the package location:
.. code:: bash
pip install . --user
To test the integrity of the installed package, check out this repository and run:
.. code:: bash
pytest
Documentation
Refer to the online documentation <https://keurfonluu.github.io/fteikpy/>
__ for detailed description of the API and examples.
Alternatively, the documentation can be built using Sphinx <https://www.sphinx-doc.org/en/master/>
__:
.. code:: bash
pip install -r doc/requirements.txt
sphinx-build -b html doc/source doc/build
Usage
The following example computes the traveltime grid in a 3D homogeneous velocity model:
.. code-block:: python
import numpy as np
from fteikpy import Eikonal3D
Velocity model
velocity_model = np.ones((8, 8, 8))
dz, dx, dy = 1.0, 1.0, 1.0
Solve Eikonal at source
eik = Eikonal3D(velocity_model, gridsize=(dz, dx, dy))
tt = eik.solve((0.0, 0.0, 0.0))
Get traveltime at specific grid point
t1 = tt[0, 1, 2]
Or get traveltime at any point in the grid
t2 = tt(np.random.rand(3) * 7.0)
Contributing
Please refer to the Contributing Guidelines <https://github.com/keurfonluu/fteikpy/blob/master/CONTRIBUTING.rst>
__ to see how you can help. This project is released with a Code of Conduct <https://github.com/keurfonluu/fteikpy/blob/master/CODE_OF_CONDUCT.rst>
__ which you agree to abide by when contributing.
.. |License| image:: https://img.shields.io/github/license/keurfonluu/fteikpy
:target: https://github.com/keurfonluu/fteikpy/blob/master/LICENSE
.. |Stars| image:: https://img.shields.io/github/stars/keurfonluu/fteikpy?logo=github
:target: https://github.com/keurfonluu/fteikpy
.. |Pyversions| image:: https://img.shields.io/pypi/pyversions/fteikpy.svg?style=flat
:target: https://pypi.org/pypi/fteikpy/
.. |Version| image:: https://img.shields.io/pypi/v/fteikpy.svg?style=flat
:target: https://pypi.org/project/fteikpy
.. |Downloads| image:: https://pepy.tech/badge/fteikpy
:target: https://pepy.tech/project/fteikpy
.. |Code style: black| image:: https://img.shields.io/badge/code%20style-black-000000.svg?style=flat
:target: https://github.com/psf/black
.. |Codacy Badge| image:: https://img.shields.io/codacy/grade/bec3d2ad6b8c45cf9bb0da110fe04838.svg?style=flat
:target: https://www.codacy.com/gh/keurfonluu/fteikpy/dashboard?utm_source=github.com&utm_medium=referral&utm_content=keurfonluu/fteikpy&utm_campaign=Badge_Grade
.. |Codecov| image:: https://img.shields.io/codecov/c/github/keurfonluu/fteikpy.svg?style=flat
:target: https://codecov.io/gh/keurfonluu/fteikpy
.. |DOI| image:: https://zenodo.org/badge/DOI/10.5281/zenodo.4269352.svg?style=flat
:target: https://doi.org/10.5281/zenodo.4269352
.. |Build| image:: https://img.shields.io/github/workflow/status/keurfonluu/fteikpy/Python%20package
:target: https://github.com/keurfonluu/fteikpy
.. |Docs| image:: https://img.shields.io/github/workflow/status/keurfonluu/fteikpy/Build%20documentation?label=docs
:target: https://keurfonluu.github.io/fteikpy/