用Python做SVD分解(包含EOF)，包含源代码与例子（11/06/2021 更新）

xiaofeixiazyh

学习一下

勿忘草

Yang.hahha 发表于 2019-6-14 20:56
你是怎么安装的python??

安装anaconda，安装成功后python就包含在里面。
anaconda网址：http://repo.continuum.io/archive/选择符合你系统的版本下载安装

铁狼

收藏了，请问 windows平台上和linux上都可以用么？

xiaofeixiazyh

学习一下

青山

很好啊，学习学习。。

失落的积雨云

勿忘草 发表于 2019-6-12 14:44
python安装成功了， xMCA 这个包安装一直有问题

这个不能直接import那个.py的程序，必须得先安装包？

失落的积雨云

请问一个nc文件里有很多变量，这种该怎么操作啊？可以用netcdf4或者numpy打开的数据来分析吗？

大果子

感谢 大神的分享！ 不过我按照这个代码运行的时候提示 ：  File "D:\Anaconda3\lib\site-packages\xmca-0.1-py3.7.egg\xMCA\core\xMCA.py", line 44, in __init__raise  TypeError('Left and Right field must be xarray DataArray.')TypeError: Left and Right field must be xarray DataArray.

请问知道是什么原因吗？


以下是我的代码：
from xMCA import xMCA
import numpy as np
import xarray as xr
#读入数据
data1=np.loadtxt("SVD-X.txt")
datal=np.array(data1).reshape(2790,42)

data2=np.loadtxt("SVD-Y.txt")
datar=np.array(data2).reshape(2790,42)

#SVD分解，并获取前4模态、其时间扩张系数以及协方差解释的百分比
print('ok')
mca = xMCA(datal, datar)
mca.solver()

lp, rp = mca.patterns(n=4) #前4模态， 左场，右场
le, re = mca.expansionCoefs(n=4) #时间扩张系数
frac = mca.covFracs(n=4) # 协方差解释百分比
print(frac)

left = open('left.txt', 'w')       #左场特征向量（空间模态）
np.savetxt(left,lp)
right = open('right.txt', 'w')     #右场特征向量（空间模态）
np.savetxt(right,rp)

t_l=open('t_l.txt', 'w')              #左场时间系数
np.savetxt(t_l,le)
t_r=open('t_r.txt', 'w')               #右场时间系数
np.savetxt(t_r,re)

superyoung16

大果子 发表于 2019-9-30 11:48
感谢 大神的分享！ 不过我按照这个代码运行的时候提示 ：  File "D:\Anaconda3\lib\site-packages\xmca-0.1 ...

数据类型问题，你load进来转成array是不行的，建议转成nc文件，用xr.open_dataarray打开

蝶恋花

今天更新了统计检验，希望对大家有所帮助。