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按照多个字符对字符串进行分割的方法
Mar 17 Python
Python中的特殊语法:filter、map、reduce、lambda介绍
Apr 14 Python
Python使用BeautifulSoup库解析HTML基本使用教程
Mar 31 Python
python中执行shell的两种方法总结
Jan 10 Python
Django Admin 实现外键过滤的方法
Sep 29 Python
pandas series序列转化为星期几的实例
Apr 11 Python
在Python中pandas.DataFrame重置索引名称的实例
Nov 06 Python
浅析python 中大括号中括号小括号的区分
Jul 29 Python
python 实现return返回多个值
Nov 19 Python
python中文分词库jieba使用方法详解
Feb 11 Python
python 等差数列末项计算方式
May 03 Python
python爬虫智能翻页批量下载文件的实例详解
Feb 02 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
我的论坛源代码(九)
2006/10/09 PHP
php var_export与var_dump 输出的不同
2013/08/09 PHP
php实现与erlang的二进制通讯实例解析
2014/07/23 PHP
php抓取网站图片并保存的实现方法
2015/10/29 PHP
php检测mysql表是否存在的方法小结
2017/07/20 PHP
JavaScript中“+=”的应用
2007/02/02 Javascript
JavaScript中this关键字使用方法详解
2007/03/08 Javascript
JS 有名函数表达式全面解析
2010/03/19 Javascript
预加载css或javascript的js代码
2010/04/23 Javascript
JavaScript 设计模式 安全沙箱模式
2010/09/24 Javascript
提高javascript效率 一次判断,而不要次次判断
2012/03/30 Javascript
javascript实现鼠标放上后下边对应内容变换的效果
2015/08/06 Javascript
AngularJS实现表单手动验证和表单自动验证
2015/12/09 Javascript
深入理解JavaScript中的对象复制(Object Clone)
2016/05/18 Javascript
vue的基本用法与常见指令
2017/08/15 Javascript
Vue中之nextTick函数源码分析详解
2017/10/17 Javascript
浅析webpack-bundle-analyzer在vue-cli3中的使用
2019/10/23 Javascript
javascript实现鼠标点击生成文字特效
2019/12/24 Javascript
[46:55]完美世界DOTA2联赛决赛 FTD vs Phoenix 第三场 11.08
2020/11/11 DOTA
Python对两个有序列表进行合并和排序的例子
2014/06/13 Python
python脚本监控docker容器
2016/04/27 Python
python使用selenium登录QQ邮箱(附带滑动解锁)
2019/01/23 Python
python利用selenium进行浏览器爬虫
2019/04/25 Python
django如何通过类视图使用装饰器
2019/07/24 Python
Python3.6 + TensorFlow 安装配置图文教程(Windows 64 bit)
2020/02/24 Python
python 连续不等式语法糖实例
2020/04/15 Python
MxNet预训练模型到Pytorch模型的转换方式
2020/05/25 Python
汉语专业应届生求职信
2013/10/01 职场文书
学生党员的自我评价范文
2014/03/01 职场文书
党的群众路线教育实践活动宣传标语口号
2014/06/06 职场文书
乡镇四风对照检查材料
2014/08/31 职场文书
公司的门卫岗位职责
2014/09/09 职场文书
秋冬农业生产标语
2014/10/09 职场文书
人身损害赔偿协议书格式
2014/11/01 职场文书
canvas绘制折线路径动画实现
2021/05/12 Javascript
《艾尔登法环》发布最新「战技」宣传片
2022/04/03 其他游戏