MeteoInfoLab脚本示例：误差棒填色图

︶ㄣ安定■丶

参考Origin官网 Error Bars with Fill Area的例子（http://www.originlab.com/doc/Tutorials/ErrorBar-FillArea），使用MILab绘制了一下误差棒填色图，绘图数据参考上述链接。


脚本程序如下：

1. #add data
   
2. fn = r'H:\test\error bars with Fill Area\Error Bars with Fill Area.txt'
   
3. ncol = numasciicol(fn)
   
4. nrow = numasciirow(fn)
   
5. a = asciiread(fn,shape=(nrow,ncol))
   
6. a[a==-999]= nan
   
7. data = a[:,0]
   
8. error = a[:,1]
   
9. data1 = a[:,2]
   
10. error1 = a[:,3]
    
11. data2 = a[:,4]
    
12. error2 = a[:,5]
    
13. x = arange(len(data))
    
14. 
    
15. #plot
    
16. axes(tickfontsize=18)
    
17. antialias(True)
    
18. #errorbar(x,data,error,color='r',marker='D',linewidth=3)
    
19. #errorbar(x,data1,error1,color='b',marker='S',linewidth=3)
    
20. #errorbar(x,data2,error2,color='m',marker='D',linewidth=3)
    
21. line1 = plot(x,data,color='g',linewidth=2,label='Effective Force')
    
22. line2 = plot(x,data1,color='r',linewidth=2,label='ineffective Force')
    
23. line3 = plot(x,data2,color='b',linewidth=2,label='Total Force')
    
24. fill_between(x,data+error,data-error,color='g', alpha=0.3)
    
25. fill_between(x,data1+error1,data1-error1,color='r', alpha=0.3)
    
26. fill_between(x,data2+error2,data2-error2,color='b', alpha=0.3)
    
27. line = [line1,line2,line3]
    
28. legend(line,fontsize=18,markerscale=2,framon=True)
    
29. xlim(0,350)
    
30. ylim(-300,450)
    
31. 
    
32. 
    

#add data<br /> fn = r'H:\test\error bars with Fill Area\Error Bars with Fill Area.txt'<br /> ncol = numasciicol(fn)<br /> nrow = numasciirow(fn)<br /> a = asciiread(fn,shape=(nrow,ncol))<br /> a[a==-999]= nan<br /> data = a[:,0]<br /> error = a[:,1]<br /> data1 = a[:,2]<br /> error1 = a[:,3]<br /> data2 = a[:,4]<br /> error2 = a[:,5]<br /> x = arange(len(data))<br /> <br /> #plot<br /> axes(tickfontsize=18)<br /> antialias(True)<br /> #errorbar(x,data,error,color='r',marker='D',linewidth=3)<br /> #errorbar(x,data1,error1,color='b',marker='S',linewidth=3)<br /> #errorbar(x,data2,error2,color='m',marker='D',linewidth=3)<br /> line1 = plot(x,data,color='g',linewidth=2,label='Effective Force')<br /> line2 = plot(x,data1,color='r',linewidth=2,label='ineffective Force')<br /> line3 = plot(x,data2,color='b',linewidth=2,label='Total Force')<br /> fill_between(x,data+error,data-error,color='g', alpha=0.3)<br /> fill_between(x,data1+error1,data1-error1,color='r', alpha=0.3)<br /> fill_between(x,data2+error2,data2-error2,color='b', alpha=0.3)<br /> line = [line1,line2,line3]<br /> legend(line,fontsize=18,markerscale=2,framon=True)<br /> xlim(0,350)<br /> ylim(-300,450)<br /> <br /> <br />

river

赞一个，很漂亮，感谢分享

Nancy_fq

{:5_213:}{:5_213:}

Persistence.9

桂师兄真厉害

爱公爵

谢楼主分享，很给力。

dada678

mark下学习

bianpuxi

为什么我画出来的在fill_between区域的边缘有一条黑色的线，要怎么设置才能取消呢

MeteoInfo

bianpuxi 发表于 2020-3-25 11:07
为什么我画出来的在fill_bwtween区域的边缘有一条黑色的线，要怎么设置才能取消呢

fill_between 函数里增加 edgecolor=None 参数。

bianpuxi

MeteoInfo 发表于 2020-3-25 12:24
fill_between 函数里增加 edgecolor=None 参数。

好的，谢谢王老师指导   

ansion

有没有数据文件下载呀