爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
12
返回列表 发新帖
楼主: MeteoInfo

MeteoInfoLab脚本示例:subplot

[复制链接]

新浪微博达人勋

发表于 2016-10-21 14:29:12 | 显示全部楼层
MeteoInfo 发表于 2016-10-21 11:44
bar图没有给出x值,软件会给缺省值。而线条图的x值是时间,这样两个图的x值不匹配。需要在绘制bar图的时 ...


王老师,您好!开始也考虑过这个问题,但是加上bar(times, b, color='b'),脚本运行就会出错了,而且右边的y轴范围发生变化,ylim(0,60)以及ylabel('rain', color='k')不起任何作用,这是怎么回事呢?如下图:
捕获.PNG 捕获2.PNG
#Read lighting data
fn = 'f:/shandianshujuchuli/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')
bar(times,b,0.3,color='b')
ylim(0,60)
#xlim(times[0], times[-1])
ylabel('rain', color='k')
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-10-21 15:01:05 | 显示全部楼层
wuwei2163 发表于 2016-10-21 14:29
王老师,您好!开始也考虑过这个问题,但是加上bar(times, b, color='b'),脚本运行就会出错了,而且右 ...

那是因为你的数组b的个数和times不一致造成的。改为(给b加了两个元素使其和times长度一致):

b=[10,15,18,20,50,8,8,20,50,8,20,50,8,8,20,10,15,18,20,50,8,9,15,30]
bar(times,b,color='b')
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-10-21 15:40:52 | 显示全部楼层
嗯,就是这个问题,两个数组不一致导致的。我太粗心了,调试的时候没把编写的样本数据弄好,谢谢王老师提醒~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2019-7-16 21:21:05 | 显示全部楼层

RE: MeteoInfoLab脚本示例

Nested pie添加Legend问题:王老师,您好,请问Nested pie图如何添加Legend的每个颜色的对应名称?目前只能显示出起始颜色
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-7-17 08:54:46 | 显示全部楼层
a.png 王老师,meteoinfo1.8.6用shp文件裁剪shaded图的时候总报错,试了以前的版本没有这个问题,向您求助。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-7-17 10:55:40 | 显示全部楼层
梦纯 发表于 2019-7-17 08:54
王老师,meteoinfo1.8.6用shp文件裁剪shaded图的时候总报错,试了以前的版本没有这个问题,向您求助。

你能不能提供一个我能够重现这个bug的例子,以便今后改进。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-7-17 17:26:18 | 显示全部楼层
王老师,附件是数据和我用的shp文件。

bs_t.csv

10.9 KB, 下载次数: 1, 下载积分: 金钱 -5

jjj_merge_shp.rar

31.16 KB, 下载次数: 1, 下载积分: 金钱 -5

密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-7-17 18:18:14 | 显示全部楼层
梦纯 发表于 2019-7-17 17:26
王老师,附件是数据和我用的shp文件。

你是怎么操作的?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-7-19 08:41:10 | 显示全部楼层
MeteoInfo 发表于 2019-7-17 18:18
你是怎么操作的?

先把csv数据shaded画图,然后对shaded的图层用shp文件进行裁剪。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表