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

147

148

149

150

151

152

153

154

# This file is part of obs_lsst. 

# 

# 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/>. 

 

"""Gen3 Butler Formatters for LSST raw data. 

""" 

 

__all__ = ("LsstCamRawFormatter", "LatissRawFormatter") 

 

from astro_metadata_translator import fix_header, merge_headers 

 

import lsst.afw.fits 

from lsst.obs.base.fitsRawFormatterBase import FitsRawFormatterBase 

 

from .instrument import LsstCamInstrument, LatissInstrument, \ 

ImsimInstrument, PhosimInstrument, Ts8Instrument, \ 

Ts3Instrument, UcdCamInstrument, LsstComCamInstrument 

from ..translators import LsstLatissTranslator, LsstCamTranslator, \ 

LsstUCDCamTranslator, LsstTS3Translator, LsstComCamTranslator, \ 

PhosimTranslator, LsstTS8Translator, ImsimTranslator 

from ..assembly import fixAmpsAndAssemble, readRawAmps 

 

 

class LsstCamRawFormatter(FitsRawFormatterBase): 

translatorClass = LsstCamTranslator 

filterDefinitions = LsstCamInstrument.filterDefinitions 

_instrument = LsstCamInstrument 

 

def readMetadata(self): 

"""Read all header metadata directly into a PropertyList. 

 

Specialist version since some of our data does not 

set INHERIT=T so we have to merge the headers manually. 

 

Returns 

------- 

metadata : `~lsst.daf.base.PropertyList` 

Header metadata. 

""" 

file = self.fileDescriptor.location.path 

phdu = lsst.afw.fits.readMetadata(file, 0) 

59 ↛ 61line 59 didn't jump to line 61, because the condition on line 59 was never true if "INHERIT" in phdu: 

# Trust the inheritance flag 

return super().readMetadata() 

 

# Merge ourselves 

md = merge_headers([phdu, lsst.afw.fits.readMetadata(file)], 

mode="overwrite") 

fix_header(md) 

return md 

 

def getDetector(self, id): 

return self._instrument.getCamera()[id] 

 

def readImage(self): 

"""Read just the image component of the Exposure. 

 

Returns 

------- 

image : `~lsst.afw.image.Image` 

In-memory image component. 

""" 

return self.readFull().getImage() 

 

def readFull(self, parameters=None): 

"""Read the complete exposure. 

 

This correctly fixes amplifier bounding box deviations from 

the camera definitions, and so should provide the safest 

interface to the data. 

 

Parameters 

---------- 

parameters : `dict`, optional 

No parameters are currently used. 

 

Returns 

------- 

exposure : `~lsst.afw.image.Exposure` 

Complete in-memory exposure representation. 

 

""" 

rawFile = self.fileDescriptor.location.path 

ccd = self.getDetector(self.observationInfo.detector_num) 

ampExps = readRawAmps(rawFile, ccd) 

exposure = fixAmpsAndAssemble(ampExps, rawFile) 

 

mask = self.readMask() 

106 ↛ 107line 106 didn't jump to line 107, because the condition on line 106 was never true if mask is not None: 

exposure.setMask(mask) 

variance = self.readVariance() 

109 ↛ 110line 109 didn't jump to line 110, because the condition on line 109 was never true if variance is not None: 

exposure.setVariance(variance) 

 

info = exposure.getInfo() 

info.setFilter(self.makeFilter()) 

info.setVisitInfo(self.makeVisitInfo()) 

info.setWcs(self.makeWcs(info.getVisitInfo(), info.getDetector())) 

 

exposure.setMetadata(self.metadata) 

 

return exposure 

 

 

class LatissRawFormatter(LsstCamRawFormatter): 

translatorClass = LsstLatissTranslator 

_instrument = LatissInstrument 

 

 

class ImsimRawFormatter(LsstCamRawFormatter): 

translatorClass = ImsimTranslator 

_instrument = ImsimInstrument 

 

 

class PhosimRawFormatter(LsstCamRawFormatter): 

translatorClass = PhosimTranslator 

_instrument = PhosimInstrument 

 

 

class Ts8RawFormatter(LsstCamRawFormatter): 

translatorClass = LsstTS8Translator 

_instrument = Ts8Instrument 

 

 

class Ts3RawFormatter(LsstCamRawFormatter): 

translatorClass = LsstTS3Translator 

_instrument = Ts3Instrument 

 

 

class LsstComCamRawFormatter(LsstCamRawFormatter): 

translatorClass = LsstComCamTranslator 

_instrument = LsstComCamInstrument 

 

 

class UcdCamRawFormatter(LsstCamRawFormatter): 

translatorClass = LsstUCDCamTranslator 

_instrument = UcdCamInstrument