def _saveImage(self, imageRequest, imageData, extra):
(mode, binning, top, left, width,
height) = self._getReadoutModeInfo(imageRequest["binning"],
imageRequest["window"])
binFactor = extra.get("binning_factor", 1.0)
pix_w, pix_h = self.getPixelSize()
focal_length = self["telescope_focal_length"]
scale_x = binFactor * (((180 / pi) / focal_length) * (pix_w * 0.001))
scale_y = binFactor * (((180 / pi) / focal_length) * (pix_h * 0.001))
full_width, full_height = self.getPhysicalSize()
CRPIX1 = ((int(full_width / 2.0)) - left) - 1
CRPIX2 = ((int(full_height / 2.0)) - top) - 1
t0 = time.time()
img = Image.create(imageData, imageRequest)
img += [
('DATE-OBS',
ImageUtil.formatDate(
extra.get("frame_start_time",
dt.datetime.utcnow())), 'Date exposure started'),
('CCD-TEMP', extra.get("frame_temperature", -275.0),
'CCD Temperature at Exposure Start [deg. C]'),
("EXPTIME", float(imageRequest['exptime'])
or -1, "exposure time in seconds"),
('IMAGETYP', imageRequest['type'].strip(), 'Image type'),
('SHUTTER', str(imageRequest['shutter']),
'Requested shutter state'),
("CRPIX1", CRPIX1, "coordinate system reference pixel"),
("CRPIX2", CRPIX2, "coordinate system reference pixel"),
("CD1_1", scale_x, "transformation matrix element (1,1)"),
("CD1_2", 0.0, "transformation matrix element (1,2)"),
("CD2_1", 0.0, "transformation matrix element (2,1)"),
("CD2_2", scale_y, "transformation matrix element (2,2)"),
('CAMERA', str(self['camera_model']), 'Camera Model'),
('CCD', str(self['ccd_model']), 'CCD Model'),
('CCD_DIMX', self.getPhysicalSize()[0], 'CCD X Dimension Size'),
('CCD_DIMY', self.getPhysicalSize()[1], 'CCD Y Dimension Size'),
('CCDPXSZX', self.getPixelSize()[0],
'CCD X Pixel Size [micrometer]'),
('CCDPXSZY', self.getPixelSize()[1],
'CCD Y Pixel Size [micrometer]')
]
# register image on ImageServer
server = getImageServer(self.getManager())
proxy = server.register(img)
# and finally compress the image
img.compress(multiprocess=True)
return proxy
def _saveImage(self, imageRequest, imageData, extra):
(mode, binning, top, left, width, height) = self._getReadoutModeInfo(
imageRequest["binning"], imageRequest["window"]
)
binFactor = extra.get("binning_factor", 1.0)
pix_w, pix_h = self.getPixelSize()
focal_length = self["telescope_focal_length"]
scale_x = binFactor * (((180 / pi) / focal_length) * (pix_w * 0.001))
scale_y = binFactor * (((180 / pi) / focal_length) * (pix_h * 0.001))
full_width, full_height = self.getPhysicalSize()
CRPIX1 = ((int(full_width / 2.0)) - left) - 1
CRPIX2 = ((int(full_height / 2.0)) - top) - 1
t0 = time.time()
img = Image.create(imageData, imageRequest)
img += [
(
"DATE-OBS",
ImageUtil.formatDate(extra.get("frame_start_time", dt.datetime.utcnow())),
"Date exposure started",
),
("CCD-TEMP", extra.get("frame_temperature", -275.0), "CCD Temperature at Exposure Start [deg. C]"),
("EXPTIME", float(imageRequest["exptime"]) or -1, "exposure time in seconds"),
("IMAGETYP", imageRequest["type"].strip(), "Image type"),
("SHUTTER", str(imageRequest["shutter"]), "Requested shutter state"),
("CRPIX1", CRPIX1, "coordinate system reference pixel"),
("CRPIX2", CRPIX2, "coordinate system reference pixel"),
("CD1_1", scale_x, "transformation matrix element (1,1)"),
("CD1_2", 0.0, "transformation matrix element (1,2)"),
("CD2_1", 0.0, "transformation matrix element (2,1)"),
("CD2_2", scale_y, "transformation matrix element (2,2)"),
("CAMERA", str(self["camera_model"]), "Camera Model"),
("CCD", str(self["ccd_model"]), "CCD Model"),
("CCD_DIMX", self.getPhysicalSize()[0], "CCD X Dimension Size"),
("CCD_DIMY", self.getPhysicalSize()[1], "CCD Y Dimension Size"),
("CCDPXSZX", self.getPixelSize()[0], "CCD X Pixel Size [micrometer]"),
("CCDPXSZY", self.getPixelSize()[1], "CCD Y Pixel Size [micrometer]"),
]
# register image on ImageServer
server = getImageServer(self.getManager())
proxy = server.register(img)
# and finally compress the image
img.compress(multiprocess=True)
return proxy
def getMetadata(self, request):
try:
return [('ENVMOD', str(self['model']), 'Weather station Model'),
('ENVTEM', self.temperature(unit_out=units.deg_C).value, '[degC] Weather station temperature'),
('ENVDEW', self.dew_point(unit_out=units.deg_C).value, '[degC] Weather station dew point temperature'),
('ENVHUM', self.humidity(unit_out=units.pct).value, '[%] Weather station relative humidity'),
('ENVWIN', self.wind_speed(unit_out=units.m / units.s).value, '[m/s] Weather station wind speed'),
('ENVDIR', self.wind_direction(unit_out=units.deg).value, '[deg] Weather station wind direction'),
('ENVPRE', self.pressure(unit_out=units.cds.mmHg).value, '[mmHg] Weather station air pressure'),
('ENVDAT', ImageUtil.formatDate(self._results['utctime']), 'UT time of the meteo observation')
]
except AttributeError:
return []
def getMetadata(self, request):
# Check first if there is metadata from an metadata override method.
md = self.getMetadataOverride(request)
if md is not None:
return md
# If not, just go on with the instrument's default metadata.
md = [('DATE-OBS', ImageUtil.formatDate(self.extra_header_info.get("frame_start_time", dt.datetime.utcnow())),
'Date exposure started'),
("EXPTIME", float(request['exptime']), "exposure time in seconds"),
('IMAGETYP', request['type'].strip(), 'Image type'),
('SHUTTER', str(request['shutter']), 'Requested shutter state'),
('INSTRUME', str(self['camera_model']), 'Name of instrument'),
('CCD', str(self['ccd_model']), 'CCD Model'),
('CCD_DIMX', self.getPhysicalSize()[0], 'CCD X Dimension Size'),
('CCD_DIMY', self.getPhysicalSize()[1], 'CCD Y Dimension Size'),
('CCDPXSZX', self.getPixelSize()[0], 'CCD X Pixel Size [micrometer]'),
('CCDPXSZY', self.getPixelSize()[1], 'CCD Y Pixel Size [micrometer]')]
if "frame_temperature" in self.extra_header_info.keys():
md += [('CCD-TEMP', self.extra_header_info["frame_temperature"],
'CCD Temperature at Exposure Start [deg. C]')]
focal_length = self["telescope_focal_length"]
if focal_length is not None: # If there is no telescope_focal_length defined, don't store WCS
mode, binning, top, left, width, height = self._getReadoutModeInfo(request["binning"], request["window"])
binFactor = self.extra_header_info.get("binning_factor", 1.0)
pix_w, pix_h = self.getPixelSize()
focal_length = self["telescope_focal_length"]
scale_x = binFactor * (((180 / pi) / focal_length) * (pix_w * 0.001))
scale_y = binFactor * (((180 / pi) / focal_length) * (pix_h * 0.001))
full_width, full_height = self.getPhysicalSize()
CRPIX1 = ((int(full_width / 2.0)) - left) - 1
CRPIX2 = ((int(full_height / 2.0)) - top) - 1
# Adding WCS coordinates according to FITS standard.
# Quick sheet: http://www.astro.iag.usp.br/~moser/notes/GAi_FITSimgs.html
# http://adsabs.harvard.edu/abs/2002A%26A...395.1061G
# http://adsabs.harvard.edu/abs/2002A%26A...395.1077C
md += [("CRPIX1", CRPIX1, "coordinate system reference pixel"),
("CRPIX2", CRPIX2, "coordinate system reference pixel"),
("CD1_1", scale_x * cos(self["rotation"]*pi/180.), "transformation matrix element (1,1)"),
("CD1_2", -scale_y * sin(self["rotation"]*pi/180.), "transformation matrix element (1,2)"),
("CD2_1", scale_x * sin(self["rotation"]*pi/180.), "transformation matrix element (2,1)"),
("CD2_2", scale_y * cos(self["rotation"]*pi/180.), "transformation matrix element (2,2)")]
return md
def getMetadata(self, request):
if not self._check():
return []
return [('ENVMOD', str(self['model']), 'Weather station Model'),
('ENVTEM', self.temperature(unit_out=units.deg_C).value,
'[degC] Weather station temperature'),
('ENVDEW', self.dew_point(unit_out=units.deg_C).value,
'[degC] Weather station dew point temperature'),
('ENVHUM', self.humidity(unit_out=units.pct).value,
'[%] Weather station relative humidity'),
('ENVWIN', self.wind_speed(unit_out=units.m / units.s).value,
'[m/s] Weather station wind speed'),
('ENVDIR', self.wind_direction(unit_out=units.deg).value,
'[deg] Weather station wind direction'),
('ENVPRE', self.pressure(unit_out=units.cds.mmHg).value,
'[mmHg] Weather station air pressure'),
('ENVDAT', ImageUtil.formatDate(self.obs_time()),
'UT time of the meteo observation')]
imageRequest.addPostHeaders(self.getManager())
binFactor = self._binning_factors[binning]
scale_x = binFactor * (((180 / pi) / cam["telescope_focal_length"]) *
(self.getPixelSize()[0] * 0.001))
scale_y = binFactor * (((180 / pi) / cam["telescope_focal_length"]) *
(self.getPixelSize()[1] * 0.001))
fullsize = self.drv.readoutModes[self.ccd][0]
CRPIX1 = ((int(fullsize.width / 2)) - left) - 1
CRPIX2 = ((int(fullsize.height / 2)) - top) - 1
img = Image.create(img, imageRequest)
img += [('DATE-OBS', ImageUtil.formatDate(self.lastFrameStartTime),
'Date exposure started'),
('CCD-TEMP', self.lastFrameTemp,
'CCD Temperature at Exposure Start [deg. C]'),
("EXPTIME", float(imageRequest['exptime'])
or -1, "exposure time in seconds"),
('IMAGETYP', imageRequest['type'].strip(), 'Image type'),
('SHUTTER', str(imageRequest['shutter']),
'Requested shutter state'),
("CRPIX1", CRPIX1, "coordinate system reference pixel"),
("CRPIX2", CRPIX2, "coordinate system reference pixel"),
("CD1_1", scale_x, "transformation matrix element (1,1)"),
("CD1_2", 0.0, "transformation matrix element (1,2)"),
("CD2_1", 0.0, "transformation matrix element (2,1)"),
("CD2_2", scale_y, "transformation matrix element (2,2)")]
def getMetadata(self, request):
# Check first if there is metadata from an metadata override method.
md = self.getMetadataOverride(request)
if md is not None:
return md
# If not, just go on with the instrument's default metadata.
md = [("EXPTIME", float(request['exptime']), "exposure time in seconds"),
('IMAGETYP', request['type'].strip(), 'Image type'),
('SHUTTER', str(request['shutter']), 'Requested shutter state'),
('INSTRUME', str(self['camera_model']), 'Name of instrument'),
('CCD', str(self['ccd_model']), 'CCD Model'),
('CCD_DIMX', self.getPhysicalSize()[0], 'CCD X Dimension Size'),
('CCD_DIMY', self.getPhysicalSize()[1], 'CCD Y Dimension Size'),
('CCDPXSZX', self.getPixelSize()[0], 'CCD X Pixel Size [micrometer]'),
('CCDPXSZY', self.getPixelSize()[1], 'CCD Y Pixel Size [micrometer]')]
if request['window'] is not None:
md += [('DETSEC', request['window'],
'Detector coodinates of the image')]
if "frame_temperature" in self.extra_header_info.keys():
md += [('CCD-TEMP', self.extra_header_info["frame_temperature"],
'CCD Temperature at Exposure Start [deg. C]')]
if "frame_start_time" in self.extra_header_info.keys():
md += [('DATE-OBS', ImageUtil.formatDate(self.extra_header_info.get("frame_start_time")),
'Date exposure started')]
mode, binning, top, left, width, height = self._getReadoutModeInfo(request["binning"], request["window"])
# Binning keyword: http://iraf.noao.edu/projects/ccdmosaic/imagedef/mosaic/MosaicV1.html
# CCD on-chip summing given as two or four integer numbers. These define
# the summing of CCD pixels in the amplifier readout order. The first
# two numbers give the number of pixels summed in the serial and parallel
# directions respectively. If the first pixel read out consists of fewer
# unbinned pixels along either direction the next two numbers give the
# number of pixels summed for the first serial and parallel pixels. From
# this it is implicit how many pixels are summed for the last pixels
# given the size of the CCD section (CCDSEC). It is highly recommended
# that controllers read out all pixels with the same summing in which
# case the size of the CCD section will be the summing factors times the
# size of the data section.
md += [("CCDSUM", binning.replace("x", " "), "CCD on-chip summing")]
focal_length = self["telescope_focal_length"]
if focal_length is not None: # If there is no telescope_focal_length defined, don't store WCS
binFactor = self.extra_header_info.get("binning_factor", 1.0)
pix_w, pix_h = self.getPixelSize()
focal_length = self["telescope_focal_length"]
scale_x = binFactor * (((180 / pi) / focal_length) * (pix_w * 0.001))
scale_y = binFactor * (((180 / pi) / focal_length) * (pix_h * 0.001))
full_width, full_height = self.getPhysicalSize()
CRPIX1 = ((int(full_width / 2.0)) - left) - 1
CRPIX2 = ((int(full_height / 2.0)) - top) - 1
# Adding WCS coordinates according to FITS standard.
# Quick sheet: http://www.astro.iag.usp.br/~moser/notes/GAi_FITSimgs.html
# http://adsabs.harvard.edu/abs/2002A%26A...395.1061G
# http://adsabs.harvard.edu/abs/2002A%26A...395.1077C
md += [("CRPIX1", CRPIX1, "coordinate system reference pixel"),
("CRPIX2", CRPIX2, "coordinate system reference pixel"),
("CD1_1", scale_x * cos(self["rotation"]*pi/180.), "transformation matrix element (1,1)"),
("CD1_2", -scale_y * sin(self["rotation"]*pi/180.), "transformation matrix element (1,2)"),
("CD2_1", scale_x * sin(self["rotation"]*pi/180.), "transformation matrix element (2,1)"),
("CD2_2", scale_y * cos(self["rotation"]*pi/180.), "transformation matrix element (2,2)")]
return md
except Exception, e:
self.log.warning("MakekFlat error: " + str(e))
#Last resort if nothing else could make a picture
if (pix == None):
pix = N.zeros((100, 100), dtype=N.int32)
imageRequest.fetchPostHeaders(self.getManager())
img = Image.create(pix, imageRequest)
# update image request
imageRequest["filename"] = img.filename()
img += [('DATE-OBS',
ImageUtil.formatDate(
dt.datetime.fromtimestamp(self.__lastFrameStart)))]
server = getImageServer(self.getManager())
proxy = server.register(img)
self.readoutComplete(proxy)
return proxy
def abortExposure(self):
if not self.isExposing(): return
self.__abort.set()
# busy waiting for exposure/readout stops
def _saveImage(self, imageRequest, imageData, extra):
(mode, binning, top, left,
width, height) = self._getReadoutModeInfo(imageRequest["binning"],
imageRequest["window"])
binFactor = extra.get("binning_factor", 1.0)
pix_w, pix_h = self.getPixelSize()
img = Image.create(imageData, imageRequest)
if self["telescope_focal_length"] is not None: # If there is no telescope_focal_length defined, don't store WCS
focal_length = self["telescope_focal_length"]
scale_x = binFactor * (((180 / pi) / focal_length) * (pix_w * 0.001))
scale_y = binFactor * (((180 / pi) / focal_length) * (pix_h * 0.001))
full_width, full_height = self.getPhysicalSize()
CRPIX1 = ((int(full_width / 2.0)) - left) - 1
CRPIX2 = ((int(full_height / 2.0)) - top) - 1
# Adding WCS coordinates according to FITS standard.
# Quick sheet: http://www.astro.iag.usp.br/~moser/notes/GAi_FITSimgs.html
# http://adsabs.harvard.edu/abs/2002A%26A...395.1061G
# http://adsabs.harvard.edu/abs/2002A%26A...395.1077C
img += [("CRPIX1", CRPIX1, "coordinate system reference pixel"),
("CRPIX2", CRPIX2, "coordinate system reference pixel"),
("CD1_1", scale_x * cos(self["rotation"]*pi/180.), "transformation matrix element (1,1)"),
("CD1_2", -scale_y * sin(self["rotation"]*pi/180.), "transformation matrix element (1,2)"),
("CD2_1", scale_x * sin(self["rotation"]*pi/180.), "transformation matrix element (2,1)"),
("CD2_2", scale_y * cos(self["rotation"]*pi/180.), "transformation matrix element (2,2)")]
img += [('DATE-OBS',
ImageUtil.formatDate(
extra.get("frame_start_time", dt.datetime.utcnow())),
'Date exposure started'),
('CCD-TEMP', extra.get("frame_temperature", -275.0),
'CCD Temperature at Exposure Start [deg. C]'),
("EXPTIME", float(imageRequest['exptime']) or -1,
"exposure time in seconds"),
('IMAGETYP', imageRequest['type'].strip(),
'Image type'),
('SHUTTER', str(imageRequest['shutter']),
'Requested shutter state'),
('INSTRUME', str(self['camera_model']), 'Name of instrument'),
('CCD', str(self['ccd_model']), 'CCD Model'),
('CCD_DIMX', self.getPhysicalSize()
[0], 'CCD X Dimension Size'),
('CCD_DIMY', self.getPhysicalSize()
[1], 'CCD Y Dimension Size'),
('CCDPXSZX', self.getPixelSize()[0],
'CCD X Pixel Size [micrometer]'),
('CCDPXSZY', self.getPixelSize()[1],
'CCD Y Pixel Size [micrometer]')]
# register image on ImageServer
server = getImageServer(self.getManager())
proxy = server.register(img)
# and finally compress the image if asked
if imageRequest['compress_format'].lower() != 'no':
img.compress(format=imageRequest['compress_format'], multiprocess=True)
return proxy
def getMetadata(self, request):
# Check first if there is metadata from an metadata override method.
md = self.getMetadataOverride(request)
if md is not None:
return md
# If not, just go on with the instrument's default metadata.
md = [('DATE-OBS',
ImageUtil.formatDate(
self.extra_header_info.get("frame_start_time",
dt.datetime.utcnow())),
'Date exposure started'),
("EXPTIME", float(request['exptime']),
"exposure time in seconds"),
('IMAGETYP', request['type'].strip(), 'Image type'),
('SHUTTER', str(request['shutter']), 'Requested shutter state'),
('INSTRUME', str(self['camera_model']), 'Name of instrument'),
('CCD', str(self['ccd_model']), 'CCD Model'),
('CCD_DIMX', self.getPhysicalSize()[0], 'CCD X Dimension Size'),
('CCD_DIMY', self.getPhysicalSize()[1], 'CCD Y Dimension Size'),
('CCDPXSZX', self.getPixelSize()[0],
'CCD X Pixel Size [micrometer]'),
('CCDPXSZY', self.getPixelSize()[1],
'CCD Y Pixel Size [micrometer]')]
if "frame_temperature" in self.extra_header_info.keys():
md += [('CCD-TEMP', self.extra_header_info["frame_temperature"],
'CCD Temperature at Exposure Start [deg. C]')]
focal_length = self["telescope_focal_length"]
if focal_length is not None: # If there is no telescope_focal_length defined, don't store WCS
mode, binning, top, left, width, height = self._getReadoutModeInfo(
request["binning"], request["window"])
binFactor = self.extra_header_info.get("binning_factor", 1.0)
pix_w, pix_h = self.getPixelSize()
focal_length = self["telescope_focal_length"]
scale_x = binFactor * (((180 / pi) / focal_length) *
(pix_w * 0.001))
scale_y = binFactor * (((180 / pi) / focal_length) *
(pix_h * 0.001))
full_width, full_height = self.getPhysicalSize()
CRPIX1 = ((int(full_width / 2.0)) - left) - 1
CRPIX2 = ((int(full_height / 2.0)) - top) - 1
# Adding WCS coordinates according to FITS standard.
# Quick sheet: http://www.astro.iag.usp.br/~moser/notes/GAi_FITSimgs.html
# http://adsabs.harvard.edu/abs/2002A%26A...395.1061G
# http://adsabs.harvard.edu/abs/2002A%26A...395.1077C
md += [("CRPIX1", CRPIX1, "coordinate system reference pixel"),
("CRPIX2", CRPIX2, "coordinate system reference pixel"),
("CD1_1", scale_x * cos(self["rotation"] * pi / 180.),
"transformation matrix element (1,1)"),
("CD1_2", -scale_y * sin(self["rotation"] * pi / 180.),
"transformation matrix element (1,2)"),
("CD2_1", scale_x * sin(self["rotation"] * pi / 180.),
"transformation matrix element (2,1)"),
("CD2_2", scale_y * cos(self["rotation"] * pi / 180.),
"transformation matrix element (2,2)")]
return md
