python 导入数据及作图的实现


Posted in Python onDecember 03, 2019

我们经常需要导入数据,按列提取 XY作图

方法一、

filename='/home/res/user/csluo/test.txt'          #将文件名赋值为变量
 
X,Y,Z=[ ],[ ],[ ]                      #给三个空列表
 
with open(filename, 'r') as f:               #打开文档
  lines=f.readlines()                   #按行读取数据
 
  for i in lines:
 
     value=[ float(s) for s in line.split( ) ]   #将 行数据 以空格分开,并以浮点型赋给value 
 
     X=append(value[0])              #将value 的第一个数值 添加到 列表X
 
     Y=append(value[1])
 
     Z=append(value[2])

方法二、

a = numpy.loadtxt('a.txt')
 
a[:,0]

方法三、

>>> import pandas as pd
>>> data = pd.read_csv("aaa.csv")

附录:IRAS16293 的数据 绘图代码

import matplotlib
import matplotlib.pyplot as plt
import numpy as np
 
##################################################################
######################################################################################### ######## 1
##################################################################
 
 
#data_a
filename = '/home/luoluo/Documents/python_scripts/Adata_cut/A1_145122-145580.txt'
X10,Y10= [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]     #以空格分开 行 并输出 浮点数
    X10.append(value[0])
    Y10.append(value[1])
 
#c2h5oh
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A1/xclass_spectrum_output.dat'
X11,Y11 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X11.append(value[0])
    Y11.append(value[1])
 
#ch3coch3
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A1/xclass_spectrum_output.dat'
X12,Y12 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X12.append(value[0])
    Y12.append(value[1])
 
 
#ch3ocho_v18=1
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A1/xclass_spectrum_output.dat'
X13,Y13 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X13.append(value[0])
    Y13.append(value[1])
 
 
#ch3ocho_v=0
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A1/xclass_spectrum_output.dat'
X14,Y14 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X14.append(value[0])
    Y14.append(value[1])
 
###################################################  data_b
 
filename = '/home/luoluo/Documents/python_scripts/Bdata_cut/B1_145122-145580.txt'
X15,Y15= [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X15.append(value[0])
    Y15.append(value[1])
 
#c2h5oh
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B1/xclass_spectrum_output.dat'
X16,Y16 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X16.append(value[0])
    Y16.append(value[1])
 
#ch3coch3
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B1/xclass_spectrum_output.dat'
X17,Y17 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X17.append(value[0])
    Y17.append(value[1])
 
#ch3ocho_v18=1
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B1/xclass_spectrum_output.dat'
X18,Y18 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X18.append(value[0])
    Y18.append(value[1])
 
#ch3ocho_v=0
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B1/xclass_spectrum_output.dat'
X19,Y19 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X19.append(value[0])
    Y19.append(value[1])
 
 
 
###################################################################  #2 plot
 
 
 
plt.figure(1)
plt.subplot(211)
plt.axis([145122,145580, -1, 15])
 
new_ticks = np.linspace(145122,145580,10)
plt.xticks(new_ticks)
 
matplotlib.rcParams['xtick.direction'] = 'in'
matplotlib.rcParams['ytick.direction'] = 'in'
plt.plot(X10,Y10,color='k',lw=1.0,ls='-',label='DATA')  
plt.plot(X11,Y11,color='r',lw=1.0,ls='-',label='C2H5OH')
plt.plot(X12,Y12,color='g',lw=1.0,ls='-',label='CH3COCH3')
plt.plot(X13,Y13,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')
plt.plot(X14,Y14,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')
plt.legend(loc='upper left')
#plt.xlabel('Freq(MHZ)')
plt.ylabel('T(K)')
plt.text(145535,13,'IRAS16293A')
 
 
#plt.show()
#plt.savefig('8.png')
 
 
plt.subplot(212)
plt.axis([145122,145580, -1, 15])
 
new_ticks = np.linspace(145122,145580,10)
plt.xticks(new_ticks)
 
plt.plot(X15,Y15,color='k',lw=1.0,ls='-',label='DATA')  
plt.plot(X16,Y16,color='r',lw=1.0,ls='-',label='C2H5OH')
plt.plot(X17,Y17,color='g',lw=1.0,ls='-',label='CH3COCH3')
plt.plot(X18,Y18,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')
plt.plot(X19,Y19,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')
plt.legend(loc='upper left')
plt.xlabel('Freq(MHZ)')
plt.ylabel('T(K)')
plt.text(145535,13,'IRAS16293B')
 
plt.plot([145263,145263],[1,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("C2H5OH", xy = (145263,1), xytext = (145280,9),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.show()
 
##################################################################
######################################################################################### ######## 2
##################################################################
 
 
#data_a
filename = '/home/luoluo/Documents/python_scripts/Adata_cut/A2_145580-146050.txt'
X20,Y20= [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X20.append(value[0])
    Y20.append(value[1])
 
#c2h5oh
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A2/xclass_spectrum_output.dat'
X21,Y21 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X21.append(value[0])
    Y21.append(value[1])
 
#ch3coch3
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A2/xclass_spectrum_output.dat'
X22,Y22 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X22.append(value[0])
    Y22.append(value[1])
 
 
#ch3ocho_v18=1
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A2/xclass_spectrum_output.dat'
X23,Y23 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X23.append(value[0])
    Y23.append(value[1])
 
 
#ch3ocho_v=0
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A2/xclass_spectrum_output.dat'
X24,Y24 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X24.append(value[0])
    Y24.append(value[1])
 
###################################################  data_b
 
filename = '/home/luoluo/Documents/python_scripts/Bdata_cut/B2.txt'
X25,Y25= [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X25.append(value[0])
    Y25.append(value[1])
 
#c2h5oh
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B2/xclass_spectrum_output.dat'
X26,Y26 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X26.append(value[0])
    Y26.append(value[1])
 
#ch3coch3
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B2/xclass_spectrum_output.dat'
X27,Y27 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X27.append(value[0])
    Y27.append(value[1])
 
#ch3ocho_v18=1
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B2/xclass_spectrum_output.dat'
X28,Y28 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X28.append(value[0])
    Y28.append(value[1])
 
#ch3ocho_v=0
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B2/xclass_spectrum_output.dat'
X29,Y29 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X29.append(value[0])
    Y29.append(value[1])
 
 
 
###################################################################  #2 plot
 
 
 
plt.figure(2)
plt.subplot(211)
plt.axis([145580,146050, -1, 15])
 
new_ticks = np.linspace(145580,146050,10)
plt.xticks(new_ticks)
 
matplotlib.rcParams['xtick.direction'] = 'in'
matplotlib.rcParams['ytick.direction'] = 'in'
plt.plot(X20,Y20,color='k',lw=1.0,ls='-',label='DATA')  
plt.plot(X21,Y21,color='r',lw=1.0,ls='-',label='C2H5OH')
plt.plot(X22,Y22,color='g',lw=1.0,ls='-',label='CH3COCH3')
plt.plot(X23,Y23,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')
plt.plot(X24,Y24,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')
plt.legend(loc='upper left')
#plt.xlabel('Freq(MHZ)')
plt.ylabel('T(K)')
plt.text(146000,13,'IRAS16293A')
 
plt.plot([145618.5,145618.5],[2.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (145618,7), xytext = (145650,11),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([146020,146020],[7,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (146020,6), xytext = (146005,10),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
 
 
#plt.show()
#plt.savefig('8.png')
 
 
plt.subplot(212)
plt.axis([145580,146050, -1, 15])
 
new_ticks = np.linspace(145580,146050,10)
plt.xticks(new_ticks)
 
plt.plot(X25,Y25,color='k',lw=1.0,ls='-',label='DATA')  
plt.plot(X26,Y26,color='r',lw=1.0,ls='-',label='C2H5OH')
plt.plot(X27,Y27,color='g',lw=1.0,ls='-',label='CH3COCH3')
plt.plot(X28,Y28,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')
plt.plot(X29,Y29,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')
plt.legend(loc='upper left')
plt.xlabel('Freq(MHZ)')
plt.ylabel('T(K)')
plt.text(146000,13,'IRAS16293B')
 
plt.plot([145618.5,145618.5],[2.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (145618,7), xytext = (145650,11),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([146020,146020],[7,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (146020,6), xytext = (145990,10),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.show()
 
 
##################################################################
######################################################################################### ######## 3
##################################################################
 
 
#data_a
filename = '/home/luoluo/Documents/python_scripts/Adata_cut/A3.txt'
X30,Y30= [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X30.append(value[0])
    Y30.append(value[1])
 
#c2h5oh
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A3/xclass_spectrum_output.dat'
X31,Y31 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X31.append(value[0])
    Y31.append(value[1])
 
#ch3coch3
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A3/xclass_spectrum_output.dat'
X32,Y32 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X32.append(value[0])
    Y32.append(value[1])
 
 
#ch3ocho_v18=1
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A3/xclass_spectrum_output.dat'
X33,Y33 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X33.append(value[0])
    Y33.append(value[1])
 
 
#ch3ocho_v=0
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A3/xclass_spectrum_output.dat'
X34,Y34 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X34.append(value[0])
    Y34.append(value[1])
 
###################################################  data_b
 
filename = '/home/luoluo/Documents/python_scripts/Bdata_cut/B3_146050-146507.txt'
X35,Y35= [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X35.append(value[0])
    Y35.append(value[1])
 
#c2h5oh
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B3/xclass_spectrum_output.dat'
X36,Y36 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X36.append(value[0])
    Y36.append(value[1])
 
#ch3coch3
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B3/xclass_spectrum_output.dat'
X37,Y37 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X37.append(value[0])
    Y37.append(value[1])
 
#ch3ocho_v18=1
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B3/xclass_spectrum_output.dat'
X38,Y38 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X38.append(value[0])
    Y38.append(value[1])
 
#ch3ocho_v=0
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B3/xclass_spectrum_output.dat'
X39,Y39 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X39.append(value[0])
    Y39.append(value[1])
 
 
 
###################################################################  #2 plot
 
 
 
plt.figure(3)
plt.subplot(211)
plt.axis([146050,146507, -1, 15])
 
new_ticks = np.linspace(146050,146507,10)
plt.xticks(new_ticks)
 
matplotlib.rcParams['xtick.direction'] = 'in'
matplotlib.rcParams['ytick.direction'] = 'in'
plt.plot(X30,Y30,color='k',lw=1.0,ls='-',label='DATA')  
plt.plot(X31,Y31,color='r',lw=1.0,ls='-',label='C2H5OH')
plt.plot(X32,Y32,color='g',lw=1.0,ls='-',label='CH3COCH3')
plt.plot(X33,Y33,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')
plt.plot(X34,Y34,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')
plt.legend(loc='upper left')
#plt.xlabel('Freq(MHZ)')
plt.ylabel('T(K)')
plt.text(146465,13,'IRAS16293A')
 
plt.plot([146172,146172],[2.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (146172,7), xytext = (146190,11),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([146232,146232],[3,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (146232,10), xytext = (146210,5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([146343,146343],[2.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (146343,7), xytext = (146320,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([146355,146355],[2.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (146355,6), xytext = (146380,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
 
#plt.show()
#plt.savefig('8.png')
 
 
plt.subplot(212)
plt.axis([146050,146507, -1, 15])
 
new_ticks = np.linspace(146050,146507,10)
plt.xticks(new_ticks)
 
plt.plot(X35,Y35,color='k',lw=1.0,ls='-',label='DATA')  
plt.plot(X36,Y36,color='r',lw=1.0,ls='-',label='C2H5OH')
plt.plot(X37,Y37,color='g',lw=1.0,ls='-',label='CH3COCH3')
plt.plot(X38,Y38,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')
plt.plot(X39,Y39,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')
plt.legend(loc='upper left')
plt.xlabel('Freq(MHZ)')
plt.ylabel('T(K)')
plt.text(146465,13,'IRAS16293B')
 
plt.plot([146343,146343],[5.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (146343,7), xytext = (146320,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([146355,146355],[5.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (146355,6), xytext = (146380,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
#plt.plot([157324,157324],[2.5,15], color ='grey', linewidth=1, linestyle="--")
#plt.annotate("CH3COCH3", xy = (157324,3), xytext = (157347,4.5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.show()
 
 
##################################################################
######################################################################################### ######## 4
##################################################################
 
 
#data_a
filename = '/home/luoluo/Documents/python_scripts/Adata_cut/A4'
X40,Y40= [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X40.append(value[0])
    Y40.append(value[1])
 
#c2h5oh
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A4/xclass_spectrum_output.dat'
X41,Y41 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X41.append(value[0])
    Y41.append(value[1])
 
#ch3coch3
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A4/xclass_spectrum_output.dat'
X42,Y42 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X42.append(value[0])
    Y42.append(value[1])
 
 
#ch3ocho_v18=1
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A4/xclass_spectrum_output.dat'
X43,Y43 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X43.append(value[0])
    Y43.append(value[1])
 
 
#ch3ocho_v=0
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A4/xclass_spectrum_output.dat'
X44,Y44 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X44.append(value[0])
    Y44.append(value[1])
 
###################################################  data_b
 
filename = '/home/luoluo/Documents/python_scripts/Bdata_cut/B4_146507-146976.txt'
X45,Y45= [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X45.append(value[0])
    Y45.append(value[1])
 
#c2h5oh
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B4/xclass_spectrum_output.dat'
X46,Y46 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X46.append(value[0])
    Y46.append(value[1])
 
#ch3coch3
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B4/xclass_spectrum_output.dat'
X47,Y47 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X47.append(value[0])
    Y47.append(value[1])
 
#ch3ocho_v18=1
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B4/xclass_spectrum_output.dat'
X48,Y48 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X48.append(value[0])
    Y48.append(value[1])
 
#ch3ocho_v=0
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B4/xclass_spectrum_output.dat'
X49,Y49 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X49.append(value[0])
    Y49.append(value[1])
 
 
 
###################################################################  #2 plot
 
 
 
plt.figure(4)
plt.subplot(211)
plt.axis([146507,146976, -1, 15])
 
new_ticks = np.linspace(146507,146976,10)
plt.xticks(new_ticks)
 
matplotlib.rcParams['xtick.direction'] = 'in'
matplotlib.rcParams['ytick.direction'] = 'in'
plt.plot(X40,Y40,color='k',lw=1.0,ls='-',label='DATA')  
plt.plot(X41,Y41,color='r',lw=1.0,ls='-',label='C2H5OH')
plt.plot(X42,Y42,color='g',lw=1.0,ls='-',label='CH3COCH3')
plt.plot(X43,Y43,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')
plt.plot(X44,Y44,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')
plt.legend(loc='upper left')
#plt.xlabel('Freq(MHZ)')
plt.ylabel('T(K)')
plt.text(146930,13,'IRAS16293A')
 
#plt.plot([157323,157323],[2.5,15], color ='grey', linewidth=1, linestyle="--")
#plt.annotate("CH3COCH3", xy = (157323,3), xytext = (157345,4.5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
 
 
#plt.show()
#plt.savefig('8.png')
 
 
plt.subplot(212)
plt.axis([146507,146976, -1, 15])
 
new_ticks = np.linspace(146507,146976,10)
plt.xticks(new_ticks)
 
plt.plot(X45,Y45,color='k',lw=1.0,ls='-',label='DATA')  
plt.plot(X46,Y46,color='r',lw=1.0,ls='-',label='C2H5OH')
plt.plot(X47,Y47,color='g',lw=1.0,ls='-',label='CH3COCH3')
plt.plot(X48,Y48,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')
plt.plot(X49,Y49,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')
plt.legend(loc='upper left')
plt.xlabel('Freq(MHZ)')
plt.ylabel('T(K)')
plt.text(146930,13,'IRAS16293B')
 
#plt.plot([157324,157324],[2.5,15], color ='grey', linewidth=1, linestyle="--")
#plt.annotate("CH3COCH3", xy = (157324,3), xytext = (157347,4.5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.show()
 
##################################################################
######################################################################################### ######## 5
##################################################################
 
#data_a
filename = '/home/luoluo/Documents/python_scripts/Adata_cut/A5_156840-157303.txt'
X50,Y50= [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X50.append(value[0])
    Y50.append(value[1])
 
 
#c2h5oh
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A5/xclass_spectrum_output.dat'
X51,Y51 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X51.append(value[0])
    Y51.append(value[1])
 
#ch3coch3
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A5/xclass_spectrum_output.dat'
X52,Y52 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X52.append(value[0])
    Y52.append(value[1])
 
 
#ch3ocho_v18=1
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A5/xclass_spectrum_output.dat'
X53,Y53 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X53.append(value[0])
    Y53.append(value[1])
 
 
#ch3ocho_v=0
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A5/xclass_spectrum_output.dat'
X54,Y54 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X54.append(value[0])
    Y54.append(value[1])
 
###################################################  data_b
 
filename = '/home/luoluo/Documents/python_scripts/Bdata_cut/B5_156840-157303.txt'
X55,Y55= [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X55.append(value[0])
    Y55.append(value[1])
 
#c2h5oh
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B5/xclass_spectrum_output.dat'
X56,Y56 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X56.append(value[0])
    Y56.append(value[1])
 
#ch3coch3
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B5/xclass_spectrum_output.dat'
X57,Y57 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X57.append(value[0])
    Y57.append(value[1])
 
#ch3ocho_v18=1
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B5/xclass_spectrum_output.dat'
X58,Y58 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X58.append(value[0])
    Y58.append(value[1])
 
#ch3ocho_v=0
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B5/xclass_spectrum_output.dat'
X59,Y59 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X59.append(value[0])
    Y59.append(value[1])
 
 
 
###################################################################  #2 plot
 
 
 
plt.figure(5)
plt.subplot(211)
plt.axis([156840, 157303, -1, 15])
 
new_ticks = np.linspace(156840,157303,10)
plt.xticks(new_ticks)
 
matplotlib.rcParams['xtick.direction'] = 'in'
matplotlib.rcParams['ytick.direction'] = 'in'
plt.plot(X50,Y50,color='k',lw=1.0,ls='-',label='DATA')  
plt.plot(X51,Y51,color='r',lw=1.0,ls='-',label='C2H5OH')
plt.plot(X52,Y52,color='g',lw=1.0,ls='-',label='CH3COCH3')
plt.plot(X53,Y53,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')
plt.plot(X54,Y54,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')
plt.legend(loc='upper left')
#plt.xlabel('Freq(MHZ)')
plt.ylabel('T(K)')
plt.text(157256,13,'IRAS16293A')
 
#plt.show()
#plt.savefig('8.png')
 
 
plt.subplot(212)
plt.axis([156820,157303, -1, 15])
 
new_ticks = np.linspace(156820,157303,10)
plt.xticks(new_ticks)
 
plt.plot(X55,Y55,color='k',lw=1.0,ls='-',label='DATA')  
plt.plot(X56,Y56,color='r',lw=1.0,ls='-',label='C2H5OH')
plt.plot(X57,Y57,color='g',lw=1.0,ls='-',label='CH3COCH3')
plt.plot(X58,Y58,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')
plt.plot(X59,Y59,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')
plt.legend(loc='upper left')
plt.xlabel('Freq(MHZ)')
plt.ylabel('T(K)')
plt.text(157255,13,'IRAS16293B')
 
#plt.plot([157324,157324],[2.5,15], color ='grey', linewidth=1, linestyle="--")
#plt.annotate("CH3COCH3", xy = (157324,3), xytext = (157347,4.5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.show()
 
##################################################################
######################################################################################### ######## 6
##################################################################
 
 
#data_a
filename = '/home/luoluo/Documents/python_scripts/Adata_cut/A6_157303-157820.txt'
X60,Y60= [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X60.append(value[0])
    Y60.append(value[1])
 
 
#c2h5oh
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A6/xclass_spectrum_output.dat'
X61,Y61 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X61.append(value[0])
    Y61.append(value[1])
 
#ch3coch3
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A6/xclass_spectrum_output.dat'
X62,Y62 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X62.append(value[0])
    Y62.append(value[1])
 
 
#ch3ocho_v18=1
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A6/xclass_spectrum_output.dat'
X63,Y63 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X63.append(value[0])
    Y63.append(value[1])
 
 
#ch3ocho_v=0
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A6/xclass_spectrum_output.dat'
X64,Y64 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X64.append(value[0])
    Y64.append(value[1])
 
###################################################  data_b
 
filename = '/home/luoluo/Documents/python_scripts/Bdata_cut/B6_157303-157820.txt'
X65,Y65= [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X65.append(value[0])
    Y65.append(value[1])
 
#c2h5oh
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B6/xclass_spectrum_output.dat'
X66,Y66 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X66.append(value[0])
    Y66.append(value[1])
 
#ch3coch3
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B6/xclass_spectrum_output.dat'
X67,Y67 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X67.append(value[0])
    Y67.append(value[1])
 
#ch3ocho_v18=1
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B6/xclass_spectrum_output.dat'
X68,Y68 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X68.append(value[0])
    Y68.append(value[1])
 
#ch3ocho_v=0
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B6/xclass_spectrum_output.dat'
X69,Y69 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X69.append(value[0])
    Y69.append(value[1])
 
 
 
###################################################################  #2 plot
 
 
 
plt.figure(6)
plt.subplot(211)
plt.axis([157303, 157820, -1, 15])
 
new_ticks = np.linspace(157303, 157820,10)
plt.xticks(new_ticks)
 
matplotlib.rcParams['xtick.direction'] = 'in'
matplotlib.rcParams['ytick.direction'] = 'in'
plt.plot(X60,Y60,color='k',lw=1.0,ls='-',label='DATA')  
plt.plot(X61,Y61,color='r',lw=1.0,ls='-',label='C2H5OH')
plt.plot(X62,Y62,color='g',lw=1.0,ls='-',label='CH3COCH3')
plt.plot(X63,Y63,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')
plt.plot(X64,Y64,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')
plt.legend(loc='upper left')
#plt.xlabel('Freq(MHZ)')
plt.ylabel('T(K)')
plt.text(157765,13,'IRAS16293A')
 
plt.plot([157323,157323],[2.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3COCH3", xy = (157323,3), xytext = (157345,4.5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
 
 
#plt.show()
#plt.savefig('8.png')
 
 
plt.subplot(212)
plt.axis([157303, 157820, -1, 15])
 
new_ticks = np.linspace(157303, 157820,10)
plt.xticks(new_ticks)
 
plt.plot(X65,Y65,color='k',lw=1.0,ls='-',label='DATA')  
plt.plot(X66,Y66,color='r',lw=1.0,ls='-',label='C2H5OH')
plt.plot(X67,Y67,color='g',lw=1.0,ls='-',label='CH3COCH3')
plt.plot(X68,Y68,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')
plt.plot(X69,Y69,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')
plt.legend(loc='upper left')
plt.xlabel('Freq(MHZ)')
plt.ylabel('T(K)')
plt.text(157765,13,'IRAS16293B')
 
plt.plot([157324,157324],[2.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3COCH3", xy = (157324,3), xytext = (157347,4.5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([157505.5,157505.5],[5.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("C2H5OH", xy = (157505,6), xytext = (157450,11),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.show()
 
###############################################################################################
################################################################################################################ 7
###############################################################################################
 
 
#data_a
filename = '/home/luoluo/Documents/python_scripts/Adata_cut/A7_157820-158228.txt'
X70,Y70= [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X70.append(value[0])
    Y70.append(value[1])
 
 
#c2h5oh
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A7/xclass_spectrum_output.dat'
X71,Y71 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X71.append(value[0])
    Y71.append(value[1])
 
#ch3coch3
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A7/xclass_spectrum_output.dat'
X72,Y72 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X72.append(value[0])
    Y72.append(value[1])
 
 
#ch3ocho_v18=1
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A7/xclass_spectrum_output.dat'
X73,Y73 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X73.append(value[0])
    Y73.append(value[1])
 
 
#ch3ocho_v=0
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A7/xclass_spectrum_output.dat'
X74,Y74 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X74.append(value[0])
    Y74.append(value[1])
#########################################################
#data_b
filename = '/home/luoluo/Documents/python_scripts/Bdata_cut/B7_157820-158228.txt'
X75,Y75= [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X75.append(value[0])
    Y75.append(value[1])
 
#c2h5oh
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B7/xclass_spectrum_output.dat'
X76,Y76 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X76.append(value[0])
    Y76.append(value[1])
 
#ch3coch3
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B7/xclass_spectrum_output.dat'
X77,Y77 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X77.append(value[0])
    Y77.append(value[1])
 
#ch3ocho_v18=1
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B7/xclass_spectrum_output.dat'
X78,Y78 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X78.append(value[0])
    Y78.append(value[1])
 
#ch3ocho_v=0
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B7/xclass_spectrum_output.dat'
X79,Y79 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X79.append(value[0])
    Y79.append(value[1])
#############################################################################################################################################
#2 plot
plt.figure(7)
plt.subplot(211)
plt.axis([157820, 158228, -1, 15])
 
new_ticks = np.linspace(157820,158228,10)
plt.xticks(new_ticks)
 
matplotlib.rcParams['xtick.direction'] = 'in'
matplotlib.rcParams['ytick.direction'] = 'in'
plt.plot(X70,Y70,color='k',lw=1.0,ls='-',label='DATA')  
plt.plot(X71,Y71,color='r',lw=1.0,ls='-',label='C2H5OH')
plt.plot(X72,Y72,color='g',lw=1.0,ls='-',label='CH3COCH3')
plt.plot(X73,Y73,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')
plt.plot(X74,Y74,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')
plt.legend(loc='upper left')
#plt.xlabel('Freq(MHZ)')
plt.ylabel('T(K)')
plt.text(158183,13,'IRAS16293A')
 
plt.plot([157894,157894],[5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (157895,7), xytext = (157880,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([158120,158120],[5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3COCH3", xy = (158120,5), xytext = (158110,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([158172,158172],[7,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3COCH3", xy = (158171,8), xytext = (158150,10),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([158225,158225],[4,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3COCH3", xy = (158225,5), xytext = (158200,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
#plt.show()
#plt.savefig('8.png')
 
 
plt.subplot(212)
plt.axis([157820, 158228, -1, 15])
 
new_ticks = np.linspace(157820, 158228,10)
plt.xticks(new_ticks)
 
plt.plot(X75,Y75,color='k',lw=1.0,ls='-',label='DATA')  
plt.plot(X76,Y76,color='r',lw=1.0,ls='-',label='C2H5OH')
plt.plot(X77,Y77,color='g',lw=1.0,ls='-',label='CH3COCH3')
plt.plot(X78,Y78,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')
plt.plot(X79,Y79,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')
plt.legend(loc='upper left')
plt.xlabel('Freq(MHZ)')
plt.ylabel('T(K)')
plt.text(158183,13,'IRAS16293B')
 
plt.plot([157895,157895],[5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (157896,7), xytext = (157880,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([158121,158121],[5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3COCH3", xy = (158121,7), xytext = (158110,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([158173,158173],[7,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3COCH3", xy = (158172,8), xytext = (158150,5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([158226,158226],[4,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3COCH3", xy = (158225,5), xytext = (158205,10),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
#plt.plot([158389,158389],[3.5,15], color ='grey', linewidth=1, linestyle="--")
#plt.annotate("CH3OCHO v18=1", xy = (158390,5), xytext = (158400,13),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
 
plt.show()
######################################################################################
###################################################################################### ###### 8
######################################################################################
 
#
#1 import data
 
#data_a
filename = 'Adata_cut/A8_158228-158680.txt'
X80,Y80= [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X80.append(value[0])
    Y80.append(value[1])
 
 
#c2h5oh
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A8/xclass_spectrum_output.dat'
X81,Y81 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X81.append(value[0])
    Y81.append(value[1])
 
#ch3coch3
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A8/xclass_spectrum_output.dat'
X82,Y82 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X82.append(value[0])
    Y82.append(value[1])
 
 
#ch3ocho_v18=1
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A8/xclass_spectrum_output.dat'
X83,Y83 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X83.append(value[0])
    Y83.append(value[1])
 
 
#ch3ocho_v=0
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A8/xclass_spectrum_output.dat'
X84,Y84 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X84.append(value[0])
    Y84.append(value[1])
#########################################################
#data_b
filename = 'Bdata_cut/B8_158228-158680.txt'
X85,Y85= [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X85.append(value[0])
    Y85.append(value[1])
 
#c2h5oh
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B8/xclass_spectrum_output.dat'
X86,Y86 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X86.append(value[0])
    Y86.append(value[1])
 
#ch3coch3
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B8/xclass_spectrum_output.dat'
X87,Y87 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X87.append(value[0])
    Y87.append(value[1])
 
#ch3ocho_v18=1
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B8/xclass_spectrum_output.dat'
X88,Y88 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X88.append(value[0])
    Y88.append(value[1])
 
#ch3ocho_v=0
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B8/xclass_spectrum_output.dat'
X89,Y89 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X89.append(value[0])
    Y89.append(value[1])
#############################################################
#2 plot
plt.figure(8)
plt.subplot(211)
plt.axis([158228, 158680, -1, 15])
 
new_ticks = np.linspace(158228,158680,10)
plt.xticks(new_ticks)
 
matplotlib.rcParams['xtick.direction'] = 'in'
matplotlib.rcParams['ytick.direction'] = 'in'
plt.plot(X80,Y80,color='k',lw=1.0,ls='-',label='DATA')  
plt.plot(X81,Y81,color='r',lw=1.0,ls='-',label='C2H5OH')
plt.plot(X82,Y82,color='g',lw=1.0,ls='-',label='CH3COCH3')
plt.plot(X83,Y83,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')
plt.plot(X84,Y84,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')
plt.legend(loc='upper left')
#plt.xlabel('Freq(MHZ)')
plt.ylabel('T(K)')
plt.text(158635,13,'IRAS16293A')
 
plt.plot([158378,158378],[1.9,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158378, 2.0), xytext = (158300,5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([158388,158388],[3.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158388, 4.5), xytext = (158370,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([158430,158430],[3.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158430, 3.5), xytext = (158400,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([158484,158484],[3.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158485, 4), xytext = (158450,9),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([158519,158519],[3.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158518,4), xytext = (158490,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([158540,158540],[3.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158541,4), xytext = (158550,11),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([158675,158675],[3.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158674,3), xytext = (158610,8),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
 
#plt.show()
#plt.savefig('8.png')
 
 
plt.subplot(212)
plt.axis([158228, 158680, -1, 15])
 
new_ticks = np.linspace(158228,158680,10)
plt.xticks(new_ticks)
 
plt.plot(X85,Y85,color='k',lw=1.0,ls='-',label='DATA')  
plt.plot(X86,Y86,color='r',lw=1.0,ls='-',label='C2H5OH')
plt.plot(X87,Y87,color='g',lw=1.0,ls='-',label='CH3COCH3')
plt.plot(X88,Y88,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')
plt.plot(X89,Y89,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')
plt.legend(loc='upper left')
plt.xlabel('Freq(MHZ)')
plt.ylabel('T(K)')
plt.text(158635,13,'IRAS16293B')
plt.plot([158389,158389],[3.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158390,5), xytext = (158400,13),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([158484,158484],[3.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158484,5), xytext = (158430,4),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([158520,158520],[3.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158520,4), xytext = (158500,6),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([158540,158540],[3.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158540,4), xytext = (158570,11),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([158675,158675],[3.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158675,4), xytext = (158640,8),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.show()
 
####################################################################################################
#################################################################################################### ######## 9
####################################################################################################
 
#data_a
 
filename = '/home/luoluo/Documents/python_scripts/Adata_cut/A9_158680-159155.txt'
X90,Y90= [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X90.append(value[0])
    Y90.append(value[1])
 
 
#c2h5oh
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A9/xclass_spectrum_output.dat'
X91,Y91 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X91.append(value[0])
    Y91.append(value[1])
 
#ch3coch3
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A9/xclass_spectrum_output.dat'
X92,Y92 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X92.append(value[0])
    Y92.append(value[1])
 
 
#ch3ocho_v18=1
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A9/xclass_spectrum_output.dat'
X93,Y93 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X93.append(value[0])
    Y93.append(value[1])
 
 
#ch3ocho_v=0
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A9/xclass_spectrum_output.dat'
X94,Y94 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X94.append(value[0])
    Y94.append(value[1])
#########################################################
#data_b
filename = '/home/luoluo/Documents/python_scripts/Bdata_cut/B9_158680-159155.txt'
X95,Y95= [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X95.append(value[0])
    Y95.append(value[1])
 
#c2h5oh
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B9/xclass_spectrum_output.dat'
X96,Y96 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X96.append(value[0])
    Y96.append(value[1])
 
#ch3coch3
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B9/xclass_spectrum_output.dat'
X97,Y97 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X97.append(value[0])
    Y97.append(value[1])
 
#ch3ocho_v18=1
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B9/xclass_spectrum_output.dat'
X98,Y98 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X98.append(value[0])
    Y98.append(value[1])
 
#ch3ocho_v=0
filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B9/xclass_spectrum_output.dat'
X99,Y99 = [],[]
with open(filename, 'r') as f:
  lines = f.readlines()
  for line in lines:
    value = [float(s) for s in line.split()]
    X99.append(value[0])
    Y99.append(value[1])
 
 
 
 
 
###2 plot
plt.figure(9)
#A
plt.subplot(211)
plt.axis([158680, 159155, -1, 15])
 
new_ticks = np.linspace(158680,159155,10)
plt.xticks(new_ticks)
 
matplotlib.rcParams['xtick.direction'] = 'in'
matplotlib.rcParams['ytick.direction'] = 'in'
plt.plot(X90,Y90,color='k',lw=1.0,ls='-',label='DATA')  
plt.plot(X91,Y91,color='r',lw=1.0,ls='-',label='C2H5OH')
plt.plot(X92,Y92,color='g',lw=1.0,ls='-',label='CH3COCH3')
plt.plot(X93,Y93,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')
plt.plot(X94,Y94,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')
plt.legend(loc='upper left')
#plt.xlabel('Freq(MHZ)')
plt.ylabel('T(K)')
plt.text(159110,13,'IRAS16293A')
 
plt.plot([158691,158691],[8,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158691, 5), xytext = (158730,3),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([158701.9,158701.9],[8,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158701, 10), xytext = (158770,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([159003.5,159003.5],[3.4,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (159005, 3), xytext = (159020,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([159106,159106],[3.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (159107, 3.5), xytext = (159110,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
#B
plt.subplot(212)
plt.axis([158680, 159155, -1, 15])
 
new_ticks = np.linspace(158680,159122,10)
plt.xticks(new_ticks)
 
plt.plot(X95,Y95,color='k',lw=1.0,ls='-',label='DATA')  
plt.plot(X96,Y96,color='r',lw=1.0,ls='-',label='C2H5OH')
plt.plot(X97,Y97,color='g',lw=1.0,ls='-',label='CH3COCH3')
plt.plot(X98,Y98,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')
plt.plot(X99,Y99,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')
plt.legend(loc='upper left')
plt.xlabel('Freq(MHZ)')
plt.ylabel('T(K)')
plt.text(159110,13,'IRAS16293B')
 
plt.plot([158689,158689],[7,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (158689, 3.5), xytext = (158700,4),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
#plt.plot([158963,158963],[2.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v=0", xy = (158694,6), xytext = (158760,9),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
#plt.plot([158963,158963],[2.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v=0", xy = (158704,9.5), xytext = (158770,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([158963,158963],[2.5,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3COCH3", xy = (158963, 3.5), xytext = (158980,6),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.plot([159030,159030],[2.4,15], color ='grey', linewidth=1, linestyle="--")
plt.annotate("CH3OCHO v18=1", xy = (159030,4), xytext = (159000,9),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))
 
plt.show()

以上这篇python 导入数据及作图的实现就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
python中string模块各属性以及函数的用法介绍
May 30 Python
Python matplotlib画图实例之绘制拥有彩条的图表
Dec 28 Python
TensorFlow实现Batch Normalization
Mar 08 Python
Flask教程之重定向与错误处理实例分析
Aug 01 Python
win10下安装Anaconda的教程(python环境+jupyter_notebook)
Oct 23 Python
python pyenv多版本管理工具的使用
Dec 23 Python
TensorFlow基本的常量、变量和运算操作详解
Feb 03 Python
基于pygame实现童年掌机打砖块游戏
Feb 25 Python
Python通过4种方式实现进程数据通信
Mar 12 Python
python实现按键精灵找色点击功能教程,使用pywin32和Pillow库
Jun 04 Python
Python读取二进制文件代码方法解析
Jun 22 Python
python通用数据库操作工具 pydbclib的使用简介
Dec 21 Python
python 读取数据库并绘图的实例
Dec 03 #Python
JupyterNotebook设置Python环境的方法步骤
Dec 03 #Python
Anaconda 查看、创建、管理和使用python环境的方法
Dec 03 #Python
使用pygame写一个古诗词填空通关游戏
Dec 03 #Python
使用python和pygame制作挡板弹球游戏
Dec 03 #Python
win10系统Anaconda和Pycharm的Tensorflow2.0之CPU和GPU版本安装教程
Dec 03 #Python
在Python中使用MongoEngine操作数据库教程实例
Dec 03 #Python
You might like
PHP 学习路线与时间表
2010/02/21 PHP
php笔记之:文章中图片处理的使用
2013/04/26 PHP
php实现的替换敏感字符串类实例
2014/09/22 PHP
浅谈PHP中单引号和双引号到底有啥区别呢?
2015/03/04 PHP
ThinkPHP3.2.3数据库设置新特性
2015/03/05 PHP
PHP中常用的魔术方法
2017/04/28 PHP
深入理解PHP的远程多会话调试
2017/09/21 PHP
TP5框架请求响应参数实例分析
2019/10/17 PHP
基于mootools 1.3框架下的图片滑动效果代码
2011/04/22 Javascript
嵌入式iframe子页面与父页面js通信的方法
2015/01/20 Javascript
Javascript验证Visa和MasterCard信用卡号的方法
2015/07/27 Javascript
jQuery学习笔记之回调函数
2016/08/15 Javascript
jQuery ui autocomplete选择列表被Bootstrap模态窗遮挡的完美解决方法
2016/09/23 Javascript
jQuery 插件封装的方法
2016/11/16 Javascript
JavaScript Ajax实现异步通信
2016/12/14 Javascript
vue.js 使用v-if v-else发现没有执行解决办法
2017/05/15 Javascript
JavaScript函数中的this四种绑定形式
2017/08/15 Javascript
基于vue.js路由参数的实例讲解——简单易懂
2017/09/07 Javascript
php中and 和 &&出坑指南
2018/07/13 Javascript
vue-router判断页面未登录自动跳转到登录页的方法示例
2018/11/04 Javascript
详解如何使用node.js的开发框架express创建一个web应用
2018/12/20 Javascript
koa大型web项目中使用路由装饰器的方法示例
2019/04/02 Javascript
axios实现简单文件上传功能
2019/09/25 Javascript
AngularJs的$http发送POST请求,php无法接收Post的数据问题及解决方案
2020/08/13 Javascript
pyqt5自定义信号实例解析
2018/01/31 Python
python 执行文件时额外参数获取的实例
2018/12/18 Python
python 判断三个数字中的最大值实例代码
2019/07/24 Python
利用python实现平稳时间序列的建模方式
2020/06/03 Python
HTML5表单验证特性(知识点小结)
2020/03/10 HTML / CSS
党员公开承诺书和承诺事项
2014/03/25 职场文书
乡镇保密工作责任书
2014/07/28 职场文书
大学毕业生管理学求职信
2014/09/01 职场文书
民主评议党员自我鉴定
2014/10/21 职场文书
2014年学生会部门工作总结
2014/11/07 职场文书
2019年怎样写好导游词?
2019/07/02 职场文书
vue中 this.$set的使用详解
2021/11/17 Vue.js