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

CDH · 发表于 2015-8-26 14:16:24

mofangbao 发表于 2015-8-26 14:13
嗯弄出来了可以分享下组织站点经纬的方式可以参考 http://bbs.06climate.com/forum.php?mod=viewthre ...

好的，结果出来了再来汇报一下。谢谢清风。

CDH · 发表于 2015-8-26 16:24:03

mofangbao 发表于 2015-8-26 14:13
嗯弄出来了可以分享下组织站点经纬的方式可以参考 http://bbs.06climate.com/forum.php?mod=viewthre ...

清风版主，我测试了一下，我打开一个格点资料，然后我想获得某个站点的资料，然后set lon 站点经度，set lat 站点纬度，然后d var的时候可以给出距离站点最近的这个格点的值，但是反过来我打开的是站点的数据，想获得临近格点的值，就提示我0 reports。不能获得这样的数据，而我现在恰恰是想用距离格点最近的站点数据代替格点数据。
不知道我表述清楚我的问题了没有，您还能再给些建议和指导吗？谢谢清风。

CDH · 发表于 2015-8-26 16:33:33

mofangbao 发表于 2015-8-26 14:13
嗯弄出来了可以分享下组织站点经纬的方式可以参考 http://bbs.06climate.com/forum.php?mod=viewthre ...

如果grads的函数实现不了，估计就只能用Fortran程序，来判断哪个站距离格点最近，然后用这个值来作为格点值了。

宁小才 · 发表于 2015-8-26 16:43:31

如果grads的函数实现不了，估计就只能用Fortran程序@CDH

CDH · 发表于 2015-8-26 22:08:47

这个问题最后通过Fortran程序解决了，核心思想就是利用经纬度，计算格点和站点之间的距离，dis(j)=sqrt(abs(lati(i)-lat(j))**2+abs(long(i)-lon(j))**2)，找到距离格点最近的站点，然后把站点值赋值给格点。

袜子女巫 · 发表于 2018-1-11 14:14:38

版主，程序还在吗，能发我一份吗

		自动登录	找回密码
密码			立即注册

[脚本编辑] 临近点插值问题(已解决)