自编双线性插值fortran程序（Bilinear Interpolation）

小虹 · 发表于 2012-9-8 11:01:20

谢谢

小虹 · 发表于 2012-9-9 21:08:23

以前自己编过，今天花了点时间看楼主的程序，觉得学到了很多，但是我有一个问题，在subroutine getNewVal(oldVal,newVal,newPoint,isFirst) 子程序的开始插值过程中，注释‘先做纬向插值’后面的if语句，
if(pos(i,1)==undef .or. pos(i,2)==undef .or. pos(i,3)==undef .or. pos(i,4)==undef)then
如果我没有理解错的话，pos(i,1)，pos(i,2)，pos(i,3)，pos(i,4)是描述所插点位于原格点的位置的，那么何来他等于undef一说，是不是楼主写错了，应该是
if(oldVal(pos(i,1),pos(i,3))==undef .or. oldVal(pos(i,2),pos(i,3))==undef .or.oldVal(pos(i,1),pos(i,4))==undef.or.oldVal(pos(i,2),pos(i,4))==undef)then
这个样子呢？

mofangbao · 发表于 2012-9-9 21:47:54

小虹发表于 2012-9-9 21:08

登录/注册后可看大图

以前自己编过，今天花了点时间看楼主的程序，觉得学到了很多，但是我有一个问题，在subroutine getNewVal(o ...

应该是你说的那样，写的有点久了，不太记得了，你可以用缺测值试试看，谢谢指出~

langyie · 发表于 2012-9-20 15:52:46

支持一下工作

gqp1221 · 发表于 2012-9-28 17:08:02

非常好的东西啊

znk1981 · 发表于 2012-10-25 17:23:44

非常感谢，刚好要用到

吕丽 · 发表于 2012-10-29 12:54:08

版主辛苦~

{:eb513:}

stopfire · 发表于 2012-10-29 14:34:34

顶一个~~~~~~~~~~~~

deecoct · 发表于 2012-11-1 22:37:32

mofangbao 发表于 2012-2-28 09:02

登录/注册后可看大图

我做的就是全球的...而且源数据固定为全球格点，做小范围的话需要稍作修改

清风君，问个问题哈
就是双线性插值的条件需不需要插值点周围的四个点在一个矩形上的？
因为实际情况是经纬度，还不是矩形呢。。。。。
我看了你的程序后唯一纠结的地方就是这里。。。。求解答

PS
我是想用etopo5.nc资料插出更小精度的地形高度值，不知道你的程序可不可以用。

mofangbao · 发表于 2012-11-1 22:52:04

deecoct 发表于 2012-11-1 22:37

登录/注册后可看大图

清风君，问个问题哈
就是双线性插值的条件需不需要插值点周围的四个点在一个矩形上的？
因为实际情况是 ...

我还没仔细考虑过，不过我直接用来插值欧洲中心的温度资料，没考虑你说的东西，画出来的图和插值之前的原图几乎是完全吻合的~

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

[源代码] 自编双线性插值fortran程序（Bilinear Interpolation）