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

mofangbao

    这两天写了一个fortran版的双线性插值程序，原来找的一个写的太乱，还是决定自己写一个，自己用插值出来的数据和源数据画图对比后发现是吻合的（当然精度有变化），把源程序发出来大家可以讨论交流一下算法上有没有什么问题。

    可用于格点到格点（如降低和提高分辨率），格点到站点的插值
    插值的原理很好懂，就是把线性插值分为横向和纵向来完成，借用一张示意图简单说明：
   


P点是目标点，四周的Q就是用于插值的源数据，插值算法的一个重要内容就是找出目标点的四周点，另一个重要内容就是找出点的权重，由于时间比较急，我在计算权重的时候是分为横向和纵向计算的，也就是先计找出两条竖线的序号，然后计算Q11和Q21横向的权重（根据到R1的距离占Q11Q21的比重来算，比重大权重小），Q12和Q22在横向上的权重是相同的，这样就可以插值出R2和R1这两个点的值，随后计算R2和R1在纵向的权重，这样就能算出P的值了。
好吧，源程序在此，收一个贡献，希望大家能把自己的东西多多上传分享，发帖分享是可以获得贡献的。
（感谢M版王老师的鼓励和指导）
bi_interp.f90 (3.95 KB, 下载次数: 1206, 售价: 1 贡献)

2012-2-27 16:34 上传

点击文件名下载附件

兰溪之水

沙发~

mofangbao

@残月哪里出问题了？贡献只扣一次的，如果扣了，你可以重复下载，不会重复扣除贡献的，但是金钱会...，不要用迅雷下载啊，嘿嘿

残月

mofangbao 发表于 2012-2-27 16:56

                               
登录/注册后可看大图

@残月哪里出问题了？贡献只扣一次的，如果扣了，你可以重复下载，不会重复扣除贡献的，但是金钱会...，不要 ...

应该是我网速的问题，已下载成功

topmad

记得深深原来貌似写过一个，不知道两个一不一样...下来学习一下

言深深

topmad 发表于 2012-2-27 19:19

                               
登录/注册后可看大图

记得深深原来貌似写过一个，不知道两个一不一样...下来学习一下

我也记得貌似我写过一个的想不起来了

topmad

言深深 发表于 2012-2-27 19:30

                               
登录/注册后可看大图

我也记得貌似我写过一个的想不起来了

咱们还讨论了一下两个点之间数据变化的斜率....

言深深

topmad 发表于 2012-2-27 19:31

                               
登录/注册后可看大图

咱们还讨论了一下两个点之间数据变化的斜率....

对呀，难道只是讨论了？没有写出来么？图看起来蛮熟悉的，公式也推过......5555555555555年代久远了

topmad

言深深 发表于 2012-2-27 19:32

                               
登录/注册后可看大图

对呀，难道只是讨论了？没有写出来么？图看起来蛮熟悉的，公式也推过......5555555555555年代久 ...

@言深深

看这里深深，自己写的东西啊，我可是找到了....累死了

言深深

topmad 发表于 2012-2-27 19:38

                               
登录/注册后可看大图

@言深深

看这里深深，自己写的东西啊，我可是找到了....累死了

哈哈，辛苦辛苦