怎么将1*1的格点数据，转化为2.5*2.5的格点数据

推开世界的门

大佬们，插值问题已解决感谢各位大佬

bc1012

例如数据是分辨率1*1的东经20-120，北纬10-60 是100*50的场，如要插值成分辨率2.5*2.5
olon1 = new(40,"float");
  olat1 = new(20,"float");
  fbl1 = new((/20,40/),"float")
  do i=0,39
     olon1(i) =20+i*2.5
  end do
  do l=0,19
     olat1(l) =15+2.5*l
  end do
lon1!0="lon"
olon1@long_name="lon"
olon1@units="degrees-east"
olon1&lon=olon1
olat1!0="lat"
olat1@long_name="lat"
olat1@units="degrees_north"
olat1&lat=olat1

lon = new(40,"float");
lat = new(20,"float");
lon!0="lon"
lon@long_name="lon"
lon@units="degrees-east"
lat!0="lat"
lat@long_name="lat"
lat@units="degrees_north"

fb5@_fillvalue=9.9621e+36
rscan=(/1/)
T1=obj_anal_ic_Wrap(zlon5, zlat5, T, olon1, olat1, rscan, False);zlon,zlat为原数据的经纬
好像有三种插值法  哪种插值方法应该看具体问题的   
  最后可以给T1定义下，给T1D的经纬上单位和赋值

Tobyright

请教一下楼主，最后问题是怎么解决的呢，能分享一下相关代码吗

推开世界的门

Tobyright 发表于 2021-7-1 10:26
请教一下楼主，最后问题是怎么解决的呢，能分享一下相关代码吗

begin

  modeltime="222.nc"
  path="F:/Piggy_NCL/Cygwin64/home/15741/"

f=addfile(path+modeltime,"r")

seat=f->a(:,{-90:90},:)

newlat = fspan(-90.,90,361)              ;格点数=361/分辨率0.5
  newlon = fspan(0.,359.5,720)             ;格点数=720/分辨率0.5
  newlat@units = "degrees_north"
  newlon@units = "degrees_east"
  newp=new((/dimsizes(seat(:,0,0)),dimsizes(newlat),dimsizes(newlon)/), float)
  newp!0="time"
  newp!1="lat"
  newp!2="lon"
  newp&time=seat&time
  newp&lat=newlat
  newp&lon=newlon
  newp@_FillValue = 1e+20
  printVarSummary(newlat)
  printVarSummary(newlon)
  newp= linint2(seat&lon,seat&lat,seat,False,newlon,newlat,0);linint2_Wrap lat,lon会变成LAT，LON
  
  
  
  ;print(seat)

  setfileoption("nc","Format","LargeFile") ;生成文件的属性(>2Gb)
  printVarSummary(newp)

  pathout="F:/Piggy_NCL/Cygwin64/home/15741/"
  outfilename = "111.nc"
  ncdf = addfile(pathout+outfilename,"c") ;输出插值后文件
  ncdf->seat = newp

end

Tobyright

推开世界的门 发表于 2021-7-1 11:07
begin

  modeltime="222.nc"

感谢楼主123456