import paths from paths import dpath from spectral_cube import SpectralCube from astropy import units as u from fnames import cube_names,rrl_cube_names from astropy.table import Table,Column from latex_info import latexdict, format_float cube_names.update(rrl_cube_names) errors = {} cubes = {} beams = {} dv = {} for name, fn in cube_names.items(): cube = SpectralCube.read(dpath(fn)).with_spectral_unit(u.km/u.s, velocity_convention='radio') print(cube) err = cube.std().to(u.mJy) print("{0} std: {1}".format(name, err)) errors[name] = err cubes[name] = cube beams[name] = cube.beam dv[name] = np.diff(cube.spectral_axis).mean() name_mapping = {'11_natural': '1-1 Natural', '11_uniform': '1-1 Uniform', '22_briggs0': '2-2 Briggs 0', '22_natural': '2-2 Natural', 'H77a': 'H$77\\alpha$ Briggs 0'} tbl = Table([Column(data=[name_mapping[name] for name in cube_names],
import numpy as np from paths import dpath from spectral_cube import SpectralCube from astropy import units as u from fnames import cube_names errors = {} for name, fn in cube_names.items(): cube = SpectralCube.read(dpath(fn)) print(cube) err = cube.std().to(u.mJy) print("{0} std: {1}".format(name, err)) errors[name] = err
from astropy import units as u from astropy.io import fits from spectral_cube import SpectralCube from paths import dpath from fnames import cube_names, continua import radio_beam import ds9 dd = ds9.ds9() cubes = {name: SpectralCube.read(dpath(fn)) for name,fn in cube_names.items()} vcubes = {name: cube.with_spectral_unit(u.km/u.s, velocity_convention='radio') for name,cube in cubes.items()} beams = {name: radio_beam.Beam.from_fits_header(cube.header) for name,cube in vcubes.items()} for cube in vcubes.values(): cube._unit = u.Jy vkcubes = {name: cube.with_flux_unit(u.K, u.brightness_temperature(beams[name], cube.wcs.wcs.restfrq*u.Hz)) for name,cube in vcubes.items()} for name,vc in vkcubes.items(): vc._header['OBJECT'] = name vc.to_ds9(dd.id, newframe=True) dd.set('tile yes')
from astropy import units as u from astropy.io import fits from spectral_cube import SpectralCube from paths import dpath from fnames import cube_names, continua import radio_beam import ds9 dd = ds9.ds9() cubes = {name: SpectralCube.read(dpath(fn)) for name, fn in cube_names.items()} vcubes = { name: cube.with_spectral_unit(u.km / u.s, velocity_convention='radio') for name, cube in cubes.items() } beams = { name: radio_beam.Beam.from_fits_header(cube.header) for name, cube in vcubes.items() } for cube in vcubes.values(): cube._unit = u.Jy vkcubes = { name: cube.with_flux_unit( u.K, u.brightness_temperature(beams[name], cube.wcs.wcs.restfrq * u.Hz)) for name, cube in vcubes.items() }