software_tools/calc_expected_visits.py at main · LSST-TVSSC/software_tools · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
# -*- coding: utf-8 -*-
"""
Created on Tue Sep 25 17:11:20 2018

@author: rstreet
"""

import numpy as np
import calculate_lsst_field_visibility_astropy

def calc_expected_visits(pointings,cadence,start_date,end_date,verbose=False):
    """Function to calculate the maximum possible number of visits to a
    given pointing, given the expected cadence of observation and within
    the date ranges given, taking target visibility into account.

    Input:
    :param array ra:            RAs, J2000.0, sexigesimal format
    :param array dec:           Decs, J2000.0, sexigesimal format
    :param float cadence:       Interval between successive visits in the
                                same single filter in hours
    :param string start_date:   Start of observing window YYYY-MM-DD
    :param string start_date:   End of observation window YYYY-MM-DD

    Output:
    :param list of arrays n_visits:       Number of visits possible per night
                                          for each pointing
    :param list of arrays hrs_visibility: Hours of visibility per night
                                          for each pointing
    """

    n_visits = []
    hrs_visibility = []

    if verbose:
        print('Calculating visbility for '+str(len(pointings))+' fields')

    for i in range(0,len(pointings),1):

        (ra, dec) = pointings[i]

        if verbose:
            print(' -> RA '+str(ra)+', Dec '+str(dec))

        (total_time_visible, hrs_visible_per_night) = calculate_lsst_field_visibility_astropy.calculate_lsst_field_visibility(ra,dec,start_date,end_date,verbose=False)

        n_visits.append( (np.array(hrs_visible_per_night) / cadence).astype(int) )
        hrs_visibility.append( np.array(hrs_visible_per_night) )

    return n_visits,hrs_visibility