@@ -27,6 +27,13 @@ def estimate_exposure_time(timestamps):
27
27
28
28
@u .quantity_input (ra = u .hourangle , dec = u .deg , fov = u .deg )
29
29
def build_exposure_regions (pointing_coords , fov = 4.5 * u .deg ):
30
+ '''
31
+ Takes a list of pointing positions and a field of view and returns
32
+ the unique pointing positions and the astropy.regions.
33
+
34
+ For an observation with N wobble positions this will return N unique
35
+ pointing positions and N circular regions.
36
+ '''
30
37
unique_pointing_positions = SkyCoord (
31
38
ra = np .unique (pointing_coords .ra ),
32
39
dec = np .unique (pointing_coords .dec )
@@ -57,7 +64,11 @@ def _build_standard_wcs(image_center, shape, naxis=2, fov=9 * u.deg):
57
64
58
65
@u .quantity_input (event_ra = u .hourangle , event_dec = u .deg , fov = u .deg )
59
66
def build_exposure_map (pointing_coords , event_time , fov = 4.5 * u .deg , wcs = None , shape = (1000 , 1000 )):
60
-
67
+ '''
68
+ Takes pointing coordinates for each event and the corresponding timestamp.
69
+ Returns a masked array containing the estimated exposure time in hours
70
+ and a WCS object so the mask can be plotted.
71
+ '''
61
72
if not wcs :
62
73
image_center = SkyCoord (ra = pointing_coords .ra .mean (), dec = pointing_coords .dec .mean ())
63
74
wcs = _build_standard_wcs (image_center , shape , fov = 2 * fov )
0 commit comments