给大家分享一个MATLAB版的cressman插值函数程序及示例

lgz

谢谢楼主分享~

jiajian0512

请问楼主，三次调用后分别得到的结果再进行平均吗？不太懂原理，谢谢啦

kc121

跟同事一起用surfer画图比较了一下。发现有时一致，有时差异还挺大的。且其中一个站所在区域明明是在大值区，这个程序画出来是小值区。楼主讲下调用这个插值程序的那几个要素嘛。还正确使用您的程序。谢谢！

xiaoxiangzi94

好东西，感谢楼主分享！

zhaohaoran999

学习学习！！！

mapleleaves99

感谢分享！！

白兰

kc121 发表于 2017-6-21 17:18
跟同事一起用surfer画图比较了一下。发现有时一致，有时差异还挺大的。且其中一个站所在区域明明是在大值区 ...

是这样的，每个格点的最终插值数值由周围半径R内的距离平方反比约束，这个请查cressman插值的原理。也可能是因为这个原因，会导致个别站点观测明明是大值，最后可能被平滑成了中等值...（极端异常站点可能就很容易被平滑掉）。。。同时，受使用者观测数据站点空间分辨率的影响，需要多次调用插值函数，选择适当的影响半径，（如调用三次，分别为R1，R2，R3），这也是由插值方法的原理决定的。之所以给三个不同的影响半径，是为了提供两次背景猜值场，最终起关键作用的是最后的小半径R3，若R3取的很小，比如小于了你观测站点空间分辨率，那画出来的肯定就是很多个圆圈，类似雷达图那种。  如果我的解释大家没有看明白的话，可以使用grads的cressman插值函数自己画图测试一下，ncl也有相应的函数，我只是移植 到了matlab上而已。

白兰

jiajian0512 发表于 2017-6-15 14:06
请问楼主，三次调用后分别得到的结果再进行平均吗？不太懂原理，谢谢啦

是这样的，每个格点的最终插值数值由周围半径R内的距离平方反比约束，这个请查cressman插值的原理。也可能是因为这个原因，会导致个别站点观测明明是大值，最后可能被平滑成了中等值...（极端异常站点可能就很容易被平滑掉）。。。同时，受使用者观测数据站点空间分辨率的影响，需要多次调用插值函数，选择适当的影响半径，（如调用三次，分别为R1，R2，R3），这也是由插值方法的原理决定的。之所以给三个不同的影响半径，是为了提供两次背景猜值场，最终起关键作用的是最后的小半径R3，若R3取的很小，比如小于了你观测站点空间分辨率，那画出来的肯定就是很多个圆圈，类似雷达图那种。  如果我的解释大家没有看明白的话，可以使用grads的cressman插值函数自己画图测试一下，ncl也有相应的函数，我只是移植 到了matlab上而已。

白兰

jiajian0512 发表于 2017-6-15 14:06
请问楼主，三次调用后分别得到的结果再进行平均吗？不太懂原理，谢谢啦

多次调用是为了使第一猜值场变得更加接近实际值，这是由插值方法的原理决定的，即首先每个格点需要有个背景值，默认都是0。那么如果先使用R1的半径对全是0的背景场做一次插值，那么生成的插值场就更接近实际观测的。基于这个得到的插值场，将其作为第二次调用函数的第一猜值场，就会使得插值场更加接近实际。同理，第三次半径R3调用之前R2的插值场，继而生成最终的插值场。因此，并非要“三次调用后分别得到的结果再进行平均”，而是递进关系，后面的调用要基于前面生成的插值场作为背景。

longlong4123

谢谢分享{:5_213:}