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

weather3699

感谢楼主分享！

shenchan

有没有人用过，程序是否有问题？

sparrow

感谢分享，学习中

Luo丹王道

亲爱的楼主，您这个程序应该就只是子程序吧，就缺个读入格点数据、站点数据和输出插值后的数据了是吧？

走路睡觉的猫WD

太好了，正好需要

mofangbao

Luo丹王道 发表于 2015-7-9 23:24
亲爱的楼主，您这个程序应该就只是子程序吧，就缺个读入格点数据、站点数据和输出插值后的数据了是吧？

那些都是外围工作，每个人读取的数据都不同，输出的要求都不同，当然需要自己实现，如果连如何调用子程序都不会或者自己都不能写数据读取和输出部分的话。。。只能说该先学习基础东西。

范海燕

很有用
谢谢分享，好好学习了

哇咔pig

楼主，你好！谢谢你的程序，我想请问下，在这个程序里面，原始数据的起始经纬度都是从0 开始的，才可以直接用，是么？如果我的数据不是从0开始的，需要在哪里进行修改？

rookie

您好，有个问题向问下您。原始网格为80*80大小 ，插值为200*200，由于80*80为lcc投影，插值后再图片上有数值的区域为扇形，也就是原来80*80网格的范围。
参考GrADS说明文档，插值采用双线性插值：
Bi-linear interpolation is used within this grid box to get down to the         position 31.24,67.88. The interpolation is linear within the i,j grid.


那么，如果我想获取 （106.23,38.44）点的值，GrADS是先在80*80网格中找到离（106.23,38.44）最近的周围四个点，然后用这四个点采用双线性插值获取（106.23,38.44）点的值。
不知道我的理解对不对？
还有插值时位置用网格点的行列号还是经纬度值？


PDEF  80  80  LCC  36.000  105.000  40.500  40.500  46.00000  26.00000  105.00000  30000.000  30000.000
XDEF  200  LINEAR  88.55157  0.12162162
YDEF  200  LINEAR  23.49021  0.12162162

mofangbao

rookie 发表于 2015-9-10 18:05
您好，有个问题向问下您。原始网格为80*80大小 ，插值为200*200，由于80*80为lcc投影，插值后再图片上有数 ...

GrADSZ中如何实现的你要去看源码啦