MeteoInfo顯著性檢驗(yàn)：

MeteoInfo官網(wǎng)：

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')levs = arange(-0.02, 0.021, 0.002)layer = contourfm(lon, lat, eof1, levs,cmap='WhBlGrYeRe',smooth=False)scatter(lon, lat, eof2, facecolor='k',size=2,edgecolor=None)title('EOF mode 1')xlim(0, 360)ylim(-90, 90)xticks(arange(0, 361, 30),bold=True,fontsize=15)yticks(arange(-90, 91, 30),bold=True,fontsize=15)colorbar(layer, orientation='horizontal',ticklen=0,extendrect=False, shrink=1, aspect=50)yaxis(tickvisible=True,location='left',tickwidth=2,linewidth=2,ticklength=3) #ticklength刻度線長度，tickwidth刻度線寬度，linewidth邊框?qū)挾葃axis(tickvisible=False,location='right',tickwidth=2,linewidth=2,ticklength=4) #分別調(diào)試tick的寬度，邊框線寬和tick的長度xaxis(tickvisible=False,location='top',tickwidth=2,linewidth=2,ticklength=4) xaxis(tickvisible=True,location='bottom',tickwidth=2,linewidth=2,ticklength=3)savefig('F:/RMeteoInfo/plot58.3.png', dpi=1200)

---

Python顯著性檢驗(yàn)：

往期推文超鏈接1：

《R語言、MeteoInfo、Python和ArcGis的Kriging、IDW空間插值結(jié)果的對比分析》

往期推文超鏈接2：

1 《Python基礎(chǔ)地圖構(gòu)建(九)》

2 《Python基礎(chǔ)地圖構(gòu)建(33)》

3 《Python基礎(chǔ)地圖構(gòu)建(28)》

4 《基于Python的NCEP再分析數(shù)據(jù)的中國區(qū)域白化(含南海小地圖)》

總結(jié)

以上是生活随笔為你收集整理的python 白化_MeteoInfo和Python显著性检验的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。