Python items Examples

Example #1

File: cube_obs_table.py Project: adamginsburg/paper_w51_evla

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],

Example #2

File: cube_noise_estimates.py Project: hongliliu/paper_w51_evla

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

Example #3

File: load_cubes_K.py Project: adamginsburg/paper_w51_evla

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')

Example #4

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()
}