临近点插值问题(已解决)

CDH

各位，我在做低能见度预报的相关工作，我现在遇到的问题是这样的，我用0.5625*0.5625°的T639数据做未来预报，我针对研究区域的每个站点确定了几个指标，我现在想用距离格点最近的站点值代替该格点值，而不是用插值的办法将站点插值到格点，因为插值后的结果降低了我选择的这几个指标的标准。
我看了清风的帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=1780中关于格点插值到站点的函数gr2stn(grid_expr, stn_expr, <-n>)后面添加上-n，就可以不做插值而是用最相邻的格点值来替代，查了一下oacres这个函数好像没有这个功能。


让各位帮忙看的就是，像我这样不做插值，而是用距离最近的站点值作为某格点值的想法能通过grads某个函数实现吗？findstn?
另外就是假如这个想法实现起来比较麻烦，而需要做插值的话，那oacres这个函数的半径序列我该如何确定呢？


关于oacres函数半径序列该如何确定的问题，我自己也在论坛里搜索学习一下，现在主要就是第一个问号的问题，谢谢各位了。

CDH

自己顶一下，谢谢各位

做个霸气的木头

帮顶喔，@mofangbao来看看

CDH

做个霸气的木头 发表于 2015-8-25 23:18
帮顶喔，@mofangbao来看看

谢谢支持，如果这个想法实现不了就只能做插值了

传说中的谁

呃。。。举个栗子，站点经纬度为a,b，那么在grads里set lon a ,set lat b，然后直接d var，出来就是最接近站点的某个格点值

CDH

传说中的谁 发表于 2015-8-26 09:01
呃。。。举个栗子，站点经纬度为a,b，那么在grads里set lon a ,set lat b，然后直接d var，出来就是最接近 ...

哎呀，就是的，您这个思路很好，我整理一下，看能不能行。谢谢传说。

CDH

传说中的谁 发表于 2015-8-26 09:01
呃。。。举个栗子，站点经纬度为a,b，那么在grads里set lon a ,set lat b，然后直接d var，出来就是最接近 ...

我现在把数据准备一下，不知道oacres这个函数能不能像gr2stn(grid_expr, stn_expr, <-n>)这个命令这样，后面加一个-n就可以不插值，用最近的值代替了

mofangbao

CDH 发表于 2015-8-26 10:00
我现在把数据准备一下，不知道oacres这个函数能不能像gr2stn(grid_expr, stn_expr, )这个命令这样，后面 ...

oacres是cressman插值，你要的是邻近插值，怎么可能是一个参数能解决的呢。笨笨说的没错，你直接设置好经纬度d一下就行了。批量做的话只要准备好经纬度信息表，然后读取，输出就可以了

CDH

mofangbao 发表于 2015-8-26 14:08
oacres是cressman插值，你要的是邻近插值，怎么可能是一个参数能解决的呢。笨笨说的没错，你直接设置好经 ...

谢谢清风的回复，我按照这个思路做一下，然后看一下结果。

mofangbao

CDH 发表于 2015-8-26 14:12
谢谢清风的回复，我按照这个思路做一下，然后看一下结果。

嗯 弄出来了可以分享下  组织站点经纬的方式可以参考  http://bbs.06climate.com/forum.php?mod=viewthread&tid=18335  这个帖子