uycvbl的博客

由于ALMA数据是4维的，不能直接用aplpy去读。 把下面的脚本放在数据所在目录下，调用一下就能用aplpy画ALMA的图了 #主要用于连续谱图，不要用输出的hdu覆盖原本的fitsfile，会丢失header里的重要信息 from astropy.io import fits from astropy.wcs import WCS def Get2DHDU(fitsfile): fits_file = fits.open(fitsfile) data = fits_file[0].data hd

##### make moment 0 map from spectral_cube import SpectralCube as sc from astropy.io import fits import astropy.units as u import astropy.constants as c import matplotlib.pyplot as plt plt.ion() import os spw = ['spw35','spw37'] def change_xaxis_in_velo.

import numpy as np import astropy as a import astropy.units as u import astropy.constants as c from scipy import optimize import matplotlib.pyplot as plt plt.ion() B = 24690.2 h = 6.626e-27 #erg/s k = 1.38e-16 #erg/K Tbg = 2.73 #K nu = [35

import matplotlib.pyplot as plt import os import numpy as np from astropy.coordinates import SkyCoord from astropy.table import Table from astropy.io import fits import warnings warnings.filterwarnings('ignore') T = Table.read('atoms_B3_12M_hotcore_candi.