登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
大家好,弱弱请教natgrid是自然邻近插值函数吧,我想把格点插值到站点上,但能力有限,搞不出来QAQ,请大家帮忙看看是哪里出问题了呢?非常谢谢!
f=addfiles("F:/test_out/CRA40_TEM_2020.nc","r") Tem1=f[:]->Tair_f_inst(:,::-1,0:719) ;这里我看说明要求x,y一样长,所以经纬度截取范围都是720 time=Tem1&time lon=Tem1&lon lat= Tem1&lat Tem1@_FillValue=-1e+10 ;站点的经纬度 xy=asciiread("F:/station/91lat_lon.txt",(/91,2/),"float") LAT=xy(:,0) LON=xy(:,1) fo= natgrid_Wrap(lon, lat, Tem1,LON, LAT) ;目的是格点插值到站点上,ncl计算耗时将近5minQAQ,结果明显有问题(见图),之前用双线性插值,得出的结果是一个时间对应一个站点值,但natgrid不是这样的,所以我在想自己可能理解错这个函数的用法了,请大家指教QAQ,ncl该如何用邻近插值办法,将格点插值到站点上,感谢感谢感谢QAQ! printVarSummary(fo)
sleep(100)
|