spherical-geometry

Python based tools for spherical geometry

These details have not been verified by PyPI

Project links

Project description

Spherical Geometry's Github Actions CI Status

https://codecov.io/gh/spacetelescope/spherical_geometry/branch/master/graph/badge.svg

The spherical_geometry library is a Python package for handling spherical polygons that represent arbitrary regions of the sky.

Installation

On PyPI:

pip install spherical-geometry

On conda:

conda install -c conda-forge spherical-geometry

Requirements

Python 3.9 or later
Numpy 1.20 or later
astropy 5.0.4 or later
qd-library 2.3.7 or later (optional: if not available, the bundled version will be used). To force using the system-installed version, build with USE_SYSTEM_QD=1 setup.py build.

Bundled qd-library

Origin: https://www.davidhbailey.com/dhbsoftware/qd-2.3.22.tar.gz

A custom libqd/include/qd/qd_config.h is provided to circumvent the need to run any configuration scripts. This generalized configuration may not be optimized for your system.
The spherical_geometry test suite fail when linked to a system-optimized qd library, because the tests are written for the general case.

Coordinate representation

Coordinates in world space are traditionally represented by right ascension and declination (ra and dec), or longitude and latitude. While these representations are convenient, they have discontinuities at the poles, making operations on them trickier at arbitrary locations on the sky sphere. Therefore, all internal operations of this library are done in 3D vector space, where coordinates are represented as (x, y, z) vectors. The spherical_geometry.vector module contains functions to convert between (ra, dec) and (x, y, z) representations.

While any (x, y, z) triple represents a vector and therefore a location on the sky sphere, a distinction must be made between normalized coordinates that fall exactly on the unit sphere, and unnormalized coordinates which do not. A normalized coordinate is defined as a vector whose length is 1, i.e.:

\begin{equation*} \sqrt{x^2 + y^2 + z^2} = 1 \end{equation*}

To prevent unnecessary recomputation, many methods in this library assume that the vectors passed in are already normalized. If this is not the case, spherical_geometry.vector.normalize_vector can be used to normalize an array of vectors.

When not working in Cartesian vectors, the library allows the user to work in either degrees or radians. All methods that require or return an angular value have a degrees keyword argument. When degrees is True, these measurements are in degrees, otherwise they are in radians.

Spherical polygons

Spherical polygons are arbitrary areas on the sky sphere enclosed by great circle arcs. They are represented by the ~spherical_geometry.polygon.SphericalPolygon class.

Representation

The points defining the polygon are available from the ~polygon.SphericalPolygon.points property. It is a Nx3 array where each row is an (x, y, z) vector, normalized. The polygon points are explicitly closed, i.e., the first and last points are the same.

Where is the inside?

The edges of a polygon serve to separate the “inside” from the “outside” area. On a traditional 2D planar surface, the “inside” is defined as the finite area and the “outside” is the infinite area. However, since the surface of a sphere is cyclical, i.e., it wraps around on itself, the a spherical polygon actually defines two finite areas. To specify which should be considered the “inside” vs. the “outside”, the definition of the polygon also has an “inside point” which is just any point that should be considered inside of the polygon.

In the following image, the inside point (marked with the red dot) declares that the area of the polygon is the green region, and not the white region.

The inside point of the the polygon can be obtained from the ~polygon.SphericalPolygon.inside property.

What is the orientation?

The correctness of several of the algorithms using polygons depends on a consistent orientation of the points defining it. That is, the points should have a clockwise order. When creating a new spherical polygon, the order of the points defining a polygon will be reversed if they are not in clockwise order. The method SphericalPolygon.is_clockwise is used to est if the points are in clockwise order. It takes two successive sides and computes the normal vector to the sides. If the normal vector points inward towards the center of the sphere, the two sides are counter clockwise. If the normal vector points outward, the two sides are clockwise. The code determines the orientation by computing the triple product of the two sides with the vertex of the the two sides. Summing the triple product over all the sides gives the predominant orientation of the points in the polygon.

Disjoint Polygons

If a polygon is the result of the intersection of polygons, it may be disjoint. Disjoint polygons are represented as a list of spherical polygons. The library handles the details of this internally. However, the user must be aware that several of the properties of polygons are generators and return the value for a single polygon at a time. To access all the values of a proeprty, either use a for loop, or coerce the property to a list. The properties which are generators are:

SphericalPolygon.points: The points defining each polygon

SphericalPolygon.inside : The inside point of each polygon

If the intersection of two polygons generates disjoint polygons the code computes a new interior point for the disjoint polygons.

Creating spherical polygons

SphericalPolygon objects have 5 different constructors:

SphericalPolygon: Takes an array of (x, y, z) points, or a list of disjoint SphericalPolygon instances.

SphericalPolygon.from_radec: Takes an array of (ra, dec) points and an inside point.

SphericalPolygon.from_cone: Creates a polygon from a cone on the sky shere. Takes (ra, dec, radius).

SphericalPolygon.from_wcs: Creates a polygon from the footprint of a FITS image using its WCS header keywords. Takes a FITS filename or a astropy.io.fits.Header object.

SphericalPolygon.convex_hull: Creates a polygon that is the convex hull of a list of points.

Operations on Spherical Polygons

Once one has a SphericalPolygon object, there are a number of operations available:

~SphericalPolygon.contains_point: Determines if the given point is inside the polygon.

~SphericalPolygon.intersects_poly: Determines if one polygon intersects with another.

~SphericalPolygon.area: Determine the area of a polygon.

~SphericalPolygon.union and ~SphericalPolygon.multi_union: Return a new polygon that is the union of two or more polygons.

~SphericalPolygon.intersection and ~SphericalPolygon.multi_intersection: Return a new polygon that is the intersection of two or more polygons.

~SphericalPolygon.overlap: Determine how much a given polygon overlaps another.

~SphericalPolygon.to_radec: Convert (x, y, z) points in the polygon to (ra, dec) points.

~SphericalPolygon.draw: Plots the polygon using matplotlib’s Basemap toolkit. This feature is rather bare and intended primarily for debugging purposes.

Great circle arcs

As seen above, great circle arcs are used to define the edges of the polygon. The spherical_geometry.great_circle_arc module contains a number of functions that are useful for dealing with them.

length: Returns the angular distance between two points on the sphere.
intersection: Returns the intersection point between two great circle arcs.
intersects: Determines if two great circle arcs intersect.
intersects_point: Determines if a point is along the great circle arc.
angle: Calculate the angle between two great circle arcs.
midpoint: Calculate the midpoint along a great circle arc.

Project details

These details have not been verified by PyPI

Project links

Release history Release notifications | RSS feed

1.3.2

Jun 13, 2024

1.3.1

Nov 5, 2023

1.3.0

Nov 2, 2023

This version

1.3.0rc1 pre-release

Oct 18, 2023

1.2.23

Oct 21, 2022

1.2.22

Jan 5, 2022

1.2.21

Dec 31, 2021

1.2.20

May 13, 2021

1.2.19

Sep 18, 2020

1.2.18

Oct 14, 2019

1.2.17

Aug 15, 2019

1.2.16

Aug 14, 2019

1.2.15

Jul 31, 2019

1.2.14

Jul 30, 2019

1.2.11

Jun 25, 2019

1.2.10

Jul 24, 2019

1.2.9

Feb 22, 2019

1.2.8

Feb 22, 2019

1.2.7

Feb 22, 2019

1.2.6

Feb 22, 2019

1.2.5

Jul 10, 2018

1.2.4

Jun 28, 2018

1.2.3

Jun 20, 2018

1.2.2

May 4, 2018

1.1.0

Jan 19, 2018

1.0.12

Aug 2, 2017

1.0.11

Jul 28, 2017

1.0.6

Mar 22, 2016

1.0.4

Aug 29, 2014

Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

spherical_geometry-1.3.0rc1.tar.gz (9.2 MB view details)

Uploaded Oct 18, 2023 Source

Built Distributions

spherical_geometry-1.3.0rc1-cp312-cp312-win_amd64.whl (7.7 MB view details)

Uploaded Oct 18, 2023 CPython 3.12 Windows x86-64

spherical_geometry-1.3.0rc1-cp312-cp312-win32.whl (7.7 MB view details)

Uploaded Oct 18, 2023 CPython 3.12 Windows x86

spherical_geometry-1.3.0rc1-cp312-cp312-musllinux_1_1_x86_64.whl (9.1 MB view details)

Uploaded Oct 18, 2023 CPython 3.12 musllinux: musl 1.1+ x86-64

spherical_geometry-1.3.0rc1-cp312-cp312-manylinux_2_17_x86_64.manylinux2014_x86_64.whl (8.6 MB view details)

Uploaded Oct 18, 2023 CPython 3.12 manylinux: glibc 2.17+ x86-64

spherical_geometry-1.3.0rc1-cp312-cp312-macosx_11_0_arm64.whl (7.7 MB view details)

Uploaded Oct 18, 2023 CPython 3.12 macOS 11.0+ ARM64

spherical_geometry-1.3.0rc1-cp312-cp312-macosx_10_9_x86_64.whl (7.7 MB view details)

Uploaded Oct 18, 2023 CPython 3.12 macOS 10.9+ x86-64

spherical_geometry-1.3.0rc1-cp311-cp311-win_amd64.whl (7.7 MB view details)

Uploaded Oct 18, 2023 CPython 3.11 Windows x86-64

spherical_geometry-1.3.0rc1-cp311-cp311-win32.whl (7.7 MB view details)

Uploaded Oct 18, 2023 CPython 3.11 Windows x86

spherical_geometry-1.3.0rc1-cp311-cp311-musllinux_1_1_x86_64.whl (9.1 MB view details)

Uploaded Oct 18, 2023 CPython 3.11 musllinux: musl 1.1+ x86-64

spherical_geometry-1.3.0rc1-cp311-cp311-manylinux_2_17_x86_64.manylinux2014_x86_64.whl (8.6 MB view details)

Uploaded Oct 18, 2023 CPython 3.11 manylinux: glibc 2.17+ x86-64

spherical_geometry-1.3.0rc1-cp311-cp311-macosx_11_0_arm64.whl (7.7 MB view details)

Uploaded Oct 18, 2023 CPython 3.11 macOS 11.0+ ARM64

spherical_geometry-1.3.0rc1-cp311-cp311-macosx_10_9_x86_64.whl (7.7 MB view details)

Uploaded Oct 18, 2023 CPython 3.11 macOS 10.9+ x86-64

spherical_geometry-1.3.0rc1-cp310-cp310-win_amd64.whl (7.7 MB view details)

Uploaded Oct 18, 2023 CPython 3.10 Windows x86-64

spherical_geometry-1.3.0rc1-cp310-cp310-win32.whl (7.7 MB view details)

Uploaded Oct 18, 2023 CPython 3.10 Windows x86

spherical_geometry-1.3.0rc1-cp310-cp310-musllinux_1_1_x86_64.whl (9.1 MB view details)

Uploaded Oct 18, 2023 CPython 3.10 musllinux: musl 1.1+ x86-64

spherical_geometry-1.3.0rc1-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl (8.6 MB view details)

Uploaded Oct 18, 2023 CPython 3.10 manylinux: glibc 2.17+ x86-64

spherical_geometry-1.3.0rc1-cp310-cp310-macosx_11_0_arm64.whl (7.7 MB view details)

Uploaded Oct 18, 2023 CPython 3.10 macOS 11.0+ ARM64

spherical_geometry-1.3.0rc1-cp310-cp310-macosx_10_9_x86_64.whl (7.7 MB view details)

Uploaded Oct 18, 2023 CPython 3.10 macOS 10.9+ x86-64

spherical_geometry-1.3.0rc1-cp39-cp39-win_amd64.whl (7.7 MB view details)

Uploaded Oct 18, 2023 CPython 3.9 Windows x86-64

spherical_geometry-1.3.0rc1-cp39-cp39-win32.whl (7.7 MB view details)

Uploaded Oct 18, 2023 CPython 3.9 Windows x86

spherical_geometry-1.3.0rc1-cp39-cp39-musllinux_1_1_x86_64.whl (9.1 MB view details)

Uploaded Oct 18, 2023 CPython 3.9 musllinux: musl 1.1+ x86-64

spherical_geometry-1.3.0rc1-cp39-cp39-manylinux_2_17_x86_64.manylinux2014_x86_64.whl (8.6 MB view details)

Uploaded Oct 18, 2023 CPython 3.9 manylinux: glibc 2.17+ x86-64

spherical_geometry-1.3.0rc1-cp39-cp39-macosx_11_0_arm64.whl (7.7 MB view details)

Uploaded Oct 18, 2023 CPython 3.9 macOS 11.0+ ARM64

spherical_geometry-1.3.0rc1-cp39-cp39-macosx_10_9_x86_64.whl (7.7 MB view details)

Uploaded Oct 18, 2023 CPython 3.9 macOS 10.9+ x86-64

File details

Details for the file spherical_geometry-1.3.0rc1.tar.gz.

File metadata

Download URL: spherical_geometry-1.3.0rc1.tar.gz
Upload date: Oct 18, 2023
Size: 9.2 MB
Tags: Source
Uploaded using Trusted Publishing? No
Uploaded via: twine/4.0.2 CPython/3.11.6

File hashes

Hashes for spherical_geometry-1.3.0rc1.tar.gz
Algorithm	Hash digest
SHA256	`d7955c12be53e8e9ee57150ec5c33b1bc281f669fa523de7760d1f9de81ba467`
MD5	`524b9265e690450c76eb7a81ab0e645d`
BLAKE2b-256	`29d00fe777031941cd5dd197f6d80646f14d2a89c666d76ffd922932a02a5805`

Algorithm	Hash digest
SHA256	`23e5fd5c58f2ec6548a3f36564a61dd07d85c1b3977f0cc5e02376b33402dca8`
MD5	`936ce9918372d4e51639da33bdd14d11`
BLAKE2b-256	`bf5bb176b9dbe3b7c5affd31f569bb8ab641ff2bdec05c824c2fb7bc5f7d270d`

Algorithm	Hash digest
SHA256	`40f040c6d54b1a083c3268fadb0553704e0b45e6aac45eaa44e91ab3add58fc4`
MD5	`4ecaf79ce6155cc0837be8f5bc11aea2`
BLAKE2b-256	`fde2d2ef39970fff251bd7c33294f2abfa7f2964bae79d8c48b13345abc45e5f`

Algorithm	Hash digest
SHA256	`0a8bc9a092568ee9020d2f69f7eb10ae9bce039222aeda3487152c874e6ce4e7`
MD5	`c820bf62b86d2c58c48d9a0ed1c19ad9`
BLAKE2b-256	`a8207e6ee088c3555568269fb858629a0e4ea9d318a6eb7d48f38fad93aae10d`

Algorithm	Hash digest
SHA256	`00f31855cfc0e8784bd02c62077d61f4ff97a7aa1dcde136949d10e78b38cf4e`
MD5	`68007629541dbd98c35cb3017687026d`
BLAKE2b-256	`cd6bd4e841baea98006758cf0cc2f5dd318e217d69b6635b0fb9d43a79cb3dbb`

Algorithm	Hash digest
SHA256	`9a0d2fb5e46fde573a8385144006138a845d58e7e5609eda1240f61cac5d3abb`
MD5	`0f20b9f1697ef90ce7c58df63d65d924`
BLAKE2b-256	`fccf313b7667a5d6665e2ef4942874372b7eb45b80a494b54bfebc74446d4b2c`

Algorithm	Hash digest
SHA256	`5ea9bcb127cbcbf082ad7dad26ec20b535b68fd7f2b99b97f8bcf970db7f089c`
MD5	`e6375d299c0b9ba606a3c850fafc693d`
BLAKE2b-256	`f7a6ab5dd8e26fc0f82f189527fbd188bd18237504e7b32975b91c5a30d8b5d6`

Algorithm	Hash digest
SHA256	`708160704a4c8b0ee01acc89f3d31bcbc19132c2202338f0224fd78f6ecd8ace`
MD5	`2958ac4c27eabc67fb41cca097976074`
BLAKE2b-256	`09c644a3d631f61442f7308a26df4e5d01aa5ad17fda44ab81aa868cc967cd13`

Algorithm	Hash digest
SHA256	`282fe222cd16e454803108c544138b9bd62809a0c0c9cebd41a16d5bf02a9079`
MD5	`de8fcdeaa3afda5bc7d4f98c599e02fb`
BLAKE2b-256	`0567ca4d5da84e847d55dcb128aa769c90fe2b09f64a23aa655d20aca06a0fec`

Algorithm	Hash digest
SHA256	`17a1ee168008a88634af0712e3408f5ee4e47ee88c752f1664d234626f7ff7c4`
MD5	`69311eb76969ec614d52992b6be8ffd0`
BLAKE2b-256	`c83a22bd626b481a932d5ff5bbe1e6413e70a159a08c805bdfdc4ac6e150a504`

Algorithm	Hash digest
SHA256	`dfce7dc21f2e3f0e89c7719670a5f286bdcd627434c0573f327bab034261be03`
MD5	`189a2296b6d53c450065971356eb04c4`
BLAKE2b-256	`e8159d9a4c059df57e44f166cc9727d69aa770e877148d74f55abee00447a982`

Algorithm	Hash digest
SHA256	`a7dfc84030f0735e2cebac710403d50759e113b947af949f8ff3fb7a2872a619`
MD5	`587fa8b947f329f7ca20a17ea8831049`
BLAKE2b-256	`9cd451c4a9ec38fb6bd9f2aeff2c53b0757f9e16a91e72554d24f870a563ce90`

Algorithm	Hash digest
SHA256	`1f4af16c8c28bd6f0bf20e6d21d875d60eadf4ef6d980973a3b2ae4255e7ef41`
MD5	`19d80d5fe98ae220fc0f0f756b7a31b0`
BLAKE2b-256	`b63663bbb1d3ca192e05d8f81098c6ac2325c921e2c32cc3765ee467457e33a2`

Algorithm	Hash digest
SHA256	`96e3bd52e4c838afd78cb38b29cef1b28b3b9fbf69db75d6d4e0dfe7987ef3a7`
MD5	`8ee7cbf75fd12a9e3ea8c88a51e27ee2`
BLAKE2b-256	`00a1b344927463b7508c7cf6022f848db034443377445a1e74a59d174bd0d740`

Algorithm	Hash digest
SHA256	`de33e15f92818987ba8d9e0341f953318d86288efa3672079ddbdd33e9891601`
MD5	`be1a13e922a0221ffbf11a75c5307500`
BLAKE2b-256	`52c21bbbc9ba98a8be2c539888b289602eb92be808c2528585540088d483b880`

Algorithm	Hash digest
SHA256	`156cfcfe1f90e35b4c602d0b5d36bf00eacadb97916104e0d643a287438bfb3d`
MD5	`6cdb020c47b440694c167852e7166258`
BLAKE2b-256	`5d576f18d56e98cb073a987ffbaf1bd23eae52801b5e48b7ef135e41aab531ef`

Algorithm	Hash digest
SHA256	`57f51e783d799048118606c0a2c8deb7dbea3757cb5a2e2fcb836e3344989332`
MD5	`2e59b79bd65fcb0af5352b7c693fb3b5`
BLAKE2b-256	`59699fa60b2bf42a46bb9acb7b84855b1b167b5ae03676633dcaa9d1a3830df9`

Algorithm	Hash digest
SHA256	`007df5d8a5273084b828d19ad7525669e29aa4dc37f8118cacdedf5e33d74e58`
MD5	`5b3fbb223734bfef5134dc949bb8eaa6`
BLAKE2b-256	`b54874c47b72945e1c4b8743aaea8aa506a4c361bc9630e63282f789b2191f1c`

Algorithm	Hash digest
SHA256	`8b92c9d3fb84a7f14dcec79a2e2e92a7aee0ca721b8239868df62f9990a8420c`
MD5	`ca64fd94b4238f053062ab7dbe6bc898`
BLAKE2b-256	`71a7801d7e2019ff41e734a057c94a9f595b8c79e3857101b6e143a667b678e2`

Algorithm	Hash digest
SHA256	`f76fd588cca919525d634121b6b995f1d2a4590bac6d49fc65f5576cc851ed54`
MD5	`ad8a57ddf4e1f5e1d2084ae9f2e00e5b`
BLAKE2b-256	`8eb0a67e807e5b815cc6c611d7b53897e85db8c02a3805e6466e956074592087`

Algorithm	Hash digest
SHA256	`d271249cfeaafdcf41661dd27848b4fca8d035ccceb7c730bfa5e28a1b3b6b8d`
MD5	`39b9242d66f0876568c6c2be093064c5`
BLAKE2b-256	`93aef8d0aa27c6caf71e435c0d6cd068ebbcefbc83eda819341ebe8e982aba05`

spherical-geometry 1.3.0rc1

Navigation

Verified details

Maintainers

Unverified details

Project links

Meta

Classifiers

Project description

Installation

Requirements

Bundled qd-library

Coordinate representation

Spherical polygons

Representation

Where is the inside?

What is the orientation?

Disjoint Polygons

Creating spherical polygons

Operations on Spherical Polygons

Great circle arcs

Project details

Verified details

Maintainers

Unverified details

Project links

Meta

Classifiers

Release history Release notifications | RSS feed

Download files

Source Distribution

Built Distributions

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

Algorithm	Hash digest
SHA256	`cc9bbe3de73afed7dd792a96b2dcaef99ba165603f6ce50dcee16b508257d8c0`
MD5	`99b02fc41b654ee24cb3927d332c3742`
BLAKE2b-256	`f14a189c7aaff846a593e3a35d81baf216ea818cccb84e2efc7537e60b12a6c0`

Algorithm	Hash digest
SHA256	`7528fad365c04b9a08b8aff3dcbb742f289d6c9824075ec7b8e57f2343c74c30`
MD5	`e82116737fde3749d6fd810ee0664616`
BLAKE2b-256	`36766e0b3a0c36e76511cbbd3119e32c6cbde287f35591c5e2004ce2f4e08e5c`

Algorithm	Hash digest
SHA256	`c0d8a8c1df6d866f88ca7df31e7ee00db83cfe9c06fdea147fb261e05ed0064d`
MD5	`6ba0b45c3f2fc1a9407712b20df1cc09`
BLAKE2b-256	`7024f741fd1f2afad40dd9dc8bf255e2c1cfec06290cfe808ec0cc4fe686ce87`

Algorithm	Hash digest
SHA256	`c2b71850bb99180ada42a72ebaf18763b06df9ca976a1ed9a54d1361aeb10d84`
MD5	`fe64f437634268d501f86d736a112281`
BLAKE2b-256	`cc146958dfe67dbe45dbc14eaeae4add67a56a342e501da9dd44f2c749e2a8e3`