How to use the pysynphot.Observation function in pysynphot
To help you get started, we’ve selected a few pysynphot examples, based on popular ways it is used in public projects.
Secure your code as it's written. Use Snyk Code to scan source code in minutes - no build needed - and fix issues immediately.
gbrammer / grizli / grizli / utils.py View on Github 
>>> import pysynphot as S
>>> bandpasses = [S.ObsBandpass('wfc3,ir,f140w')]
Returns
-------
ab_mags : dict
Dictionary with keys from `bandpasses` and the integrated magnitudes
"""
import pysynphot as S
flat = S.FlatSpectrum(0, fluxunits='ABMag')
ab_mags = OrderedDict()
for bp in bandpasses:
flat_obs = S.Observation(flat, bp)
spec_obs = S.Observation(spectrum, bp)
ab_mags[bp.name] = -2.5*np.log10(spec_obs.countrate()/flat_obs.countrate())
return ab_mags
def calc_fluxes():
return np.array([S.Observation(sp,i).effstim("flam") for i in pfilts])if False:
import pysynphot as S
n = 1.e-20
spec = S.FlatSpectrum(n, fluxunits='flam')
photflam = {}
photplam = {}
for filter in ['F098M', 'F105W', 'F110W', 'F125W', 'F140W', 'F160W']:
bp = S.ObsBandpass('wfc3,ir,%s' %(filter.lower()))
photplam[filter] = bp.pivot()
obs = S.Observation(spec, bp)
photflam[filter] = n/obs.countrate()
for filter in ['F435W', 'F606W', 'F775W', 'F814W']:
bp = S.ObsBandpass('acs,wfc1,%s' %(filter.lower()))
photplam[filter] = bp.pivot()
obs = S.Observation(spec, bp)
photflam[filter] = n/obs.countrate()
class GrismFLT(object):
"""
Scripts for simple modeling of individual grism FLT images
tbd:
o helper functions for extracting 2D spectra
o lots of book-keeping for handling SExtractor objects & catalogs
...
"""
def __init__(self, flt_file='ico205lwq_flt.fits', sci_ext=('SCI',1),
direct_image=None, refimage=None, segimage=None, refext=0,
verbose=True, pad=100, shrink_segimage=True,>>> bandpasses = [S.ObsBandpass('wfc3,ir,f140w')]
Returns
-------
ab_mags : dict
Dictionary with keys from `bandpasses` and the integrated magnitudes
"""
import pysynphot as S
flat = S.FlatSpectrum(0, fluxunits='ABMag')
ab_mags = OrderedDict()
for bp in bandpasses:
flat_obs = S.Observation(flat, bp)
spec_obs = S.Observation(spectrum, bp)
ab_mags[bp.name] = -2.5*np.log10(spec_obs.countrate()/flat_obs.countrate())
return ab_mags
'VISTAH':1.6433e+04,
'GRISM': 1.6e4}
no_newline = '\x1b[1A\x1b[1M' # character to skip clearing line on STDOUT printing
### Demo for computing photflam and photplam with pysynphot
if False:
import pysynphot as S
n = 1.e-20
spec = S.FlatSpectrum(n, fluxunits='flam')
photflam = {}
photplam = {}
for filter in ['F098M', 'F105W', 'F110W', 'F125W', 'F140W', 'F160W']:
bp = S.ObsBandpass('wfc3,ir,%s' %(filter.lower()))
photplam[filter] = bp.pivot()
obs = S.Observation(spec, bp)
photflam[filter] = n/obs.countrate()
for filter in ['F435W', 'F606W', 'F775W', 'F814W']:
bp = S.ObsBandpass('acs,wfc1,%s' %(filter.lower()))
photplam[filter] = bp.pivot()
obs = S.Observation(spec, bp)
photflam[filter] = n/obs.countrate()
class GrismFLT(object):
"""
Scripts for simple modeling of individual grism FLT images
tbd:
o helper functions for extracting 2D spectra
o lots of book-keeping for handling SExtractor objects & catalogs
