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 13 Python
python制作花瓣网美女图片爬虫
Oct 28 Python
浅谈python中scipy.misc.logsumexp函数的运用场景
Jun 23 Python
Python极简代码实现杨辉三角示例代码
Nov 15 Python
Python数据结构与算法之二叉树结构定义与遍历方法详解
Dec 12 Python
Python实现模拟浏览器请求及会话保持操作示例
Jul 30 Python
python调用百度语音识别api
Aug 30 Python
Python中最大递归深度值的探讨
Mar 05 Python
Django如何将URL映射到视图
Jul 29 Python
你还在@微信官方?聊聊Python生成你想要的微信头像
Sep 25 Python
Python实现井字棋小游戏
Mar 09 Python
Django对接elasticsearch实现全文检索的示例代码
Aug 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
PHP下载生成的csv文件及问题总结
2015/08/06 PHP
php实现图片按比例截取的方法
2017/02/06 PHP
Js获取事件对象代码
2010/08/05 Javascript
可兼容IE的获取及设置cookie的jquery.cookie函数方法
2013/09/02 Javascript
JS设置获取cookies的方法
2014/01/26 Javascript
jQuery学习笔记之 Ajax操作篇(二) - 数据传递
2014/06/23 Javascript
javascript函数中参数传递问题示例探讨
2014/07/31 Javascript
5个数组Array方法: indexOf、filter、forEach、map、reduce使用实例
2015/01/29 Javascript
jQuery解决input超多的表单提交
2015/08/10 Javascript
JS鼠标拖拽实例分析
2015/11/23 Javascript
JavaScript中eval函数的问题
2016/01/31 Javascript
JavaScript函数节流和函数防抖之间的区别
2017/02/15 Javascript
用纯Node.JS弹出Windows系统消息提示框实例(MessageBox)
2017/05/17 Javascript
js+css实现红包雨效果
2018/07/12 Javascript
JS中DOM元素的attribute与property属性示例详解
2018/09/04 Javascript
微信小程序控制台提示warning:Now you can provide attr "wx:key" for a "wx:for" to improve performance解决方法
2019/02/21 Javascript
通过JavaScript下载文件到本地的方法(单文件)
2019/03/17 Javascript
通过jQuery学习js类型判断的技巧
2019/05/27 jQuery
vue 组件中使用 transition 和 transition-group实现过渡动画
2019/07/09 Javascript
layui表格数据重载
2019/07/27 Javascript
NodeJS有难度的面试题(能答对几个)
2019/10/09 NodeJs
微信小程序手动添加收货地址省市区联动
2020/05/18 Javascript
python网络编程学习笔记(五):socket的一些补充
2014/06/09 Python
Django中利用filter与simple_tag为前端自定义函数的实现方法
2017/06/15 Python
如何在sae中设置django,让sae的工作环境跟本地python环境一致
2017/11/21 Python
Golang GBK转UTF-8的例子
2019/08/26 Python
基于python+selenium自动健康打卡的实现代码
2021/01/13 Python
adidas官方旗舰店:德国运动用品制造商
2017/11/25 全球购物
德国前卫设计师时装在线商店:Luxury Loft
2019/11/04 全球购物
PHP面试题及答案一
2012/06/18 面试题
有关九一八事变的演讲稿
2014/09/14 职场文书
导游词欢迎词
2015/02/02 职场文书
产品调价通知函
2015/04/20 职场文书
入党转正介绍人意见
2015/06/03 职场文书
小组组名及励志口号
2015/12/24 职场文书
用python开发一款操作MySQL的小工具
2021/05/12 Python