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/) 。