keflavich

import numpy as np
from astropy import units as u
from astropy.convolution import Gaussian1DKernel, Gaussian2DKernel
import os
import glob
from spectral_cube import SpectralCube


for fn in [
    "G001_13CO21_SEDIGISM_DR1c.fits",

keflavich

import numpy as np
from astropy.io import fits
import reproject
from astropy import wcs

base_header = fits.getheader('moments/G0_13CO21_downsampled_mom0.fits')
base_wcs = wcs.WCS(base_header).celestial

new_header = base_wcs.to_header()

keflavich

import numpy as np
from astropy.io import fits
import reproject
from astropy import wcs

base_header = fits.getheader('G000_13CO21_SEDIGISM_DR1c_spatialdownsample_4x_downsampled_2kms.fits')
base_wcs = wcs.WCS(base_header)

new_header = base_wcs.to_header()

keflavich

import os
from astropy import constants, units as u, table, stats, coordinates, wcs, log, coordinates as coord, convolution, modeling, time; from astropy.io import fits, ascii
import reproject
import glob
from reproject.mosaicking import find_optimal_celestial_wcs
from reproject.mosaicking import reproject_and_coadd
from reproject import reproject_interp

import requests
from bs4 import BeautifulSoup

keflavich

keflavich

import glob
from astropy import table
import requests
import tarfile
import os

url = "http://specmodels.iag.usp.br/fits_search/compress/s_coelho14_sed.tgz"
fn = os.path.basename(url)
if not os.path.exists(fn):
    resp = requests.get(url)

keflavich

import numpy as np

from glue.core import Data, DataCollection
from glue.app.qt.application import GlueApplication
try:
    from glue.viewers.image.qt.data_viewer import ImageViewer
except ImportError:
    from glue.viewers.image.qt.viewer_widget import ImageWidget as ImageViewer

from glue.core.coordinates import coordinates_from_wcs

keflavich

import numpy as np
from astroquery.vizier import Vizier
import pylab as pl
from astropy import constants, units as u, table, stats, coordinates, wcs, log, coordinates as coord, convolution, modeling, time; from astropy.io import fits, ascii

Vizier.ROW_LIMIT = 3e5
rslt = Vizier.query_constraints(catalog="II/246", GLAT='<5 & >-5', Kmag='<7', GLON='<45 | >315')[0]
crds = coordinates.SkyCoord(rslt['RAJ2000'], rslt['DEJ2000'], frame='fk5', unit=(u.deg, u.deg)).galactic

gridy,gridx = np.mgrid[-5:5:0.333, -45:45:0.333]

keflavich

keflavich