# 运行信息，rebinning和interval的乘积就是选定的时间区间 $ Name of the window file ('-'for default window) ##默认回车可看到一些时间等信息 $ Newbin Time or negative rebinning 20 ##输入重新并合的时间bin大小 $ Number of Newbins/Interval 50 ##输入时间重新bin后可用的通道数 $ Name of output file ##输入输出文件名或默认回车 $ Do you want to plot your results? ##输入yes或回车 $ Enter PGPLOT device ##默认回车或输入 /xw
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
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 #生成文件
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 #强度分布栏
import pandas as pd import numpy as np from sklearn.neighbors import KNeighborsClassifier from sklearn.model_selection import train_test_split data = pd.read_csv('./iris.csv') print(data.head())
