爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 11853|回复: 1

[作图] 请教各位大神,关于ncl里natgrid插值函数的使用,求进~~~

[复制链接]
发表于 2021-7-7 09:53:54 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
大家好,弱弱请教natgrid是自然邻近插值函数吧,我想把格点插值到站点上,但能力有限,搞不出来QAQ,请大家帮忙看看是哪里出问题了呢?非常谢谢!

f=addfiles("F:/test_out/CRA40_TEM_2020.nc","r")
Tem1=f[:]->Tair_f_inst(:,::-1,0:719)  ;这里我看说明要求x,y一样长,所以经纬度截取范围都是720
time=Tem1&time
  lon=Tem1&lon
  lat= Tem1&lat
Tem1@_FillValue=-1e+10
;站点的经纬度
  xy=asciiread("F:/station/91lat_lon.txt",(/91,2/),"float")
LAT=xy(:,0)
  LON=xy(:,1)
  fo= natgrid_Wrap(lon, lat, Tem1,LON, LAT) ;目的是格点插值到站点上,ncl计算耗时将近5minQAQ,结果明显有问题(见图),之前用双线性插值,得出的结果是一个时间对应一个站点值,但natgrid不是这样的,所以我在想自己可能理解错这个函数的用法了,请大家指教QAQ,ncl该如何用邻近插值办法,将格点插值到站点上,感谢感谢感谢QAQ!
printVarSummary(fo)

sleep(100)





气象家园.png
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2021-7-8 11:01:35 | 显示全部楼层
泪流满面,自己给自己顶顶
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表