爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

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

MeteoInfoLab脚本示例:EOF分析

[复制链接]

新浪微博达人勋

 楼主| 发表于 2022-11-16 08:30:02 | 显示全部楼层
guo066600 发表于 2022-11-15 21:41
王老师,EOF时报错信息如下
\MeteoInfo_3.3.8\MeteoInfo\pylib\mipylib\numeric\linalg\linalg.py", line  ...

参考这里:http://www.meteothink.org/docs/i ... nfomap-meteoinfolab
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-11-16 14:52:21 | 显示全部楼层
本帖最后由 guo066600 于 2022-11-16 14:58 编辑

王老师,您的示例txt文件就能正常运行出图,是不是说明Java版本没有问题?
但是换成其他数据就报错\MeteoInfo\pylib\mipylib\numeric\linalg\linalg.py", line 207, in svd
    r = LinalgUtil.svd(a.asarray())
        at org.meteoinfo.math.matrix.Matrix.<init>(Matrix.java:81)
        at org.meteoinfo.math.matrix.Matrix.<init>(Matrix.java:63)
        at org.meteoinfo.math.matrix.Matrix.svd(Matrix.java:1862)
        at org.meteoinfo.math.matrix.Matrix.svd(Matrix.java:1835)
        at org.meteoinfo.math.linalg.LinalgUtil.svd(LinalgUtil.java:160)
        at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
        at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)
        at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)
        at java.lang.reflect.Method.invoke(Unknown Source)
因为SVD分解时对数据解析出了问题?数据中有一些为0或者nan,是研究区外没有数据的地方,会有影响吗?示例的txt文件里没有0值
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2022-11-16 17:05:48 | 显示全部楼层
guo066600 发表于 2022-11-16 14:52
王老师,您的示例txt文件就能正常运行出图,是不是说明Java版本没有问题?
但是换成其他数据就报错\MeteoI ...

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

使用道具 举报

新浪微博达人勋

发表于 2022-11-25 21:39:36 | 显示全部楼层
王老师,请问您MI软件可以绘制站点的EOF图吗?还是需要插值成格点再进行EOF?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2022-11-25 23:17:22 | 显示全部楼层
guo066600 发表于 2022-11-25 21:39
王老师,请问您MI软件可以绘制站点的EOF图吗?还是需要插值成格点再进行EOF?

应该是可以的
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-11-26 20:00:04 | 显示全部楼层
王老师,如果数据的范围是某个区域,比如全国的降水,0.5*0.5格点。代码只在出图阶段改变范围吗?
#Plot
lon = linspace(0, 360, nx)
lat = linspace(-90, 90, ny)
还是需要在数据处理
ss1 = asciiread(fn, shape=(71,144,n))
ss1 = ss1[::-1,::1,:]
X = ss1.reshape(ny * nx, n)
以后加什么语句限定?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2022-11-27 11:21:47 | 显示全部楼层
guo066600 发表于 2022-11-26 20:00
王老师,如果数据的范围是某个区域,比如全国的降水,0.5*0.5格点。代码只在出图阶段改变范围吗?
#Plot
...

没看明白你的问题
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

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

本版积分规则

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

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

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