nondefined-CSDN博客

原创射电观测中的流量密度单位Jy/beam转换成亮温度K

from astropy.io import fitsimport numpy as npimport astropy.constants as cimport astropy.units as ucont_file = 'member.uid___A001_X1465_X284f.I08448-4343_sci.spw25_27_29_31_33_35_37_39.cont.I.tt0.pbcor.fits'hdu = fits.open(cont_file)hdr = hdu[0].hea

2021-05-20 11:32:35 996

原创 line survey后续二：标分子名

为了把标完分子后的图不至于很杂乱，尽量做到清晰简洁。思路：窗口很宽，先把窗口拆分为二；标分子是个大工程。我们可以先把跃迁线较少的分子标上；对于很多跃迁线的分子，如果每一个跃迁位置都标分子就会显得很杂乱，我的做法是把靠的比较近的跃迁统一只标一个分子名；调整标记的上下左右位置。代码如下：import matplotlib.pyplot as pltimport numpy as npplt.ion()spw = 'spw35'data_file='I17233-3606_%s_17h

2021-01-08 13:37:53 205

原创备忘--identify.py

myXCLASS的拟合非常耗时间，如果不小心关掉面板再想重新看图就特别耗时间。这里提供了一个简单的命令，在不改变molfit文件，只重复前一次拟合的结果时，is_molfit_file_changed = ‘no’再执行就能不拟合，直接画图。##remove the old runing dataimport osimport shutilpath = '/home/zhang/software/XCLASS-Interface/run/myXCLASS/'try: files =

2020-12-22 15:55:35 160

原创备忘--myXCLASSFit.py

备忘录把拟合结果移动到当前文件夹##remove the old runing dataimport osimport shutilpath = '/home/zhang/software/XCLASS-Interface/run/myXCLASSFit/'files = os.listdir(path)if len(files)>0: shutil.rmtree(path)MolfitsFileName = "mol.molfit"experimentalDa

2020-12-11 16:38:15 121

原创获取aplpy可直接读取画图的2Dhdu

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

2020-11-25 20:19:12 311

原创从ALMA的data cube导出moment0

##### make moment 0 mapfrom spectral_cube import SpectralCube as scfrom astropy.io import fitsimport astropy.units as uimport astropy.constants as cimport matplotlib.pyplot as pltplt.ion()import osspw = ['spw35','spw37']def change_xaxis_in_velo.

2020-11-23 19:40:36 255

原创 line survey后续一：找到对辐射谱有贡献的跃迁

有两个功能没写好。XCLASS的分子线证认已经写得很好，但是后续还要自己做一些比较繁杂的工作，这里放了部分的操作步骤。根据CDMS的database文件和已经证认好分子的molfit文件，这里可以搜寻一下，得出拟合结果大于3sigma的跃迁。这里有个task文件，定义了一些功能。#!/usr/bin/python# -*- coding: utf-8 -*-##********************************************************************

2020-11-23 19:31:07 245

原创画rotational-diagram

import numpy as npimport astropy as aimport astropy.units as uimport astropy.constants as cfrom scipy import optimizeimport matplotlib.pyplot as pltplt.ion()B = 24690.2h = 6.626e-27 #erg/sk = 1.38e-16 #erg/KTbg = 2.73 #Knu = [35

2020-11-15 15:32:23 163

原创画图，subplot多子图

import matplotlib.pyplot as pltimport osimport numpy as npfrom astropy.coordinates import SkyCoordfrom astropy.table import Tablefrom astropy.io import fitsimport warningswarnings.filterwarnings('ignore')T = Table.read('atoms_B3_12M_hotcore_candi.

2020-07-24 09:21:50 257

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

from astropy.io import fitsfrom spectral_cube import SpectralCube as specimport numpy as npimport matplotlib.pyplot as pltimport osimport sysfrom astropy.coordinates import SkyCoordfrom astropy.wcs import WCSfrom astropy.table import Tableimport .

2020-07-17 12:58:36 661

原创批量压缩多个文件夹下的多个文件

移动硬盘坏了，所以还没有验证是否可用import osimport sysimport gzippath="/home/redhat/ATOMS/2019.1.00685.S/science_goal.uid___A001_X1465_X2848/"list1 = os.listdir(path)#便历两级目录dirs = []for i in range(len(list1)): list2 = os.listdir(list1[i]) for j in range(len(l

2020-07-14 18:59:15 553

原创运用shell命令对超大的数据/压缩文件的拆分和合并

使用split命令拆分文件。split -b 4000 -d -a 2 split_file split_file.-b：以文件大小拆分，默认单位M，可以通过修改后数字以及换单位来自由选择拆分后文件大小，比如-b 200K-l：以文件行数拆分-d：拆分后的的文件命令以数字为后缀，split_file.00，split_file.01等等-a：后缀的长度，默认为两位数，可以更改数字。使用cat命令来还原拆分的文件。cat split_file.* split_file.

2020-06-29 14:29:43 918

原创 Herschel-sed-fitter代码

翻到去年写的sed点到点拟合的代码。乱七八糟还没写完，还有很多写得不够智能。建议直接用王科老师的higal_sedfitter.######! /usr/bin python3#This progress is wrote by Lixia Yuan#Version v2.0 modified by Chao Zhang at Dec. 26, 2019#This script requires dust_emissivity which can be installed by "pip in.

2020-06-22 10:31:33 840 4

uycvbl的博客