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

# 

# LSST Data Management System 

# Copyright 2016 LSST Corporation. 

# 

# This product includes software developed by the 

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

# 

# 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 LSST License Statement and 

# the GNU General Public License along with this program. If not, 

# see <http://www.lsstcorp.org/LegalNotices/>. 

# 

 

import astropy.units 

 

from lsst.afw.coord import Observatory, Weather 

from lsst.afw.geom import SpherePoint 

from lsst.obs.base import MakeRawVisitInfo 

 

__all__ = ["MakeMegacamRawVisitInfo"] 

 

 

class MakeMegacamRawVisitInfo(MakeRawVisitInfo): 

 

"""Make a VisitInfo from the FITS header of a raw Megacam image 

""" 

 

def setArgDict(self, md, argDict): 

"""Set an argument dict for VisitInfo and pop associated metadata 

 

Parameters 

---------- 

md : `lsst.daf.base.PropertyList` or `lsst.daf.base.PropertySet` 

Metadata to extract from. Extracted values are removed. 

argdict : `dict` 

A dict of arguments to add to values to. 

""" 

MakeRawVisitInfo.setArgDict(self, md, argDict) 

argDict["darkTime"] = self.popFloat(md, "DARKTIME") 

argDict["boresightAzAlt"] = SpherePoint( 

self.popAngle(md, "TELAZ"), 

self.popAngle(md, "TELALT"), 

) 

argDict["boresightRaDec"] = SpherePoint( 

self.popAngle(md, "RA_DEG",), 

self.popAngle(md, "DEC_DEG"), 

) 

argDict["boresightAirmass"] = self.popFloat(md, "AIRMASS") 

argDict["observatory"] = Observatory( 

self.popAngle(md, "LONGITUD"), 

self.popAngle(md, "LATITUDE"), 

4204, # from Wikipedia 

) 

argDict["weather"] = Weather( 

self.popFloat(md, "TEMPERAT"), 

self.popFloat(md, "PRESSURE")*100.0, # 100 Pascal per millibar 

self.popFloat(md, "RELHUMID"), 

) 

# Using LST to compute ERA until we get UT1 (see: DM-8053) 

LST = self.popAngle(md, "LST-OBS", units=astropy.units.h) 

argDict['era'] = self.eraFromLstAndLongitude(LST, argDict["observatory"].getLongitude()) 

 

def getDateAvg(self, md, exposureTime): 

"""Return date at the middle of the exposure 

 

Parameters 

---------- 

md : `lsst.daf.base.PropertyList` or `lsst.daf.base.PropertySet` 

FITS metadata; changed in place 

exposureTime : `float` 

exposure time in sec 

""" 

dateObs = self.popMjdDate(md, "MJD-OBS") 

return self.offsetDate(dateObs, 0.5*exposureTime)