• Advertisement
Sign in to follow this  

OpenGL Aligning an object to a path

This topic is 2508 days old which is more than the 365 day threshold we allow for new replies. Please post a new topic.

If you intended to correct an error in the post then please contact us.

Recommended Posts

Hey guys,

Long time since I last posted around here.

I'm trying to make a object follow a certain path using opengl and C++.

The path I created is an 8 figure and as we speak, my object can follow the path through translation but I'm having a very hard time figuring out how to make sure the front of my object faces the direction of the path as it moves forward.

Basically, what I'm trying to figure out is how to rotate my object so that it faces the direction of the tangent of the position at which it is at.

The only thing I have ended up with so far is trying to figure out the angle between the 2 points by making those 2 points the ending points of vectors both having the origin as a starting point and calculating the angle between those and rotating my object. It doesn't well at all and I believe I do not have the proper mathematical knowledge to achieve what I want to do.

What should I look for in order to make my object face the proper direction when moving ?

Thanks!

Share this post


Link to post
Share on other sites
Advertisement
Don't compute an angle and rotate by that angle, but use the tangent as forward vector directly. Use the standard method of constructing a orientation matrix from the forward vector and a given up vector (e.g. one that is perpendicular to the figure, or the world up vector, or something else).

How to compute the tangent ... depends on what kind of path you have. If you have a curve that can be differentiated, then you may use the 1st derivative. The mathematical approach can be read in the [url="http://en.wikipedia.org/wiki/Frenet–Serret_formulas"]Frenet Formulas article[/url] on wikipedia. Of course, you can ever step a little bit forward and a little bit backward on the curve and use the difference vector after normalization as approximation.

Share this post


Link to post
Share on other sites
I expected my method was wrong so that is why I decided to start this thread.

I'm currently using a [url="http://en.wikipedia.org/wiki/Lemniscate_of_Bernoulli"]Lemniscate of Bernoulli [/url]for my path. The link also provides both the first and second derivatives but other than that, I have no clue how to obtain the exact information I want or what kind of information.

I find the Frenet-Serret Formulas extremely confusing. The issue also arise from the fact that I have no clue what is my T N and B information from the what I have in my scene or how to obtain it and what to do with it.

I understand the logic behind it from the formula but I can't relate this to the code with my path and objects.

The only thing I understand right now is that my parameter t will provide me with the proper X and Z position on the curve. After that I have no clue numbers I must get from what I have and where to plug them to get the rest of the data I need.

Where should I start?

Share this post


Link to post
Share on other sites
The linked article does show derivatives dx/dy and dy/dx, but you're using the parametric form of the curve. Hence you need dx/dt and dy/dt. Setting in t into the both derivatives gives you the x and y components of the unnormalized tangent (which you are using as x and z components, if I understood your latest post correctly). The Frenet Frame article shows this in the section [url="http://en.wikipedia.org/wiki/Frenet–Serret_formulas#Other_expressions_of_the_frame"]"Other expressions of the frame"[/url], where
T(t) = r'(t) / ||r'(t)||
is mentioned (the r' means the derivative dr/dt).

The T is the normalized tangent vector along the curve, the N is the normal (i.e. it is perpendicular to the curve at location t) with a relation to the curvature of the curve at location t, and B is the bi-normal (i.e. also perpendicular to the curve at t). All together, the point at location t, T, N, B define a Cartesian axis, i.e. a local co-ordinate system similar to one used for e.g. model vertices.

Assume that you have a unit-length forward vector [b]f[/b] and a unit-length up vector [b]u[/b] and both are perpendicular (so [b]f[/b] . [b]u[/b] = [b]0[/b]). Then you can compute a side vector [b]s[/b] using the cross-product
[b]s[/b] := [b]f[/b] x [b]u[/b]
These 3 vectors are pairwise perpendicular and all of unit length. So they build a basis, written as the column vectors (or row vectors, depending on whether you use column or row vector math) of a matrix
[ [b]s f u[/b] ]
and you'll have the rotational part of a local space transformation. Append a position [b]p[/b] and you have the full transformation matrix:
[ [b]s f u p[/b] ]

If the forward vector [b]f[/b] and the up vector [b]u'[/b] are not perpendicular by itself, you can compute the side vector as usual but with a additional normalization
[b]s'[/b] := [b]f[/b] x [b]u'[/b]
[b]s [/b]:=[b] s'[/b] / ||[b]s'[/b]||
and compute the actual up vector from them
[b]u[/b] := [b]s[/b] x [b]f[/b]
[b]
[/b]
So the question is how to get a pair of those vectors.

One way (the simple one) is to compute the forward vector as difference vector between two points on the curve that are both close to the current point [b]p[/b] at location t[sub]0[/sub],
[b]p[/b] := [b]r[/b]( t0 )
so that
[b]d'[/b] := [b]r[/b]( t[sub]0[/sub] + 0.05 ) - [b]r[/b]( t0 - 0.05 )
[b]d[/b] := [b]d'[/b] / ||[b]d'[/b]||
and choosing an up vector [b]u'[/b], e.g. [ 0 1 0 ] in your case where the curve is in x-z. Then follow the steps above. Of course, picking a delta (here 0.05) is a bit of voodoo.

Another, more exact way would be to compute dx/dt and dy/dt from your curve to determine the tangent T and use it as forward vector. Choose an up vector like above and complete the matrix as usual.

The 3rd way would be to compute the 2nd derivatives as well to yield in N, and use that as side vector.


Now it's up to you to decide how precise you want the tangent to be, and whether the up vector should be global or curvature dependent or what else. Its easy to mix up the vectors in the above, so check what I wrote if you want to use it.

Share this post


Link to post
Share on other sites
[quote name='haegarr' timestamp='1299439313' post='4782506']
The linked article does show derivatives dx/dy and dy/dx, but you're using the parametric form of the curve. Hence you need dx/dt and dy/dt. Setting in t into the both derivatives gives you the x and y components of the unnormalized tangent (which you are using as x and z components, if I understood your latest post correctly). The Frenet Frame article shows this in the section [url="http://en.wikipedia.org/wiki/Frenet%E2%80%93Serret_formulas#Other_expressions_of_the_frame"]"Other expressions of the frame"[/url], where
T(t) = r'(t) / ||r'(t)||
is mentioned (the r' means the derivative dr/dt).

The T is the normalized tangent vector along the curve, the N is the normal (i.e. it is perpendicular to the curve at location t) with a relation to the curvature of the curve at location t, and B is the bi-normal (i.e. also perpendicular to the curve at t). All together, the point at location t, T, N, B define a Cartesian axis, i.e. a local co-ordinate system similar to one used for e.g. model vertices.

Assume that you have a unit-length forward vector [b]f[/b] and a unit-length up vector [b]u[/b] and both are perpendicular (so [b]f[/b] . [b]u[/b] = [b]0[/b]). Then you can compute a side vector [b]s[/b] using the cross-product
[b]s[/b] := [b]f[/b] x [b]u[/b]
These 3 vectors are pairwise perpendicular and all of unit length. So they build a basis, written as the column vectors (or row vectors, depending on whether you use column or row vector math) of a matrix
[ [b]s f u[/b] ]
and you'll have the rotational part of a local space transformation. Append a position [b]p[/b] and you have the full transformation matrix:
[ [b]s f u p[/b] ]

If the forward vector [b]f[/b] and the up vector [b]u'[/b] are not perpendicular by itself, you can compute the side vector as usual but with a additional normalization
[b]s'[/b] := [b]f[/b] x [b]u'[/b]
[b]s [/b]:=[b] s'[/b] / ||[b]s'[/b]||
and compute the actual up vector from them
[b]u[/b] := [b]s[/b] x [b]f[/b]
[b]
[/b]
So the question is how to get a pair of those vectors.

One way (the simple one) is to compute the forward vector as difference vector between two points on the curve that are both close to the current point [b]p[/b] at location t[sub]0[/sub],
[b]p[/b] := [b]r[/b]( t0 )
so that
[b]d'[/b] := [b]r[/b]( t[sub]0[/sub] + 0.05 ) - [b]r[/b]( t0 - 0.05 )
[b]d[/b] := [b]d'[/b] / ||[b]d'[/b]||
and choosing an up vector [b]u'[/b], e.g. [ 0 1 0 ] in your case where the curve is in x-z. Then follow the steps above. Of course, picking a delta (here 0.05) is a bit of voodoo.

Another, more exact way would be to compute dx/dt and dy/dt from your curve to determine the tangent T and use it as forward vector. Choose an up vector like above and complete the matrix as usual.

The 3rd way would be to compute the 2nd derivatives as well to yield in N, and use that as side vector.


Now it's up to you to decide how precise you want the tangent to be, and whether the up vector should be global or curvature dependent or what else. Its easy to mix up the vectors in the above, so check what I wrote if you want to use it.
[/quote]

Thanks for the in depth information. I wanted to have tried something before posting back. I have to admit the math isn't too easy for me but manageable on paper but relating this to OpenGL is confusing as hell.

So far my train is following the curve but it's a bit crappy at points. At 4 spots on my curve, the train suddenly flips over 180degrees and still follows the curve backwards and then will flip back again when it reached another one of those points.

I'll explain what I'm doing. Perhaps you could see what is wrong.

At the moment. I keep the position X and Z as global variables by obtaining the value from T in my parametric equations( Y is always 0). Then, using both points from the position, I pass them in their respective derivative to obtain a set of points that is a vector tangent to the X Z position on my curve.

I then normalize the tangent vector by its length but I'm not sure I'm doing this part right. I know the magnitude of the vector is sqrt(x^2+y^2+z^2) but isn't that only we have a point (X,Y,Z) and the other at (0,0,0) ? At the moment, I have that my length is sqrt(X^2+Z^2) since Y is always 0.

Following from this, I calculate the dot product between my new derived points and the vector (-1,0,0) which is the direction vector in which my train is pointing at initially or pretty much all the time before being translated and rotated.

With the dot product result, I calculate the rotation angle = acos(dotproduct)*(180/pi) and use this angle to rotate my train after I translated it at the proper position. It's not perfect though and no matter how long I try to figure out what is wrong I have no clue what else to do.

Also, For the train flipping issue. I just managed to do a kind of hackjob that detects if the angle change is large and toogle a boolean that I check to rotate the train back 180 degree while it is on. I would like to get rid of this too with proper calculations.

I've also tried to have my camera follow the train from the back but I have no idea how to achieve this. The camera can follow the train from it's origin by changing the eye position but if I try to move the camera as well it becomes erratic and goes crazy.

Share this post


Link to post
Share on other sites
Sign in to follow this  

  • Advertisement
  • Advertisement
  • Popular Now

  • Advertisement
  • Similar Content

    • By Balma Alparisi
      i got error 1282 in my code.
      sf::ContextSettings settings; settings.majorVersion = 4; settings.minorVersion = 5; settings.attributeFlags = settings.Core; sf::Window window; window.create(sf::VideoMode(1600, 900), "Texture Unit Rectangle", sf::Style::Close, settings); window.setActive(true); window.setVerticalSyncEnabled(true); glewInit(); GLuint shaderProgram = createShaderProgram("FX/Rectangle.vss", "FX/Rectangle.fss"); float vertex[] = { -0.5f,0.5f,0.0f, 0.0f,0.0f, -0.5f,-0.5f,0.0f, 0.0f,1.0f, 0.5f,0.5f,0.0f, 1.0f,0.0f, 0.5,-0.5f,0.0f, 1.0f,1.0f, }; GLuint indices[] = { 0,1,2, 1,2,3, }; GLuint vao; glGenVertexArrays(1, &vao); glBindVertexArray(vao); GLuint vbo; glGenBuffers(1, &vbo); glBindBuffer(GL_ARRAY_BUFFER, vbo); glBufferData(GL_ARRAY_BUFFER, sizeof(vertex), vertex, GL_STATIC_DRAW); GLuint ebo; glGenBuffers(1, &ebo); glBindBuffer(GL_ELEMENT_ARRAY_BUFFER, ebo); glBufferData(GL_ELEMENT_ARRAY_BUFFER, sizeof(indices), indices,GL_STATIC_DRAW); glVertexAttribPointer(0, 3, GL_FLOAT, false, sizeof(float) * 5, (void*)0); glEnableVertexAttribArray(0); glVertexAttribPointer(1, 2, GL_FLOAT, false, sizeof(float) * 5, (void*)(sizeof(float) * 3)); glEnableVertexAttribArray(1); GLuint texture[2]; glGenTextures(2, texture); glActiveTexture(GL_TEXTURE0); glBindTexture(GL_TEXTURE_2D, texture[0]); glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_CLAMP_TO_EDGE); glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_CLAMP_TO_EDGE); glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_LINEAR); glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_LINEAR); sf::Image* imageOne = new sf::Image; bool isImageOneLoaded = imageOne->loadFromFile("Texture/container.jpg"); if (isImageOneLoaded) { glTexImage2D(GL_TEXTURE_2D, 0, GL_RGBA, imageOne->getSize().x, imageOne->getSize().y, 0, GL_RGBA, GL_UNSIGNED_BYTE, imageOne->getPixelsPtr()); glGenerateMipmap(GL_TEXTURE_2D); } delete imageOne; glActiveTexture(GL_TEXTURE1); glBindTexture(GL_TEXTURE_2D, texture[1]); glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_CLAMP_TO_EDGE); glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_CLAMP_TO_EDGE); glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_LINEAR); glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_LINEAR); sf::Image* imageTwo = new sf::Image; bool isImageTwoLoaded = imageTwo->loadFromFile("Texture/awesomeface.png"); if (isImageTwoLoaded) { glTexImage2D(GL_TEXTURE_2D, 0, GL_RGBA, imageTwo->getSize().x, imageTwo->getSize().y, 0, GL_RGBA, GL_UNSIGNED_BYTE, imageTwo->getPixelsPtr()); glGenerateMipmap(GL_TEXTURE_2D); } delete imageTwo; glUniform1i(glGetUniformLocation(shaderProgram, "inTextureOne"), 0); glUniform1i(glGetUniformLocation(shaderProgram, "inTextureTwo"), 1); GLenum error = glGetError(); std::cout << error << std::endl; sf::Event event; bool isRunning = true; while (isRunning) { while (window.pollEvent(event)) { if (event.type == event.Closed) { isRunning = false; } } glClear(GL_COLOR_BUFFER_BIT); if (isImageOneLoaded && isImageTwoLoaded) { glActiveTexture(GL_TEXTURE0); glBindTexture(GL_TEXTURE_2D, texture[0]); glActiveTexture(GL_TEXTURE1); glBindTexture(GL_TEXTURE_2D, texture[1]); glUseProgram(shaderProgram); } glBindVertexArray(vao); glDrawElements(GL_TRIANGLES, 6, GL_UNSIGNED_INT, nullptr); glBindVertexArray(0); window.display(); } glDeleteVertexArrays(1, &vao); glDeleteBuffers(1, &vbo); glDeleteBuffers(1, &ebo); glDeleteProgram(shaderProgram); glDeleteTextures(2,texture); return 0; } and this is the vertex shader
      #version 450 core layout(location=0) in vec3 inPos; layout(location=1) in vec2 inTexCoord; out vec2 TexCoord; void main() { gl_Position=vec4(inPos,1.0); TexCoord=inTexCoord; } and the fragment shader
      #version 450 core in vec2 TexCoord; uniform sampler2D inTextureOne; uniform sampler2D inTextureTwo; out vec4 FragmentColor; void main() { FragmentColor=mix(texture(inTextureOne,TexCoord),texture(inTextureTwo,TexCoord),0.2); } I was expecting awesomeface.png on top of container.jpg

    • By khawk
      We've just released all of the source code for the NeHe OpenGL lessons on our Github page at https://github.com/gamedev-net/nehe-opengl. code - 43 total platforms, configurations, and languages are included.
      Now operated by GameDev.net, NeHe is located at http://nehe.gamedev.net where it has been a valuable resource for developers wanting to learn OpenGL and graphics programming.

      View full story
    • By TheChubu
      The Khronos™ Group, an open consortium of leading hardware and software companies, announces from the SIGGRAPH 2017 Conference the immediate public availability of the OpenGL® 4.6 specification. OpenGL 4.6 integrates the functionality of numerous ARB and EXT extensions created by Khronos members AMD, Intel, and NVIDIA into core, including the capability to ingest SPIR-V™ shaders.
      SPIR-V is a Khronos-defined standard intermediate language for parallel compute and graphics, which enables content creators to simplify their shader authoring and management pipelines while providing significant source shading language flexibility. OpenGL 4.6 adds support for ingesting SPIR-V shaders to the core specification, guaranteeing that SPIR-V shaders will be widely supported by OpenGL implementations.
      OpenGL 4.6 adds the functionality of these ARB extensions to OpenGL’s core specification:
      GL_ARB_gl_spirv and GL_ARB_spirv_extensions to standardize SPIR-V support for OpenGL GL_ARB_indirect_parameters and GL_ARB_shader_draw_parameters for reducing the CPU overhead associated with rendering batches of geometry GL_ARB_pipeline_statistics_query and GL_ARB_transform_feedback_overflow_querystandardize OpenGL support for features available in Direct3D GL_ARB_texture_filter_anisotropic (based on GL_EXT_texture_filter_anisotropic) brings previously IP encumbered functionality into OpenGL to improve the visual quality of textured scenes GL_ARB_polygon_offset_clamp (based on GL_EXT_polygon_offset_clamp) suppresses a common visual artifact known as a “light leak” associated with rendering shadows GL_ARB_shader_atomic_counter_ops and GL_ARB_shader_group_vote add shader intrinsics supported by all desktop vendors to improve functionality and performance GL_KHR_no_error reduces driver overhead by allowing the application to indicate that it expects error-free operation so errors need not be generated In addition to the above features being added to OpenGL 4.6, the following are being released as extensions:
      GL_KHR_parallel_shader_compile allows applications to launch multiple shader compile threads to improve shader compile throughput WGL_ARB_create_context_no_error and GXL_ARB_create_context_no_error allow no error contexts to be created with WGL or GLX that support the GL_KHR_no_error extension “I’m proud to announce OpenGL 4.6 as the most feature-rich version of OpenGL yet. We've brought together the most popular, widely-supported extensions into a new core specification to give OpenGL developers and end users an improved baseline feature set. This includes resolving previous intellectual property roadblocks to bringing anisotropic texture filtering and polygon offset clamping into the core specification to enable widespread implementation and usage,” said Piers Daniell, chair of the OpenGL Working Group at Khronos. “The OpenGL working group will continue to respond to market needs and work with GPU vendors to ensure OpenGL remains a viable and evolving graphics API for all its customers and users across many vital industries.“
      The OpenGL 4.6 specification can be found at https://khronos.org/registry/OpenGL/index_gl.php. The GLSL to SPIR-V compiler glslang has been updated with GLSL 4.60 support, and can be found at https://github.com/KhronosGroup/glslang.
      Sophisticated graphics applications will also benefit from a set of newly released extensions for both OpenGL and OpenGL ES to enable interoperability with Vulkan and Direct3D. These extensions are named:
      GL_EXT_memory_object GL_EXT_memory_object_fd GL_EXT_memory_object_win32 GL_EXT_semaphore GL_EXT_semaphore_fd GL_EXT_semaphore_win32 GL_EXT_win32_keyed_mutex They can be found at: https://khronos.org/registry/OpenGL/index_gl.php
      Industry Support for OpenGL 4.6
      “With OpenGL 4.6 our customers have an improved set of core features available on our full range of OpenGL 4.x capable GPUs. These features provide improved rendering quality, performance and functionality. As the graphics industry’s most popular API, we fully support OpenGL and will continue to work closely with the Khronos Group on the development of new OpenGL specifications and extensions for our customers. NVIDIA has released beta OpenGL 4.6 drivers today at https://developer.nvidia.com/opengl-driver so developers can use these new features right away,” said Bob Pette, vice president, Professional Graphics at NVIDIA.
      "OpenGL 4.6 will be the first OpenGL release where conformant open source implementations based on the Mesa project will be deliverable in a reasonable timeframe after release. The open sourcing of the OpenGL conformance test suite and ongoing work between Khronos and X.org will also allow for non-vendor led open source implementations to achieve conformance in the near future," said David Airlie, senior principal engineer at Red Hat, and developer on Mesa/X.org projects.

      View full story
    • By _OskaR
      Hi,
      I have an OpenGL application but without possibility to wite own shaders.
      I need to perform small VS modification - is possible to do it in an alternative way? Do we have apps or driver modifictions which will catch the shader sent to GPU and override it?
    • By xhcao
      Does sync be needed to read texture content after access texture image in compute shader?
      My simple code is as below,
      glUseProgram(program.get());
      glBindImageTexture(0, texture[0], 0, GL_FALSE, 3, GL_READ_ONLY, GL_R32UI);
      glBindImageTexture(1, texture[1], 0, GL_FALSE, 4, GL_WRITE_ONLY, GL_R32UI);
      glDispatchCompute(1, 1, 1);
      // Does sync be needed here?
      glUseProgram(0);
      glBindFramebuffer(GL_READ_FRAMEBUFFER, framebuffer);
      glFramebufferTexture2D(GL_READ_FRAMEBUFFER, GL_COLOR_ATTACHMENT0,
                                     GL_TEXTURE_CUBE_MAP_POSITIVE_X + face, texture[1], 0);
      glReadPixels(0, 0, kWidth, kHeight, GL_RED_INTEGER, GL_UNSIGNED_INT, outputValues);
       
      Compute shader is very simple, imageLoad content from texture[0], and imageStore content to texture[1]. Does need to sync after dispatchCompute?
  • Advertisement