MeteoInfo显著性检验:
MeteoInfo官网:
http://www.meteothink.org
源代码:
https://github.com/meteoinfo
fn = 'F:/RMeteoInfo/data/test.txt'ny = 71nx = 144m = ny * nxn = 62ss1 = asciiread(fn, shape=(71,144,n))ss1 = ss1[::-1,::1,:]X = ss1.reshape(ny * nx, n)#Time-space transformC = dot(X.T, X)E1, EOF1 = linalg.eig(C)EOF1 = EOF1.T[:,::-1]E = diag(E1[::-1])EOFa = dot(X, EOF1)EOF = zeros((m,n))for i in range(n): EOF[:,i] = EOFa[:,i]/sqrt(E[i,i])PC = dot(EOF.T, X)PC = PC[::-1,:]eof1 = EOF[:,0].reshape(ny, nx)eof2=eof1.copy()eof2[eof2>-0.014]=nan#Plotaxesm()lon = linspace(0, 360, nx)lat = linspace(-90, 90, ny)geoshow('country', edgecolor='k&