NCL--插值

1. 將二維網(wǎng)格數(shù)據(jù)插值到網(wǎng)格。

就是為了將NCEP月資料2.5度網(wǎng)格插值到1度網(wǎng)格,與日資料進(jìn)行相應(yīng)計算。

linint2

;二維插值 lat and lon from 2.5 to 1

u_avg_2=linint2 (lon_0,lat_0,u_avg_1, True, lon_1,lat_1, 0)

lon_0,lat_0:原始經(jīng)緯度夜牡,單調(diào)遞增(劃重點(diǎn),否則會出現(xiàn)缺省值)

lon_1,lat_1:插值經(jīng)緯度,單調(diào)遞增

http://www.ncl.ucar.edu/Document/Functions/Built-in/linint2.shtml

2. 一維插值 p?

linint1_n_Wrap

u_avg=linint1_n_Wrap (p_0,u_avg_2, False, p_1, 0, 0)

p_0:原始?xì)鈮禾磷埃瑔握{(diào)遞增或遞減

p_1:插值氣壓急迂,單調(diào)遞增或遞減(遞增或遞減兩者保持一致)

http://www.ncl.ucar.edu/Document/Functions/Built-in/linint1_n.shtml

3.?曲線網(wǎng)格插值到直線網(wǎng)格

rcm2rgrid和rcm2rgrid_Wrap函數(shù)可以用于將WRF等模式產(chǎn)品插值到標(biāo)準(zhǔn)網(wǎng)格上

官網(wǎng)說明

http://www.ncl.ucar.edu/Document/Functions/Built-in/rcm2rgrid.shtml

emi_grid = rcm2rgrid_Wrap(lat_emi,lon_emi,emi_mean,lat_grid,lon_grid,0)

需要注意的點(diǎn):

標(biāo)準(zhǔn)網(wǎng)格的經(jīng)緯度都是一維的,且為單調(diào)遞增

如果插值前的數(shù)據(jù)中有缺省值蹦肴,必須給出定義emi_mean@_FillValue=-9999

此外關(guān)于插值還有系列函數(shù)僚碎,可以去官網(wǎng)仔細(xì)查看

rgrid2rcm:從直線網(wǎng)格插值到曲線網(wǎng)格

rcm2points:從曲線網(wǎng)格插值到不規(guī)則站點(diǎn)

triple2grid:將不規(guī)則站點(diǎn)數(shù)據(jù)填充到最近的直線網(wǎng)格

最近使用了一下triple2grid,將站點(diǎn)數(shù)據(jù)插值到與再分析資料匹配的網(wǎng)格阴幌。由于我的站點(diǎn)資料非常細(xì)勺阐,option選的false,計算速度很快矛双。

triple2grid2d:將不規(guī)則站點(diǎn)數(shù)據(jù)填充到最近的曲線網(wǎng)格

4. 站點(diǎn)數(shù)據(jù)插值成格點(diǎn)

渊抽;讀取原始數(shù)據(jù)及其經(jīng)緯度

a = asciiread("/public/home/ncl/wrfout/pm25-2013-11-china.txt",(/row,col/),"float")

lat = a(:,1)

lon = a(:,2)

pm25 = a(:,3)

olon=fspan(72,138,132);設(shè)置插值范圍及格點(diǎn)數(shù)议忽,即空間分辨率

olat=fspan(17,57,80)

olon!0???= "lon"

olon@long_name?= "lon"

olon@units?= "degrees-east"

olon&lon??= olon

olat!0???= "lat"

olat@long_name?= "lat"

olat@units?= "degrees_north"

olat&lat??= olat

rscan = (/10,5,3/);連續(xù)的有效半徑大小懒闷,最大為10,依次遞減

pm25_reshape = obj_anal_ic_deprecated(lon,lat,pm25,olon,olat,rscan, False)?;Creanm插值

Creanm插值會對高值有所損耗徙瓶,但空間分布結(jié)果最為合理毛雇。

通過調(diào)整rsan來控制選取的影響范圍嫉称,rsacn數(shù)值越大侦镇,影響范圍越大,選取影響點(diǎn)越多织阅,曲線越平滑壳繁,但是相對來說更容易失真。rsacn的選取單位為維度荔棉。例如闹炉,上述例子所使用的插值分辨率為0.5deg,影響半徑為10润樱,其影響范圍為10個緯度渣触,20個格點(diǎn)*20個格點(diǎn),導(dǎo)致其插值結(jié)果失真比較嚴(yán)重壹若。

下圖為插值前后的對比嗅钻,圖1為未插值,直接用contour函數(shù)畫出的圖店展;圖二rscan為 (/10,5,3/) 养篓;圖三rscan為 (/2,1,0.5/) 。


?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末赂蕴,一起剝皮案震驚了整個濱河市柳弄,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌概说,老刑警劉巖碧注,帶你破解...
    沈念sama閱讀 221,430評論 6 515
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件嚣伐,死亡現(xiàn)場離奇詭異,居然都是意外死亡应闯,警方通過查閱死者的電腦和手機(jī)纤控,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 94,406評論 3 398
  • 文/潘曉璐 我一進(jìn)店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來碉纺,“玉大人船万,你說我怎么就攤上這事」翘铮” “怎么了耿导?”我有些...
    開封第一講書人閱讀 167,834評論 0 360
  • 文/不壞的土叔 我叫張陵,是天一觀的道長态贤。 經(jīng)常有香客問我舱呻,道長,這世上最難降的妖魔是什么悠汽? 我笑而不...
    開封第一講書人閱讀 59,543評論 1 296
  • 正文 為了忘掉前任箱吕,我火速辦了婚禮,結(jié)果婚禮上柿冲,老公的妹妹穿的比我還像新娘茬高。我一直安慰自己,他們只是感情好假抄,可當(dāng)我...
    茶點(diǎn)故事閱讀 68,547評論 6 397
  • 文/花漫 我一把揭開白布怎栽。 她就那樣靜靜地躺著,像睡著了一般宿饱。 火紅的嫁衣襯著肌膚如雪熏瞄。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 52,196評論 1 308
  • 那天谬以,我揣著相機(jī)與錄音强饮,去河邊找鬼。 笑死为黎,一個胖子當(dāng)著我的面吹牛邮丰,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播碍舍,決...
    沈念sama閱讀 40,776評論 3 421
  • 文/蒼蘭香墨 我猛地睜開眼柠座,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了片橡?” 一聲冷哼從身側(cè)響起妈经,我...
    開封第一講書人閱讀 39,671評論 0 276
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后吹泡,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體骤星,經(jīng)...
    沈念sama閱讀 46,221評論 1 320
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 38,303評論 3 340
  • 正文 我和宋清朗相戀三年爆哑,在試婚紗的時候發(fā)現(xiàn)自己被綠了洞难。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點(diǎn)故事閱讀 40,444評論 1 352
  • 序言:一個原本活蹦亂跳的男人離奇死亡揭朝,死狀恐怖队贱,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情潭袱,我是刑警寧澤柱嫌,帶...
    沈念sama閱讀 36,134評論 5 350
  • 正文 年R本政府宣布,位于F島的核電站屯换,受9級特大地震影響编丘,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜彤悔,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,810評論 3 333
  • 文/蒙蒙 一嘉抓、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧晕窑,春花似錦抑片、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,285評論 0 24
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽级遭。三九已至望拖,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間挫鸽,已是汗流浹背说敏。 一陣腳步聲響...
    開封第一講書人閱讀 33,399評論 1 272
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留丢郊,地道東北人盔沫。 一個月前我還...
    沈念sama閱讀 48,837評論 3 376
  • 正文 我出身青樓,卻偏偏與公主長得像枫匾,于是被迫代替她去往敵國和親架诞。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 45,455評論 2 359

推薦閱讀更多精彩內(nèi)容