首页
所有分类
TAGS
文字工具 EMOJI BIBLE
编程 Python

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 导入数据及作图的实现

- Author -

csluo19550418

声明：登载此文出于传递更多信息之目的，并不意味着赞同其观点或证实其描述。

Python 相关文章推荐
Python实现list反转实例汇总
Nov 11 Python
python集合用法实例分析
May 30 Python
13个最常用的Python深度学习库介绍
Oct 28 Python
如何通过python画loss曲线的方法
Jun 26 Python
python使用 __init__初始化操作简单示例
Sep 26 Python
python 求10个数的平均数实例
Dec 16 Python
Tensorflow 多线程设置方式
Feb 06 Python
python读取当前目录下的CSV文件数据
Mar 11 Python
python小技巧——将变量保存在本地及读取
Nov 13 Python
golang特有程序结构入门教程
Jun 02 Python
Python中的套接字编程是什么?
Jun 21 Python
python数字图像处理：图像的绘制
Jun 28 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
HAM-2000(1) 闭幕词(9) 慰问信(53) 通知(159) 导游词(322) 观后感(479) 实用文书(952) 申请书(21) 开幕词(63) 文书(291)
You might like
php中目录,文件操作详谈
2007/03/19 PHP
基于simple_html_dom的使用小结
2013/07/01 PHP
php检测mysql表是否存在的方法小结
2017/07/20 PHP
php 判断IP为有效IP地址的方法
2018/01/28 PHP
laravel7学习之无限级分类的最新实现方法
2020/09/30 PHP
区分JS中的undefined,null,"",0和false
2007/03/08 Javascript
jquery 注意事项与常用语法小结
2010/06/07 Javascript
Javascript 面向对象编程(coolshell)
2012/03/18 Javascript
Extjs4 GridPanel 的几种样式使用介绍
2013/04/18 Javascript
AngularJS控制器之间的数据共享及通信详解
2016/08/01 Javascript
30分钟快速实现小程序语音识别功能
2018/11/27 Javascript
JavaScript设计模式之命令模式实例分析
2019/01/16 Javascript
[39:07]LGD vs VP 2018国际邀请赛淘汰赛BO3 第二场 8.21
2018/08/22 DOTA
编写Python脚本来获取Google搜索结果的示例
2015/05/04 Python
python实现随机梯度下降（SGD）
2020/03/24 Python
简述Python2与Python3的不同点
2018/01/21 Python
为什么选择python编程语言入门黑客攻防 给你几个理由!
2018/02/02 Python
pandas将numpy数组写入到csv的实例
2018/07/04 Python
运用PyTorch动手搭建一个共享单车预测器
2019/08/06 Python
关于pandas的离散化,面元划分详解
2019/11/22 Python
详解pandas绘制矩阵散点图(scatter_matrix)的方法
2020/04/23 Python
Ubuntu20下的Django安装的方法步骤
2021/01/24 Python
python 多线程爬取壁纸网站的示例
2021/02/20 Python
一款纯css3实现的竖形二级导航的实例教程
2014/12/11 HTML / CSS
浅谈CSS3中display属性的Flex布局的方法
2017/08/14 HTML / CSS
Columbia美国官网：美国著名的户外服装品牌
2016/11/24 全球购物
地球上最先进的胡子和头发修剪器：Bevel
2018/01/23 全球购物
意大利比基尼品牌：MISS BIKINI
2019/11/02 全球购物
村官学习十八大感想
2014/01/15 职场文书
茶叶店创业计划书范文
2014/01/19 职场文书
病人写给医生的感谢信
2015/01/23 职场文书
幼儿园小班个人工作总结
2015/02/12 职场文书
六年级数学教学反思
2016/02/16 职场文书
2020年元旦晚会策划书模板
2019/12/30 职场文书
抖音短视频(douyin)去水印工具的实现代码
2021/03/30 Javascript
日元符号 ¥
2022/02/17 杂记