Source code for pydl.photoop.photoobj

# Licensed under a 3-clause BSD style license - see LICENSE.rst
# -*- coding: utf-8 -*-
"""This module corresponds to the photoobj directory in photoop.
"""
import numpy as np
import astropy.units as u


[docs] def sdss_calibv(): """Return calibration for velocities from pix/frame to deg/day. Returns ------- :class:`~astropy.units.quanity.Quantity` The conversion from pixels per frame to degrees per day Notes ----- Assumes frame time difference of 71.72 seconds and pixel scale of 0.396 arcsec, both fixed. Also note that observations of the same part of sky from adjacent bands are separated by *two* frame numbers, so we multiply by a factor two. """ pixscale = 0.396 * u.arcsec ftime = 71.72 * u.s return (2.0*pixscale/ftime).to('deg / d')
[docs] def unwrap_objid(objid): """Unwrap CAS-style objID into run, camcol, field, id, rerun. See :func:`~pydl.pydlutils.sdss.sdss_objid` for details on how the bits within an objID are assigned. Parameters ---------- objid : :class:`numpy.ndarray` An array containing 64-bit integers or strings. If strings are passed, they will be converted to integers internally. Returns ------- :class:`numpy.recarray` A record array with the same length as `objid`, with the columns 'skyversion', 'rerun', 'run', 'camcol', 'firstfield', 'frame', 'id'. Raises ------ :exc:`ValueError` If the input objID has a type that can't be converted into 64-bit integer. Notes ----- For historical reasons, the inverse of this function, :func:`~pydl.pydlutils.sdss.sdss_objid` is not in the same namespace as this function. 'frame' is used instead of 'field' because record arrays have a method of the same name. Examples -------- >>> from numpy import array >>> from pydl.photoop.photoobj import unwrap_objid >>> unwrap_objid(array([1237661382772195474])) rec.array([(2, 301, 3704, 3, 0, 91, 146)], dtype=[('skyversion', '<i4'), ('rerun', '<i4'), ('run', '<i4'), ('camcol', '<i4'), ('firstfield', '<i4'), ('frame', '<i4'), ('id', '<i4')]) """ try: np_string = np.string_ # Numpy 1.x np_unicode = np.unicode_ except AttributeError: np_string = np.bytes_ # Numpy 2.x np_unicode = np.str_ if objid.dtype.type is np_string or objid.dtype.type is np_unicode: tempobjid = objid.astype(np.int64) elif objid.dtype.type is np.int64: tempobjid = objid.copy() else: raise ValueError('Unrecognized type for objid!') unwrap = np.recarray(objid.shape, dtype=[('skyversion', 'i4'), ('rerun', 'i4'), ('run', 'i4'), ('camcol', 'i4'), ('firstfield', 'i4'), ('frame', 'i4'), ('id', 'i4')]) unwrap.skyversion = np.bitwise_and(tempobjid >> 59, 2**4 - 1) unwrap.rerun = np.bitwise_and(tempobjid >> 48, 2**11 - 1) unwrap.run = np.bitwise_and(tempobjid >> 32, 2**16 - 1) unwrap.camcol = np.bitwise_and(tempobjid >> 29, 2**3 - 1) unwrap.firstfield = np.bitwise_and(tempobjid >> 28, 2**1 - 1) unwrap.frame = np.bitwise_and(tempobjid >> 16, 2**12 - 1) unwrap.id = np.bitwise_and(tempobjid, 2**16 - 1) return unwrap