bicubic spline

http://en.wikipedia.org/wiki/Bicubic_interpolation

blockmedian は実装して，
GMTのsurfaceコマンド代わりにbicubic interpolation を実装してみたんだけど．．．
どうもグリッド上の値は保存されて，
近傍16点だと0は0のまま保存されて，周りが高くなっちゃうんだけど，，，
うーーん，どこか間違えている気がする．

Wikipediaにはα＝-0.5って書いてあるけど，他のサイトでは-1.0を使っていた．