OpenGL数据可视化（2）绘制二维地震剖面

上篇文章实现二维 RGB 数组转为纹理，本文展示如何将二维地震剖面转为纹理。由于代码比较多，此处只讲实现思路，本系列完结后会将代码开源。

实现思路

1. 读取二维地震剖面数据 seis，尺寸为 imgh*imgw；
2. 将地震数据的范围线性缩放至 [0, 255]，即
3. 构造二维 colormap_promax 数组，尺寸为 256*3；
4. 以 seis 就近取整的值为索引，利用 colormap 获得纹理数组 data，即

   unsigned char* data = (unsigned char*)calloc(imgh * imgw * 3, sizeof(unsigned char));
   int now_id, now_v;
   for (j = 0; j < imgh; j++) {
       for (i = 0; i < imgw; i++) {
           now_id = (j * imgw + i) * 3;
           now_v = round(data_raw[i][j]);
           if (now_v < 0)
               now_v = 0;
           data[now_id] = colormap_promax[now_v][0];
           data[now_id + 1] = colormap_promax[now_v][1];
           data[now_id + 2] = colormap_promax[now_v][2];
       }
   }

结果展示

注意

数据的尺寸不能超过 GL_MAX_TEXTURE_SIZE，否则绘制结果将是一片黑色。在我的电脑上 GL_MAX_TEXTURE_SIZE=128*128=16384，该最大尺寸限制的是数据单个方向的尺寸，即二维数据的宽和高均不能超过16384。（我最开始还以为是面积）

对于二维地震剖面来说，这已经足够了。单炮数据的道数不会超过这个值，时间采样点数如果超过的话，可以抽稀一下。