卫星数据变量坐标变量如何定义

tiantianyixin

所选卫星数据为AQUA/MODIS 1B 数据（250M），其中关于卫星数据相关维数信息为
dimensions:
      fakeDim0 = 2
      fakeDim1 = 15920
      fakeDim2 = 5416
      fakeDim6 = 3980
      fakeDim7 = 1354
      fakeDim8 = 3980
      fakeDim9 = 1354
; 变量空间分辨率为250M，对应的经纬度格点数为（15920,5416）
   variables:
      ushort EV_250_RefSB ( fakeDim0, fakeDim1, fakeDim2 )
         long_name : Earth View 250M Reflective Solar Bands Scaled Integers
         units : none
         valid_range : ( 0, 32767 )
         _FillValue : 65535
         band_names : 1,2
         radiance_scales : ( 0.02943603, 0.009811864 )
         radiance_offsets : ( 316.9722, 316.9722 )
         radiance_units : Watts/m^2/micrometer/steradian
         reflectance_scales : ( 5.872472e-05, 3.175246e-05 )
         reflectance_offsets : ( 316.9722, 316.9722 )
         reflectance_units : none
         corrected_counts_scales : ( 0.126194, 0.126194 )
         corrected_counts_offsets : ( 316.9722, 316.9722 )
         corrected_counts_units : counts
         hdf_name : EV_250_RefSB
; 坐标变量的经纬度格点数为（3980,1354）
    float Latitude ( fakeDim6, fakeDim7 )
         units : degrees
         valid_range : ( -90, 90 )
         _FillValue : -999
         hdf_name : Latitude
      float Longitude ( fakeDim8, fakeDim9 )
         units : degrees
         valid_range : ( -180, 180 )
         _FillValue : -999
         hdf_name : Longitude

如上所描述，变量空间分辨率为250M，对应的经纬度格点数为（15920,5416），坐标变量的经纬度格点数为（3980,1354），变量经纬度格点数是坐标格点数的4倍。我需要对变量的坐标进行定义，选取begin 后的脚本：
begin
  fi = addfile("./AQUA_2015_08_26_17_10_HY.MOD02QKM.hdf","r")
  data = new(dimsizes(fi->EV_250_RefSB),float,-99999)
;  由于Latitude,Lontitude格点数与变量data格点数不一致，坐标变量该如何设置呢？谢谢各位。
  lat2d = fi->Latitude
  lat2d@units ="degrees_north"
  lon2d = fi->Longitude
  lon2d@units ="degrees_east"
  data@lat2d = lat2d
  data@lon2d = lon2d

szqx9

szqx9

longlivehj

看产品说明

tiantianyixin

longlivehj 发表于 2015-8-27 13:27
看产品说明

好的，我在找找产品说明文件，谢谢

勇敢的悲伤

请问你这个问题解决了吗？