【NCL】反距离权重插值问题

s543098687

写了个用dsgrid2函数做插值的脚本：
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"

begin

filepath = "/data2/public1/shengjw/work/shice.txt"
lines = asciiread(filepath,-1,"string")
lon = stringtofloat(str_get_field(lines(0:),2," "))
lat = stringtofloat(str_get_field(lines(0:),1," "))
R  = stringtofloat(str_get_field(lines(0:),3,"  "))
;print(R)
;print(lat)
;print(lon)
olon = new((/15/),float)
olat = new((/9/),float)

do i=0,14
   olon(i) = 89+i
end do
do j=0,8
   olat(j) = 31+j
end do

dssetp("exp",0.5)
data1 = dsgrid2(lon,lat,R,olon,olat)
print(data1)
;system("/bin/rm -f "+"out_ds.txt")

;write_table("./out_ds.txt","a",[/"LON","LAT","R"/],"%s %s  %s")
;do x=0,1400
;   do y=0,800
;       write_table("./out_ds.txt","a",[/olon(x),olat(y),data1(x,y)/],"%0.2f    %0.2f   %f")
;   end do
;end do

end

shice.txt如下格式：
38.25    90.85    44
38.74    93.34    35.4
38.81    98.42    315.3
38.43    99.59    472.8
38.18    100.24    433.3
36.8    93.68    24.8
37.85    95.35    65.1
37.37    97.38    165.7
37.3    99.02    428.9
37.33    100.14    423.7
37.38    101.61    504.2
...

为什么插值完成后的同格点预测值和实测值会差别巨大呢？是参数定义的原因么？

s543098687

这个是1°的精度，我插到0.01°后所有格点预测值都为380左右，但实际值有很多站点值都是两位数或者800多，不知道是哪里出了问题？

ymy_00000

我想请问下在数据文件中的R是如何生成的

yeah...

ymy_00000 发表于 2018-12-12 10:43
我想请问下在数据文件中的R是如何生成的

请问楼主知道了吗？我也想问