python 导入数据及作图的实现

2023-08-05 01:59:04 64

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

方法一、

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

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

withopen(filename,'r')asf:#打开文档
lines=f.readlines()#按行读取数据

foriinlines:

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

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

Y=append(value[1])

Z=append(value[2])

方法二、

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

a[:,0]

方法三、

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

附录：IRAS16293的数据绘图代码

importmatplotlib
importmatplotlib.pyplotasplt
importnumpyasnp

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


#data_a
filename='/home/luoluo/Documents/python_scripts/Adata_cut/A1_145122-145580.txt'
X10,Y10=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.split()]
X19.append(value[0])
Y19.append(value[1])



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



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='CH3OCHOv18=1')
plt.plot(X14,Y14,color='b',lw=1.0,ls='-',label='CH3OCHOv=0')
plt.legend(loc='upperleft')
#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='CH3OCHOv18=1')
plt.plot(X19,Y19,color='b',lw=1.0,ls='-',label='CH3OCHOv=0')
plt.legend(loc='upperleft')
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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.split()]
X24.append(value[0])
Y24.append(value[1])

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

filename='/home/luoluo/Documents/python_scripts/Bdata_cut/B2.txt'
X25,Y25=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.split()]
X29.append(value[0])
Y29.append(value[1])



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



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='CH3OCHOv18=1')
plt.plot(X24,Y24,color='b',lw=1.0,ls='-',label='CH3OCHOv=0')
plt.legend(loc='upperleft')
#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("CH3OCHOv18=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("CH3OCHOv18=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='CH3OCHOv18=1')
plt.plot(X29,Y29,color='b',lw=1.0,ls='-',label='CH3OCHOv=0')
plt.legend(loc='upperleft')
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("CH3OCHOv18=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("CH3OCHOv18=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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.split()]
X39.append(value[0])
Y39.append(value[1])



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



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='CH3OCHOv18=1')
plt.plot(X34,Y34,color='b',lw=1.0,ls='-',label='CH3OCHOv=0')
plt.legend(loc='upperleft')
#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("CH3OCHOv18=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("CH3OCHOv18=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("CH3OCHOv18=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("CH3OCHOv18=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='CH3OCHOv18=1')
plt.plot(X39,Y39,color='b',lw=1.0,ls='-',label='CH3OCHOv=0')
plt.legend(loc='upperleft')
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("CH3OCHOv18=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("CH3OCHOv18=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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.split()]
X49.append(value[0])
Y49.append(value[1])



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



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='CH3OCHOv18=1')
plt.plot(X44,Y44,color='b',lw=1.0,ls='-',label='CH3OCHOv=0')
plt.legend(loc='upperleft')
#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='CH3OCHOv18=1')
plt.plot(X49,Y49,color='b',lw=1.0,ls='-',label='CH3OCHOv=0')
plt.legend(loc='upperleft')
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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.split()]
X59.append(value[0])
Y59.append(value[1])



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



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='CH3OCHOv18=1')
plt.plot(X54,Y54,color='b',lw=1.0,ls='-',label='CH3OCHOv=0')
plt.legend(loc='upperleft')
#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='CH3OCHOv18=1')
plt.plot(X59,Y59,color='b',lw=1.0,ls='-',label='CH3OCHOv=0')
plt.legend(loc='upperleft')
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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.split()]
X69.append(value[0])
Y69.append(value[1])



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



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='CH3OCHOv18=1')
plt.plot(X64,Y64,color='b',lw=1.0,ls='-',label='CH3OCHOv=0')
plt.legend(loc='upperleft')
#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='CH3OCHOv18=1')
plt.plot(X69,Y69,color='b',lw=1.0,ls='-',label='CH3OCHOv=0')
plt.legend(loc='upperleft')
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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.split()]
X79.append(value[0])
Y79.append(value[1])
#############################################################################################################################################
#2plot
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='CH3OCHOv18=1')
plt.plot(X74,Y74,color='b',lw=1.0,ls='-',label='CH3OCHOv=0')
plt.legend(loc='upperleft')
#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("CH3OCHOv18=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='CH3OCHOv18=1')
plt.plot(X79,Y79,color='b',lw=1.0,ls='-',label='CH3OCHOv=0')
plt.legend(loc='upperleft')
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("CH3OCHOv18=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("CH3OCHOv18=1",xy=(158390,5),xytext=(158400,13),arrowprops=dict(facecolor="y",headlength=10,headwidth=3,width=0.1))


plt.show()
######################################################################################
############################################################################################8
######################################################################################

#
#1importdata

#data_a
filename='Adata_cut/A8_158228-158680.txt'
X80,Y80=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.split()]
X84.append(value[0])
Y84.append(value[1])
#########################################################
#data_b
filename='Bdata_cut/B8_158228-158680.txt'
X85,Y85=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.split()]
X89.append(value[0])
Y89.append(value[1])
#############################################################
#2plot
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='CH3OCHOv18=1')
plt.plot(X84,Y84,color='b',lw=1.0,ls='-',label='CH3OCHOv=0')
plt.legend(loc='upperleft')
#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("CH3OCHOv18=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("CH3OCHOv18=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("CH3OCHOv18=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("CH3OCHOv18=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("CH3OCHOv18=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("CH3OCHOv18=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("CH3OCHOv18=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='CH3OCHOv18=1')
plt.plot(X89,Y89,color='b',lw=1.0,ls='-',label='CH3OCHOv=0')
plt.legend(loc='upperleft')
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("CH3OCHOv18=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("CH3OCHOv18=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("CH3OCHOv18=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("CH3OCHOv18=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("CH3OCHOv18=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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.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=[],[]
withopen(filename,'r')asf:
lines=f.readlines()
forlineinlines:
value=[float(s)forsinline.split()]
X99.append(value[0])
Y99.append(value[1])





###2plot
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='CH3OCHOv18=1')
plt.plot(X94,Y94,color='b',lw=1.0,ls='-',label='CH3OCHOv=0')
plt.legend(loc='upperleft')
#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("CH3OCHOv18=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("CH3OCHOv18=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("CH3OCHOv18=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("CH3OCHOv18=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='CH3OCHOv18=1')
plt.plot(X99,Y99,color='b',lw=1.0,ls='-',label='CH3OCHOv=0')
plt.legend(loc='upperleft')
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("CH3OCHOv18=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("CH3OCHOv=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("CH3OCHOv=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("CH3OCHOv18=1",xy=(159030,4),xytext=(159000,9),arrowprops=dict(facecolor="y",headlength=10,headwidth=3,width=0.1))

plt.show()

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

声明：本文内容来源于网络，版权归原作者所有，内容由互联网用户自发贡献自行上传，本网站不拥有所有权，未作人工编辑处理，也不承担相关法律责任。如果您发现有涉嫌版权的内容，欢迎发送邮件至：czq8825#qq.com（发邮件时，请将#更换为@）进行举报，并提供相关证据，一经查实，本站将立刻删除涉嫌侵权内容。

python 导入数据及作图的实现

热门推荐

随机推荐