linint2插值后经纬度范围出错【已解决】

Hurricane_Hu

刚刚开始研究双线性插值，打算用NCL做处理，但是不是用NCL出图。
信息如下
原数据：  格点数据，维度： [time | 1872] x [lat | 96] x [lon | 96]；
打算插值到纬度为（-90,90），精度为（0，360）的1*1分辨率的均匀格点；
linint2 语句如下：
mlon=360
mlat=180

klon=ispan(0,mlon-1,1)
klat=ispan(-90,89,1）
tas_new=linint2(tas&lon,tas&lat,tas,True,klon,klat,0)

插值之后发现 tas_new 的纬度是（0,180），并不是（-90,90），试了很多次都是这样，不知道是不是NCL自身插值函数的问题，后面用UV-CDAT画出的图就变成了只有半个地球。。。（用NCL画无地图的等值线图的纬度为0-180的时候没有问题，所以应该不是插值错误）
请问这个问题应该如何解决，把纬度调成（-90,90），谢谢！

解决方法：使用 linint_Wrap 函数（附上解决之后的图)

金玉姻缘

我用linint2这个函数差值，参照的是http://www.ncl.ucar.edu/Applications/regrid.shtml，，
这里面的regrid_1.ncl即这个，http://www.ncl.ucar.edu/Applications/Scripts/regrid_1.ncl，
我模仿者变换了一下输入文件，改为air_mon_mean.nc,具体情况是air（time,lat,lon）,lat=360,lon=720,是0.5*0.5度的数据，我自己设定的插值网格为
newlat=fspan(50,35,16),newlon=fspan(75,95,21),是1*1度的均匀网格，然后newair=linint2(lon,lat,air,True,newlon,newlat,0)差值完成后，发现经纬度显示是正确的，newlon=75,76,77,78,79,......94,95,newlat=50,49,48,47,.....36,35,和我预期的一致，但是newair这个数组里面没数据，全是空的，显示为newair= _, _, _, _,
   等等，这应该是没插值成功吧，想了好久，想不通，求解释，具体的ncl脚本就是参考上面的链接regrdi_1.ncl的，只改了以上提到的几处，

longlivehj

用linint2_Wrap

Hurricane_Hu

longlivehj 发表于 2014-8-17 16:02
用linint2_Wrap

成功了！！太谢谢了！！！！
研究研究这个函数！！感谢感谢！！

zht6707

太谢谢了！！

柚柚子

我也遇到了这个问题，请问你后来是怎么解决的？

柚柚子

金玉姻缘 发表于 2014-11-29 21:10
我用linint2这个函数差值，参照的是http://www.ncl.ucar.edu/Applications/regrid.shtml，，
这里面的regr ...

我也遇到了这个问题，请问你后来是怎么解决的？

clannand

你好我的文件是nc格式，但没有经纬度信息，文件名称后缀不是.nc，可以插值吗