请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
楼主: MeteoInfo

格点插值为站点数据批量处理

  [复制链接]

新浪微博达人勋

 楼主| 发表于 2016-2-17 10:13:45 | 显示全部楼层
zyfws 发表于 2016-2-17 10:01
嗯,解决了,是不是方法名称变了,建议把一楼的示例改一下

嗯,目前addcol方法只是增加一个列,addcoldata增加列及列中的数据。也许应该将两个方法合并为一个addcol方法更简便些。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-2-17 17:33:30 | 显示全部楼层
MeteoInfo 发表于 2016-2-17 10:13
嗯,目前addcol方法只是增加一个列,addcoldata增加列及列中的数据。也许应该将两个方法合并为一个addcol ...

请问一下,如果数据文件缺失,以上脚本会抛出异常提示但是能执行完成,插值结果是站点都为零,我想这种情况下让插值生成无效值(比如NaN或-999),怎么实现呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-2-18 10:32:51 | 显示全部楼层
zyfws 发表于 2016-2-17 17:33
请问一下,如果数据文件缺失,以上脚本会抛出异常提示但是能执行完成,插值结果是站点都为零,我想这种情 ...

如果数据文件不存在,可以创建一个空的一维数组(数组长度为站点数),然后将数组的值赋为nan,再用addcoldata()函数将该数组加入结果表格中即可:

>>> a = zeros(20)
>>> a
array([0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0])
>>> a[:] = nan
>>> a
array([NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN])
>>>
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-2-18 10:38:34 | 显示全部楼层
MeteoInfo 发表于 2016-2-18 10:32
如果数据文件不存在,可以创建一个空的一维数组(数组长度为站点数),然后将数组的值赋为nan,再用addco ...

问题是要自己去判断文件不存在是吧,那么怎么从f=addfile(...)得到的f实例获取文件的路径名称呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-2-18 11:02:18 | 显示全部楼层
zyfws 发表于 2016-2-18 10:38
问题是要自己去判断文件不存在是吧,那么怎么从f=addfile(...)得到的f实例获取文件的路径名称呢?

addfile()里面的参数不就是文件名吗,你在网上搜一下Python判断文件是否存在的例子,如果不存在就不要用addfile()语句了,这样才不会报错。

fn = 'D:/Temp/test.ctl'
if os.path.exists(fn):
    f = addfile(fn)
else:
    ...
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-2-18 11:11:02 | 显示全部楼层
MeteoInfo 发表于 2016-2-18 11:02
addfile()里面的参数不就是文件名吗,你在网上搜一下Python判断文件是否存在的例子,如果不存在就不要用a ...

我知道python如何判断文件存在,但是addfile添加的是ctl文件,我用另一个子程序批量生成ctl文件,但里面指向的dat数据文件不一定存在,所以想在插值子程序里面进行判断。。。如果插值函数能够自动判断无资料就生成缺测值就好了。。。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-2-18 11:14:27 | 显示全部楼层
zyfws 发表于 2016-2-18 11:11
我知道python如何判断文件存在,但是addfile添加的是ctl文件,我用另一个子程序批量生成ctl文件,但里面 ...

批量生成ctl文件时可以判断一下dat文件是否存在,不存在就不要生成ctl文件。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-2-28 14:20:04 | 显示全部楼层
谢谢王老师~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-3-2 21:39:38 | 显示全部楼层
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2016-3-18 20:13:52 | 显示全部楼层
sinohack 发表于 2015-5-11 21:05
再次感谢 老师指教。但是老师,我运行的时候为什么老是报tn=psv.dimlen(0)这一行的错呢?AttributeError: D ...

你好  请问 我知道网格4个点的数据 想根据距离4个格点的距离 将4个点的数据插值到网格中的点上 怎么实现呢?
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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