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