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.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)
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):
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):
# 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'
