NCL插值

cpp803520

程序如下：
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
a = addfile("AVHRR-Land_v004_AVH15C1_NOAA-18_20080101.nc","r")
lai = short2flt(a->LAI)
LON=fspan(73,135,1008)
LAT=fspan(14,53,640)
lai1=area_hi2lores_Wrap(lai&longitude,lai&latitude,lai,True,1.0,LON,LAT,0)
b = addfile("2008-01-01.nc","c")
b->LAI=lai1
end
报错信息如下： warning:Argument 7 of the current function or procedure was coerced to the appropriate type and thus will not change if the function or procedure modifies its value
        file: NclNewFile.c, line: 6839
        Could not get varnode->name: <time>
        file: NclNewFile.c, line: 6839
        Could not get varnode->name: <LATITUDE>
        file: NclNewFile.c, line: 6839
        Could not get varnode->name: <LONGITUDE>
这是只有6839这一行的数据格式不符合还是所有的数据格式都不符合啊？还有这个数据我是下载后将nan设置为0后就插值为什么会说数据格式不符合呢？求大神们解答，谢谢啦~~~

低调高姿态

area_hi2lores_Wrap,师姐函数用错了，你做的LAI，应该用area_hi2lores。

cpp803520

低调高姿态 发表于 2016-7-9 09:25
area_hi2lores_Wrap,师姐函数用错了，你做的LAI，应该用area_hi2lores。

可是改成lai1=area_hi2lores(lai&longitude,lai&latitude,lai,True,1.0,LON,LAT,0)
出错信息是Segmentation fault (core dumped)

低调高姿态

cpp803520 发表于 2016-7-9 09:47
可是改成lai1=area_hi2lores(lai&longitude,lai&latitude,lai,True,1.0,LON,LAT,0)
出错信息是Segmentat ...

师姐，现在解决了吗？

cpp803520

低调高姿态 发表于 2016-7-9 23:26
师姐，现在解决了吗？

没有呢，总是说是块错误，我猜测是没有坐标系的原因，在考虑要不要给它加个坐标系，不过也有可能是元数据有问题~

低调高姿态

cpp803520 发表于 2016-7-11 09:02
没有呢，总是说是块错误，我猜测是没有坐标系的原因，在考虑要不要给它加个坐标系，不过也有可能是元数据 ...

元数据可以用软件看看。昨天开会了，不好意思，刚看到。。