ncl 做双线性插值的问题

橙子鲜

求助，请问一下，想要对数据做双线性插值，原本数据为txt，东亚区域，因而读入之后对其重新赋lat，lon.之后对数值做插值，为什么得出来的结果大部分为缺测值呢
脚本
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
;************************************************
; open file and read in variable
;***********************************************
f2 =addfile("/cygdrive/e/LMDZ/data/LMDZ_day_t2m1961.nc","r")
lat=f2->lat
lon=f2->lon

f1=new((/80,77/),float)
filelat="e:/LMDZ/data/1.txt"
f1 = asciiread(filelat,(/80,77/),"float")
  ;print(f1)
f1!0="lat"
  f1&lat=lat
f1!1="lon"
  f1&lon=lon;;;;;;;;;;;;;;;;;;;;;;;;将f2的经纬度赋值给f1
; printVarSummary(f1)

    ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;做双线性插值
  LON = fspan(86.15,134.15,49);;;;0.6
  LAT = fspan(53.81723,5.81723,81);;;0.6
;print(LAT)

  fo = linint2_Wrap (lon,lat,f3, True, LON,LAT,0)
   print(fo)
结果如图

niuda

双线性插值的边界是没有值的

橙子鲜

niuda 发表于 2015-6-28 21:51
双线性插值的边界是没有值的

请问这个意思是？我哪里编写有问题呢？

subtropical

楼主解决了吗，我也是插值后全是缺测值

lmyyky

subtropical 发表于 2017-3-15 10:27
楼主解决了吗，我也是插值后全是缺测值

请问解决了吗？我现在插值后也全部为缺测，用的linint2_points_Wrap这个函数

悠然一梦123

楼主，想问下，你这个插值的数据分辨率是0.6*0.6度的吗， LON = fspan(86.15,134.15,49)这里的49应该是81吧，求解答，目前我也在做插值这方面的，不太懂，求楼主解答下

机智的大气宝宝

你试一下把维数lat的方向反过来，我之前lat由北到南出来也是全是缺测值，反向就好了

Nalanalanaa

机智的大气宝宝 发表于 2018-1-30 13:17
你试一下把维数lat的方向反过来，我之前lat由北到南出来也是全是缺测值，反向就好了

你好 我也遇到这个问题 请问怎么反过来呢

liweier123

Nalanalanaa 发表于 2020-3-9 13:45
你好 我也遇到这个问题 请问怎么反过来呢

lat=f2->lat(::-1)

dengruzhang

有可能输入或者输出的纬度数组是递减的，读取的时候设置成：：-1