Hide keyboard shortcuts

Hot-keys on this page

r m x p   toggle line displays

j k   next/prev highlighted chunk

0   (zero) top of page

1   (one) first highlighted chunk

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

50

51

52

53

54

55

56

57

58

59

60

61

62

63

64

65

66

67

68

69

70

71

72

73

74

75

76

77

78

79

80

81

82

83

84

85

86

87

88

89

90

91

92

93

94

95

96

97

98

99

100

101

102

103

104

105

106

107

108

109

110

111

112

113

114

115

116

117

118

119

120

121

122

123

124

125

126

127

128

129

130

131

132

133

134

135

136

137

138

139

140

141

142

143

144

145

146

# This file is part of daf_butler. 

# 

# Developed for the LSST Data Management System. 

# This product includes software developed by the LSST Project 

# (http://www.lsst.org). 

# See the COPYRIGHT file at the top-level directory of this distribution 

# for details of code ownership. 

# 

# This program is free software: you can redistribute it and/or modify 

# it under the terms of the GNU General Public License as published by 

# the Free Software Foundation, either version 3 of the License, or 

# (at your option) any later version. 

# 

# This program is distributed in the hope that it will be useful, 

# but WITHOUT ANY WARRANTY; without even the implied warranty of 

# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 

# GNU General Public License for more details. 

# 

# You should have received a copy of the GNU General Public License 

# along with this program. If not, see <http://www.gnu.org/licenses/>. 

 

__all__ = ("Instrument", "makeExposureEntryFromVisitInfo", "makeVisitEntryFromVisitInfo") 

 

from lsst.daf.base import DateTime 

 

 

# TODO: all code in this module probably needs to be moved to a higher-level 

# package (perhaps obs_base), but it is needed by the gen2convert subpackage. 

# We should probably move that as well. 

 

 

class Instrument: 

"""A template method class that can register itself with a 

`Registry. 

 

This class should be subclassed by various implementations. 

Subclasses should provide all relevant attributes, as documented 

below. 

 

Attributes 

---------- 

camera : `str` 

Name of the camera. Must be provided by subclass. 

physicalFilters : `list` 

List of PhysicalFilter entries (each entry being a dict). 

sensors : `list` 

List of Sensor entries (each entry being a dict). 

""" 

camera = None 

physicalFilters = [] 

sensors = [] 

 

def register(self, registry): 

"""Register an instance of this `Instrument` with a `Registry`. 

 

Creates all relevant `DataUnit` entries. 

""" 

assert self.camera is not None 

self._addCamera(registry) 

self._addPhysicalFilters(registry) 

self._addSensors(registry) 

 

def _addCamera(self, registry): 

registry.addDataUnitEntry("Camera", {"camera": self.camera}) 

 

def _addPhysicalFilters(self, registry): 

for entry in self.physicalFilters: 

68 ↛ 70line 68 didn't jump to line 70, because the condition on line 68 was never false if "camera" not in entry: 

entry["camera"] = self.camera 

registry.addDataUnitEntry("PhysicalFilter", entry) 

 

def _addSensors(self, registry): 

for entry in self.sensors: 

74 ↛ 76line 74 didn't jump to line 76, because the condition on line 74 was never false if 'camera' not in entry: 

entry['camera'] = self.camera 

registry.addDataUnitEntry('Sensor', entry) 

 

 

def makeExposureEntryFromVisitInfo(dataId, visitInfo, snap=0): 

"""Construct an Exposure DataUnit entry from `afw.image.VisitInfo`. 

 

Parameters 

---------- 

dataId : `dict` 

Dictionary of DataUnit primary/foreign key values for Exposure 

("camera", "exposure", optionally "visit" and "physical_filter"). 

visitInfo : `afw.image.VisitInfo` 

A `VisitInfo` object corresponding to the Exposure. 

snap : `int` 

Snap index of the Exposure. 

 

Returns 

------- 

entry : `dict` 

A dictionary containing all fields in the Exposure table. 

""" 

avg = visitInfo.getDate() 

begin = DateTime(int(avg.nsecs(DateTime.TAI) - 0.5E9*visitInfo.getExposureTime()), DateTime.TAI) 

result = { 

"datetime_begin": begin.toPython(), 

"exposure_time": visitInfo.getExposureTime(), 

"boresight_az": visitInfo.getBoresightAzAlt().getLongitude().asDegrees(), 

"boresight_alt": visitInfo.getBoresightAzAlt().getLatitude().asDegrees(), 

"rot_angle": visitInfo.getBoresightRotAngle().asDegrees(), 

"snap": snap, 

"dark_time": visitInfo.getDarkTime() 

} 

result.update(dataId) 

return result 

 

 

def makeVisitEntryFromVisitInfo(dataId, visitInfo): 

"""Construct a Visit DataUnit entry from `afw.image.VisitInfo`. 

 

Parameters 

---------- 

dataId : `dict` 

Dictionary of DataUnit primary/foreign key values for Visit ("camera", 

"visit", optionally "physical_filter"). 

visitInfo : `afw.image.VisitInfo` 

A `VisitInfo` object corresponding to the Visit. 

 

Returns 

------- 

entry : `dict` 

A dictionary containing all fields in the Visit table aside from 

"region". 

""" 

avg = visitInfo.getDate() 

begin = DateTime(int(avg.nsecs(DateTime.TAI) - 0.5E9*visitInfo.getExposureTime()), DateTime.TAI) 

end = DateTime(int(avg.nsecs(DateTime.TAI) + 0.5E9*visitInfo.getExposureTime()), DateTime.TAI) 

result = { 

"datetime_begin": begin.toPython(), 

"datetime_end": end.toPython(), 

"exposure_time": visitInfo.getExposureTime(), 

"boresight_az": visitInfo.getBoresightAzAlt().getLongitude().asDegrees(), 

"boresight_alt": visitInfo.getBoresightAzAlt().getLatitude().asDegrees(), 

"rot_angle": visitInfo.getBoresightRotAngle().asDegrees(), 

"earth_rotation_angle": visitInfo.getEra().asDegrees(), 

"boresight_ra": visitInfo.getBoresightRaDec().getLongitude().asDegrees(), 

"boresight_dec": visitInfo.getBoresightRaDec().getLatitude().asDegrees(), 

"boresight_parallactic_angle": visitInfo.getBoresightParAngle().asDegrees(), 

"local_era": visitInfo.getLocalEra().asDegrees(), 

} 

result.update(dataId) 

return result