MeteoInfoLab脚本示例：EOF分析

MeteoInfo

guo066600 发表于 2022-11-15 16:16
sample在哪里？
这又这些index modules |next |home| news| examples| downloads| docs| Introduction

MeteoInfo软件所在的目录中，在你的电脑硬盘里，不在MeteoInfo网站。

guo066600

MeteoInfo 发表于 2022-11-15 18:46
MeteoInfo软件所在的目录中，在你的电脑硬盘里，不在MeteoInfo网站。

找到了，谢谢王老师

guo066600

王老师，EOF时报错信息如下
\MeteoInfo_3.3.8\MeteoInfo\pylib\mipylib\numeric\linalg\linalg.py", line 207, in svd
    r = LinalgUtil.svd(a.asarray())
java.lang.OutOfMemoryError: Java heap space

java.lang.OutOfMemoryError: java.lang.OutOfMemoryError: Java heap space

是计算空间问题吗？还是我的数据经纬度格点数设置不对？我调换nx和ny都试过了

MeteoInfo

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

guo066600

王老师，您的示例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值

MeteoInfo

guo066600 发表于 2022-11-16 14:52
王老师，您的示例txt文件就能正常运行出图，是不是说明Java版本没有问题？
但是换成其他数据就报错\MeteoI ...

信息太少

guo066600

王老师，请问您MI软件可以绘制站点的EOF图吗？还是需要插值成格点再进行EOF？

MeteoInfo

guo066600 发表于 2022-11-25 21:39
王老师，请问您MI软件可以绘制站点的EOF图吗？还是需要插值成格点再进行EOF？

应该是可以的

guo066600

王老师，如果数据的范围是某个区域，比如全国的降水，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)
以后加什么语句限定？

MeteoInfo

guo066600 发表于 2022-11-26 20:00
王老师，如果数据的范围是某个区域，比如全国的降水，0.5*0.5格点。代码只在出图阶段改变范围吗？
#Plot
...

没看明白你的问题