phn

PyFITS: FITS files in Python

In this article, we provide examples of using the python module PyFITS for working with FITS data. We first go through a brief

phn

import jdcal
from jdcal import ipart, fpart, is_leap, gcal2jd, jd2gcal

CJ = 36525.0


def days_in_year(year, c="g"):
    """Number of days in a calendar year (Gregorian/Julian)."""
    if c.lower() == "g":
        return 366 if is_leap(year) else 365

phn

"""Conversion between 6D Cartesian and spherical coordinates.

Cartesian6D represents the coordinates (x, y, z, xdot, ydot, zdot).
Spherical6D represents the coordinates (r, alpha, delta, rdot,
alphadot, deltadot). Alpha is the longitudinal angle and delta is the
latitudinal angle. The latter goes from 90 degrees at the +ve z-axis to
-90 at the -ve z-axis.

The Cartesian coordinates can have any units. The radial spherical
coordinate has the some unit as the input Cartesian

phn

import matplotlib as mpl
import lineid_plot


def color_text_boxes(ax, labels, colors, color_arrow=True):
    assert len(labels) == len(colors), \
        "Equal no. of colors and lables must be given"
    boxes = ax.findobj(mpl.text.Annotation)
    box_lables = lineid_plot.unique_labels(labels)
    for box in boxes:

phn

# Code moved to http://github.com/phn/jdcal.

phn

# Code moved to http://github.com/phn/angles.