Jump to content
  • Advertisement
Sign in to follow this  

OpenGL marching cube in opengl

This topic is 3753 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

hi guys... i have a big problem...i spent a lot of time on this topic...but i didn't get the right path... i want to draw triangular marching cube...(a famous marching cube algorithm) but i don't know how...which type of data is used in marching cube algorithm... i know...first of all we need to convert my model into slices...then draw voxels... but how to convert slices...and what type of data should be exported from slices for draw the voxels...means slices should be PNG(image formate) or ASCII data...???? or what kind of data should be exported for march the cube... i have a triangular model and i want to convert that model into voxel model... how can i convert my model into slices..??? i would appreciate you kind help... much thanks in advance...

Share this post

Link to post
Share on other sites

To be honest I'm not really sure what you mean, but when I was experimenting with marching cubes, I got a lot of help here.

kind regards

Share this post

Link to post
Share on other sites
yes...i already read that article and got a good concept...

please only tell me what kind of data is used in that article for marching cube...

is that ASCII data...means (XYZ vertices and indices) or other type...

i found a one file is there name(cell.gz) when i extract this file (9mb file size) exported. what type of that file...there is no extension of that file...

and also brother i couldn't compile this (http://local.wasp.uwa.edu.au/~pbourke/geometry/polygonise/source/) opengl program...

it needs another files...

please if you have this compilable complete marching cube program...please give me or some references...

UncleRemus, you were working on marching cube...please sir i want to contact you...you can help me...if you have some time...please tell me your Email...

much thanks...

Share this post

Link to post
Share on other sites

Marching Cubes uses a scalar field to denote the relative density across a 3D space.

Converting a mesh to a scalar field (aka voxel field) is not straight-forward. The most basic way would be to establish a regular grid, then step through it cell by cell. If the cell is "inside" of the mesh, you can mark it as being dense (a scalar value of 1). If the cell is "outside" of the mesh, you can mark it as being sparse (a scalar value of 0). The key thing with Marching Cubes is that the scalar field can have any value within the range of 0 to 1. So when you are converting your mesh to a grid, you will eventually also have to take into consideration that certain cells will be both "inside" and "outside" of the mesh, and that you will have to calculate the percentage that the cell is inside of the mesh, so that you can come up with a value that is appropriately in-between 0 and 1.

The main function of Marching Cubes is then to take this scalar field and convert it to a mesh.

My page has an example program:

I am also wondering why you would want to take a mesh, convert it to a scalar field, then convert it back to a mesh using Marching Cubes. It seems like a lot of work, where the result is simply a mesh that looks bad compared to the original. The main purpose of Marching Cubes it to convert density data to a mesh, ex: information gained from sonography, tomography, seismology, implicit functions, etc.

Share this post

Link to post
Share on other sites
much Thanks taby..

yes i can get a lot of concept from your example code...

my problem is this...

I have a triangular surface model data...(triangular vertices + indices)

i imported this data and draw in OpenGL...

now i want to draw that model with marching cube algorithm...

i read your volume text file...but i couldn't understand...from where did you get that file...and what is that data...i just understand header of that data file...(size of grids, resolution, data type etc)..

but what is lower data...from where did you export that data...

Please tell me ...how can i convert my triangular surface model into volume data...

means if i use your example program with my model data...how can i convert my model data as same as your (input.txt)...grid data..

please help me...


Share this post

Link to post
Share on other sites
Hi Lucky,

The data beyond the header is the density information for the cells. I believe the header says that the grid is 10x10x10, so there should be 1000 density numbers after the header (one per cell). So the first density number should represent cell {0, 0, 0}, the second one should represent cell {0, 0, 1}, and the last one should represent cell {9, 9, 9}.

I exported this density data from a fluid simulation. Basically you are seeing a portion of the top of a cloud.

The first step to converting your model to a scalar field is finding out how to detect whether a point {x, y, z} is "inside" the model, or "outside". I'm not an expert at this, so I recommend that you create a new topic called something like "detecting whether a point is inside/outside a mesh". Once you have figured this out, come back to this thread and we can continue.

Share this post

Link to post
Share on other sites
thanks you very much sir...

now i am clear...

actually sometime i thought about this kind of thing...and i studied some papers related to detection of points inside or outside of mesh...

i believe, i will do this in some coming days...then again...i will come back...

please sir guide me...otherwise it would be very difficult for me...

thanks a lot...

Share this post

Link to post
Share on other sites
Sounds good. If you make a reply and I don't respond within a day or two, just send me a private message to remind me. :)

Share this post

Link to post
Share on other sites
thanks taby...

yes... i will remind you...

i am trying to do your approach...( triangular surface mesh data(IDs, x,y,z) into 3D scalar volume data(density, dimensions, type of data, 0-1 array, etc)

but can you clear me one thing..???

i spent a lot of time for searching this kind of program or paper..in which anybody draw voxel with surface data(ID, x,y,z)...
but i did not find any material about this...

i saw, many author used marching cube for displaying 3D model...i saw a lot of images on the net...my question is this...what kind of data they use in marching cube algorithm...if they use scalar data as you described...from where they get that data...i think they must would have exported from somewhere/from some kind of software/or algorithm as you did in your program file...

like this... below link...

in above link three model you can see...these model drawn with marching cube algorithm...how they draw these model...???from where they did get scalar data...if they use some kind of scanner then which scanner.???..if they produce scalar data themself..then which method they choose for displaying voxels...???
are they use your method as you explained above...detect inside or outside points...or anyother...

here is a one example... if you want draw a tooth model for some kind of simulation, not you anybody...and you have triangular surface data of that tooth, scanned from some kind of 3D scanner...

now you want to convert that model into voxel model(volume model)...

what will you do.??...means in your case...which method you choose for display your model as voxels...

now i think you can fully understand my problem...means we don't want to use any machine for exporting scalar data of tooth model...

in case of MetaBalls...its easy to create scalar field of sphere...many author did this...

but in my case ...i have a triangular surface(face) data of tooth...
please firstly you suggest me which method i need to choose for display surface model into voxel model...which method is suitable for tooth model...every tooth model have millions of vertices...

one is another method in which we convert surface model into slices and then draw into voxel model...

please suggest me...which method is suitable for me...


[Edited by - luckyyyyyy on March 10, 2008 7:21:52 AM]

Share this post

Link to post
Share on other sites

A common source of density data is tomography and imaging:

Rendering volumetric information (density data) is a fairly common topic these days, especially with fluid dynamics becoming so popular. Jos Stam and so many others have written good papers about raytracing semi-opaque materials.

Marching Cubes is used if you want to go one step further and compute the triangulated shell of the density. This can be used for things like collision detection, volume calculation, etc. Most importantly, having a set of triangles allows you to render the object using standard GLSL shaders (or even the fixed-function pipeline if you really must).

If a particular type of scanner generates a triangular mesh by default, then Marching Cubes is no longer necessary.

Share this post

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

  • Advertisement
  • Advertisement
  • Popular Tags

  • Similar Content

    • By McGrane
      My laptop recently decided to die, so Ive been transferring my project to my work laptop just to get it up to date, and commit it. I was banging my head against the wall all day, as my textures where not displaying in my program- I was getting no errors and no indication of why it was occurring so I have been just trying to figure it out- I know the image loading was working ok, as im using image data elsewhere, I was pretty confident that the code was fine also, as ive never had an issue with displaying textures before, so I thought it might be the drivers on this laptop, (my old one was just using the built in IntelHD, while this laptop has a NVIDIA graphics card) but all seems to be up to date.
      Below are my basic shaders:
      Vertex Shader
      #version 330 core layout(location = 0) in vec3 position; layout(location = 1) in vec3 color; layout(location = 2) in vec3 normal; layout(location = 3) in vec2 texCoord; uniform mat4 Projection; uniform mat4 Model; out vec3 Color; out vec3 Normal; out vec2 TexCoord; void main() { gl_Position = Projection * Model * vec4( position, 1.0 ); Color = color; Normal = normal; TexCoord = vec2( texCoord.x, texCoord.y); } Fragment Shader
      #version 330 core in vec3 Color; in vec3 Normal; in vec2 TexCoord; uniform sampler2D textureData; void main() { vec4 textureColor = texture( textureData, TexCoord ); vec4 finalColor = textureColor * vec4( Color, 1.0f); gl_FragColor = finalColor; } Calling Code
      glActiveTexture(GL_TEXTURE0); glBindTexture(GL_TEXTURE_2D, textureID); glUniform1i(glGetUniformLocation(shaderID, "textureData"), textureID); Now this is the part i dont understand, I worked through my program, until I got to the above 'Calling Code'. This just displays a black texture.. my original issue. Out of desperation, I just tried changing the name in glGetUniformLocation from "textureData" to "textureData_invalid" to see if my error checks would through up something, but in actual fact, it is now displaying the texture as expected. Can anyone fathom a guess as too why this is occurring.. im assuming the random text is just picking up the correct location by c++ witchcraft, but why is the original one not getting picked up correctly and/or not working as expected
      I realize more code is probably needed to see how it all hangs together.. but it seems to come down to this as the issue
    • By QQemka
      Hello. So far i got decently looking 3d scene. I also managed to render a truetype font, on my way to implementing gui (windows, buttons and textboxes). There are several issues i am facing, would love to hear your feedback.
      1) I render text using atlas with VBO containing x/y/u/v of every digit in the atlas (calculated basing on x/y/z/width/height/xoffset/yoffset/xadvance data in binary .fnt format file, screenshot 1). I generated a Comic Sans MS with 32 size and times new roman with size 12 (screenshot 2 and 3). The first issue is the font looks horrible when rescaling. I guess it is because i am using fixed -1 to 1 screen space coords. This is where ortho matrix should be used, right?
      2) Rendering GUI. Situation is similar to above. I guess the widgets should NOT scale when scaling window, am i right? So what am i looking for is saying "this should be always in the middle, 200x200 size no matter the display window xy", and "this should stick to the bottom left corner". Is ortho matrix the cure for all such problems?
      3) The game is 3D but i have to go 2D to render static gui elements over the scene - and i want to do it properly! At the moment i am using matrix 3x3 for 2d transformations and vec3 for all kinds of coordinates. In shaders tho i technically still IS 3D. I have to set all 4 x y z w of the gl_Position while it would be much much more conventient to... just do the maths in 2d space. Can i achieve it somehow?
      4) Text again. I am kind of confused what is the reason of artifacts in Times New Roman font displaying (screenshot 1). I render from left to right, letter after letter. You can clearly see that letters on the right (so the ones rendered after ones on the left are covered by the previous one). I was toying around with blending options but no luck. I do not support kerning at the moment but that's definitely not the cause of error. The display of the small font looks dirty aliased too. I am knd of confused how to interpret the integer data and how should be scaled/adapted to the screen view. Is it just store the data as constant size and again - use ortho matrix?
      Thanks in advance for all your ideas and suggestions!
    • By plz717
      Hello, everyone! I hope my problem isn't too 'beginnerish'. I'm doing research on motion synthesis now, trying to implement the Deep Mimic paper (DeepMimic) by BINPENG XUE, in this paper, I need to first retarget character A's motion to another character B to make the reference motion clips for character B, since we don't have character B‘s reference motion. The most important thing is that in the paper, the author copied character A's joint's rotation with respective to joint's local coordinate system (not the parent) to character B. In my personal understanding, the joint's rotation with respective to joint's local coordinate system is something like that in the attached photo, where for the Elbow joint, i need to get the Elbow's rotation in the elbow's local coordinate system (i'm very grateful for you to share your ideas if i have misunderstanding about it 🙂)
      I have searched many materials on the internet about how to extract the local joint's information from FBX, the most relative one i found is the pivot rotation( and geometric transformation, object offset transformation). I'm a beginner in computer graphics, and i'm confused about whether the pivot rotation( or geometric transformation, object offset transformation) is exactly the joint's local rotation i'm seeking? I hope someone that have any ideas can help me, I'd be very grateful for any pointers in the right direction. Thanks in advance! 

    • By nOoNEE
      hello guys , i have some questions  what does glLinkProgram  and  glBindAttribLocation do?  i searched but there wasnt any good resource 
    • By owenjr
      Hi, I'm a Multimedia Engineering student. I am about to finish my dergree and I'm already thinking about what topic to cover in my final college project.
      I'm interested in the procedural animation with c++ and OpenGL of creatures, something like a spider for example. Can someone tell me what are the issues I should investigate to carry it out? I understand that it has some dependence on artificial intelligence but I do not know to what extent. Can someone help me to find information about it? Thank you very much.
      - Procedural multi-legged walking animation
      - Procedural Locomotion of Multi-Legged Characters in Dynamic Environments
    • By Lewa
      So, i'm still on my quest to unterstanding the intricacies of HDR and implementing this into my engine. Currently i'm at the step to implementing tonemapping. I stumbled upon this blogposts:
      and tried to implement some of those mentioned tonemapping methods into my postprocessing shader.
      The issue is that none of them creates the same results as shown in the blogpost which definitely has to do with the initial range in which the values are stored in the HDR buffer. For simplicity sake i store the values between 0 and 1 in the HDR buffer (ambient light is 0.3, directional light is 0.7)
      This is the tonemapping code:
      vec3 Uncharted2Tonemap(vec3 x) { float A = 0.15; float B = 0.50; float C = 0.10; float D = 0.20; float E = 0.02; float F = 0.30; return ((x*(A*x+C*B)+D*E)/(x*(A*x+B)+D*F))-E/F; } This is without the uncharted tonemapping:
      This is with the uncharted tonemapping:
      Which makes the image a lot darker.
      The shader code looks like this:
      void main() { vec3 color = texture2D(texture_diffuse, vTexcoord).rgb; color = Uncharted2Tonemap(color); //gamma correction (use only if not done in tonemapping code) color = gammaCorrection(color); outputF = vec4(color,1.0f); } Now, from my understanding is that tonemapping should bring the range down from HDR to 0-1.
      But the output of the tonemapping function heavily depends on the initial range of the values in the HDR buffer. (You can't expect to set the sun intensity the first time to 10 and the second time to 1000 and excpect the same result if you feed that into the tonemapper.) So i suppose that this also depends on the exposure which i have to implement?
      To check this i plotted the tonemapping curve:
      You can see that the curve goes only up to around to a value of 0.21 (while being fed a value of 1) and then basically flattens out. (which would explain why the image got darker.)
      My guestion is: In what range should the values in the HDR buffer be which then get tonemapped? Do i have to bring them down to a range of 0-1 by multiplying with the exposure?
      For example, if i increase the values of the light by 10 (directional light would be 7 and ambient light 3) then i would need to divide HDR values by 10 in order to get a value range of 0-1 which then could be fed into the tonemapping curve. Is that correct?
    • By nOoNEE
      i am reading this book : link
      in the OpenGL Rendering Pipeline section there is a picture like this: link
      but the question is this i dont really understand why it is necessary to turn pixel data in to fragment and then fragment into pixel could please give me a source or a clear Explanation that why it is necessary ? thank you so mu
    • By Inbar_xz
      I'm using the OPENGL with eclipse+JOGL.
      My goal is to create movement of the camera and the player.
      I create main class, which create some box in 3D and hold 
      an object of PlayerAxis.
      I create PlayerAxis class which hold the axis of the player.
      If we want to move the camera, then in the main class I call to 
      the func "cameraMove"(from PlayerAxis) and it update the player axis.
      That's work good.
      The problem start if I move the camera on 2 axis, 
      for example if I move with the camera right(that's on the y axis)
      and then down(on the x axis) -
      in some point the move front is not to the front anymore..
      In order to move to the front, I do
      player.playerMoving(0, 0, 1);
      And I learn that in order to keep the front move, 
      I need to convert (0, 0, 1) to the player axis, and then add this.
      I think I dont do the convert right.. 
      I will be glad for help!

      Here is part of my PlayerAxis class:
      //player coordinate float x[] = new float[3]; float y[] = new float[3]; float z[] = new float[3]; public PlayerAxis(float move_step, float angle_move) { x[0] = 1; y[1] = 1; z[2] = -1; step = move_step; angle = angle_move; setTransMatrix(); } public void cameraMoving(float angle_step, String axis) { float[] new_x = x; float[] new_y = y; float[] new_z = z; float alfa = angle_step * angle; switch(axis) { case "x": new_z = addVectors(multScalar(z, COS(alfa)), multScalar(y, SIN(alfa))); new_y = subVectors(multScalar(y, COS(alfa)), multScalar(z, SIN(alfa))); break; case "y": new_x = addVectors(multScalar(x, COS(alfa)), multScalar(z, SIN(alfa))); new_z = subVectors(multScalar(z, COS(alfa)), multScalar(x, SIN(alfa))); break; case "z": new_x = addVectors(multScalar(x, COS(alfa)), multScalar(y, SIN(alfa))); new_y = subVectors(multScalar(y, COS(alfa)), multScalar(x, SIN(alfa))); } x = new_x; y = new_y; z = new_z; normalization(); } public void playerMoving(float x_move, float y_move, float z_move) { float[] move = new float[3]; move[0] = x_move; move[1] = y_move; move[2] = z_move; setTransMatrix(); float[] trans_move = transVector(move); position[0] = position[0] + step*trans_move[0]; position[1] = position[1] + step*trans_move[1]; position[2] = position[2] + step*trans_move[2]; } public void setTransMatrix() { for (int i = 0; i < 3; i++) { coordiTrans[0][i] = x[i]; coordiTrans[1][i] = y[i]; coordiTrans[2][i] = z[i]; } } public float[] transVector(float[] v) { return multiplyMatrixInVector(coordiTrans, v); }  
      and in the main class i have this:
      public void keyPressed(KeyEvent e) { if (e.getKeyCode()== KeyEvent.VK_ESCAPE) { System.exit(0); //player move } else if (e.getKeyCode()== KeyEvent.VK_W) { //front //moveAmount[2] += -0.1f; player.playerMoving(0, 0, 1); } else if (e.getKeyCode()== KeyEvent.VK_S) { //back //moveAmount[2] += 0.1f; player.playerMoving(0, 0, -1); } else if (e.getKeyCode()== KeyEvent.VK_A) { //left //moveAmount[0] += -0.1f; player.playerMoving(-1, 0, 0); } else if (e.getKeyCode()== KeyEvent.VK_D) { //right //moveAmount[0] += 0.1f; player.playerMoving(1, 0, 0); } else if (e.getKeyCode()== KeyEvent.VK_E) { //moveAmount[0] += 0.1f; player.playerMoving(0, 1, 0); } else if (e.getKeyCode()== KeyEvent.VK_Q) { //moveAmount[0] += 0.1f; player.playerMoving(0, -1, 0); //camera move } else if (e.getKeyCode()== KeyEvent.VK_I) { //up player.cameraMoving(1, "x"); } else if (e.getKeyCode()== KeyEvent.VK_K) { //down player.cameraMoving(-1, "x"); } else if (e.getKeyCode()== KeyEvent.VK_L) { //right player.cameraMoving(-1, "y"); } else if (e.getKeyCode()== KeyEvent.VK_J) { //left player.cameraMoving(1, "y"); } else if (e.getKeyCode()== KeyEvent.VK_O) { //right round player.cameraMoving(-1, "z"); } else if (e.getKeyCode()== KeyEvent.VK_U) { //left round player.cameraMoving(1, "z"); } }  
      finallt found it.... i confused with the transformation matrix row and col. thanks anyway!
    • By Lewa
      So, i'm currently trying to implement an SSAO shader from THIS tutorial and i'm running into a few issues here.
      Now, this SSAO method requires view space positions and normals. I'm storing the normals in my deferred renderer in world-space so i had to do a conversion and reconstruct the position from the depth buffer.
      And something there goes horribly wrong (which has probably to do with worldspace to viewspace transformations).
      (here is the full shader source code if someone wants to take a look at it)
      Now, i suspect that the normals are the culprit.
      vec3 normal = ((uNormalViewMatrix*vec4(normalize(texture2D(sNormals, vTexcoord).rgb),1.0)).xyz); "sNormals" is a 2D texture which stores the normals in world space in a RGB FP16 buffer.
      Now i can't use the camera viewspace matrix to transform the normals into viewspace as the cameras position isn't set at (0,0,0), thus skewing the result.
      So what i did is to create a new viewmatrix specifically for this normal without the position at vec3(0,0,0);
      //"camera" is the camera which was used for rendering the normal buffer renderer.setUniform4m(ressources->shaderSSAO->getUniform("uNormalViewMatrix"), glmExt::createViewMatrix(glm::vec3(0,0,0),camera.getForward(),camera.getUp())//parameters are (position,forwardVector,upVector) ); Though i have the feeling this is the wrong approach. Is this right or is there a better/correct way of transforming a world space normal into viewspace?
    • By HawkDeath
      I'm trying mix two textures using own shader system, but I have a problem (I think) with uniforms.
      Code: https://github.com/HawkDeath/shader/tree/test
      To debug I use RenderDocs, but I did not receive good results. In the first attachment is my result, in the second attachment is what should be.
      PS. I base on this tutorial https://learnopengl.com/Getting-started/Textures.

  • Advertisement
  • Popular Now

  • Forum Statistics

    • Total Topics
    • Total Posts

Important Information

By using GameDev.net, you agree to our community Guidelines, Terms of Use, and Privacy Policy.

Participate in the game development conversation and more when you create an account on GameDev.net!

Sign me up!