0

我正在开发我的体积渲染应用程序（C# + OpenTK）。体积是使用光线投射渲染的，我在这个网站上找到了很多灵感：http: //graphicsrunner.blogspot.sk/2009/01/volume-rendering-101.html，即使我的应用程序使用 OpenGL，使用 3D 纹理和其他东西的主要思想是相同的。应用程序工作正常，但在我“流入体积”（意味着在边界框内）之后，一切都消失了，我想阻止这种情况。那么有一些简单的方法可以做到这一点吗？--> 我将能够在体积中流动或在体积中移动。

这是片段着色器的代码：

#version 330

in vec3 EntryPoint;
in vec4 ExitPointCoord;

uniform sampler2D ExitPoints;
uniform sampler3D VolumeTex;
uniform sampler1D TransferFunc;  
uniform float     StepSize;
uniform float     AlphaReduce;
uniform vec2      ScreenSize;
layout (location = 0) out vec4 FragColor;

void main()
{
//gl_FragCoord --> http://www.txutxi.com/?p=182
vec3 exitPoint = texture(ExitPoints, gl_FragCoord.st/ScreenSize).xyz;

//background need no raycasting
if (EntryPoint == exitPoint)
    discard;

vec3 rayDirection = normalize(exitPoint - EntryPoint);
vec4 currentPosition = vec4(EntryPoint, 0.0f);
vec4 colorSum = vec4(.0f,.0f,.0f,.0f);
vec4 color = vec4(0.0f,0.0f,0.0f,0.0f);
vec4 value = vec4(0.0f);

vec3 Step = rayDirection * StepSize;
float stepLength= length(Step);
float LengthSum = 0.0f;
float Length = length(exitPoint - EntryPoint);

for(int i=0; i < 16000; i++)
{
    currentPosition.w = 0.0f;
    value = texture(VolumeTex, currentPosition.xyz);
    color = texture(TransferFunc, value.a);

    //reduce the alpha to have a more transparent result
    color.a *= AlphaReduce;

    //Front to back blending
    color.rgb *= color.a;
    colorSum = (1.0f - colorSum.a) * color + colorSum;

    //accumulate length
    LengthSum += stepLength;

    //break from the loop when alpha gets high enough
    if(colorSum.a >= .95f)
        break;

    //advance the current position
    currentPosition.xyz += Step;

    //break if the ray is outside of the bounding box
    if(LengthSum >= Length)
        break;
}
FragColor = colorSum;
}

下面的代码基于 https://github.com/toolchainX/Volume_Rendering_Using_GLSL

显示（）函数：

    public void Display()
    {
        // the color of the vertex in the back face is also the location
        // of the vertex
        // save the back face to the user defined framebuffer bound
        // with a 2D texture named `g_bfTexObj`
        // draw the front face of the box
        // in the rendering process, i.e. the ray marching process
        // loading the volume `g_volTexObj` as well as the `g_bfTexObj`
        // after vertex shader processing we got the color as well as the location of
        // the vertex (in the object coordinates, before transformation).
        // and the vertex assemblied into primitives before entering
        // fragment shader processing stage.
        // in fragment shader processing stage. we got `g_bfTexObj`
        // (correspond to 'VolumeTex' in glsl)and `g_volTexObj`(correspond to 'ExitPoints')
        // as well as the location of primitives.

        // draw the back face of the box
        GL.Enable(EnableCap.DepthTest);

        //"vykreslim" front || back face objemu do framebuffru --> teda do 2D textury s ID bfTexID 
        //(pomocou backface.frag &.vert)
        GL.BindFramebuffer(FramebufferTarget.Framebuffer, frameBufferID);
        GL.Viewport(0, 0, width, height);
        LinkShader(spMain.GetProgramHandle(), bfVertShader.GetShaderHandle(), bfFragShader.GetShaderHandle());
        spMain.UseProgram();
        //cull front face
        Render(CullFaceMode.Front);
        spMain.UseProgram(0);
        //klasicky framebuffer --> "obrazovka"
        GL.BindFramebuffer(FramebufferTarget.Framebuffer, 0);

        GL.Viewport(0, 0, width, height);
        LinkShader(spMain.GetProgramHandle(), rcVertShader.GetShaderHandle(), rcFragShader.GetShaderHandle());
        spMain.UseProgram();
        SetUniforms();
        Render(CullFaceMode.Back);
        spMain.UseProgram(0);

        GL.Disable(EnableCap.DepthTest);
    }

    private void DrawBox(CullFaceMode mode)
    {
        // --> Face culling allows non-visible triangles of closed surfaces to be culled before expensive Rasterization and Fragment Shader operations.
        GL.Enable(EnableCap.CullFace);
        GL.CullFace(mode);
        GL.BindVertexArray(VAO);
        GL.DrawElements(PrimitiveType.Triangles, 36, DrawElementsType.UnsignedInt, 0);
        GL.BindVertexArray(0);
        GL.Disable(EnableCap.CullFace);
        spMain.UseProgram(0);//zapnuty bol v Render() ktora DrawBox zavolala
    }

    private void Render(CullFaceMode mode)
    {
        GL.ClearColor(0.0f, 0.0f, 0.0f, 1.0f);
        GL.Clear(ClearBufferMask.ColorBufferBit | ClearBufferMask.DepthBufferBit);
        spMain.UseProgram();
        spMain.SetUniform("modelViewMatrix", Current);
        spMain.SetUniform("projectionMatrix", projectionMatrix);
        DrawBox(mode);
    }

问题是（我认为）当我向音量移动时（我不移动相机，只是缩放音量），如果比例因子> 2.7，我在音量中，这意味着“之后渲染最终图片的飞机”，所以什么都看不到。我能想到的解决方案（也许）是这样的：如果我达到比例因子 = 2.7something：

1.) -> 不要缩放音量

2.) -> 以某种方式告诉片段着色器将 EntryPoint 向 RayDirection 移动一定长度（可能基于比例因子）。

现在，我尝试了这个“方法”，它似乎可以工作：

vec3 entryPoint = EntryPoint + some_value * rayDirection;

some_value 必须被限制在 [0,1[ 区间（或 [0,1]？） 之间，但也许这并不重要：

if (EntryPoint == exitPoint)
discard;

所以现在，也许（如果我的解决方案不是那么糟糕），我可以改变我的答案：如何计算 some_value（基于我发送到片段着色器的比例因子）？

if(scale_factor < 2.7something)
    work like before;
else
{
    compute some_value; //(I need help with this part)
    change entry point;
    work like before;
}

（我不是英语母语者，所以如果文本中有一些大错误并且您不理解某些内容，请告诉我，我会尝试修复这些错误）

谢谢。

4

2 回答 2

1

我解决了我的问题。它不会产生“被体积包围”的错觉，但现在，我可以在体积中流动，什么都没有消失。这是我添加到片段着色器的解决方案的代码：

vec3 entryPoint = vec3(0.0f);

if(scaleCoeff >= 2.7f)
{
    float tmp = min((scaleCoeff - 2.7f) * 0.1f, 1.0f);
    entryPoint = EntryPoint + tmp * (exitPoint - EntryPoint);
}
else
{
    entryPoint = EntryPoint;
}
//

但是，如果您知道或可以考虑更好的解决方案，使“被音量包围”效果，如果您让我知道，我会很高兴。

谢谢你。

于 2016-04-07T14:10:13.723 回答
0

如果理解正确，我认为您应该使用平面剪裁来浏览卷。（如果你附上这个解决方案，我可以根据你的代码给你一个简单的例子。将整个 C++ 项目翻译成 C# 太耗时了。）

于 2019-03-13T01:59:18.710 回答