首页
所有分类
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 smtplib模块发送SSL/TLS安全邮件实例
Apr 08 Python
Python实现随机生成手机号及正则验证手机号的方法
Apr 25 Python
mvc框架打造笔记之wsgi协议的优缺点以及接口实现
Aug 01 Python
python实现flappy bird游戏
Dec 24 Python
python实现年会抽奖程序
Jan 22 Python
python3 tkinter实现点击一个按钮跳出另一个窗口的方法
Jun 13 Python
python接口调用已训练好的caffe模型测试分类方法
Aug 26 Python
pytorch自定义二值化网络层方式
Jan 07 Python
Python读取表格类型文件代码实例
Feb 17 Python
django从后台返回html代码的实例
Mar 11 Python
python爬虫使用requests发送post请求示例详解
Aug 05 Python
python实现自动打卡的示例代码
Oct 10 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
荷兰购物网站(99) 葡萄牙购物网站(13) 奥地利购物网站(28) 南非购物网站(8) 欧洲购物网站(15) 土耳其购物网站(10) 波兰购物网站(52) 阿联酋购物网站(26) 哥伦比亚购物网站(5) 泰国购物网站(30)
You might like
thinkphp3.2嵌入百度编辑器ueditor的实例代码
2017/07/13 PHP
JavaScript this调用规则说明
2010/03/08 Javascript
设为首页加入收藏兼容360/火狐/谷歌/IE等主流浏览器的代码
2013/03/26 Javascript
利用NodeJS和PhantomJS抓取网站页面信息以及网站截图
2013/11/18 NodeJs
JavaScript采用递归算法计算阶乘实例
2015/08/04 Javascript
JQuery核心函数是什么及使用方法介绍
2016/05/03 Javascript
浅析js绑定事件的常用方法
2016/05/15 Javascript
Vue.js每天必学之指令系统与自定义指令
2016/09/07 Javascript
require.js 加载 vue组件 r.js 合并压缩的实例
2016/10/14 Javascript
实例讲解vue源码架构
2019/01/24 Javascript
JS集合set类的实现与使用方法示例
2019/02/01 Javascript
微信小程序定义和调用全局变量globalData的实现
2019/11/01 Javascript
Vuejs中的watch实例详解(监听者)
2020/01/05 Javascript
JavaScript设计模式--简单工厂模式定义与应用案例详解
2020/05/23 Javascript
element-ui点击查看大图的方法示例
2020/12/14 Javascript
Python 字符串操作实现代码(截取/替换/查找/分割)
2013/06/08 Python
Python3 socket同步通信简单示例
2017/06/07 Python
python利用dir函数查看类中所有成员函数示例代码
2017/09/08 Python
python3+PyQt5实现自定义窗口部件Counters
2018/04/20 Python
python实现本地图片转存并重命名的示例代码
2018/10/27 Python
python+opencv边缘提取与各函数参数解析
2020/03/09 Python
浅谈Python中os模块及shutil模块的常规操作
2020/04/03 Python
Jupyter notebook 启动闪退问题的解决
2020/04/13 Python
详解canvas.toDataURL()报错的解决方案全都在这了
2020/03/31 HTML / CSS
非洲NO.1网上商店：Jumia肯尼亚
2016/08/18 全球购物
狗狗玩具、零食和咀嚼物的月度送货服务：Super Chewer
2018/08/22 全球购物
北美最大的参茸药食商城：德成行
2020/12/06 全球购物
局部内部类是否可以访问非final变量？
2013/04/20 面试题
自我鉴定范文
2013/11/10 职场文书
上学迟到的检讨书
2014/01/11 职场文书
国际贸易专业个人鉴定
2014/02/22 职场文书
阳光体育活动实施方案
2014/05/25 职场文书
4s店销售经理岗位职责
2014/07/19 职场文书
2014年政协工作总结
2014/12/09 职场文书
拔河比赛队名及霸气口号
2015/12/24 职场文书
VMware虚拟机安装 Windows Server 2022的详细图文教程
2022/09/23 Servers