1671

pending completion

Build # 1671

Build Type

Pull #253

travis-ci

Committed by

web-flow

Commit Message

<a href="https://github.com/LCOGT/banzai/commit/<a class=hub.com/LCOGT/banzai/commit/0371f6d1063dce82c5fe514896440fd254b257e5">0371f6d10<a href="https://github.com/LCOGT/banzai/commit/0371f6d1063dce82c5fe514896440fd254b257e5">">Merge </a><a class="double-link" href="https://github.com/LCOGT/banzai/commit/<a class="double-link" href="https://github.com/LCOGT/banzai/commit/f52e4d0615d885e64fdee0c6aab5281dbbb2f025">f52e4d061</a>">f52e4d061</a><a href="https://github.com/LCOGT/banzai/commit/0371f6d1063dce82c5fe514896440fd254b257e5"> into 4ed7aa072">4ed7aa072</a>

Pull Request Pull Request #253: Merge dev into master

Run Details

1380 of 1380 new or added lines in 31 files covered. (100.0%)

1772 of 2819 relevant lines covered (62.86%)

0.63 hits per line

Source File
Press 'n' to go to next uncovered line, 'b' for previous

75.76

/banzai/qc/pointing.py

import logging

import astropy.units as u
from astropy.coordinates import SkyCoord

from banzai.stages import Stage
from banzai.utils import qc

logger = logging.getLogger('banzai')


class PointingTest(Stage):
    """
    A test to determine  whether or not the poiting error on the frame
    (as determined by a WCS solve) is within tolerance.
    """

    # Typical pointing is within 5" of requested pointing (decimal degrees).
    WARNING_THRESHOLD = 30.0
    SEVERE_THRESHOLD = 300.0

    def __init__(self, runtime_context):
        super(PointingTest, self).__init__(runtime_context)

    def do_stage(self, image):
        try:
            # OFST-RA/DEC is the same as CAT-RA/DEC but includes user requested offset
            requested_coords = SkyCoord(image.meta['OFST-RA'], image.meta['OFST-DEC'],
                                        unit=(u.hour, u.deg), frame='icrs')
        except ValueError as e:
            try:
                # Fallback to CAT-RA and CAT-DEC
                requested_coords = SkyCoord(image.meta['CAT-RA'], image.meta['CAT-DEC'],
                                            unit=(u.hour, u.deg), frame='icrs')
            except:
                logger.error(e, image=image)
                return image
        # This only works assuming CRPIX is at the center of the image
        solved_coords = SkyCoord(image.meta['CRVAL1'], image.meta['CRVAL2'],
                                 unit=(u.deg, u.deg), frame='icrs')

        angular_separation = solved_coords.separation(requested_coords).arcsec

        logging_tags = {'PNTOFST': angular_separation}

        pointing_severe = abs(angular_separation) > self.SEVERE_THRESHOLD
        pointing_warning = abs(angular_separation) > self.WARNING_THRESHOLD
        if pointing_severe:
            logger.error('Pointing offset exceeds threshold', image=image, extra_tags=logging_tags)
        elif pointing_warning:
            logger.warning('Pointing offset exceeds threshhold', image=image, extra_tags=logging_tags)
        qc_results = {'pointing.failed': pointing_severe,
                      'pointing.failed_threshold': self.SEVERE_THRESHOLD,
                      'pointing.warning': pointing_warning,
                      'pointing.warning_threshold': self.WARNING_THRESHOLD,
                      'pointing.offset': angular_separation}
        qc.save_qc_results(self.runtime_context, qc_results, image)

        image.meta['PNTOFST'] = (angular_separation, '[arcsec] offset of requested and solved center')

        return image

1	import logging	1✔
2
3	import astropy.units as u	1✔
4	from astropy.coordinates import SkyCoord	1✔
5
6	from banzai.stages import Stage	1✔
7	from banzai.utils import qc	1✔
8
9	logger = logging.getLogger('banzai')	1✔
10
11
12	class PointingTest(Stage):	1✔
13	"""
14	A test to determine whether or not the poiting error on the frame
15	(as determined by a WCS solve) is within tolerance.
16	"""
17
18	# Typical pointing is within 5" of requested pointing (decimal degrees).
19	WARNING_THRESHOLD = 30.0	1✔
20	SEVERE_THRESHOLD = 300.0	1✔
21
22	def __init__(self, runtime_context):	1✔
23	super(PointingTest, self).__init__(runtime_context)	1✔
24
25	def do_stage(self, image):	1✔
26	try:	1✔
27	# OFST-RA/DEC is the same as CAT-RA/DEC but includes user requested offset
28	requested_coords = SkyCoord(image.meta['OFST-RA'], image.meta['OFST-DEC'],	1✔
29	unit=(u.hour, u.deg), frame='icrs')
30	except ValueError as e:	×
31	try:	×
32	# Fallback to CAT-RA and CAT-DEC
33	requested_coords = SkyCoord(image.meta['CAT-RA'], image.meta['CAT-DEC'],	×
34	unit=(u.hour, u.deg), frame='icrs')
35	except:	×
36	logger.error(e, image=image)	×
37	return image	×
38	# This only works assuming CRPIX is at the center of the image
39	solved_coords = SkyCoord(image.meta['CRVAL1'], image.meta['CRVAL2'],	1✔
40	unit=(u.deg, u.deg), frame='icrs')
41
42	angular_separation = solved_coords.separation(requested_coords).arcsec	1✔
43
44	logging_tags = {'PNTOFST': angular_separation}	1✔
45
46	pointing_severe = abs(angular_separation) > self.SEVERE_THRESHOLD	1✔
47	pointing_warning = abs(angular_separation) > self.WARNING_THRESHOLD	1✔
48	if pointing_severe:	1✔
49	logger.error('Pointing offset exceeds threshold', image=image, extra_tags=logging_tags)	×
50	elif pointing_warning:	1✔
51	logger.warning('Pointing offset exceeds threshhold', image=image, extra_tags=logging_tags)	×
52	qc_results = {'pointing.failed': pointing_severe,	1✔
53	'pointing.failed_threshold': self.SEVERE_THRESHOLD,
54	'pointing.warning': pointing_warning,
55	'pointing.warning_threshold': self.WARNING_THRESHOLD,
56	'pointing.offset': angular_separation}
57	qc.save_qc_results(self.runtime_context, qc_results, image)	1✔
58
59	image.meta['PNTOFST'] = (angular_separation, '[arcsec] offset of requested and solved center')	1✔
60
61	return image	1✔

LCOGT / banzai / 1671

Source File Press 'n' to go to next uncovered line, 'b' for previous

Source File
Press 'n' to go to next uncovered line, 'b' for previous