用spectral_cube抽取cube文件上某个点或者面源的光谱

最新推荐文章于 2023-07-28 13:27:45 发布

原创

最新推荐文章于 2023-07-28 13:27:45 发布 · 809 阅读

5 ·

CC 4.0 BY-SA版权

文章标签：

#python

from astropy.io import fits
from spectral_cube import SpectralCube as spec
import numpy as np
import matplotlib.pyplot as plt
import os
import sys
from astropy.coordinates import SkyCoord
from astropy.wcs import WCS
from astropy.table import Table
import warnings
warnings.filterwarnings('ignore')
T = Table.read('atoms_B3_12M_allsour_cont.tbl',format='ascii')
#把需要的位置信息列成表格，表格格式如下


def seek_file(sourcename):
    files = []
    for filename in os.listdir(os.getcwd()):
        fp = os.path.join(filename)
        spw35 = sourcename+'_sci.spw35.cube.I.pbcor.fits'
        spw37 = sourcename+'_sci.spw37.cube.I.pbcor.fits'
        spw39 = sourcename+'_sci.spw39.cube.I.pbcor.fits'
        if os.path.isfile(fp) and spw35 in filename:
            files.append(fp)
        if os.path.isfile(fp) and spw37 in filename:
            files.append(fp)
        if os.path.isfile(fp) and spw39 in filename:
            files.append(fp)
    if len(files)==3:
        files.so