A Python interface to TRIPACK and STRIPACK Fortran code for (constrained) triangulation in Cartesian coordinates and on a sphere. Stripy is an object-oriented package and includes routines from SRFPACK and SSRFPACK for interpolation (nearest neighbor, linear and hermite cubic) and to evaluate derivatives (Renka 1996a,b and 1997a,b).
stripy is bundled with litho1pt0 which is a python interface to the crust 1.0 dataset and the lithospheric part of the litho 1.0 dataset (Laske et al, 2013 and Pasyanos et al, 2014) which both requires and demonstrates the triangulation / searching and interpolation on the sphere that is provided by stripy.
Sample images created with stripy illustrating the meshing capability: ocean age data can be triangulated on the sphere with no need for points on land. Once stripy ingests your data points, you can sample another dataset to your grid (bathymetry on the right), smooth, find the derivatives of your data, or interpolate to another set of points.
There are two matching sets of stripy notebooks - one set for Cartesian Triangulations and one for Spherical Triangulations. For most geographical applications, the spherical triangulations are the natural choice as they expect longitude and latitude coordinates (admittedly in radians). There are some worked examples
which use the companion package litho1pt0
- Documentation / Notebooks https://underworldcode.github.io/stripy/2.0.5b2
- API documentation https://underworldcode.github.io/stripy/2.0.5b2_api
- Documentation / Notebooks https://underworldcode.github.io/stripy/2.1.0b1
- API documentation https://underworldcode.github.io/stripy/2.1.0b1_api
For previous versions, see the changelog
Launch the demonstration
(This is the development branch)
Moresi, L. and Mather, B.R., (2019). Stripy: A Python module for (constrained) triangulation in Cartesian coordinates and on a sphere.. Journal of Open Source Software, 4(38), 1410, https://doi.org/10.21105/joss.01410
There are two matching sets of stripy notebooks - one set for Cartesian Triangulations and one for Spherical Triangulations. For most geographical applications, the spherical triangulations are the natural choice as they expect longitude and latitude coordinates (admittedly in radians).
Note: the Cartesian and Spherical notebooks can be obtained / installed from stripy itself as follows:
python -c 'import stripy; stripy.documentation.install_documentation(path="Notebooks")'- Ex1-Cartesian-Triangulations.ipynb
- Ex2-CartesianGrids.ipynb
- Ex3-Interpolation.ipynb
- Ex4-Gradients.ipynb
- Ex5-Smoothing.ipynb
- Ex6-Scattered-Data.ipynb
- Ex7-Refinement-of-Triangulations.ipynb
- Ex8-Spline-Tension.ipynb
- Ex9-Voronoi-Diagram.ipynb
- Ex1-Spherical-Triangulations.ipynb
- Ex2-SphericalGrids.ipynb
- Ex3-Interpolation.ipynb
- Ex4-Gradients.ipynb
- Ex5-Smoothing.ipynb
- Ex6-Scattered-Data.ipynb
- Ex7-Refinement-of-Triangulations.ipynb
- Ex8-Spline-Tension.ipynb
- Ex9-Voronoi-Diagram.ipynb
Note, these examples are the notebooks from litho1pt0 which are installed from the
package itself:
python -c 'import litho1pt0; litho1pt0.documentation.install_documentation(path="Notebooks")'The first three notebooks are an introduction to litho1pt0 that does not explicitly mention stripy but
the next two worked examples show how to search, interpolate and plot with the help of stripy routines.
- Ex1-Litho1Layers.ipynb
- Ex2-Litho1Properties.ipynb
- Ex3-CrustalRegionalisation.ipynb
- WorkEx1-CratonAverageProperties.ipynb
- WorkEx2-OceanDepthAge.ipynb
You will need Python 3.6+. Also, the following packages are required:
Recommended Packages for running the notebooks:
All of which should be available from pip or anaconda (conda-forge) for most platforms.
You can install stripy using the
pip package manager with either version of Python:
python3 -m pip install stripyAll the dependencies will be automatically installed by pip, except for gfortran
(or any Fortran compiler). It must be installed in your system before installing
stripy with pip.
If you change the Fortran compiler, you may have to add the
flags config_fc --fcompiler=<compiler name> when setup.py is run
(see docs for numpy.distutils).
If you use the anaconda packaging system, then you should be able to
conda install -c geo-down-under stripyTwo classes are included as part of the Stripy package:
sTriangulation(Spherical coordinates)Triangulation(Cartesian coordinates)
These classes share similar methods and can be easily interchanged. In addition, there are many helper functions provided for building meshes.
A series of tests are located in the tests subdirectory.
In order to perform these tests clone the repository and run pytest:
git checkout https://github.com/underworldcode/stripy.git
cd stripy
pytest -v-
Laske, G., G. Masters, and Z. Ma (2013), Update on CRUST1. 0—A 1-degree global model of Earth's crust, Geophys Research Abstracts, 15, EGU2013–2658.
-
Pasyanos, M. E., T. G. Masters, G. Laske, and Z. Ma (2014), LITHO1.0: An updated crust and lithospheric model of the Earth, Journal of Geophysical Research-Solid Earth, 119(3), 2153–2173, doi:10.1002/2013JB010626.
-
R. J. Renka, "ALGORITHM 751: TRIPACK: A Constrained Two- Dimensional Delaunay Triangulation Package" ACM Trans. Math. Software, Vol. 22, No. 1, 1996, pp. 1-8.
-
R. J. Renka, "ALGORITHM 752: SRFPACK: Software for Scattered Data Fitting with a Constrained Surface under Tension", ACM Trans. Math. Software, Vol. 22, No. 1, 1996, pp. 9-17.
-
R. J. Renka, "ALGORITHM 772: STRIPACK: Delaunay Triangulation and Voronoi Diagram on the Surface of a Sphere" ACM Trans. Math. Software, Vol. 23, No. 3, 1997, pp. 416-434.
-
R. J. Renka, "ALGORITHM 773: SSRFPACK: Interpolation of Scattered Data on the Surface of a Sphere with a Surface under Tension", ACM Trans. Math. Software, Vol. 23, No. 3, 1997, pp. 437-439.