ncl利用linint2_points进行插值结果全为缺省值

tiantianyixin

各位大神，我利用linint2_points将全球的格点数据插值到一些站点上，插值后的结果全为缺省值，现将脚本等信息黏贴如下，还望各位指点：

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

  all_files = systemfunc("ls /home/dih/file/gfs.*.grb2")
  a  = addfiles(all_files,"r")
  ListSetType(a, "join")

  latori = a[0]->lat_0   
  lon = a[0]->lon_0
; 将latori变量按递增的方式赋值给变量lat
; 将原始纬度数组按从小到大的顺序重新排序
  lat = new(721, float)    ;创建一个容量为121，元素类型为float的数组
do i = 0,720
   lat(i) = latori(720-i)
end do

  u10  = a[:]->UGRD_P0_L103_GLL0(:,:,::-1,:)
  v10  = a[:]->VGRD_P0_L103_GLL0(:,:,::-1,:)
printVarSummary(u10)
; List of stations
  data=asciiread("station.txt",(/119,3/),"float")
  stat=data(:,0)
  lon1=data(:,1)
  lat1=data(:,2)
  num_st = dimsizes(stat)
; 进行双线性插值计算
  n_u = linint2_points(lon,lat,u10(:,0,:,:),True,lon1,lat1,0)
  n_v = linint2_points(lon,lat,v10(:,0,:,:),True,lon1,lat1,0)

其中u10的变量信息为：
Variable: u10
Type: float
Total Size: 37376640 bytes
            9344160 values
Number of Dimensions: 4
Dimensions and sizes: [ncl_join | 3] x [lv_HTGL7 | 3] x [lat_0 | 721] x [lon_0 | 1440]
Coordinates:
            lv_HTGL7: [10..100]
            lat_0: [-90..90]
            lon_0: [ 0..359.75]
Number Of Attributes: 12
  center : US National Weather Service - NCEP (WMC)
  production_status : Operational products
  long_name : U-component of wind
  units : m s-1
  _FillValue : 1e+20
  grid_type : Latitude/longitude
  parameter_discipline_and_category : Meteorological products, Momentum
  parameter_template_discipline_category_number : ( 0, 0, 2, 2 )
  level_type : Specified height level above ground (m)
  forecast_time : 0
  forecast_time_units : hours
  initial_time : 04/01/2015 (12:00)

chen--

怎么解决原插值数据有缺测值这一问题啊？

哈库拉玛塔塔000

请问楼主最后是怎么解决的呢

权威

检测经纬度是不是递增的

liuxiaoyue

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

Tobyright

liuxiaoyue 发表于 2017-11-24 23:21
数据中不能有缺省值，如果有一个格点有缺省值，那么所得的结果都是缺省值。

那如果数据中有缺测值的情况该怎么解决呢？

这一天

楼主，请问有解决的办法吗