Astrometry Database Interface
This module provides an interface for accessing a database which contains updated WCS solutions for many images in the HST archive. The exact set of images with new solutions in this database will change over time, with the initial release containing updates only for ACS and WFC3 direct images. This database will be designed to provide multiple WCS solutions for each exposure depending on what astrometric frame was used to define the new WCS solution; for example, an image may be matched to Pan-STARRS sources in one solution and to GAIA sources in another solution.
Each new WCS solution will be appended to the exposure by this module as a headerlet extension, taking advantage of the headerlet definition provided by the STWCS package, found in Headerlet documentation. The headerlet library allows a user to manage all appended solutions and determine which one gets used by the exposure at anytime, since only 1 WCS can be used at a time.
Usage
This interface defines an object that provides the services for updating an exposure with new WCS solutions. Only 2 calls are typically necessary to update any exposure or set of exposures.
0. Import the package
>>> from stwcs.updatewcs import astrometry_utils as stwcsau
1. Establish a working connection with the database
>>> db = stwcsau.AstrometryDB()
2. For each exposure, query the database and update with any new WCS solutions
>>> db.updateObs(filename)
The environment variables used by this module can be defined by the user to change what database is being accessed, whether to raise exceptions instead of warnings when problems arise in connecting to the database or otherwise updating an exposure, or to completely turn on/off use of this module. These allow external code, such as updatewcs (see updatewcs), to include this code and control it as necessary.
API
Interface for Astrometry database service.
This module contains interface functions for the AstrometryDB Restful service based on interfaces/code provided by B. McLean 11-Oct-2017.
The code checks for the existence of environmental variables to help control the operation of this interface; namely,
- RAISE_PIPELINE_ERRORS - boolean to specify whether to raise exceptions
during processing or simply log errors and quit gracefully. If not set, default behavior will be to log errors and quit gracefully.
- ASTROMETRY_SERVICE_URL - URL pointing to user-specified web service that
will provide updated astrometry solutions for observations being processed. This will replace the built-in URL included in the base class. This value will also be replaced by any URL provided by the user as an input parameter
url
.- ASTROMETRY_STEP_CONTROL - String specifying whether or not to perform the
astrometry update processing at all. Valid Values: “ON”, “On”, “on”, “OFF”, “Off”, “off” If not set, default value is “ON”.
- GSSS_WEBSERVICES_URL - URL point to user-specified web service which provides
information on the guide stars used for taking HST observations. This value will replace the default URL included in this module as the
gsss_url
variable.
- class stwcs.updatewcs.astrometry_utils.AstrometryDB(url=None, raise_errors=None, perform_step=True, write_log=False)
Bases:
object
Base class for astrometry database interface.
- apply_new_apriori(obsname)
Compute and apply a new a priori WCS based on offsets from astrometry database.
- Parameters:
obsname (str) – Full filename or
astropy.io.fits.HDUList
object for the observation to be corrected- Returns:
wcsname – Value of WCSNAME keyword for this new WCS
- Return type:
- findObservation(observationID)
Find whether there are any entries in the AstrometryDB for the observation with
observationID
.- Parameters:
observationID (str) – base rootname for observation to be updated (eg.,
iab001a1q
)- Returns:
entry – Database entry for this observation, if found. It will return None if there was an error in accessing the database and
self.raise_errors
was not set to True.- Return type:
obj
- getObservation(observationID)
Get solutions for observation from AstrometryDB.
- Parameters:
observationID (str) – base rootname for observation to be updated (eg.,
iab001a1q
)- Returns:
headerlets (dict) – Dictionary containing all solutions found for exposure in the form of headerlets labelled by the name given to the solution in the database.
best_solution_id (str) – WCSNAME of the WCS solution flagged as ‘best’ in the astrometry database for the observation. The ‘best’ solution should be the one that aligns the observation as close to the GAIA frame as possible.
- isAvailable()
Test availability of astrometryDB web-service.
- updateObs(obsname, all_wcs=False, remove_duplicates=True)
Update observation with any available solutions.
- Parameters:
obsname (str) – Filename for observation to be updated
all_wcs (bool) – If True, all solutions from the Astrometry database are appended to the input file as separate FITS extensions. If False, only those solutions based on the same IDCTAB will be appended.
remove_duplicates (bool) – If True, any headerlet extensions with the same HDRNAME are found, the copies will be deleted until only the first version added remains.
- stwcs.updatewcs.astrometry_utils.apply_astrometric_updates(obsnames, **pars)
Apply new astrometric solutions to observation.
Functional stand-alone interface for applying new astrometric solutions found in the astrometry dB to the given observation(s).
- Parameters:
obsnames (str, list) – Filename or list of filenames of observation(s) to be updated
url (str, optional) – URL of astrometry database web-interface to use. If None (default), it will use built-in URL for STScI web-interface
raise_errors (bool, optional) – Specify whether or not to raise Exceptions and stop processing when an error in either accessing the database, retrieving a solution from the database or applying the new solution to the observation. If None, it will look to see whether the environmental variable
RAISE_PIPELINE_ERRORS
was set, otherwise, it will default to ‘False’.
Version
- stwcs.__version__ = '1.7.4.dev28+g9b8269b.d20240918'
str(object=’’) -> str str(bytes_or_buffer[, encoding[, errors]]) -> str
Create a new string object from the given object. If encoding or errors is specified, then the object must expose a data buffer that will be decoded using the given encoding and error handler. Otherwise, returns the result of object.__str__() (if defined) or repr(object). encoding defaults to sys.getdefaultencoding(). errors defaults to ‘strict’.