MeteoInfoLab脚本示例：subplot

MeteoInfo

subplot命令可以自动等间距分配多个坐标系（Axes），命令中有三个参数，前两个定义了行数和列数，第三个指定了当前的坐标系（Axes），绘图命令是作用在当前坐标系中的。


脚本程序：

1. def f(t):
   
2.   return exp(-t) * cos(2*pi*t)
   
3. 
   
4. t1 = arange(0., 5., 0.1)
   
5. t2 = arange(0., 5., 0.02)
   
6. 
   
7. subplot(2,1,1)
   
8. plot(t1, f(t1), 'bo', t2, f(t2), 'k')
   
9. 
   
10. subplot(2,1,2)
    
11. plot(t2, cos(2*pi*t2), 'r--')
    

def f(t):<br />   return exp(-t) * cos(2*pi*t)<br /> <br /> t1 = arange(0., 5., 0.1)<br /> t2 = arange(0., 5., 0.02)<br /> <br /> subplot(2,1,1)<br /> plot(t1, f(t1), 'bo', t2, f(t2), 'k')<br /> <br /> subplot(2,1,2)<br /> plot(t2, cos(2*pi*t2), 'r--')<br />

wuwei2163

王老师，您好！今天在学习教程时，二维坐标轴是实线完全包围封闭，算是一个轴箱（axes box），怎么不能用box off关闭呢

MeteoInfo

wuwei2163 发表于 2016-10-20 13:06
王老师，您好！今天在学习教程时，二维坐标轴是实线完全包围封闭，算是一个轴箱（axes box），怎么不能用bo ...

MeteoInfoLab目前还没有做这样的函数，可以考虑在后续的版本中加上。大家有什么好的建议都可以提出来，帮着MeteoInfoLab不断提高。

wuwei2163

MeteoInfo 发表于 2016-10-20 13:38
MeteoInfoLab目前还没有做这样的函数，可以考虑在后续的版本中加上。大家有什么好的建议都可以提出来，帮 ...

这一点按照Matlab中的box命令就可以吧，二维坐标轴这个轴箱，可以用box off关闭，也可以box on（默认）打开。

MeteoInfo

wuwei2163 发表于 2016-10-20 13:45
这一点按照Matlab中的box命令就可以吧，二维坐标轴这个轴箱，可以用box off关闭，也可以box on（默认）打 ...

已经加了一个box函数，应该比较好用，你如果着急用的话我可以尽快上传一个更新版本。

wuwei2163

MeteoInfo 发表于 2016-10-20 13:53
已经加了一个box函数，应该比较好用，你如果着急用的话我可以尽快上传一个更新版本。

谢谢王老师，我不急着用，现在正根据教程和测试数据不断学习。这个也不是什么大问题，等到大家碰到的问题多了，一并解决再更新版吧~

wuwei2163

MeteoInfo 发表于 2016-10-20 13:38
MeteoInfoLab目前还没有做这样的函数，可以考虑在后续的版本中加上。大家有什么好的建议都可以提出来，帮 ...

王老师，还有个问题想请教您，比如研究站点上的温度以及降水情况，如下图
温度为折线，降水为柱形图，横坐标为站点，图中表示雨量的柱形图跟温度或者站点对应不上？
温度和站点没关系，折线能否变成的平滑的曲线？
ax1 = axes()
xaxis(ax1, color='k')
yaxis(ax1, color='k')
t = arange(1,10)
#s1 = t**3
#plot(t, s1, '-ro')
s2 =[14,15,18,17,12,18,10,11,12,13]
plot(t, s2, '-ro')
xlabel('station')
ylabel('temperature', color='k')
ax2 = twinx(ax1)
yaxis(ax2, color='k')
b=[10,15,18,20,50,8,20,30,16,40]
bar(b,0.2,color='b')
ylim(0,60)
ylabel('rain', color='k')
title('Two y axis example')

MeteoInfo

wuwei2163 发表于 2016-10-20 15:23
王老师，还有个问题想请教您，比如研究站点上的温度以及降水情况，如下图
温度为折线，降水为柱形图， ...

给两个axes都加上xlim语句使其具有相同的x轴范围。

1. ax1 = axes()
   
2. xaxis(ax1, color='k')
   
3. yaxis(ax1, color='k')
   
4. t = arange(1,11)
   
5. #s1 = t**3
   
6. #plot(t, s1, '-ro')
   
7. s2 =[14,15,18,17,12,18,10,11,12,13]
   
8. plot(t, s2, '-ro')
   
9. xlim(0, 11)
   
10. xlabel('station')
    
11. ylabel('temperature', color='k')
    
12. ax2 = twinx(ax1)
    
13. yaxis(ax2, color='k')
    
14. b=[10,15,18,20,50,8,20,30,16,40]
    
15. bar(b,0.2,color='b')
    
16. ylim(0,60)
    
17. xlim(0,11)
    
18. ylabel('rain', color='k')
    
19. title('Two y axis example')

ax1 = axes()<br /> xaxis(ax1, color='k')<br /> yaxis(ax1, color='k')<br /> t = arange(1,11)<br /> #s1 = t**3<br /> #plot(t, s1, '-ro')<br /> s2 =[14,15,18,17,12,18,10,11,12,13]<br /> plot(t, s2, '-ro')<br /> xlim(0, 11)<br /> xlabel('station')<br /> ylabel('temperature', color='k')<br /> ax2 = twinx(ax1)<br /> yaxis(ax2, color='k')<br /> b=[10,15,18,20,50,8,20,30,16,40]<br /> bar(b,0.2,color='b')<br /> ylim(0,60)<br /> xlim(0,11)<br /> ylabel('rain', color='k')<br /> title('Two y axis example')

wuwei2163

MeteoInfo 发表于 2016-10-20 17:28
给两个axes都加上xlim语句使其具有相同的x轴范围。

王老师，可以出来了，我原以为只要对应好站点，第二个axes默认跟第一个axes相同呢，越来越对MeteoInfolab感兴趣了。为了进一步巩固这种不同要素多Y轴画法，我又将前几天的闪电数据跟降水或者温度等要素画一起，发现这种不同要素多Y轴方法还没完全学会，给两个axes都加上xlim语句使其相同的x轴，则不出Bar图；第二个axes不加上xlim语句，两个又对应不上，这是什么原因啊？如下图

#Read lighting data
fn = 'd:/shujulianxi/79.txt'
tf = open(fn)
times = []
nums = []
for h in range(0, 24):
    times.append(datetime.datetime(2016,7,9,h))
    nums.append(0)
for aline in tf:   
    datalist = aline.split()
    h = int(datalist[2].split(':')[0])
    nums[h] = nums[h] + 1
#Plot
ax1 = axes()
xaxis(ax1, color='k')
yaxis(ax1, color='k')
plot(times, nums, '-ro')
xlim(times[0], times[-1])
xaxis(axistype='time')
ylabel('Lighting number')
title('Lighting number hourly variation in ' + times[0].strftime('%Y/%m/%d'))
#bar
ax2 = twinx(ax1)
yaxis(ax2, color='k')
b=[10,15,18,20,50,8,8,20,50,8,20,50,8,8,20,10,15,18,20,50,8,9]
bar(b,0.3,color='b')
ylim(0,60)
xlim(times[0], times[-1])
ylabel('rain', color='k')

MeteoInfo

wuwei2163 发表于 2016-10-21 09:14
王老师，可以出来了，我原以为只要对应好站点，第二个axes默认跟第一个axes相同呢，越来越对MeteoInfolab ...

bar图没有给出x值，软件会给缺省值。而线条图的x值是时间，这样两个图的x值不匹配。需要在绘制bar图的时候也给出x值。

bar(times, b, color='b')

bar函数的帮助文档在这里：http://www.meteothinker.com/docs ... s/plot/bar.html#bar

在Console中输入 >>> help(bar)   也可以看到。