這一部分涉及到了常用的操作，比如調(diào)換維度的位置，給數(shù)據(jù)重新reshape換形等等，建議大家可以認(rèn)真閱讀這部分。老樣子，先新建一個(gè)數(shù)組 ex1調(diào)換維度順序 比如說(shuō)在求某個(gè)東西時(shí)...
發(fā)簡(jiǎn)信
IP屬地：內(nèi)蒙古
-
Python氣象數(shù)據(jù)處理進(jìn)階之Xarray(6)：數(shù)據(jù)重組與換形
-
Python氣象數(shù)據(jù)處理進(jìn)階之Xarray(4)：計(jì)算（完結(jié)）
最近開(kāi)始忙了。更新的速度會(huì)放緩一些，見(jiàn)諒。這部分主要講一下Xarray內(nèi)部的計(jì)算函數(shù)，可以直接對(duì)dataarray進(jìn)行操作。還是以舉例子的方式來(lái)展示。先創(chuàng)建一個(gè)dataarr...
-
Python氣象數(shù)據(jù)處理進(jìn)階之Xarray(3)：插值Xarray庫(kù)關(guān)于插值的操作內(nèi)容不多，但是很實(shí)用，需要安裝scipy庫(kù)。首先我們還是創(chuàng)建一個(gè)DataArray: exp1 我們先對(duì)時(shí)間維插值，比如說(shuō)這里有第0天，1天，2天...
-
Python氣象數(shù)據(jù)處理進(jìn)階之Xarray(2)：如何優(yōu)雅的索引數(shù)據(jù)在了解了數(shù)據(jù)的結(jié)構(gòu)之后，我們便可以開(kāi)始著手處理數(shù)據(jù)了，那么如何才能從一個(gè)數(shù)據(jù)或數(shù)據(jù)集中拿到我們所需要的那部分?jǐn)?shù)據(jù)呢？Xarray官方提供了三種方法用來(lái)索引數(shù)據(jù)：1.利用下標(biāo)索...
-
Python氣象數(shù)據(jù)處理進(jìn)階之Xarray(1)：Xarray的數(shù)據(jù)結(jié)構(gòu)首先聲明的是前一系列并沒(méi)有完結(jié)，只是暫時(shí)新開(kāi)一個(gè)新系列。關(guān)于畫(huà)圖部分目前也沒(méi)有什么太好的想法，所以暫時(shí)擱置，等有想法了會(huì)繼續(xù)寫(xiě)下去的。考慮到目前大家最感興趣的還是數(shù)據(jù)讀取和處...
-
Python氣象數(shù)據(jù)處理與繪圖(17)：加快循環(huán)的運(yùn)算速度
最近比較忙好幾天沒(méi)有更新，今天就更新一個(gè)簡(jiǎn)單且實(shí)用的數(shù)據(jù)處理方法吧。在使用python時(shí)大家會(huì)發(fā)現(xiàn)有時(shí)候程序運(yùn)行的異常緩慢。數(shù)據(jù)讀取環(huán)節(jié)是我們不能控制的，那么只有在數(shù)據(jù)處理環(huán)...
-
Python氣象數(shù)據(jù)處理與繪圖(16)：去趨勢(shì)與濾波在最開(kāi)始，我首先給出一個(gè)原始序列，序列并沒(méi)有任何意義，藍(lán)線(xiàn)為原始序列，黃線(xiàn)為其線(xiàn)性回歸，可以看到序列有一個(gè)明顯的線(xiàn)性上升的趨勢(shì)，并且序列有很強(qiáng)的年際變化和年代際變化特征。 去...
-
Python氣象數(shù)據(jù)處理與繪圖(15)：兩種波作用通量計(jì)算的python實(shí)現(xiàn)及對(duì)比(Plumb & T-N) （已更正）大氣動(dòng)力學(xué)中通常用波作用通量來(lái)診斷 Rossby波的傳播。常用的三種波作用通量分別為局地E-P 通量，Plumb 波作用通量和T-N 波作用通量。局地E-P 通量可以診斷一段...
-
通過(guò)Python語(yǔ)言實(shí)現(xiàn)大氣科學(xué)研究中常用的統(tǒng)計(jì)方法與氣象數(shù)據(jù)可視化
-
Python氣象數(shù)據(jù)處理與繪圖(14)：封裝繪圖函數(shù)(簡(jiǎn)潔繪制帶地圖子圖的代碼)
Cartopy雖然對(duì)地理圖形的繪制提供了極大的方便，但是個(gè)人感覺(jué)還是有很多反人類(lèi)的地方，尤其是在同一個(gè)代碼中繪制多副圖形的時(shí)候，比如說(shuō)地理坐標(biāo)軸的設(shè)置，每幅子圖都要設(shè)置一遍，...
-
Python氣象數(shù)據(jù)處理與繪圖(13)：聚類(lèi)算法(K-means軌跡聚類(lèi))前一篇文章講到了軌跡的繪制，那么今天就順著講下軌跡的聚類(lèi)，比如說(shuō)我們常見(jiàn)的，寒潮研究把寒潮的冷空氣路徑分為三類(lèi)(西路，北路，西北路)，或者臺(tái)風(fēng)研究中也可以根據(jù)臺(tái)風(fēng)的速度，或者...
-
Python氣象數(shù)據(jù)處理與繪圖(12)：軌跡(臺(tái)風(fēng)路徑，寒潮路徑，水汽軌跡)繪制寒潮是筆者主要的研究方向，寒潮路徑作為寒潮重要的特征，是寒潮預(yù)報(bào)的重點(diǎn)之一，同樣的道理也適用在臺(tái)風(fēng)研究以及降水的水汽來(lái)源研究中。關(guān)于路徑的計(jì)算以及獲取方法(比如軌跡倒推，模型...
-
Python氣象數(shù)據(jù)處理與繪圖(11)：矢量箭頭圖(風(fēng)場(chǎng)，通量場(chǎng))
最近在計(jì)算波作用通量，正好就畫(huà)到了矢量箭頭圖，畫(huà)的過(guò)程中發(fā)現(xiàn)還是有很多細(xì)節(jié)需要注意的，那就直接進(jìn)入正題吧。首先，矢量箭頭圖是分為兩種的，一種類(lèi)似我們常見(jiàn)的風(fēng)場(chǎng)（左圖），另一種...
-
Python氣象數(shù)據(jù)處理與繪圖(10)：SVD分析(奇異值分解)猶豫了很久要不要講SVD的python實(shí)現(xiàn)，今天還是寫(xiě)了吧，糾結(jié)的原因在于我對(duì)SVD也是新手理解，很多東西也是一知半解，怕誤導(dǎo)大家。所以這篇文章也只是我個(gè)人的理解，不對(duì)之處還...
-
Python氣象數(shù)據(jù)處理與繪圖(9)：更自由的多子圖組圖繪制今天分享一個(gè)比較實(shí)用的多子圖繪制方法，由于氣象科研中經(jīng)常需要將多圖拼成組圖，并且需要將多種類(lèi)型的圖拼在一起(比如EOF，將填色圖和折線(xiàn)圖拼在一起)，雖說(shuō)在前幾篇文章均有子圖的...
-
Python氣象數(shù)據(jù)處理與繪圖(8)：緯高圖(對(duì)數(shù)坐標(biāo)軸)鴿了好幾天，一直沒(méi)想好更新什么，翻了下以前畫(huà)過(guò)的一些圖，突然想到緯高圖也是氣象研究中心很常見(jiàn)的一類(lèi)圖型，通常用于分析高低空的環(huán)流配置，垂直運(yùn)動(dòng)等等，那么今天就更新一下。緯高圖...
-
-
Python氣象數(shù)據(jù)處理與繪圖(7)：站點(diǎn)氣象數(shù)據(jù)繪圖(散點(diǎn)圖)站點(diǎn)數(shù)據(jù)也是氣象數(shù)據(jù)很重要的一部分，尤其是在觀測(cè)數(shù)據(jù)方面，站點(diǎn)數(shù)據(jù)一般存放于.txt(.csv)文件中，在"Python氣象數(shù)據(jù)處理與繪圖(1)：數(shù)據(jù)讀取"中提到過(guò)這類(lèi)文件的...
-
Python氣象數(shù)據(jù)處理與繪圖(6)：極地投影繪圖
當(dāng)我們研究高緯(極地)區(qū)域時(shí)，往往需要極地投影。在前文中通過(guò)EOF介紹了圓柱投影，實(shí)際上極地投影的繪制僅需改變其中幾個(gè)參數(shù)便可以實(shí)現(xiàn)。還是先上圖： 填色圖的意義并不重要(我之...
-
Python氣象數(shù)據(jù)處理與繪圖(5)：氣候突變檢驗(yàn)(年代際突變檢驗(yàn))長(zhǎng)時(shí)間尺度的氣候研究中通常涉及到年際變化和年代際變化。文獻(xiàn)中也經(jīng)常出現(xiàn)所謂年代際突變這類(lèi)的描述，這次便介紹兩個(gè)用來(lái)檢驗(yàn)?zāi)甏H變化的方法， 其一是滑動(dòng)t檢驗(yàn)(sliding t...
