ncl用linint2_Wrap函数插值后得到的都是缺省值

最后的秦羽 · 发表于 2017-4-16 18:38:40

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

goal = addfile("CN05.1_Tmax_1961_2014_daily_025x025.nc","r")
t = goal->tmax
lon = goal->lon ；lat: [14.75..55.25]
lat = goal->lat ；lon: [69.75..140.25]

intmax = addfile("NOSN_out04_198812_201011_TMAX_2M_daily_rg.nc","r")
tmax = intmax->TMAX_2M
rlon = intmax->rlon ；rlat: [-22.1..32.9]
rlat = intmax->rlat ；rlon: [-43.23..43.67000000000001]

TMAX_CCLM_IN = linint2(rlon,rlat,tmax,False,lon,lat,0)

我是想把CCLM的251*396的数据插值为163*283的0.25*0.25网格上

结果显示都是缺省值，是不是插值函数选用有问题？求大神们解答

lgz · 发表于 2017-4-16 19:00:35

原来的数据有缺测值

最后的秦羽 · 发表于 2017-4-16 19:14:48

lgz 发表于 2017-4-16 19:00
原来的数据有缺测值

应该不是这个原因吧，我打印下来全部都是缺省值

最后的秦羽 · 发表于 2017-4-16 21:28:54

问下大家原数据有缺省值，双线性插值会使缺省值扩散吗

爱吃爱吃 · 发表于 2017-4-27 18:29:56

要把经纬度具体数值标出来，例如 lon = goal->lon({0:360}),lat = goal->lat({-90:90}),理由我不是清楚

权威 · 发表于 2017-5-2 18:41:04

检测经纬度是否都是递增的

梁静茹93 · 发表于 2017-6-4 10:35:54

仔细看一下 nitmax 的lat/lon

liuxiaoyue · 发表于 2017-11-24 23:19:03

请问楼主解决了吗？

liuxiaoyue · 发表于 2017-11-24 23:21:43

数据中不能有缺省值，如果有一个格点有缺省值，那么所得的结果都是缺省值。

ransijingscu · 发表于 2017-11-25 12:33:06

经纬是不是反了

		自动登录	找回密码
密码			立即注册

[其他] ncl用linint2_Wrap函数插值后得到的都是缺省值

登录后查看更多精彩内容~

浏览过的版块