接着上一篇文章，试图画出速度模型时，中间会有用到调用C库函数--gsl库；记一下怎么使用的：

和C语言中类似，但要使用到条件编译：

#ifdef  __cplusplus
extern "C" {
#endif
#include <gsl/gsl_math.h>
#include <gsl/gsl_interp2d.h>
#include <gsl/gsl_spline2d.h>
#ifdef  __cplusplus
}
#endif

将这部分包含到相应的文件中，一般是.h中；

关于GSL手册，可以参考官网：

GNU Scientific Library — GSL 2.7 documentationhttps://www.gnu.org/software/gsl/doc/html/index.html

然后对于插值部分，需要确定插值类型，插值工作区，加速区，插值设置，初始化： 

const gsl_interp2d_type *T = gsl_interp2d_bilinear;//插值类型（方法）spline = gsl_spline2d_alloc(T, s->nx, s->ny);//gsl_spline2d，插值工作区xacc = gsl_interp_accel_alloc();//gsl_interp_acc，加速工作区yacc = gsl_interp_accel_alloc();for(int i=0;i<s->nx;i++){for(int j=0;j<s->ny;j++){gsl_spline2d_set(spline, s->za, i, j, s->array.at(i*(s->ny)+j) );}}    
gsl_spline2d_init(spline, s->xa, s->ya, s->za, s->nx, s->ny);

 以及调用

QwtRasterData类的value函数

返回某坐标下的值。 

double vv = gsl_spline2d_eval(spline, x, y, xacc, yacc);       
return vv;