Skip to main content
Join the official 2019 Python Developers SurveyStart the survey!

minor planet checking

Project description

pympc

perform checks for the presence of minor bodies at astronomical locations for a given epoch.

installation

pip install .

or

python setup.py install

examples

define where we want the orbit catalogue to be downloaded to (see next section)

$ export MPCORB_JSON_PATH = /my/file/path/mpcorb_extended.json

import the package and (optionally) grab the catalogue - this must be done at least once prior to doing any searches and can be run to overwrite the catalogue with the latest version from the Minor Planet Center

>>> import pympc
>>> pympc.update_catalogue()
example 1

define our search location, epoch and radius and run the check

>>> import astropy.units as u
>>> from astropy.time import Time
>>> ra = 230.028 * u.deg
>>> dec = -11.774 * u.deg
>>> epoch = Time("2019-01-01T00:00")
>>> search_radius = 0.5 * u.arcmin
>>> max_mag = 20
>>> pympc.minor_planet_check(ra, dec, epoch, search_radius, max_mag)
example 2

here we use float arguments, and the program assumes the units (see pympc.minor_planet_check() docstring)

>>> ra = 230.028  # assumed degrees
>>> dec = -11.774  # assumed degrees
>>> epoch = 58484.  # assumed MJD
>>> search_radius = 30  # assumed arcseconds
>>> max_mag = 20
>>> pympc.minor_planet_check(ra, dec, epoch, search_radius, max_mag)

setting the orbit catalogue location

the location to download the MPCORB catalogue is set via the environment variable MPCORB_JSON_PATH (e.g. export MPCORB_JSON_PATH = /my/file/path/mpcorb.json). If this is not found the catalogue will download to the default location of /tmp/mpcorb_extended.json.

speed and multiprocessing

the check should take of order a second or two, depending on multiprocessing capabilities.

the private function which actually performs the calculation is _minor_planet_check() (note leading underscore). this can be called to avoid the overhead associated with converting input arguments to minor_planet_check(), if you can provide them directly as required (see _minor_planet_check() docstring).

by default the program will calculate positions of bodies in the catalogue multiprocessed. to switch this off set chunk_size = 0, i.e.:

>>> pympc.minor_planet_check(ra, dec, epoch, search_radius, chunk_size=0)

limitations

the orbits are propagated by xephem, and that package does not account for perturbations of the orbits. thus the accuracy of the position is dependant on the time difference between the epoch of the orbit elements and the epoch at which the search is being performed. epoch differences of around a month or less should be fine for typical positional accuracies of ground-based observations (~arcsecond), but time differences of several months can accrue to inaccuracies of around ~arcminute for certain bodies.

currently the epoch of the orbit elements is visible in the xephem db strings returned by minor_planet_check() as a decimal year format (e.g. ..,2019.317808,..). some diagnostic information and warning when using large time differences is to be implemented.

the filtering of matches based on magnitude via max_mag is limited by the accuracy of the magnitude information in the database so some buffer should be applied to the desired magnitude cutoff to allow for this.

acknowledgments

this package makes use of orbit elements provided by the International Astronomical Union's Minor Planet Center.

based from a package developed by Chris Klein and Duncan Galloway

Project details


Release history Release notifications

This version

0.2

Download files

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

Files for pympc, version 0.2
Filename, size File type Python version Upload date Hashes
Filename, size pympc-0.2.tar.gz (5.5 kB) File type Source Python version None Upload date Hashes View hashes

Supported by

Elastic Elastic Search Pingdom Pingdom Monitoring Google Google BigQuery Sentry Sentry Error logging AWS AWS Cloud computing DataDog DataDog Monitoring Fastly Fastly CDN SignalFx SignalFx Supporter DigiCert DigiCert EV certificate StatusPage StatusPage Status page