【转】OpenGL 显示DICOM医学图像

OpenGL 显示DICOM医学图像

                                                                                                                                                关键词：OpenGL、GLSL、DICOM

一、数据读取

        DICOM文件的读取使用DCMTK库，DcmFileFormat 成员函数findAndGetUint16Array(DCM_PixelData, pPixels)；pPixels是指向数据的指针。

二、创建查找表并保存为纹理

glGenTextures(1, &lutTexture);
glBindTexture(GL_TEXTURE_1D,lutTexture);
glTexParameteri(GL_TEXTURE_1D, GL_TEXTURE_MAG_FILTER, GL_LINEAR);
glTexParameteri(GL_TEXTURE_1D, GL_TEXTURE_MIN_FILTER, GL_LINEAR);
        glTexParameteri(GL_TEXTURE_1D, GL_TEXTURE_WRAP_S, GL_CLAMP_TO_EDGE);
glTexImage1D(GL_TEXTURE_1D, 0, GL_LUMINANCE16, texSize, 0, GL_LUMINANCE, GL_UNSIGNED_SHORT, lut16U);

                     lut16U是查找表指针。

三、将图像数据保存为纹理

          glPixelStorei(GL_UNPACK_ALIGNMENT,2);
glGenTextures(1, &imgTexture);
glBindTexture(GL_TEXTURE_2D, imgTexture);
glTexEnvi(GL_TEXTURE_ENV, GL_TEXTURE_ENV_MODE, GL_REPLACE);
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_LINEAR);
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_LINEAR);
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_CLAMP_TO_EDGE);
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_CLAMP_TO_EDGE);
glTexImage2D(GL_TEXTURE_2D, 0, GL_LUMINANCE16, img_width, img_height,0, GL_LUMINANCE, GL_UNSIGNED_SHORT,pData);

                     pData是数据指针

三、着色器代码

void main()
{
    vec4 finalCol;
    float tmpColor;

    vec4 color = vec4(texture2D(mainTexture, gl_TexCoord[0].st)) * multiplier;
    color.rgb -= minVal;

    if (invertGray == true) {
        finalCol.r = 1-texture1D(lutTexture, color.r).r;
        finalCol.g = 1-texture1D(lutTexture, color.g).g;
        finalCol.b = 1-texture1D(lutTexture, color.b).b;
        finalCol.a = color.a;
    }
    else{
        finalCol.r = texture1D(lutTexture, color.r).r;
        finalCol.g = texture1D(lutTexture, color.g).g;
        finalCol.b = texture1D(lutTexture, color.b).b;
        finalCol.a = color.a;
    }
    gl_FragColor = finalCol;
}

 

已实现功能：右键调节窗宽窗位、反色、缩放、平移；

待实现功能：测量工具、汉字显示、高级应用。

 

 

 

本文来源：http://blog.csdn.net/yunzhaoqiang/article/details/8603463