Astro

Gildas_soft

Gildas是一个具有强大功能的天文数据处理软件包，可以处理单张光谱，批量处理光谱，抑或是输出FITS图像文件

开始前，我们首先进入文件的文件夹作为工作目录

1
cd radio/ #your_folder_name/

然后启动gildas

1
class #直接输入class即可直接进入Gildas环境

Step1.单个谱线数据处理

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
file in G000.31.30m #输入谱线文件，由于G000.31.30m就在当前工作目录，相对路径即可
find #将G000.31.30m文件内所存储的数据输入缓存
list #将缓存中的内容列表输出

set tele *e0*li* #只读取带有e0的字符串，这是读取特定望远镜的数据
find #存入新筛选的数据，替换原来缓存的数据
list #查看
ave #将新选出的数据做平均处理

mod f 97980.953 #选定中心频率

# 设置图标坐标显示及画图
set unit v f #设置x轴显示速度和频率
set unit v #设置x轴仅显示速度
# set unit y t 
set mod x -25 75 #设置坐标起始点和终点
set win -20 70
ba #拟合
plot #绘制谱线，可简写p
mod line C32S(2-1) #更改图像名称
label "T\dM\dB(K)" /Y #给予y轴标签

# 拟合曲线
line 2 #选定成分数量，然后在谱线图上选择对应成分分布范围
min #进行Gauss拟合
view #查看拟合曲线，简写v

# 保存图像
hardcopy C32S21.eps #存储图像为C32S21.eps

效果如图

Step2.批量处理谱线

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
file in G000.31.30m #打开数据文件
find #读入
set tele #设置所有望远镜数据
find #读入
set character 0.6 #设置字体大小
set unit v f #同Step1设置坐标标签
set unit v
set mod y t 
set mod x -50 50 #同Step1
set win -45 45
set tele *e0*li* #选取数据
find
ave #同Step1

# 谱线97980.953
mod f 97980.953 #选定中心谱线
ba
# clear 若之前的图未清除可以使用clear清除
set view 0.3 0.7 0.7 0.9 #设置图像放置的位置，前两项数字为大小，后两项数字为位置
box n o #让box不绘制标签
spec #作图
line 1 #同Step1
min #做Gauss拟合
v /pen 1 #查看同时改变画笔颜色为红色

#谱线97172.063
## 步骤与97980.953相同
mod f 97172.063
ba 
set view 0.3 0.7 0.5 0.7 #注意要将第二幅图往下画
box n o
spec
line 1 
min
v /pen 1

#谱线96412.95
## 步骤与97980.953相同
mod f 96412.95
ba
set view 0.3 0.7 0.3 0.5 #注意要将第三幅图往下画
box n o
spec
line 1
min
ba /pen 1

#谱线92494.308
set tele *e0*lo* #谱线92494.308并不在上述的望远镜数据中，在另一个数据里，要进行切换
find
ave
list
ave
mod f 92494.308
ba
set view 0.3 0.7 0.1 0.3 #最后一幅图像
box #由于是最后一幅图，所以不需要隐藏坐标标签
line 1
min
v /pen 1

#补充说明
set character 0.6
draw text 10 18 G000.31 # 总谱线名称
draw text 17 18 "CS(2-1)" # 谱线97980.953名称
draw text 17 13.7 "C\u3\u3S(2-1)" # 谱线97172.063名称
draw text 17 9.5 "C\u3\u4S(2-1)" # 谱线96412.95名称
draw text 17 5.2 "\u1\u3CS(2-1)" # 谱线92494.308名称
set view 0.3 0.35 0.48 0.58
label "T\dM\dB(K)" /Y #y轴坐标标签
\sic dele g000.31.eps #若先前已经有这个文件则删除
hardcopy g000.31.eps #生成文件

效果如图

Step3.FITS彩图绘制

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
fil out test.bur single #创建一个test.bur
lmv OrionA_20171031L.fits #输入fits文件并转换
file in test.bur #将转换后的fits文件数据读入test.bur
list
find
list
set off 0 0 #设置中心位置
find
set unit v f
get f #选取中心位置的第一个scan
p
set mod x -50 50
p
set win -35 35
ba 1;p #基线拟合，数字表示基线拟合阶数
set off
set ang min
find
file out test.bas s #创建批量保存的数据文件
find
get n
for i 1 to found #对数据进行遍历
get n
bas 1 #即ba 1
write #写入数据
next #退出循环
# las\get n

#谱线绘制
file in test.bas
find
list
set mod y -0.2 1.5 #看谱线定值，X或Y必须给一个值，一般选择Y
clear #清除先前绘制查看的全谱线图
map /gr #画图

#输出data文件
file in test.bas
find
list
sic dele test.dat
print area -35 35 /output test.dat
exit #退出class

greg #进入greg环境，角分画图，使用相对路径
col x 2 y 3 z 4 /file test.dat #选定坐标轴
set box m #定义画框
lev 1 to 10 #强度划分
lim /rev x
rand 100 /blank -1 -1
box
lab "\gDL (arcmin)" /x #x轴标签
lab "\gDB (arcmin)" /y #y轴标签
set mar 3 2 0.2
plot /scalin lin -2 10
rgmap /per 10 #上色
wedge #强度分布栏

hardcopy H(41).eps /dev ps color #输出图像文件

效果如图