def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) self.meta['CUNIT1'] = self.meta.get('CUNIT1', 'arcsec').lower() self.meta['CUNIT2'] = self.meta.get('CUNIT2', 'arcsec').lower() # Fill in some missing or broken info # Test if change has already been applied if 'T' not in self.meta['date-obs']: datestr = "{date}T{time}".format(date=self.meta.get('date-obs', self.meta.get('date_obs') ), time=self.meta.get('time-obs', self.meta.get('time_obs') ) ) self.meta['date-obs'] = datestr # If non-standard Keyword is present, correct it too, for compatibility. if 'date_obs' in self.meta: self.meta['date_obs'] = self.meta['date-obs'] self._nickname = self.instrument + "-" + self.detector self.plot_settings['cmap'] = plt.get_cmap('soholasco{det!s}'.format(det=self.detector[1])) self.plot_settings['norm'] = ImageNormalize(stretch=source_stretch(self.meta, PowerStretch(0.5)))
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) self.meta['CUNIT1'] = self.meta['CUNIT1'].lower() self.meta['CUNIT2'] = self.meta['CUNIT2'].lower() # Fill in some missing or broken info datestr = "{date}T{time}".format(date=self.meta.get('date-obs', self.meta.get('date_obs') ), time=self.meta.get('time-obs', self.meta.get('time_obs') ) ) self.meta['date-obs'] = datestr # If non-standard Keyword is present, correct it too, for compatibility. if 'date_obs' in self.meta: self.meta['date_obs'] = self.meta['date-obs'] self.meta['wavelnth'] = np.nan self.meta['waveunit'] = 'nm' self._nickname = self.instrument + "-" + self.detector self.plot_settings['cmap'] = cm.get_cmap('soholasco{det!s}'.format(det=self.detector[1])) self.plot_settings['norm'] = ImageNormalize(stretch=PowerStretch(0.5))
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) self._name = self.observatory + " " + self.detector + " " + str(self.measurement) self._nickname = "{0}-{1}".format(self.detector, self.observatory[-1]) self.cmap = cm.get_cmap('stereocor%s' % self.detector[-1])
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) self.meta['detector'] = "HMI" # self.meta['instrme'] = "HMI" # self.meta['obsrvtry'] = "SDO" self._nickname = self.detector
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) # Fill in some missing info self.meta['detector'] = "AIA" self._nickname = self.detector self.plot_settings['cmap'] = plt.get_cmap(self._get_cmap_name()) self.plot_settings['norm'] = ImageNormalize(stretch=source_stretch(self.meta, AsinhStretch(0.01)))
def __init__(self, data, header, **kwargs): # Assume pixel units are arcesc if not given header['cunit1'] = header.get('cunit1', 'arcsec') header['cunit2'] = header.get('cunit2', 'arcsec') GenericMap.__init__(self, data, header, **kwargs) self.meta['detector'] = "SJI" self.meta['waveunit'] = "Angstrom" self.meta['wavelnth'] = header['twave1']
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) # Fill in some missing or broken info self.meta['detector'] = "MDI" self._fix_dsun() self._name = self.detector + " " + self.measurement self._nickname = self.detector + " " + self.measurement
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) # Fill in some missing or broken info self.meta['detector'] = "MDI" self._fix_dsun() self._name = self.observatory + " " + self.detector self._nickname = self.detector + " " + self.measurement
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) # Fill in some missing info self.meta['detector'] = "EIT" self.meta['waveunit'] = "Angstrom" self._fix_dsun() self._nickname = self.detector self.plot_settings['cmap'] = plt.get_cmap(self._get_cmap_name()) self.plot_settings['norm'] = ImageNormalize(stretch=source_stretch(self.meta, PowerStretch(0.5)))
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) self._nickname = "{0}-{1}".format(self.detector, self.observatory[-1]) self.plot_settings['cmap'] = cm.get_cmap('stereohi{det!s}'.format(det=self.detector[-1])) self.plot_settings['norm'] = ImageNormalize(stretch=source_stretch(self.meta, PowerStretch(0.25))) # Try to identify when the FITS meta data does not have the correct # date FITS keyword if ('date_obs' in self.meta) and not('date-obs' in self.meta): self.meta['date-obs'] = self.meta['date_obs']
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) # Fill in some missing info self.meta['detector'] = "AIA" # self.meta['instrme'] = "AIA" self._nickname = self.detector self._name = self.detector + " " + str(self.measurement) self.cmap = cm.get_cmap('sdoaia{wl:d}'.format(wl=self.wavelength))
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) # Fill in some missing info self.meta['detector'] = "EIT" self._fix_dsun() self._name = self.detector + " " + str(self.measurement) self._nickname = self.detector self.cmap = cm.get_cmap('sohoeit{wl:d}'.format(wl=self.wavelength))
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) # It needs to be verified that these must actually be set and are not # already in the header. self.meta['detector'] = "TRACE" self.meta['obsrvtry'] = "TRACE" self._nickname = self.detector # Colour maps self.plot_settings['cmap'] = cm.get_cmap('trace' + str(self.meta['WAVE_LEN'])) self.plot_settings['norm'] = ImageNormalize(stretch=source_stretch(self.meta, LogStretch()))
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) self._nickname = "{0}-{1}".format(self.detector, self.observatory[-1]) self.plot_settings['cmap'] = cm.get_cmap('sohoeit{wl:d}'.format(wl=int(self.wavelength.value))) self.plot_settings['norm'] = ImageNormalize(stretch=source_stretch(self.meta, PowerStretch(0.25))) self.meta['waveunit'] = 'Angstrom' # Try to identify when the FITS meta data does not have the correct # date FITS keyword if ('date_obs' in self.meta) and not('date-obs' in self.meta): self.meta['date-obs'] = self.meta['date_obs']
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) # It needs to be verified that these must actually be set and are not # already in the header. self.meta['detector'] = "TRACE" self.meta['obsrvtry'] = "TRACE" self._nickname = self.detector # Colour maps self.plot_settings['cmap'] = cm.get_cmap('trace' + self.measurement) self.plot_settings['norm'] = colors.LogNorm()
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) # It needs to be verified that these must actually be set and # are not already in the header. self.meta['detector'] = "SWAP" # self.meta['instrme'] = "SWAP" self.meta['obsrvtry'] = "PROBA2" self._nickname = self.detector self.plot_settings['cmap'] = plt.get_cmap(name='sdoaia171')
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) self._name = self.observatory + " " + self.detector + " " + str(self.measurement) self._nickname = "{0}-{1}".format(self.detector, self.observatory[-1]) self.cmap = cm.get_cmap('stereohi%s' % self.detector[-1]) # Try to identify when the FITS meta data does not have the correct # date FITS keyword if ('date_obs' in self.meta) and not('date-obs' in self.meta): self.meta['date-obs'] = self.meta['date_obs']
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) self._name = "RHESSI {measure[0]:.0f} - {measure[1]:.0f} keV".format(measure=self.measurement) self._nickname = self.detector # Fix some broken/misapplied keywords if self.meta['ctype1'] == 'arcsec': self.meta['cunit1'] = 'arcsec' self.meta['ctype1'] = 'HPLN-TAN' if self.meta['ctype2'] == 'arcsec': self.meta['cunit2'] = 'arcsec' self.meta['ctype2'] = 'HPLT-TAN'
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) # It needs to be verified that these must actually be set and are not # already in the header. self.meta['detector'] = "SWAP" # self.meta['instrme'] = "SWAP" self.meta['obsrvtry'] = "PROBA2" self._name = self.detector + " " + str(self.measurement) self._nickname = self.detector self.cmap = cm.get_cmap(name='sdoaia171')
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) # Fill in some missing or broken info self.meta['detector'] = "MDI" self._fix_dsun() self._nickname = self.detector + " " + self.measurement vmin = np.nanmin(self.data) vmax = np.nanmax(self.data) if abs(vmin) > abs(vmax): self.plot_settings['norm'] = colors.Normalize(-vmin, vmin) else: self.plot_settings['norm'] = colors.Normalize(-vmax, vmax)
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) self._name = "RHESSI %d - %d keV" % (self.measurement[0], self.measurement[1]) self._nickname = self.detector # Fix some broken/misapplied keywords if self.meta['ctype1'] == 'arcsec': self.meta['cunit1'] = 'arcsec' self.meta['ctype1'] = 'HPLN-TAN' if self.meta['ctype2'] == 'arcsec': self.meta['cunit2'] = 'arcsec' self.meta['ctype2'] = 'HPLT-TAN'
def __init__(self, data, header, **kwargs): # Assume pixel units are arcesc if not given header['cunit1'] = header.get('cunit1', 'arcsec') header['cunit2'] = header.get('cunit2', 'arcsec') GenericMap.__init__(self, data, header, **kwargs) # Fill in some missing info self.meta['detector'] = "EIT" self.meta['waveunit'] = "Angstrom" self._fix_dsun() self._nickname = self.detector self.plot_settings['cmap'] = plt.get_cmap(self._get_cmap_name()) self.plot_settings['norm'] = ImageNormalize( stretch=source_stretch(self.meta, PowerStretch(0.5)))
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) # It needs to be verified that these must actually be set and are not # already in the header. self.meta['detector'] = "TRACE" self.meta['obsrvtry'] = "TRACE" # Name that will appear at the top of a TRACE image plot self._name = self.detector + " " + self.measurement self._nickname = self.detector # Colour maps self.cmap = cm.get_cmap('trace' + self.measurement)
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) # Fill in some missing or broken info datestr = "%sT%s" % (self.meta.get('date-obs', self.meta.get( 'date_obs')), self.meta.get('time-obs', self.meta.get('time_obs'))) self.meta['date-obs'] = datestr # If non-standard Keyword is present, correct it too, for compatibility. if 'date_obs' in self.meta: self.meta['date_obs'] = self.meta['date-obs'] self._name = self.instrument + " " + self.detector + " " + self.measurement self._nickname = self.instrument + "-" + self.detector self.cmap = cm.get_cmap('soholasco%s' % self.detector[1])
def __init__(self, data, header, **kwargs): # Assume pixel units are arcesc if not given header['cunit1'] = header.get('cunit1', 'arcsec') header['cunit2'] = header.get('cunit2', 'arcsec') GenericMap.__init__(self, data, header, **kwargs) # Fill in some missing or broken info self.meta['detector'] = "MDI" self._fix_dsun() self._nickname = self.detector + " " + self.measurement vmin = np.nanmin(self.data) vmax = np.nanmax(self.data) if abs(vmin) > abs(vmax): self.plot_settings['norm'] = colors.Normalize(-vmin, vmin) else: self.plot_settings['norm'] = colors.Normalize(-vmax, vmax)
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) self._nickname = self.detector #TODO Currently (8/29/2011), cannot read fits files containing more than one image (schriste) # Fix some broken/misapplied keywords if self.meta['ctype1'] == 'arcsec': self.meta['cunit1'] = 'arcsec' self.meta['ctype1'] = 'HPLN-TAN' if self.meta['ctype2'] == 'arcsec': self.meta['cunit2'] = 'arcsec' self.meta['ctype2'] = 'HPLT-TAN' self.meta['waveunit'] = 'keV' self.meta['wavelnth'] = [self.meta['energy_l'], self.meta['energy_h']] self.plot_settings['cmap'] = cm.get_cmap('rhessi')
def __init__(self, data, header, **kwargs): # Assume pixel units are arcesc if not given header['cunit1'] = header.get('cunit1', 'arcsec') header['cunit2'] = header.get('cunit2', 'arcsec') GenericMap.__init__(self, data, header, **kwargs) # It needs to be verified that these must actually be set and are not # already in the header. self.meta['detector'] = "TRACE" self.meta['obsrvtry'] = "TRACE" self._nickname = self.detector # Colour maps self.plot_settings['cmap'] = 'trace' + str(self.meta['WAVE_LEN']) self.plot_settings['norm'] = ImageNormalize(stretch=source_stretch( self.meta, LogStretch()), clip=False)
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) if header.get('lvl_num') == 2: self.meta['wavelnth'] = header.get('twave1') self.meta['detector'] = header.get('instrume') self.meta['waveunit'] = "Angstrom" if header.get('lvl_num') == 1: self.meta['wavelnth'] = int(header.get('img_path').split('_')[1]) self.meta['waveunit'] = "Angstrom" self.meta['detector'] = "SJI" self.meta['waveunit'] = "Angstrom" palette = cm.get_cmap('irissji' + str(int(self.meta['wavelnth']))) palette.set_bad('black') self.plot_settings['cmap'] = palette self.plot_settings['norm'] = ImageNormalize( stretch=visualization.AsinhStretch(0.1))
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) self.meta['detector'] = "SOT" self.meta['telescop'] = "Hinode" self._nickname = self.detector #TODO (add other options, Now all threated as intensity. This follows Hinode SDC archive) # StokesQUV -> grey, Velocity -> EIS, Width -> EIS, Mag Field Azi -> IDL 5 (STD gamma II) # 'WB' -> red # 'NB'(0 = red); (>0 = gray), # nb has 1 stokes I, the rest quv # 'SP' (<=1 = red); (>1 = gray) #sp has 2 stokes I, the rest quv color = {'SOT/WB': 'intensity', 'SOT/NB': 'intensity', # For the 1st dimension 'SOT/SP': 'intensity', # For the 1st 2 dimensions } self.plot_settings['cmap'] = cm.get_cmap('hinodesot' + color[self.instrument])
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) self._nickname = "{}-{}".format(self.detector, self.observatory[-1]) self.plot_settings['cmap'] = 'sohoeit{wl:d}'.format( wl=int(self.wavelength.value)) self.plot_settings['norm'] = ImageNormalize(stretch=source_stretch( self.meta, PowerStretch(0.25)), clip=False) self.meta['waveunit'] = 'Angstrom' # Try to identify when the FITS meta data does not have the correct # date FITS keyword if ('date_obs' in self.meta) and not ('date-obs' in self.meta): self.meta['date-obs'] = self.meta['date_obs'] # fix CROTA to CROTAn if "crota" in self.meta and "crota2" not in self.meta: log.debug("EUVIMap: Changing the CROTA keyword to CROTA2") self.meta["crota2"] = self.meta.pop("crota")
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) self.meta['detector'] = "SXT" self.meta['telescop'] = "Yohkoh" self.plot_settings['cmap'] = cm.get_cmap(name='yohkohsxt' + self.measurement[0:2].lower()) self.plot_settings['norm'] = ImageNormalize(stretch=PowerStretch(0.5)) # 2012/12/19 - the SXT headers do not have a value of the distance from # the spacecraft to the center of the Sun. The FITS keyword 'DSUN_OBS' # appears to refer to the observed diameter of the Sun. Until such # time as that is calculated and properly included in the file, we will # use simple trigonometry to calculate the distance of the center of # the Sun from the spacecraft. Note that the small angle approximation # is used, and the solar radius stored in SXT FITS files is in arcseconds. self.meta['dsun_apparent'] = constants.au if 'solar_r' in self.meta: self.meta['dsun_apparent'] = constants.radius/(np.deg2rad(self.meta['solar_r']/3600.0))
def __init__(self, data, header, **kwargs): GenericMap.__init__(self, data, header, **kwargs) # converting data array to masked array # self.data = ma.masked_where(self.data > SATURATION_LIMIT, self.data) fw1 = header.get('EC_FW1_') if fw1.lower() not in _lower_list(self.filter_wheel1_measurements): raise ValueError('Unpexpected filter wheel 1 in header.') fw1 = fw1.replace("_", " ") fw2 = header.get('EC_FW2_') if fw2.lower() not in _lower_list(self.filter_wheel2_measurements): raise ValueError('Unpexpected filter wheel 2 in header.') fw2 = fw2.replace("_", " ") self.meta['detector'] = "XRT" # self.meta['instrume'] = "XRT" self.meta['telescop'] = "Hinode" self.plot_settings['cmap'] = 'hinodexrt'