###############
pangeo-pyinterp
###############
|Build Status| |Downloads| |Platforms|
|Latest Release Date| |License| |Binder|
Python library for optimized geo-referenced interpolation.
About
=====
The motivation of this project is to provide tools for interpolating
geo-referenced data used in the field of geosciences. Other libraries cover this
problem, but written entirely in Python, the performance of these projects was
not quite sufficient for our needs. That is why this project started.
With this library, you can interpolate `2D
<https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.grid.Grid2D.html#pyinterp.grid.Grid2D>`_,
`3D
<https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.grid.Grid3D.html#pyinterp.grid.Grid3D>`_,
or `4D
<https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.grid.Grid4D.html#pyinterp.grid.Grid4D>`_
fields using ``n-variate`` and ``bicubic`` `interpolators
<https://pangeo-pyinterp.readthedocs.io/en/latest/api.html#cartesian-interpolators>`_
and `unstructured grids
<https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.RTree.html>`_.
You can also apply for a data `binning
<https://pangeo-pyinterp.readthedocs.io/en/latest/api.html#binning>`_ on the
bivariate area by simple or linear binning.
The library core is written in C++ using the `Boost C++ Libraries
<https://www.boost.org/>`_, `Eigen3 <http://eigen.tuxfamily.org/>`_ and
`pybind11 <https://github.com/pybind/pybind11/>`_ libraries.
This software also uses `CMake <https://cmake.org/>`_ to configure the project
and `Googletest <https://github.com/google/googletest>`_ to perform unit testing
of the library kernel.
Fill undefined values
=====================
The undefined values in the grids do not allow interpolation of values located
in the neighborhood. This behavior is a concern when you need to interpolate
values near the mask of some fields. The library provides utilities to fill the
undefined values:
* `loess <https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.fill.loess.html>`_
to fill the undefined values on the boundary between the defined/undefined
values using local regression.
* `gauss_seidel <https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.fill.gauss_seidel.html>`_
to fill all undefined values in a grid using the Gauss-Seidel method by
relaxation.
Geographic indexers
===================
N-Dimensional Grids
-------------------
N-dimensional grid is a grid defined by a matrix, in a 2D space, by a cube in a
3D space, etc. Each dimension of the grid is associated with a vector
corresponding to its coordinates or axes. Axes used to locate a pixel in the
grid from the coordinates of a point. These axes are either:
* *regular*: a vector of 181 latitudes spaced a degree from -90 to 90 degrees;
* *irregular*: a vector of 109 latitudes irregularly spaced from -90 to
89.940374 degrees.
These objects are manipulated by the class `pyinterp.Axis
<https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.Axis.html>`_,
which will choose, according to Axis definition, the best implementation. This
object will allow you to find the two indexes framing a given value. This
operating mode allows better performance when searching for a regular axis (a
simple calculation will enable you to see the index of a point immediately). In
contrast, in the case of an irregular axis, the search will be performed using a
binary search.
Finally, this class can define a circular axis from a vector to correctly
locate a value on the circle. This type of Axis will is used for handling
longitudes.
Temporal Axes
-------------
The `pyinterp.TemporalAxis
<https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.TemporalAxis.html>`_
class handles temporal axes, i.e., axes defined by 64-bit integer vectors, which
is the encoding used by `numpy
<https://docs.scipy.org/doc/numpy/reference/arrays.datetime.html>`_ to control
dates. This class allows you to process dates using integer arithmetic to ensure
that no information is lost during calculations. These objects are used by
spatiotemporal grids to perform temporal interpolations.
Unstructured Grids
------------------
In the case of unstructured grids, the index used is a *R\*Tree*. These trees
have better performance than the *KDTree* generally found in Python library
implementations.
The tree used here is the implementation provided by the `C++ Boost library
<https://www.boost.org/doc/libs/1_70_0/libs/geometry/doc/html/geometry/reference/spatial_indexes/boost__geometry__index__rtree.html>`_.
An adaptation has been introduced to address spherical equatorial coordinates
effectively. Although the Boost library allows these coordinates to manipulate
natively, the performance is lower than in the case of Cartesian space. Thus, we
have chosen to implement a conversion of Longitude Latitude Altitude (*LLA*)
coordinates into Earth-Centered, Earth-Fixed (*ECEF*) coordinates transparently
for the user to ensure that we can preserve excellent performance. The
disadvantage of this implementation is that it requires fairly more memory, as
one more element gets used to index the value of the Cartesian space.
The management of the `LLA
<https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.geodetic.Coordinates.ecef_to_lla.html>`_/`ECEF
<https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.geodetic.Coordinates.lla_to_ecef.html>`_
coordinate conversion is managed to use the `Olson, D.K.
<https://ieeexplore.ieee.org/document/481290>`_ algorithm. It has an excellent
performance with the accuracy of 1e-8 meters for altitude.
Geohash
-------
Geohashing is a geocoding method used to encode geographic coordinates
(latitude and longitude) into a short string of digits and letters delineating
an area on a map, which is called a cell, with varying resolutions. The more
characters in the string, the more precise the location.
Geohashes use Base-32 alphabet encoding (characters can be ``0`` to ``9`` and
``A`` to ``Z``, excl ``A``, ``I``, ``L`` and ``O``).
The geohash is a compact way of representing a location, and is useful for
storing a location in a database, or for indexing a location in a database.
.. |Build Status| image:: https://github.com/CNES/pangeo-pyinterp/actions/workflows/ci.yaml/badge.svg
:target: https://github.com/CNES/pangeo-pyinterp/actions
.. |Downloads| image:: https://anaconda.org/conda-forge/pyinterp/badges/downloads.svg?service=github
:target: image:: https://www.anaconda.com/distribution/
.. |Platforms| image:: https://anaconda.org/conda-forge/pyinterp/badges/platforms.svg?service=github
:target: https://anaconda.org/conda-forge/pyinterp
.. |Latest Release Date| image:: https://anaconda.org/conda-forge/pyinterp/badges/latest_release_date.svg?service=github
:target: https://github.com/CNES/pangeo-pyinterp/commits/master
.. |License| image:: https://anaconda.org/conda-forge/pyinterp/badges/license.svg?service=github
:target: https://opensource.org/licenses/BSD-3-Clause
.. |Binder| image:: https://mybinder.org/badge_logo.svg
:target: https://mybinder.org/v2/gh/CNES/pangeo-pyinterp/master?filepath=notebooks/auto_examples/
Raw data
{
"_id": null,
"home_page": "https://github.com/CNES/pangeo-pyinterp",
"name": "pyinterp",
"maintainer": null,
"docs_url": null,
"requires_python": ">=3.10",
"maintainer_email": null,
"keywords": "interpolation, geospatial, geohash, geodetic",
"author": "CNES/CLS",
"author_email": "fbriol@gmail.com",
"download_url": "https://files.pythonhosted.org/packages/f3/bf/d38d8b649f9628e4490831ec503e6fd5796bb34e5701237d06bda4f60620/pyinterp-2024.11.0.tar.gz",
"platform": "POSIX",
"description": "###############\npangeo-pyinterp\n###############\n\n|Build Status| |Downloads| |Platforms|\n|Latest Release Date| |License| |Binder|\n\nPython library for optimized geo-referenced interpolation.\n\nAbout\n=====\n\nThe motivation of this project is to provide tools for interpolating\ngeo-referenced data used in the field of geosciences. Other libraries cover this\nproblem, but written entirely in Python, the performance of these projects was\nnot quite sufficient for our needs. That is why this project started.\n\nWith this library, you can interpolate `2D\n<https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.grid.Grid2D.html#pyinterp.grid.Grid2D>`_,\n`3D\n<https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.grid.Grid3D.html#pyinterp.grid.Grid3D>`_,\nor `4D\n<https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.grid.Grid4D.html#pyinterp.grid.Grid4D>`_\nfields using ``n-variate`` and ``bicubic`` `interpolators\n<https://pangeo-pyinterp.readthedocs.io/en/latest/api.html#cartesian-interpolators>`_\nand `unstructured grids\n<https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.RTree.html>`_.\nYou can also apply for a data `binning\n<https://pangeo-pyinterp.readthedocs.io/en/latest/api.html#binning>`_ on the\nbivariate area by simple or linear binning.\n\nThe library core is written in C++ using the `Boost C++ Libraries\n<https://www.boost.org/>`_, `Eigen3 <http://eigen.tuxfamily.org/>`_ and\n`pybind11 <https://github.com/pybind/pybind11/>`_ libraries.\n\nThis software also uses `CMake <https://cmake.org/>`_ to configure the project\nand `Googletest <https://github.com/google/googletest>`_ to perform unit testing\nof the library kernel.\n\nFill undefined values\n=====================\n\nThe undefined values in the grids do not allow interpolation of values located\nin the neighborhood. This behavior is a concern when you need to interpolate\nvalues near the mask of some fields. The library provides utilities to fill the\nundefined values:\n\n* `loess <https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.fill.loess.html>`_\n to fill the undefined values on the boundary between the defined/undefined\n values using local regression.\n* `gauss_seidel <https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.fill.gauss_seidel.html>`_\n to fill all undefined values in a grid using the Gauss-Seidel method by\n relaxation.\n\nGeographic indexers\n===================\n\nN-Dimensional Grids\n-------------------\n\nN-dimensional grid is a grid defined by a matrix, in a 2D space, by a cube in a\n3D space, etc. Each dimension of the grid is associated with a vector\ncorresponding to its coordinates or axes. Axes used to locate a pixel in the\ngrid from the coordinates of a point. These axes are either:\n\n* *regular*: a vector of 181 latitudes spaced a degree from -90 to 90 degrees;\n* *irregular*: a vector of 109 latitudes irregularly spaced from -90 to\n 89.940374 degrees.\n\nThese objects are manipulated by the class `pyinterp.Axis\n<https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.Axis.html>`_,\nwhich will choose, according to Axis definition, the best implementation. This\nobject will allow you to find the two indexes framing a given value. This\noperating mode allows better performance when searching for a regular axis (a\nsimple calculation will enable you to see the index of a point immediately). In\ncontrast, in the case of an irregular axis, the search will be performed using a\nbinary search.\n\nFinally, this class can define a circular axis from a vector to correctly\nlocate a value on the circle. This type of Axis will is used for handling\nlongitudes.\n\nTemporal Axes\n-------------\n\nThe `pyinterp.TemporalAxis\n<https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.TemporalAxis.html>`_\nclass handles temporal axes, i.e., axes defined by 64-bit integer vectors, which\nis the encoding used by `numpy\n<https://docs.scipy.org/doc/numpy/reference/arrays.datetime.html>`_ to control\ndates. This class allows you to process dates using integer arithmetic to ensure\nthat no information is lost during calculations. These objects are used by\nspatiotemporal grids to perform temporal interpolations.\n\nUnstructured Grids\n------------------\n\nIn the case of unstructured grids, the index used is a *R\\*Tree*. These trees\nhave better performance than the *KDTree* generally found in Python library\nimplementations.\n\nThe tree used here is the implementation provided by the `C++ Boost library\n<https://www.boost.org/doc/libs/1_70_0/libs/geometry/doc/html/geometry/reference/spatial_indexes/boost__geometry__index__rtree.html>`_.\n\nAn adaptation has been introduced to address spherical equatorial coordinates\neffectively. Although the Boost library allows these coordinates to manipulate\nnatively, the performance is lower than in the case of Cartesian space. Thus, we\nhave chosen to implement a conversion of Longitude Latitude Altitude (*LLA*)\ncoordinates into Earth-Centered, Earth-Fixed (*ECEF*) coordinates transparently\nfor the user to ensure that we can preserve excellent performance. The\ndisadvantage of this implementation is that it requires fairly more memory, as\none more element gets used to index the value of the Cartesian space.\n\nThe management of the `LLA\n<https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.geodetic.Coordinates.ecef_to_lla.html>`_/`ECEF\n<https://pangeo-pyinterp.readthedocs.io/en/latest/generated/pyinterp.geodetic.Coordinates.lla_to_ecef.html>`_\ncoordinate conversion is managed to use the `Olson, D.K.\n<https://ieeexplore.ieee.org/document/481290>`_ algorithm. It has an excellent\nperformance with the accuracy of 1e-8 meters for altitude.\n\nGeohash\n-------\n\nGeohashing is a geocoding method used to encode geographic coordinates\n(latitude and longitude) into a short string of digits and letters delineating\nan area on a map, which is called a cell, with varying resolutions. The more\ncharacters in the string, the more precise the location.\n\nGeohashes use Base-32 alphabet encoding (characters can be ``0`` to ``9`` and\n``A`` to ``Z``, excl ``A``, ``I``, ``L`` and ``O``).\n\nThe geohash is a compact way of representing a location, and is useful for\nstoring a location in a database, or for indexing a location in a database.\n\n\n.. |Build Status| image:: https://github.com/CNES/pangeo-pyinterp/actions/workflows/ci.yaml/badge.svg\n :target: https://github.com/CNES/pangeo-pyinterp/actions\n.. |Downloads| image:: https://anaconda.org/conda-forge/pyinterp/badges/downloads.svg?service=github\n :target: image:: https://www.anaconda.com/distribution/\n.. |Platforms| image:: https://anaconda.org/conda-forge/pyinterp/badges/platforms.svg?service=github\n :target: https://anaconda.org/conda-forge/pyinterp\n.. |Latest Release Date| image:: https://anaconda.org/conda-forge/pyinterp/badges/latest_release_date.svg?service=github\n :target: https://github.com/CNES/pangeo-pyinterp/commits/master\n.. |License| image:: https://anaconda.org/conda-forge/pyinterp/badges/license.svg?service=github\n :target: https://opensource.org/licenses/BSD-3-Clause\n.. |Binder| image:: https://mybinder.org/badge_logo.svg\n :target: https://mybinder.org/v2/gh/CNES/pangeo-pyinterp/master?filepath=notebooks/auto_examples/\n",
"bugtrack_url": null,
"license": "BSD License",
"summary": "Interpolation of geo-referenced data for Python.",
"version": "2024.11.0",
"project_urls": {
"Homepage": "https://github.com/CNES/pangeo-pyinterp"
},
"split_keywords": [
"interpolation",
" geospatial",
" geohash",
" geodetic"
],
"urls": [
{
"comment_text": "",
"digests": {
"blake2b_256": "f3bfd38d8b649f9628e4490831ec503e6fd5796bb34e5701237d06bda4f60620",
"md5": "59aff2eade1b5d7c5905ed805f0e4e74",
"sha256": "e1c3f684953d351cb4485b8816a098351991702d78118e1daeea8e840decb896"
},
"downloads": -1,
"filename": "pyinterp-2024.11.0.tar.gz",
"has_sig": false,
"md5_digest": "59aff2eade1b5d7c5905ed805f0e4e74",
"packagetype": "sdist",
"python_version": "source",
"requires_python": ">=3.10",
"size": 4818364,
"upload_time": "2024-11-07T14:37:38",
"upload_time_iso_8601": "2024-11-07T14:37:38.124316Z",
"url": "https://files.pythonhosted.org/packages/f3/bf/d38d8b649f9628e4490831ec503e6fd5796bb34e5701237d06bda4f60620/pyinterp-2024.11.0.tar.gz",
"yanked": false,
"yanked_reason": null
}
],
"upload_time": "2024-11-07 14:37:38",
"github": true,
"gitlab": false,
"bitbucket": false,
"codeberg": false,
"github_user": "CNES",
"github_project": "pangeo-pyinterp",
"travis_ci": false,
"coveralls": false,
"github_actions": true,
"lcname": "pyinterp"
}
JSON
