// 作業用メモリ
    std::vector<GLubyte> texture;
    
    // ノイズ関数を初期化する
    const Noise3 noise(5, 5, 5);
    
    // ノイズ関数を使ってテクスチャを作る
    for (GLint k = 0; k < depth; ++k)
    {
      const double z(double(k) / double(depth - 1));
      
      for (GLint j = 0; j < height; ++j)
      {
        const double y(double(j) / double(height - 1));
        
        for (GLint i = 0; i < width; ++i)
        {
          const double x(double(i) / double(width - 1));
          
#if VOLUMETYPE == CHECKER
          texture.push_back(((i >> 4) + (j >> 4) + (k >> 4)) & 1 ? 0 : 255);
#elif VOLUMETYPE == NOISE
          texture.push_back(GLubyte(noise.noise(x, y, z) * 255.0));
#elif VOLUMETYPE == PERLIN
          texture.push_back(GLubyte(noise.perlin(x, y, z, 4, 0.5) * 255.0));
#elif VOLUMETYPE == TURBULENCE
          texture.push_back(GLubyte(noise.turbulence(x, y, z, 4, 0.5) * 255.0));
#elif VOLUMETYPE == SPHERE
          const double px(2.0 * x - 1.0), py(2.0 * y - 1.0), pz(2.0 * z - 1.0);
          texture.push_back(GLubyte(255.0 - sqrt(px * px + py * py + pz * pz) * 127.5));
#else
          texture.push_back(255);
#endif
        }
      }
    }

    // テクスチャオブジェクトを作成して結合する
    GLuint tex;
    glGenTextures(1, &tex);
    glBindTexture(GL_TEXTURE_3D, tex);
    
    // テクスチャを割り当てる
    glTexImage3D(GL_TEXTURE_3D, 0, GL_R8, width, height, depth, 0,
      GL_RED, GL_UNSIGNED_BYTE, &texture[0]);
    
    // テクスチャの拡大・縮小に線形補間を用いる
    glTexParameteri(GL_TEXTURE_3D, GL_TEXTURE_MAG_FILTER, GL_LINEAR);
    glTexParameteri(GL_TEXTURE_3D, GL_TEXTURE_MIN_FILTER, GL_LINEAR);
    
    // テクスチャからはみ出た部分には境界色を用いる
    glTexParameteri(GL_TEXTURE_3D, GL_TEXTURE_WRAP_S, GL_CLAMP_TO_BORDER);
    glTexParameteri(GL_TEXTURE_3D, GL_TEXTURE_WRAP_T, GL_CLAMP_TO_BORDER);
    glTexParameteri(GL_TEXTURE_3D, GL_TEXTURE_WRAP_R, GL_CLAMP_TO_BORDER);
    
    // テクスチャの境界色を設定する (ボリュームの外には何もない)
    static const GLfloat black[] = { 0.0f, 0.0f, 0.0f, 0.0f };
    glTexParameterfv(GL_TEXTURE_3D, GL_TEXTURE_BORDER_COLOR, black);

    // スライスの描画
    glBindVertexArray(slice);
    glDrawArraysInstanced(GL_TRIANGLE_FAN, 0, 4, slices);

#version 150 core
#extension GL_ARB_explicit_attrib_location : enable
 
// テクスチャ座標の変換行列
uniform mat4 mt;
 
// モデルビュー変換行列
uniform mat4 mw;
 
// プロジェクション変換行列
uniform mat4 mp;
 
// スライスの間隔
uniform float spacing;
 
// [-0.5, 0.5] の正方形の頂点位置
layout (location = 0) in vec2 pv;
 
// スライスの頂点位置
out vec4 p;
 
// スライスのテクスチャ座標
out vec3 t;
 
void main()
{
  // スライスを gl_InstanceID でずらす
  p = vec4(pv, (float(gl_InstanceID) + 0.5) * spacing - 0.5, 1.0);

  // スライスのテクスチャ座標はスライスの中心を基準に √3 倍に拡大してから回転する
  t = (mat3(mt) * p.xyz) * 1.732 + 0.5;
  
  // 頂点位置を視点座標系に移す
  p = mw * p;
  
  // モデルビュープロジェクション変換をしてからラスタライザに送る
  gl_Position = mp * p;
}

... (中略) ...
 
// テクスチャのサンプラ
uniform sampler3D volume;
 
// クリッピング座標系への変換行列
uniform mat4 mc;
 
// 閾値
uniform float threshold;
 
// スライスの表面上の位置
in vec4 p;
 
// テクスチャ座標
in vec3 t;
 
// フレームバッファに出力するデータ
layout (location = 0) out vec4 fc;
 
void main()
{
  // 元のボリュームの濃度と閾値の差
  float v = texture(volume, t).r - threshold;
  
  // 濃度が閾値以下ならフラグメントを捨てる
  if (v <= 0.0) discard;

  // 濃度の勾配を求める
  vec3 g = vec3(
    textureOffset(volume, t, ivec3(-1, 0, 0)).r - textureOffset(volume, t, ivec3(1, 0, 0)).r,
    textureOffset(volume, t, ivec3(0, -1, 0)).r - textureOffset(volume, t, ivec3(0, 1, 0)).r,
    textureOffset(volume, t, ivec3(0, 0, -1)).r - textureOffset(volume, t, ivec3(0, 0, 1)).r
  );

  // 勾配をそのままフラグメントの色に使う
  fc = vec4(normalize(g) * 0.5 + 0.5, v);

  vec3 l = normalize((pl * p.w - p * pl.w).xyz);  // 光線ベクトル
  vec3 n = normalize(g * mat3(mt));               // 法線ベクトル
  vec3 h = normalize(l - normalize(p.xyz));       // 中間ベクトル
  
  // 拡散反射光＋環境光の反射光
  vec4 idiff = max(dot(n, l), 0.0) * kdiff * ldiff + kamb * lamb;
  
  // 鏡面反射光
  vec4 ispec = pow(max(dot(n, h), 0.0), kshi) * kspec * lspec;
  
  // フラグメントの色
  fc = vec4((idiff + ispec).rgb, v);
}