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

方法一、

filename='/home/res/user/csluo/test.txt' #将文件名赋值为变量

X,Y,Z=[ ],[ ],[ ] #给三个空列表

with open(filename, 'r') as f: #打开文档

lines=f.readlines() #按行读取数据

for i in lines:

value=[ float(s) for s in line.split( ) ] #将 行数据 以空格分开，并以浮点型赋给value

X=append(value[0]) #将value 的第一个数值 添加到 列表X

Y=append(value[1])

Z=append(value[2])

方法二、

a = numpy.loadtxt('a.txt')

a[:,0]

方法三、

>>> import pandas as pd

>>> data = pd.read_csv("aaa.csv")

附录：IRAS16293 的数据 绘图代码

import matplotlib

import matplotlib.pyplot as plt

import numpy as np

##################################################################

######################################################################################### ######## 1

##################################################################

#data_a

filename = '/home/luoluo/Documents/python_scripts/Adata_cut/A1_145122-145580.txt'

X10,Y10= [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()] #以空格分开 行 并输出 浮点数

X10.append(value[0])

Y10.append(value[1])

#c2h5oh

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A1/xclass_spectrum_output.dat'

X11,Y11 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X11.append(value[0])

Y11.append(value[1])

#ch3coch3

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A1/xclass_spectrum_output.dat'

X12,Y12 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X12.append(value[0])

Y12.append(value[1])

#ch3ocho_v18=1

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A1/xclass_spectrum_output.dat'

X13,Y13 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X13.append(value[0])

Y13.append(value[1])

#ch3ocho_v=0

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A1/xclass_spectrum_output.dat'

X14,Y14 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X14.append(value[0])

Y14.append(value[1])

################################################### data_b

filename = '/home/luoluo/Documents/python_scripts/Bdata_cut/B1_145122-145580.txt'

X15,Y15= [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X15.append(value[0])

Y15.append(value[1])

#c2h5oh

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B1/xclass_spectrum_output.dat'

X16,Y16 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X16.append(value[0])

Y16.append(value[1])

#ch3coch3

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B1/xclass_spectrum_output.dat'

X17,Y17 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X17.append(value[0])

Y17.append(value[1])

#ch3ocho_v18=1

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B1/xclass_spectrum_output.dat'

X18,Y18 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X18.append(value[0])

Y18.append(value[1])

#ch3ocho_v=0

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B1/xclass_spectrum_output.dat'

X19,Y19 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X19.append(value[0])

Y19.append(value[1])

################################################################### #2 plot

plt.figure(1)

plt.subplot(211)

plt.axis([145122,145580, -1, 15])

new_ticks = np.linspace(145122,145580,10)

plt.xticks(new_ticks)

matplotlib.rcParams['xtick.direction'] = 'in'

matplotlib.rcParams['ytick.direction'] = 'in'

plt.plot(X10,Y10,color='k',lw=1.0,ls='-',label='DATA')

plt.plot(X11,Y11,color='r',lw=1.0,ls='-',label='C2H5OH')

plt.plot(X12,Y12,color='g',lw=1.0,ls='-',label='CH3COCH3')

plt.plot(X13,Y13,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')

plt.plot(X14,Y14,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')

plt.legend(loc='upper left')

#plt.xlabel('Freq(MHZ)')

plt.ylabel('T(K)')

plt.text(145535,13,'IRAS16293A')

#plt.show()

#plt.savefig('8.png')

plt.subplot(212)

plt.axis([145122,145580, -1, 15])

new_ticks = np.linspace(145122,145580,10)

plt.xticks(new_ticks)

plt.plot(X15,Y15,color='k',lw=1.0,ls='-',label='DATA')

plt.plot(X16,Y16,color='r',lw=1.0,ls='-',label='C2H5OH')

plt.plot(X17,Y17,color='g',lw=1.0,ls='-',label='CH3COCH3')

plt.plot(X18,Y18,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')

plt.plot(X19,Y19,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')

plt.legend(loc='upper left')

plt.xlabel('Freq(MHZ)')

plt.ylabel('T(K)')

plt.text(145535,13,'IRAS16293B')

plt.plot([145263,145263],[1,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("C2H5OH", xy = (145263,1), xytext = (145280,9),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.show()

##################################################################

######################################################################################### ######## 2

##################################################################

#data_a

filename = '/home/luoluo/Documents/python_scripts/Adata_cut/A2_145580-146050.txt'

X20,Y20= [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X20.append(value[0])

Y20.append(value[1])

#c2h5oh

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A2/xclass_spectrum_output.dat'

X21,Y21 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X21.append(value[0])

Y21.append(value[1])

#ch3coch3

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A2/xclass_spectrum_output.dat'

X22,Y22 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X22.append(value[0])

Y22.append(value[1])

#ch3ocho_v18=1

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A2/xclass_spectrum_output.dat'

X23,Y23 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X23.append(value[0])

Y23.append(value[1])

#ch3ocho_v=0

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A2/xclass_spectrum_output.dat'

X24,Y24 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X24.append(value[0])

Y24.append(value[1])

################################################### data_b

filename = '/home/luoluo/Documents/python_scripts/Bdata_cut/B2.txt'

X25,Y25= [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X25.append(value[0])

Y25.append(value[1])

#c2h5oh

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B2/xclass_spectrum_output.dat'

X26,Y26 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X26.append(value[0])

Y26.append(value[1])

#ch3coch3

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B2/xclass_spectrum_output.dat'

X27,Y27 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X27.append(value[0])

Y27.append(value[1])

#ch3ocho_v18=1

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B2/xclass_spectrum_output.dat'

X28,Y28 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X28.append(value[0])

Y28.append(value[1])

#ch3ocho_v=0

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B2/xclass_spectrum_output.dat'

X29,Y29 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X29.append(value[0])

Y29.append(value[1])

################################################################### #2 plot

plt.figure(2)

plt.subplot(211)

plt.axis([145580,146050, -1, 15])

new_ticks = np.linspace(145580,146050,10)

plt.xticks(new_ticks)

matplotlib.rcParams['xtick.direction'] = 'in'

matplotlib.rcParams['ytick.direction'] = 'in'

plt.plot(X20,Y20,color='k',lw=1.0,ls='-',label='DATA')

plt.plot(X21,Y21,color='r',lw=1.0,ls='-',label='C2H5OH')

plt.plot(X22,Y22,color='g',lw=1.0,ls='-',label='CH3COCH3')

plt.plot(X23,Y23,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')

plt.plot(X24,Y24,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')

plt.legend(loc='upper left')

#plt.xlabel('Freq(MHZ)')

plt.ylabel('T(K)')

plt.text(146000,13,'IRAS16293A')

plt.plot([145618.5,145618.5],[2.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (145618,7), xytext = (145650,11),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([146020,146020],[7,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (146020,6), xytext = (146005,10),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.show()

#plt.savefig('8.png')

plt.subplot(212)

plt.axis([145580,146050, -1, 15])

new_ticks = np.linspace(145580,146050,10)

plt.xticks(new_ticks)

plt.plot(X25,Y25,color='k',lw=1.0,ls='-',label='DATA')

plt.plot(X26,Y26,color='r',lw=1.0,ls='-',label='C2H5OH')

plt.plot(X27,Y27,color='g',lw=1.0,ls='-',label='CH3COCH3')

plt.plot(X28,Y28,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')

plt.plot(X29,Y29,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')

plt.legend(loc='upper left')

plt.xlabel('Freq(MHZ)')

plt.ylabel('T(K)')

plt.text(146000,13,'IRAS16293B')

plt.plot([145618.5,145618.5],[2.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (145618,7), xytext = (145650,11),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([146020,146020],[7,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (146020,6), xytext = (145990,10),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.show()

##################################################################

######################################################################################### ######## 3

##################################################################

#data_a

filename = '/home/luoluo/Documents/python_scripts/Adata_cut/A3.txt'

X30,Y30= [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X30.append(value[0])

Y30.append(value[1])

#c2h5oh

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A3/xclass_spectrum_output.dat'

X31,Y31 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X31.append(value[0])

Y31.append(value[1])

#ch3coch3

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A3/xclass_spectrum_output.dat'

X32,Y32 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X32.append(value[0])

Y32.append(value[1])

#ch3ocho_v18=1

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A3/xclass_spectrum_output.dat'

X33,Y33 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X33.append(value[0])

Y33.append(value[1])

#ch3ocho_v=0

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A3/xclass_spectrum_output.dat'

X34,Y34 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X34.append(value[0])

Y34.append(value[1])

################################################### data_b

filename = '/home/luoluo/Documents/python_scripts/Bdata_cut/B3_146050-146507.txt'

X35,Y35= [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X35.append(value[0])

Y35.append(value[1])

#c2h5oh

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B3/xclass_spectrum_output.dat'

X36,Y36 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X36.append(value[0])

Y36.append(value[1])

#ch3coch3

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B3/xclass_spectrum_output.dat'

X37,Y37 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X37.append(value[0])

Y37.append(value[1])

#ch3ocho_v18=1

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B3/xclass_spectrum_output.dat'

X38,Y38 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X38.append(value[0])

Y38.append(value[1])

#ch3ocho_v=0

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B3/xclass_spectrum_output.dat'

X39,Y39 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X39.append(value[0])

Y39.append(value[1])

################################################################### #2 plot

plt.figure(3)

plt.subplot(211)

plt.axis([146050,146507, -1, 15])

new_ticks = np.linspace(146050,146507,10)

plt.xticks(new_ticks)

matplotlib.rcParams['xtick.direction'] = 'in'

matplotlib.rcParams['ytick.direction'] = 'in'

plt.plot(X30,Y30,color='k',lw=1.0,ls='-',label='DATA')

plt.plot(X31,Y31,color='r',lw=1.0,ls='-',label='C2H5OH')

plt.plot(X32,Y32,color='g',lw=1.0,ls='-',label='CH3COCH3')

plt.plot(X33,Y33,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')

plt.plot(X34,Y34,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')

plt.legend(loc='upper left')

#plt.xlabel('Freq(MHZ)')

plt.ylabel('T(K)')

plt.text(146465,13,'IRAS16293A')

plt.plot([146172,146172],[2.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (146172,7), xytext = (146190,11),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([146232,146232],[3,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (146232,10), xytext = (146210,5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([146343,146343],[2.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (146343,7), xytext = (146320,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([146355,146355],[2.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (146355,6), xytext = (146380,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.show()

#plt.savefig('8.png')

plt.subplot(212)

plt.axis([146050,146507, -1, 15])

new_ticks = np.linspace(146050,146507,10)

plt.xticks(new_ticks)

plt.plot(X35,Y35,color='k',lw=1.0,ls='-',label='DATA')

plt.plot(X36,Y36,color='r',lw=1.0,ls='-',label='C2H5OH')

plt.plot(X37,Y37,color='g',lw=1.0,ls='-',label='CH3COCH3')

plt.plot(X38,Y38,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')

plt.plot(X39,Y39,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')

plt.legend(loc='upper left')

plt.xlabel('Freq(MHZ)')

plt.ylabel('T(K)')

plt.text(146465,13,'IRAS16293B')

plt.plot([146343,146343],[5.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (146343,7), xytext = (146320,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([146355,146355],[5.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (146355,6), xytext = (146380,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.plot([157324,157324],[2.5,15], color ='grey', linewidth=1, linestyle="--")

#plt.annotate("CH3COCH3", xy = (157324,3), xytext = (157347,4.5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.show()

##################################################################

######################################################################################### ######## 4

##################################################################

#data_a

filename = '/home/luoluo/Documents/python_scripts/Adata_cut/A4'

X40,Y40= [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X40.append(value[0])

Y40.append(value[1])

#c2h5oh

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A4/xclass_spectrum_output.dat'

X41,Y41 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X41.append(value[0])

Y41.append(value[1])

#ch3coch3

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A4/xclass_spectrum_output.dat'

X42,Y42 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X42.append(value[0])

Y42.append(value[1])

#ch3ocho_v18=1

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A4/xclass_spectrum_output.dat'

X43,Y43 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X43.append(value[0])

Y43.append(value[1])

#ch3ocho_v=0

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A4/xclass_spectrum_output.dat'

X44,Y44 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X44.append(value[0])

Y44.append(value[1])

################################################### data_b

filename = '/home/luoluo/Documents/python_scripts/Bdata_cut/B4_146507-146976.txt'

X45,Y45= [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X45.append(value[0])

Y45.append(value[1])

#c2h5oh

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B4/xclass_spectrum_output.dat'

X46,Y46 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X46.append(value[0])

Y46.append(value[1])

#ch3coch3

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B4/xclass_spectrum_output.dat'

X47,Y47 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X47.append(value[0])

Y47.append(value[1])

#ch3ocho_v18=1

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B4/xclass_spectrum_output.dat'

X48,Y48 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X48.append(value[0])

Y48.append(value[1])

#ch3ocho_v=0

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B4/xclass_spectrum_output.dat'

X49,Y49 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X49.append(value[0])

Y49.append(value[1])

################################################################### #2 plot

plt.figure(4)

plt.subplot(211)

plt.axis([146507,146976, -1, 15])

new_ticks = np.linspace(146507,146976,10)

plt.xticks(new_ticks)

matplotlib.rcParams['xtick.direction'] = 'in'

matplotlib.rcParams['ytick.direction'] = 'in'

plt.plot(X40,Y40,color='k',lw=1.0,ls='-',label='DATA')

plt.plot(X41,Y41,color='r',lw=1.0,ls='-',label='C2H5OH')

plt.plot(X42,Y42,color='g',lw=1.0,ls='-',label='CH3COCH3')

plt.plot(X43,Y43,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')

plt.plot(X44,Y44,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')

plt.legend(loc='upper left')

#plt.xlabel('Freq(MHZ)')

plt.ylabel('T(K)')

plt.text(146930,13,'IRAS16293A')

#plt.plot([157323,157323],[2.5,15], color ='grey', linewidth=1, linestyle="--")

#plt.annotate("CH3COCH3", xy = (157323,3), xytext = (157345,4.5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.show()

#plt.savefig('8.png')

plt.subplot(212)

plt.axis([146507,146976, -1, 15])

new_ticks = np.linspace(146507,146976,10)

plt.xticks(new_ticks)

plt.plot(X45,Y45,color='k',lw=1.0,ls='-',label='DATA')

plt.plot(X46,Y46,color='r',lw=1.0,ls='-',label='C2H5OH')

plt.plot(X47,Y47,color='g',lw=1.0,ls='-',label='CH3COCH3')

plt.plot(X48,Y48,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')

plt.plot(X49,Y49,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')

plt.legend(loc='upper left')

plt.xlabel('Freq(MHZ)')

plt.ylabel('T(K)')

plt.text(146930,13,'IRAS16293B')

#plt.plot([157324,157324],[2.5,15], color ='grey', linewidth=1, linestyle="--")

#plt.annotate("CH3COCH3", xy = (157324,3), xytext = (157347,4.5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.show()

##################################################################

######################################################################################### ######## 5

##################################################################

#data_a

filename = '/home/luoluo/Documents/python_scripts/Adata_cut/A5_156840-157303.txt'

X50,Y50= [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X50.append(value[0])

Y50.append(value[1])

#c2h5oh

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A5/xclass_spectrum_output.dat'

X51,Y51 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X51.append(value[0])

Y51.append(value[1])

#ch3coch3

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A5/xclass_spectrum_output.dat'

X52,Y52 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X52.append(value[0])

Y52.append(value[1])

#ch3ocho_v18=1

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A5/xclass_spectrum_output.dat'

X53,Y53 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X53.append(value[0])

Y53.append(value[1])

#ch3ocho_v=0

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A5/xclass_spectrum_output.dat'

X54,Y54 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X54.append(value[0])

Y54.append(value[1])

################################################### data_b

filename = '/home/luoluo/Documents/python_scripts/Bdata_cut/B5_156840-157303.txt'

X55,Y55= [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X55.append(value[0])

Y55.append(value[1])

#c2h5oh

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B5/xclass_spectrum_output.dat'

X56,Y56 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X56.append(value[0])

Y56.append(value[1])

#ch3coch3

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B5/xclass_spectrum_output.dat'

X57,Y57 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X57.append(value[0])

Y57.append(value[1])

#ch3ocho_v18=1

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B5/xclass_spectrum_output.dat'

X58,Y58 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X58.append(value[0])

Y58.append(value[1])

#ch3ocho_v=0

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B5/xclass_spectrum_output.dat'

X59,Y59 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X59.append(value[0])

Y59.append(value[1])

################################################################### #2 plot

plt.figure(5)

plt.subplot(211)

plt.axis([156840, 157303, -1, 15])

new_ticks = np.linspace(156840,157303,10)

plt.xticks(new_ticks)

matplotlib.rcParams['xtick.direction'] = 'in'

matplotlib.rcParams['ytick.direction'] = 'in'

plt.plot(X50,Y50,color='k',lw=1.0,ls='-',label='DATA')

plt.plot(X51,Y51,color='r',lw=1.0,ls='-',label='C2H5OH')

plt.plot(X52,Y52,color='g',lw=1.0,ls='-',label='CH3COCH3')

plt.plot(X53,Y53,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')

plt.plot(X54,Y54,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')

plt.legend(loc='upper left')

#plt.xlabel('Freq(MHZ)')

plt.ylabel('T(K)')

plt.text(157256,13,'IRAS16293A')

#plt.show()

#plt.savefig('8.png')

plt.subplot(212)

plt.axis([156820,157303, -1, 15])

new_ticks = np.linspace(156820,157303,10)

plt.xticks(new_ticks)

plt.plot(X55,Y55,color='k',lw=1.0,ls='-',label='DATA')

plt.plot(X56,Y56,color='r',lw=1.0,ls='-',label='C2H5OH')

plt.plot(X57,Y57,color='g',lw=1.0,ls='-',label='CH3COCH3')

plt.plot(X58,Y58,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')

plt.plot(X59,Y59,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')

plt.legend(loc='upper left')

plt.xlabel('Freq(MHZ)')

plt.ylabel('T(K)')

plt.text(157255,13,'IRAS16293B')

#plt.plot([157324,157324],[2.5,15], color ='grey', linewidth=1, linestyle="--")

#plt.annotate("CH3COCH3", xy = (157324,3), xytext = (157347,4.5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.show()

##################################################################

######################################################################################### ######## 6

##################################################################

#data_a

filename = '/home/luoluo/Documents/python_scripts/Adata_cut/A6_157303-157820.txt'

X60,Y60= [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X60.append(value[0])

Y60.append(value[1])

#c2h5oh

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A6/xclass_spectrum_output.dat'

X61,Y61 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X61.append(value[0])

Y61.append(value[1])

#ch3coch3

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A6/xclass_spectrum_output.dat'

X62,Y62 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X62.append(value[0])

Y62.append(value[1])

#ch3ocho_v18=1

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A6/xclass_spectrum_output.dat'

X63,Y63 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X63.append(value[0])

Y63.append(value[1])

#ch3ocho_v=0

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A6/xclass_spectrum_output.dat'

X64,Y64 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X64.append(value[0])

Y64.append(value[1])

################################################### data_b

filename = '/home/luoluo/Documents/python_scripts/Bdata_cut/B6_157303-157820.txt'

X65,Y65= [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X65.append(value[0])

Y65.append(value[1])

#c2h5oh

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B6/xclass_spectrum_output.dat'

X66,Y66 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X66.append(value[0])

Y66.append(value[1])

#ch3coch3

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B6/xclass_spectrum_output.dat'

X67,Y67 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X67.append(value[0])

Y67.append(value[1])

#ch3ocho_v18=1

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B6/xclass_spectrum_output.dat'

X68,Y68 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X68.append(value[0])

Y68.append(value[1])

#ch3ocho_v=0

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B6/xclass_spectrum_output.dat'

X69,Y69 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X69.append(value[0])

Y69.append(value[1])

################################################################### #2 plot

plt.figure(6)

plt.subplot(211)

plt.axis([157303, 157820, -1, 15])

new_ticks = np.linspace(157303, 157820,10)

plt.xticks(new_ticks)

matplotlib.rcParams['xtick.direction'] = 'in'

matplotlib.rcParams['ytick.direction'] = 'in'

plt.plot(X60,Y60,color='k',lw=1.0,ls='-',label='DATA')

plt.plot(X61,Y61,color='r',lw=1.0,ls='-',label='C2H5OH')

plt.plot(X62,Y62,color='g',lw=1.0,ls='-',label='CH3COCH3')

plt.plot(X63,Y63,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')

plt.plot(X64,Y64,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')

plt.legend(loc='upper left')

#plt.xlabel('Freq(MHZ)')

plt.ylabel('T(K)')

plt.text(157765,13,'IRAS16293A')

plt.plot([157323,157323],[2.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3COCH3", xy = (157323,3), xytext = (157345,4.5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.show()

#plt.savefig('8.png')

plt.subplot(212)

plt.axis([157303, 157820, -1, 15])

new_ticks = np.linspace(157303, 157820,10)

plt.xticks(new_ticks)

plt.plot(X65,Y65,color='k',lw=1.0,ls='-',label='DATA')

plt.plot(X66,Y66,color='r',lw=1.0,ls='-',label='C2H5OH')

plt.plot(X67,Y67,color='g',lw=1.0,ls='-',label='CH3COCH3')

plt.plot(X68,Y68,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')

plt.plot(X69,Y69,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')

plt.legend(loc='upper left')

plt.xlabel('Freq(MHZ)')

plt.ylabel('T(K)')

plt.text(157765,13,'IRAS16293B')

plt.plot([157324,157324],[2.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3COCH3", xy = (157324,3), xytext = (157347,4.5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([157505.5,157505.5],[5.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("C2H5OH", xy = (157505,6), xytext = (157450,11),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.show()

###############################################################################################

################################################################################################################ 7

###############################################################################################

#data_a

filename = '/home/luoluo/Documents/python_scripts/Adata_cut/A7_157820-158228.txt'

X70,Y70= [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X70.append(value[0])

Y70.append(value[1])

#c2h5oh

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A7/xclass_spectrum_output.dat'

X71,Y71 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X71.append(value[0])

Y71.append(value[1])

#ch3coch3

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A7/xclass_spectrum_output.dat'

X72,Y72 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X72.append(value[0])

Y72.append(value[1])

#ch3ocho_v18=1

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A7/xclass_spectrum_output.dat'

X73,Y73 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X73.append(value[0])

Y73.append(value[1])

#ch3ocho_v=0

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A7/xclass_spectrum_output.dat'

X74,Y74 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X74.append(value[0])

Y74.append(value[1])

#########################################################

#data_b

filename = '/home/luoluo/Documents/python_scripts/Bdata_cut/B7_157820-158228.txt'

X75,Y75= [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X75.append(value[0])

Y75.append(value[1])

#c2h5oh

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B7/xclass_spectrum_output.dat'

X76,Y76 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X76.append(value[0])

Y76.append(value[1])

#ch3coch3

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B7/xclass_spectrum_output.dat'

X77,Y77 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X77.append(value[0])

Y77.append(value[1])

#ch3ocho_v18=1

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B7/xclass_spectrum_output.dat'

X78,Y78 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X78.append(value[0])

Y78.append(value[1])

#ch3ocho_v=0

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B7/xclass_spectrum_output.dat'

X79,Y79 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X79.append(value[0])

Y79.append(value[1])

#############################################################################################################################################

#2 plot

plt.figure(7)

plt.subplot(211)

plt.axis([157820, 158228, -1, 15])

new_ticks = np.linspace(157820,158228,10)

plt.xticks(new_ticks)

matplotlib.rcParams['xtick.direction'] = 'in'

matplotlib.rcParams['ytick.direction'] = 'in'

plt.plot(X70,Y70,color='k',lw=1.0,ls='-',label='DATA')

plt.plot(X71,Y71,color='r',lw=1.0,ls='-',label='C2H5OH')

plt.plot(X72,Y72,color='g',lw=1.0,ls='-',label='CH3COCH3')

plt.plot(X73,Y73,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')

plt.plot(X74,Y74,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')

plt.legend(loc='upper left')

#plt.xlabel('Freq(MHZ)')

plt.ylabel('T(K)')

plt.text(158183,13,'IRAS16293A')

plt.plot([157894,157894],[5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (157895,7), xytext = (157880,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158120,158120],[5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3COCH3", xy = (158120,5), xytext = (158110,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158172,158172],[7,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3COCH3", xy = (158171,8), xytext = (158150,10),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158225,158225],[4,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3COCH3", xy = (158225,5), xytext = (158200,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.show()

#plt.savefig('8.png')

plt.subplot(212)

plt.axis([157820, 158228, -1, 15])

new_ticks = np.linspace(157820, 158228,10)

plt.xticks(new_ticks)

plt.plot(X75,Y75,color='k',lw=1.0,ls='-',label='DATA')

plt.plot(X76,Y76,color='r',lw=1.0,ls='-',label='C2H5OH')

plt.plot(X77,Y77,color='g',lw=1.0,ls='-',label='CH3COCH3')

plt.plot(X78,Y78,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')

plt.plot(X79,Y79,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')

plt.legend(loc='upper left')

plt.xlabel('Freq(MHZ)')

plt.ylabel('T(K)')

plt.text(158183,13,'IRAS16293B')

plt.plot([157895,157895],[5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (157896,7), xytext = (157880,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158121,158121],[5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3COCH3", xy = (158121,7), xytext = (158110,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158173,158173],[7,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3COCH3", xy = (158172,8), xytext = (158150,5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158226,158226],[4,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3COCH3", xy = (158225,5), xytext = (158205,10),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.plot([158389,158389],[3.5,15], color ='grey', linewidth=1, linestyle="--")

#plt.annotate("CH3OCHO v18=1", xy = (158390,5), xytext = (158400,13),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.show()

######################################################################################

###################################################################################### ###### 8

######################################################################################

#

#1 import data

#data_a

filename = 'Adata_cut/A8_158228-158680.txt'

X80,Y80= [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X80.append(value[0])

Y80.append(value[1])

#c2h5oh

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A8/xclass_spectrum_output.dat'

X81,Y81 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X81.append(value[0])

Y81.append(value[1])

#ch3coch3

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A8/xclass_spectrum_output.dat'

X82,Y82 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X82.append(value[0])

Y82.append(value[1])

#ch3ocho_v18=1

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A8/xclass_spectrum_output.dat'

X83,Y83 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X83.append(value[0])

Y83.append(value[1])

#ch3ocho_v=0

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A8/xclass_spectrum_output.dat'

X84,Y84 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X84.append(value[0])

Y84.append(value[1])

#########################################################

#data_b

filename = 'Bdata_cut/B8_158228-158680.txt'

X85,Y85= [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X85.append(value[0])

Y85.append(value[1])

#c2h5oh

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B8/xclass_spectrum_output.dat'

X86,Y86 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X86.append(value[0])

Y86.append(value[1])

#ch3coch3

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B8/xclass_spectrum_output.dat'

X87,Y87 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X87.append(value[0])

Y87.append(value[1])

#ch3ocho_v18=1

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B8/xclass_spectrum_output.dat'

X88,Y88 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X88.append(value[0])

Y88.append(value[1])

#ch3ocho_v=0

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B8/xclass_spectrum_output.dat'

X89,Y89 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X89.append(value[0])

Y89.append(value[1])

#############################################################

#2 plot

plt.figure(8)

plt.subplot(211)

plt.axis([158228, 158680, -1, 15])

new_ticks = np.linspace(158228,158680,10)

plt.xticks(new_ticks)

matplotlib.rcParams['xtick.direction'] = 'in'

matplotlib.rcParams['ytick.direction'] = 'in'

plt.plot(X80,Y80,color='k',lw=1.0,ls='-',label='DATA')

plt.plot(X81,Y81,color='r',lw=1.0,ls='-',label='C2H5OH')

plt.plot(X82,Y82,color='g',lw=1.0,ls='-',label='CH3COCH3')

plt.plot(X83,Y83,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')

plt.plot(X84,Y84,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')

plt.legend(loc='upper left')

#plt.xlabel('Freq(MHZ)')

plt.ylabel('T(K)')

plt.text(158635,13,'IRAS16293A')

plt.plot([158378,158378],[1.9,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (158378, 2.0), xytext = (158300,5),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158388,158388],[3.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (158388, 4.5), xytext = (158370,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158430,158430],[3.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (158430, 3.5), xytext = (158400,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158484,158484],[3.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (158485, 4), xytext = (158450,9),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158519,158519],[3.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (158518,4), xytext = (158490,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158540,158540],[3.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (158541,4), xytext = (158550,11),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158675,158675],[3.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (158674,3), xytext = (158610,8),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.show()

#plt.savefig('8.png')

plt.subplot(212)

plt.axis([158228, 158680, -1, 15])

new_ticks = np.linspace(158228,158680,10)

plt.xticks(new_ticks)

plt.plot(X85,Y85,color='k',lw=1.0,ls='-',label='DATA')

plt.plot(X86,Y86,color='r',lw=1.0,ls='-',label='C2H5OH')

plt.plot(X87,Y87,color='g',lw=1.0,ls='-',label='CH3COCH3')

plt.plot(X88,Y88,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')

plt.plot(X89,Y89,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')

plt.legend(loc='upper left')

plt.xlabel('Freq(MHZ)')

plt.ylabel('T(K)')

plt.text(158635,13,'IRAS16293B')

plt.plot([158389,158389],[3.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (158390,5), xytext = (158400,13),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158484,158484],[3.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (158484,5), xytext = (158430,4),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158520,158520],[3.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (158520,4), xytext = (158500,6),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158540,158540],[3.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (158540,4), xytext = (158570,11),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158675,158675],[3.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (158675,4), xytext = (158640,8),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.show()

####################################################################################################

#################################################################################################### ######## 9

####################################################################################################

#data_a

filename = '/home/luoluo/Documents/python_scripts/Adata_cut/A9_158680-159155.txt'

X90,Y90= [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X90.append(value[0])

Y90.append(value[1])

#c2h5oh

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A9/xclass_spectrum_output.dat'

X91,Y91 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X91.append(value[0])

Y91.append(value[1])

#ch3coch3

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A9/xclass_spectrum_output.dat'

X92,Y92 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X92.append(value[0])

Y92.append(value[1])

#ch3ocho_v18=1

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A9/xclass_spectrum_output.dat'

X93,Y93 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X93.append(value[0])

Y93.append(value[1])

#ch3ocho_v=0

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A9/xclass_spectrum_output.dat'

X94,Y94 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X94.append(value[0])

Y94.append(value[1])

#########################################################

#data_b

filename = '/home/luoluo/Documents/python_scripts/Bdata_cut/B9_158680-159155.txt'

X95,Y95= [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X95.append(value[0])

Y95.append(value[1])

#c2h5oh

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B9/xclass_spectrum_output.dat'

X96,Y96 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X96.append(value[0])

Y96.append(value[1])

#ch3coch3

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B9/xclass_spectrum_output.dat'

X97,Y97 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X97.append(value[0])

Y97.append(value[1])

#ch3ocho_v18=1

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B9/xclass_spectrum_output.dat'

X98,Y98 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X98.append(value[0])

Y98.append(value[1])

#ch3ocho_v=0

filename = '/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B9/xclass_spectrum_output.dat'

X99,Y99 = [],[]

with open(filename, 'r') as f:

lines = f.readlines()

for line in lines:

value = [float(s) for s in line.split()]

X99.append(value[0])

Y99.append(value[1])

###2 plot

plt.figure(9)

#A

plt.subplot(211)

plt.axis([158680, 159155, -1, 15])

new_ticks = np.linspace(158680,159155,10)

plt.xticks(new_ticks)

matplotlib.rcParams['xtick.direction'] = 'in'

matplotlib.rcParams['ytick.direction'] = 'in'

plt.plot(X90,Y90,color='k',lw=1.0,ls='-',label='DATA')

plt.plot(X91,Y91,color='r',lw=1.0,ls='-',label='C2H5OH')

plt.plot(X92,Y92,color='g',lw=1.0,ls='-',label='CH3COCH3')

plt.plot(X93,Y93,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')

plt.plot(X94,Y94,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')

plt.legend(loc='upper left')

#plt.xlabel('Freq(MHZ)')

plt.ylabel('T(K)')

plt.text(159110,13,'IRAS16293A')

plt.plot([158691,158691],[8,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (158691, 5), xytext = (158730,3),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158701.9,158701.9],[8,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (158701, 10), xytext = (158770,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([159003.5,159003.5],[3.4,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (159005, 3), xytext = (159020,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([159106,159106],[3.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (159107, 3.5), xytext = (159110,7),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#B

plt.subplot(212)

plt.axis([158680, 159155, -1, 15])

new_ticks = np.linspace(158680,159122,10)

plt.xticks(new_ticks)

plt.plot(X95,Y95,color='k',lw=1.0,ls='-',label='DATA')

plt.plot(X96,Y96,color='r',lw=1.0,ls='-',label='C2H5OH')

plt.plot(X97,Y97,color='g',lw=1.0,ls='-',label='CH3COCH3')

plt.plot(X98,Y98,color='cyan',lw=1.0,ls='-',label='CH3OCHO v18=1')

plt.plot(X99,Y99,color='b',lw=1.0,ls='-',label='CH3OCHO v=0')

plt.legend(loc='upper left')

plt.xlabel('Freq(MHZ)')

plt.ylabel('T(K)')

plt.text(159110,13,'IRAS16293B')

plt.plot([158689,158689],[7,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (158689, 3.5), xytext = (158700,4),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.plot([158963,158963],[2.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v=0", xy = (158694,6), xytext = (158760,9),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

#plt.plot([158963,158963],[2.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v=0", xy = (158704,9.5), xytext = (158770,12),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([158963,158963],[2.5,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3COCH3", xy = (158963, 3.5), xytext = (158980,6),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.plot([159030,159030],[2.4,15], color ='grey', linewidth=1, linestyle="--")

plt.annotate("CH3OCHO v18=1", xy = (159030,4), xytext = (159000,9),arrowprops = dict(facecolor = "y", headlength=10, headwidth =3, width =0.1))

plt.show()

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

python画图代码的输入数据可以取出来_python 导入数据及作图的实现相关推荐

1. python画图代码的输入数据可以取出来_用Python写了个小程序：最小二乘法、读取文件、作图以及数据输出到文件...

   #最小二乘法 import numpy as np #numpy是一个基于python的基础的科学计算包,在本代码中我们会用它来实现方差和标准偏差的计算 from decimal import * # ...

2. 小猪的Python学习之旅 —— 20.抓取Gank.io所有数据存储到MySQL中

   小猪的Python学习之旅 -- 20.抓取Gank.io所有数据存储到MySQL中 标签:Python 一句话概括本文: 内容较多,建议先mark后看,讲解了一波MySQL安装,基本操作,语法速成, ...

3. 灰色模型代码GM(1,1)，从excel导入数据，亦可导出数据到excel中。

   灰色模型代码GM(1,1),从excel导入数据,亦可导出数据到excel中. 总结1:直接输出结果 %clc W= MicrosoftExcel;(从EXCEL导入数据后,matlab自动储存的数组 ...

4. 用python画图代码意思_Python科学画图代码分享

   Python画图主要用到matplotlib这个库.Matplotlib 是一个 Python 的 2D绘图库,它以各种硬拷贝格式和跨平台的交互式环境生成出版质量级别的图形. 这里有一本电子书供大家参 ...

5. 用python画图代码-Python科学画图代码分享

   Python画图主要用到matplotlib这个库.Matplotlib 是一个 Python 的 2D绘图库,它以各种硬拷贝格式和跨平台的交互式环境生成出版质量级别的图形. 具体来说是pylab和p ...

6. python画图代码-Python实战小程序利用matplotlib模块画图代码分享

   Python中的数据可视化 matplotlib 是python最著名的绘图库,它提供了一整套和matlab相似的命令API,十分适合交互式地进行制图.而且也可以方便地将它作为绘图控件. 实战小程序: ...

7. 用python画图代码-Python实战小程序利用matplotlib模块画图代码分享

   Python中的数据可视化 matplotlib 是python最著名的绘图库,它提供了一整套和matlab相似的命令API,十分适合交互式地进行制图.而且也可以方便地将它作为绘图控件. 实战小程序: ...

8. python画图代码大全-Python科学画图代码分享

   Python画图主要用到matplotlib这个库.Matplotlib 是一个 Python 的 2D绘图库,它以各种硬拷贝格式和跨平台的交互式环境生成出版质量级别的图形. 这里有一本电子书供大家参 ...

9. python 画图代码及解释_python画图(图例、网格、注释) - v2

   每件事情总有两面性,当我们看到不太好的一面时,一定要想着去寻找一下好的一面.就比如,原本计划好的五一行程被打乱之后,我脑子里便诞生了一个"伟大"的想法,立马买好了乳胶漆打算利用闲暇 ...

最新文章

1. SQLite.NET.0.17 的离奇之处, BUG??? BY DESIGN??
2. python3的float数精度_Python numpy 浮点数精度问题
3. android搜索框功能实现_Android实现滑动解锁功能
4. openpyxl删除添加excel列_Python | 如何使用Python操作Excel（二）
5. Linux6新建分区,centos6中添加一块新的硬盘并分区的方法介绍
6. MFC开发之显示网页
7. SSD1306 OLED驱动芯片 详细介绍
8. 用代码实现自反闭包，对称闭包，传递闭包
9. Android 仿微信全局字体大小调整
10. series 合并pandas_在python中pandas的series合并方法
11. antd 覆盖css样式不生效（antd避坑）
12. MySQL之函数和谓词
13. C语言基础——1024我把函数献给你
14. python swapcase用法_python中swapcase是什么意思
15. 软件功能介绍之(数据维护)3.1数据编辑（1）
16. 微信论坛交流小程序系统毕业设计毕设（6）开题答辩PPT
17. onMeasure实例分析
18. SpringMVC几个核心类(控制器核心类,加载配置文件核心类,处理url影射核心类,处理视图资源核心类,方法动态调用核心类)
19. Backup Exec安装
20. 今日头条信息流 - 创意设置详解

热门文章

1. java long初始化_Java变量的初始化
2. 紫禁城里“冷宫”在哪儿? 故宫的六大宫廷之迷
3. TinyP2P如何运行 How tinyp2p works!
4. linux操作系统安装man,在deepin linux系统终端中安装man的中文手册
5. 【笔记】5、初学python3网络爬虫——正则表达式的基本使用
6. 【CSDN竞赛第26期】赢实体图书盲盒和高级定制背包等周边
7. CATIA二次开发CAA：（一）开发环境搭建
8. Android TV Leanback (七)（详情视图）
9. 什么是硬盘？什么是分区？
10. MAC下IDEA快捷键