python 导入数据及作图的实现
作者:csluo19550418
今天小编就为大家分享一篇python 导入数据及作图的实现,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
我们经常需要导入数据,按列提取 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 导入数据及作图的实现就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持脚本之家。